Linking globular clusters (GCs) to the assembly of their host galaxies is an overarching goal in GC studies. The inference of tight scaling relations between GC system properties and the mass of both the stellar and dark halo components of their host galaxies are indicative of an intimate physical connection, yet have also raised fundamental questions about how and when GCs form. Specifically, the inferred correlation between the mass of a GC system (M_gc) and the dark matter halo mass (M_halo) of a galaxy has been posited as a consequence of a causal relation between the formation of dark matter mini-haloes and GC formation during the early epochs of galaxy assembly. We present the first results from a new simulation of a cosmological volume (L=34.4cMpc on a side) from the E-MOSAICS suite, which includes treatments of the formation and evolution of GCs within the framework of a detailed galaxy formation model. The simulated M_gc-M_halo relation is linear for halo masses $>5\times10^{11}$ Msun, and is driven by the hierarchical assembly of galaxies, in agreement with previous studies. Below this halo mass, the simulated relation features a downturn, which we show is consistent with observations, and is driven by the underlying stellar mass (M_star)-halo mass relation of galaxies. Our fiducial model reproduces the observed M_gc-M_star relation across the full mass range, which we argue is more physically relevant than the M_gc-M_halo relation. We also explore the physical processes driving the observed constant value of M_gc / M_halo $\sim 5\times10^{-5}$ and find that it is the result of a combination of cluster formation physics and cluster disruption.

Globular clusters (GCs) are typically old, with most having formed at z >~ 2. This makes understanding their birth environments difficult, as they are typically too distant to observe with sufficient angular resolution to resolve GC birth sites. Using 25 cosmological zoom-in simulations of Milky Way-like galaxies from the E-MOSAICS project, with physically-motivated models for star formation, feedback, and the formation, evolution, and disruption of GCs, we identify the birth environments of present-day GCs. We find roughly half of GCs in these galaxies formed in-situ (52.0 +/- 1.0 per cent) between z ~ 2 - 4, in turbulent, high-pressure discs fed by gas that was accreted without ever being strongly heated through a virial shock or feedback. A minority of GCs form during mergers (12.6 +/- 0.6 per cent in major mergers, and 7.2 +/- 0.5 per cent in minor mergers), but we find that mergers are important for preserving the GCs seen today by ejecting them from their natal, high density interstellar medium (ISM), where proto-GCs are rapidly destroyed due to tidal shocks from ISM substructure. This chaotic history of hierarchical galaxy assembly acts to mix the spatial and kinematic distribution of GCs formed through different channels, making it difficult to use observable GC properties to distinguish GCs formed in mergers from ones formed by smooth accretion, and similarly GCs formed in-situ from those formed ex-situ. These results suggest a simple picture of GC formation, in which GCs are a natural outcome of normal star formation in the typical, gas-rich galaxies that are the progenitors of present-day galaxies.

We present a detailed comparison of the Milky Way (MW) globular cluster (GC) kinematics with the 25 Milky Way-mass cosmological simulations from the E-MOSAICS project. While the MW falls within the kinematic distribution of GCs spanned by the simulations, the relative kinematics of its metal-rich ($[\rm{Fe/H}]>-1.2$) versus metal-poor ($[\rm{Fe/H}]<-1.2$), and inner ($r<8\rm{kpc}$) versus outer ($r>8\rm{kpc}$) populations are atypical for its mass. To understand the origins of these features, we perform a comprehensive statistical analysis of the simulations, and find 18 correlations describing the assembly of $L^*$ galaxies and their dark matter haloes based on their GC population kinematics. The correlations arise because the orbital distributions of accreted and in-situ GCs depend on the masses and accretion redshifts of accreted satellites, driven by the combined effects of dynamical fraction, tidal stripping, and dynamical heating. Because the kinematics of in-situ/accreted GCs are broadly traced by the metal-rich/metal-poor and inner/outer populations, the observed GC kinematics are a sensitive probe of galaxy assembly. We predict that relative to the population of $L^*$ galaxies, the MW assembled its dark matter and stellar mass rapidly through a combination of in-situ star formation, more than a dozen low-mass mergers, and $1.4\pm1.2$ early ($z=3.1\pm1.3$) major merger. The rapid assembly period ended early, limiting the fraction of accreted stars. We conclude by providing detailed quantitative predictions for the assembly history of the MW.

Globular clusters (GCs) formed when the Milky Way experienced a phase of rapid assembly. We use the wealth of information contained in the Galactic GC population to quantify the properties of the satellite galaxies from which the Milky Way assembled. To achieve this, we train an artificial neural network on the E-MOSAICS cosmological simulations of the co-formation and co-evolution of GCs and their host galaxies. The network uses the ages, metallicities, and orbital properties of GCs that formed in the same progenitor galaxies to predict the stellar masses and accretion redshifts of these progenitors. We apply the network to Galactic GCs associated with five progenitors: {\it Gaia}-Enceladus, the Helmi streams, Sequoia, Sagittarius, and the recently discovered, `low-energy' GCs, which provide an excellent match to the predicted properties of the enigmatic galaxy `Kraken'. The five galaxies cover a narrow stellar mass range [$M_\star=(0.6{-}4.6)\times10^8~{\rm M}_\odot$], but have widely different accretion redshifts ($z_{\rm acc}=0.57{-}2.65$). All accretion events represent minor mergers, but Kraken likely represents the most major merger ever experienced by the Milky Way, with stellar and virial mass ratios of $r_{M_\star}=1$:$31^{+34}_{-16}$ and $r_{M_{\rm h}}=1$:$7^{+4}_{-2}$, respectively. The progenitors match the $z=0$ relation between GC number and halo virial mass, but have elevated specific frequencies, suggesting an evolution with redshift. Even though these progenitors likely were the Milky Way's most massive accretion events, they contributed a total mass of only $\log{(M_{\rm \star,tot}/{\rm M}_\odot)}=9.0\pm0.1$, similar to the stellar halo. This implies that the Milky Way grew its stellar mass mostly by in-situ star formation. We conclude by organising these accretion events into the most detailed reconstruction to date of the Milky Way's merger tree.

The ages and metallicities of globular clusters (GCs) are known to be powerful tracers of the properties of their progenitor galaxies, enabling their use in determining the merger histories of galaxies. However, while useful in separating GCs into individual accretion events, the orbits of GC groups themselves have received less attention as probes of their progenitor galaxy properties. In this work, we use simulations of galaxies and their GC systems from the E-MOSAICS project to explore how the present-day orbital properties of GCs are related to the properties of their progenitor galaxies. We find that the orbits of GCs deposited by accretion events are sensitive to the mass and merger redshift of the satellite galaxy. Earlier mergers and larger galaxy masses deposit GCs at smaller median apocentres and lower total orbital energy. The orbital properties of accreted groups of GCs can therefore be used to infer the properties of their progenitor galaxy, though there exists a degeneracy between galaxy mass and accretion time. Combining GC orbits with other tracers (GC ages, metallicities) will help to break the galaxy mass/accretion time degeneracy, enabling stronger constraints on the properties of their progenitor galaxy. In situ GCs generally orbit at lower energies (small apocentres) than accreted GCs, however they exhibit a large tail to high energies and even retrograde orbits (relative to the present-day disc), showing significant overlap with accreted GCs. Applying the results to Milky Way GCs groups suggests a merger redshift $z \sim 1.5$ for the Gaia Sausage/Enceladus and $z>2$ for the `low-energy'/Kraken group, adding further evidence that the Milky Way had two significant mergers in its past.

We present a simple physical model for the minimum mass of bound stellar clusters as a function of the galactic environment. The model evaluates which parts of a hierarchically-clustered star-forming region remain bound given the time-scales for gravitational collapse, star formation, and stellar feedback. We predict the initial cluster mass functions (ICMFs) for a variety of galaxies and we show that these predictions are consistent with observations of the solar neighbourhood and nearby galaxies, including the Large Magellanic Cloud and M31. In these galaxies, the low minimum cluster mass of $\sim10^2~\rm{M}_{\odot}$ is caused by sampling statistics, representing the lowest mass at which massive (feedback-generating) stars are expected to form. At the high gas density and shear found in the Milky Way's Central Molecular Zone and the nucleus of M82, the model predicts that a mass $>10^2~\rm{M}_{\odot}$ must collapse into a single cluster prior to feedback-driven dispersal, resulting in narrow ICMFs with elevated characteristic masses. We find that the minimum cluster mass is a sensitive probe of star formation physics due to its steep dependence on the star formation efficiency per free-fall time. Finally, we provide predictions for globular cluster (GC) populations, finding a narrow ICMF for dwarf galaxy progenitors at high redshift, which can explain the high specific frequency of GCs at low metallicities observed in Local Group dwarfs like Fornax and WLM. The predicted ICMFs in high-redshift galaxies constitute a critical test of the model, ideally-suited for the upcoming generation of telescopes.

The number density of field galaxies per rotation velocity, referred to as the velocity function, is an intriguing statistical measure probing the smallest scales of structure formation. In this paper we point out that the velocity function is sensitive to small shifts in key cosmological parameters such as the amplitude of primordial perturbations ($\sigma_8$) or the total matter density ($\Omega_{\rm m}$). Using current data and applying conservative assumptions about baryonic effects, we show that the observed velocity function of the Local Volume favours cosmologies in tension with the measurements from Planck but in agreement with the latest findings from weak lensing surveys. While the current systematics regarding the relation between observed and true rotation velocities are potentially important, upcoming data from HI surveys as well as new insights from hydrodynamical simulations will dramatically improve the situation in the near future.

The dearth of dwarf galaxies in the local universe is hard to reconcile with the large number of low mass haloes expected within the concordance $\Lambda$CDM paradigm. In this paper we perform a systematic evaluation of the uncertainties affecting the measurement of DM halo abundance using galaxy kinematics. Using a large sample of dwarf galaxies with spatially resolved kinematic data we derive a correction to obtain the observed abundance of galaxies as a function of their halo maximum circular velocity from the line-of-sight velocity function in the Local Volume. This estimate provides a direct means of comparing the predictions of theoretical models and simulations (including nonstandard cosmologies and novel galaxy formation physics) to the observational constraints. The new "galactic $V_{max}$" function is steeper than the line-of-sight velocity function but still shallower than the theoretical CDM expectation, showing that some unaccounted physical process is necessary to reduce the abundance of galaxies and/or drastically modify their density profiles compared to CDM haloes. Using this new galactic $V_{max}$ function, we investigate the viability of baryonic solutions such as feedback-powered outflows and photoevaporation of gas from an ionising radiation background. At the 3-$\sigma$ confidence level neither energetic feedback nor photoevaporation are effective enough to reconcile the disagreement. In the case of maximum baryonic effects, the theoretical estimate still deviates significantly from the observations for $V_{max} < 60$ km/s. CDM predicts at least 1.8 times more galaxies with $V_{max} = 50$ km/s and 2.5 times more than observed at $30$ km/s. Recent hydrodynamic simulations seem to resolve the discrepancy but disagree with the properties of observed galaxies with resolved kinematics. (abridged)

The observed number of dwarf galaxies as a function of rotation velocity is significantly smaller than predicted by the standard model of cosmology. This discrepancy cannot be simply solved by assuming strong baryonic feedback processes, since they would violate the observed relation between maximum circular velocity ($v_{\rm max}$) and baryon mass of galaxies. A speculative but tantalising possibility is that the mismatch between observation and theory points towards the existence of non-cold or non-collisionless dark matter (DM). In this paper, we investigate the effects of warm, mixed (i.e warm plus cold), and self-interacting DM scenarios on the abundance of dwarf galaxies and the relation between observed HI line-width and maximum circular velocity. Both effects have the potential to alleviate the apparent mismatch between the observed and theoretical abundance of galaxies as a function of $v_{\rm max}$. For the case of warm and mixed DM, we show that the discrepancy disappears, even for luke-warm models that evade stringent bounds from the Lyman-$\alpha$ forest. Self-interacting DM scenarios can also provide a solution as long as they lead to extended ($\gtrsim 1.5$ kpc) dark matter cores in the density profiles of dwarf galaxies. Only models with velocity-dependent cross sections can yield such cores without violating other observational constraints at larger scales.

We present a detailed analysis of a large-scale galactic outflow in the CGM of a massive (M_h ~ 10^12.5 Msun), star-forming (6.9 Msun/yr), sub-L* (0.5 L_B*) galaxy at z=0.39853 that exhibits a wealth of metal-line absorption in the spectra of the background quasar Q 0122-003 at an impact parameter of 163 kpc. The galaxy inclination angle (i=63 degree) and the azimuthal angle (Phi=73 degree) imply that the QSO sightline is passing through the projected minor-axis of the galaxy. The absorption system shows a multiphase, multicomponent structure with ultra-strong, wide velocity spread OVI (logN = 15.16\pm0.04, V_{90} = 419 km/s) and NV (logN = 14.69\pm0.07, V_{90} = 285 km/s) lines that are extremely rare in the literature. The highly ionized absorption components are well explained as arising in a low density (10^{-4.2} cm^{-3}), diffuse (10 kpc), cool (10^4 K) photoionized gas with a super-solar metallicity ([X/H] > 0.3). From the observed narrowness of the Lyb profile, the non-detection of SIV absorption, and the presence of strong CIV absorption we rule out equilibrium/non-equilibrium collisional ionization models. The low-ionization photoionized gas with a density of 10^{-2.5} cm^{-3} and a metallicity of [X/H] > -1.4 is possibly tracing recycled halo gas. We estimate an outflow mass of ~2x10^{10} Msun, a mass-flow rate of ~54 Msun/yr, a kinetic luminosity of ~9x10^{41} erg/s, and a mass loading factor of ~8 for the outflowing high-ionization gas. These are consistent with the properties of "down-the-barrel" outflows from infrared-luminous starbursts as studied by Rupke et al. Such powerful, large-scale, metal-rich outflows are the primary means of sufficient mechanical and chemical feedback as invoked in theoretical models of galaxy formation and evolution.

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

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

We studied a sample of 14 galaxies (0.1 < z < 0.7) using HST/WFPC2 imaging and high-resolution HST/COS or HST/STIS quasar spectroscopy of Lya, Lyb, OVI1031, and OVI1037 absorption. The galaxies, having 10.8 < log(M/M_solar) < 12.2, lie within D = 300 kpc of quasar sightlines, probing out to D/R_vir = 3. When the full range of galaxy virial masses and D/R_vir of the sample are examined, 40% of the HI absorbing clouds can be inferred to be escaping their host halo. The fraction of bound clouds decreases as D/R_vir increases such that the escaping fraction is around 15% for D/R_vir < 1, around 45% for 1 < D/R_vir < 2, and around 90% for 2 < D/R_vir < 3. Adopting the median mass log(M/M_solar) = 11.5 to divide the sample into "higher" and "lower" mass galaxies, we find mass dependency for the hot CGM kinematics. To our survey limits, OVI absorption is found in only 40% of the HI clouds in and around lower mass halos as compared to 85% around higher mass halos. For D/R < 1, lower mass halos have an escape fraction of 65%, whereas higher mass halos have an escape fraction of 5%. For 1 < D/R_vir < 2, the escape fractions are 55% and 35% for lower mass and higher mass halos, respectively. For 2 < D/R_vir < 3, the escape fraction for lower mass halos is around 90%. We show that it is highly likely that the absorbing clouds reside within 4 virial radii of their host galaxies and that the kinematics are dominated by outflows. Our finding of "differential kinematics" is consistent with the scenario of "differential wind recycling" proposed by Oppenheimer et al. We discuss the implications for galaxy evolution, the stellar to halo mass function, and the mass metallicity relationship of galaxies.

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

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

In Churchill et al., we used halo abundance matching applied to 182 galaxies in the MgII Absorber-Galaxy Catalog (MAGIICAT, Nielsen et al.) and showed that the mean MgII 2796 equivalent width follows a tight inverse-square power law, W_r(2796) ~ (D/R_vir)^-2, with projected location relative to the galaxy virial radius and that the MgII absorption covering fraction is invariant with galaxy virial mass, M_h, over the range 10.7 < M_h/M_solar < 13.9. In this work, we explore multivariate relationships between W_r(2796), virial mass, impact parameter, virial radius, and the theoretical cooling radius that further elucidate self-similarity in the cool/warm (T=10^{4-4.5} K) circumgalactic medium (CGM) with virial mass. We show that virial mass determines the extent and strength of the MgII absorbing gas such that the mean W_r(2796) increases with virial mass at fixed distance while decreasing with galactocentric distance for fixed virial mass. The majority of the absorbing gas resides within D ~ 0.3 R_vir, independent of both virial mass and minimum absorption threshold; inside this region, and perhaps also in the region 0.3 < D/R_vir = 1, the mean W_r(2796) is independent of virial mass. Contrary to absorber-galaxy cross-correlation studies, we show there is no anti-correlation between W_r(2796) and virial mass. We discuss how simulations and theory constrained by observations fully support self-similarity of the cool/warm CGM via the physics governing star formation, gas-phase metal enrichment, recycling efficiency of galactic scale winds, filament and merger accretion, and overdensity of local environment as a function of virial mass.

We apply halo abundance matching to obtain galaxy virial masses, M_h, and radii, Rvir, for the 183 isolated galaxies in the "MgII Absorber-Galaxy Catalog" (MAGIICAT, Nielsen et al.). All galaxies have spectroscopic redshifts (0.1 < z < 1.1) and their circumgalactic medium (CGM) is probed in MgII absorption within projected galactocentric distances D < 200 kpc. We examine the behavior of equivalent width, W(2796), and covering fraction, f_c, as a function of D, D/Rvir, and M_h. We find: [1] systematic segregation of M_h on the W(2796)-D plane (4.0 sigma); high-mass halos are found at higher D with larger W(2796) compared to lower mass halos. On the W(2796)-D/Rvir plane, mass segregation vanishes and we find W(2796) ~ (D/Rvir)^-2 (8.9 sigma); [2] higher mass halos have larger f_c at a given D, whereas f_c is independent of M_h at all D/Rvir; [3] f_c is constant with M_h over the range 10.7 < log(M_h/Msun) < 13.8 within a given D or D/Rvir. The combined results suggest that the MgII absorbing CGM is self-similar with halo mass, even above log(M_h/Msun)~12, where cold mode accretion is theoretically predicted to be quenched. If theory is correct, either outflows or sub-halos must contribute to absorption in high-mass halos such that low- and high-mass halos are observationally indistinguishable using MgII absorption strength once impact parameter is scaled by halo mass. Alternatively, the data may indicate that predictions of a universal shut down of cold-mode accretion in high-mass halos may require revision.

Using HST/COS/STIS and HIRES/Keck high-resolution spectra, we have studied a remarkable HI absorbing complex at z=0.672 toward the quasar Q1317+277. The HI absorption has a velocity spread of 1600 km/s, comprises 21 Voigt profile components, and resides at an impact parameter of D=58 kpc from a bright, high mass [log(M_vir/M_sun) ~ 13.7] elliptical galaxy that is deduced to have a 6 Gyr old, solar metallicity stellar population. Ionization models suggest the majority of the structure is cold gas surrounding a shock heated cloud that is kinematically adjacent to a multi-phase group of clouds with detected CIII, CIV and OVI absorption, suggestive of a conductive interface near the shock. The deduced metallicities are consistent with the moderate in situ enrichment relative to the levels observed in the z ~ 3 Ly-alpha forest. We interpret the HI complex as a metal-poor filamentary structure being shock heated as it accretes into the halo of the galaxy. The data support the scenario of an early formation period (z > 4) in which the galaxy was presumably fed by cold-mode gas accretion that was later quenched via virial shocking by the hot halo such that, by intermediate redshift, the cold filamentary accreting gas is continuing to be disrupted by shock heating. Thus, continued filamentary accretion is being mixed into the hot halo, indicating that the star formation of the galaxy will likely remain quenched. To date, the galaxy and the HI absorption complex provide some of the most compelling observational data supporting the theoretical picture in which accretion is virial shocked in the hot coronal halos of high mass galaxies.

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

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