Search

Simulations are the best approximation to experimental laboratories in astrophysics and cosmology. However, the complexity, richness, and large size of their outputs severely limit the interpretability of their predictions. We describe a new, unbiased, and machine learning based approach to obtaining useful scientific insights from a broad range of simulations. The method can be used on today's largest simulations and will be essential to solve the extreme data exploration and analysis challenges posed by the Exascale era. Furthermore, this concept is so flexible, that it will also enable explorative access to observed data. Our concept is based on applying nonlinear dimensionality reduction to learn compact representations of the data in a low-dimensional space. The simulation data is projected onto this space for interactive inspection, visual interpretation, sample selection, and local analysis. We present a prototype using a rotational invariant hyperspherical variational convolutional autoencoder, utilizing a power distribution in the latent space, and trained on galaxies from IllustrisTNG simulation. Thereby, we obtain a natural Hubble tuning fork like similarity space that can be visualized interactively on the surface of a sphere by exploiting the power of HiPS tilings in Aladin Lite., Comment: 4 pages, 1 figure

Understanding what shapes the cold gas component of galaxies, which both provides the fuel for star formation and is strongly affected by the subsequent stellar feedback, is a crucial step towards a better understanding of galaxy evolution. Here, we analyse the HI properties of a sample of 46 Milky Way halo-mass galaxies, drawn from cosmological simulations (EMP-Pathfinder and FIREbox). This set of simulations comprises galaxies evolved self-consistently across cosmic time with different baryonic sub-grid physics: three different star formation models [constant star formation efficiency (SFE) with different star formation eligibility criteria, and an environmentally-dependent, turbulence-based SFE] and two different feedback prescriptions, where only one sub-sample includes early stellar feedback. We use these simulations to assess the impact of different baryonic physics on the HI content of galaxies. We find that the galaxy-wide HI properties agree with each other and with observations. However, differences appear for small-scale properties. The thin HI discs observed in the local Universe are only reproduced with a turbulence-dependent SFE and/or early stellar feedback. Furthermore, we find that the morphology of HI discs is particularly sensitive to the different physics models: galaxies simulated with a turbulence-based SFE have discs that are smoother and more rotationally symmetric, compared to those simulated with a constant SFE; galaxies simulated with early stellar feedback have more regular discs than supernova-feedback-only galaxies. We find that the rotational asymmetry of the HI discs depends most strongly on the underlying physics model, making this a promising observable for understanding the physics responsible for shaping the interstellar medium of galaxies., Comment: 13 pages, 10 figures + appendices (7 pages, 10 figures); updated to version accepted by MNRAS

We present an end-to-end description of the formation process of globular clusters (GCs) which combines a treatment for their formation and dynamical evolution within galaxy haloes with a state-of-the-art semi-analytic simulation of galaxy formation. Our approach allows us to obtain exquisite statistics to study the effect of the environment and assembly history of galaxies, while still allowing a very efficient exploration of the parameter space of star cluster physics. Our reference model, including both efficient cluster disruption during galaxy mergers and a model for the dynamical friction of GCs within the galactic potential, accurately reproduces the observed correlation between the total mass in GCs and the parent halo mass. A deviation from linearity is predicted at low halo masses, which is driven by a strong dependence on morphological type: bulge-dominated galaxies tend to host larger masses of GCs than their later-type counterparts of similar stellar mass. While the significance of the difference might be affected by resolution at the lowest halo masses considered, this is a robust prediction of our model and represents a natural consequence of the assumption that cluster migration from the disk to the halo is triggered by galaxy mergers. Our model requires an environmental dependence of GC radii to reproduce the observed mass distribution of GCs in our Galaxy at the low-mass end. At GC masses $>10^6\,{\rm M}_\odot$, our model predicts fewer GCs than observed due to an overly aggressive treatment of dynamical friction. The metallicity distribution measured for Galactic GCs is well reproduced by our model, even though it predicts systematically younger GCs than observed. We argue that this adds further evidence for an anomalously early formation of the stars in our Galaxy., Comment: 18 pages, 16 figures, accepted for publication in MNRAS

Globular clusters (GCs) are powerful tracers of the galaxy assembly process, and have already been used to obtain a detailed picture of the progenitors of the Milky Way. Using the E-MOSAICS cosmological simulation of a (34.4 Mpc)$^3$ volume that follows the formation and co-evolution of galaxies and their star cluster populations, we develop a method to link the origin of GCs to their observable properties. We capture this complex link using a supervised deep learning algorithm trained on the simulations, and predict the origin of individual GCs (whether they formed in the main progenitor or were accreted from satellites) based solely on extragalactic observables. An artificial neural network classifier trained on $\sim50,000$ GCs hosted by $\sim 700$ simulated galaxies successfully predicts the origin of GCs in the test set with a mean accuracy of $89$ per cent for the objects with [Fe/H]<-0.5 that have unambiguous classifications. The network relies mostly on the alpha-element abundances, metallicities, projected positions, and projected angular momenta of the clusters to predict their origin. A real-world test using the known progenitor associations of the Milky Way GCs achieves up to $90$ per cent accuracy, and successfully identifies as accreted most of the GCs in the inner Galaxy associated to the Kraken progenitor, as well as all the Gaia-Enceladus GCs. We demonstrate that the model is robust to observational uncertainties, and develop a method to predict the classification accuracy across observed galaxies. The classifier can be optimized for available observables (e.g. to improve the accuracy by including GC ages), making it a valuable tool to reconstruct the assembly histories of galaxies in upcoming wide-field surveys., Comment: 21 pages, 22 figures. Submitted to MNRAS. Comments welcome

The metallicity distributions of globular cluster (GC) systems in galaxies are a critical test of any GC formation scenario. In this work, we investigate the predicted GC metallicity distributions of galaxies in the MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE (E-MOSAICS) simulation of a representative cosmological volume ($L = 34.4$ comoving Mpc). We find that the predicted GC metallicity distributions and median metallicities from the fiducial E-MOSAICS GC formation model agree well the observed distributions, except for galaxies with masses $M_\ast \sim 2 \times 10^{10}$ M$_\odot$, which contain an overabundance of metal-rich GCs. The predicted fraction of galaxies with bimodal GC metallicity distributions ($37 \pm 2$ per cent in total; $45 \pm 7$ per cent for $M_\ast > 10^{10.5}$ M$_\odot$) is in good agreement with observed fractions ($44^{+10}_{-9}$ per cent), as are the mean metallicities of the metal-poor and metal-rich peaks. We show that, for massive galaxies ($M_\ast > 10^{10}$ M$_\odot$), bimodal GC distributions primarily occur as a result of cluster disruption from initially-unimodal distributions, rather than as a result of cluster formation processes. Based on the distribution of field stars with GC-like abundances in the Milky Way, we suggest that the bimodal GC metallicity distribution of Milky Way GCs also occurred as a result of cluster disruption, rather than formation processes. We conclude that separate formation processes are not required to explain metal-poor and metal-rich GCs, and that GCs can be considered as the surviving analogues of young massive star clusters that are readily observed to form in the local Universe today., Comment: 19 pages, 16 figures. Published in MNRAS

It was recently proposed that the dark matter-deficient ultra-diffuse galaxies DF2 and DF4 in the NGC1052 group could be the products of a "bullet dwarf" collision between two gas-rich progenitor galaxies. In this model DF2 and DF4 formed at the same time in the immediate aftermath of the collision, and a strong prediction is that their globular clusters should have nearly identical stellar populations. Here we test this prediction by measuring accurate F606W-F814W colors from deep HST/ACS imaging. We find that the clusters are extremely homogeneous. The mean color difference between the globular clusters in DF2 and DF4 is $-0.003\pm 0.005$ mag and the observed scatter for the combined sample of 18 clusters with $M_V<-8.6$ in both galaxies is $0.015 \pm 0.002$ mag. After accounting for observational uncertainties and stochastic cluster-to-cluster variation in the number of red giants, the remaining scatter is $0.008^{+0.005}_{-0.006}$ mag. Both the color difference and the scatter are an order of magnitude smaller than in other dwarf galaxies, and we infer that the bullet scenario passes an important test that could have falsified it. No other formation models have predicted this extreme uniformity of the globular clusters in the two galaxies. We find that the galaxies themselves are slightly redder than the clusters, consistent with a previously-measured metallicity difference. Numerical simulations have shown that such differences are expected in the bullet scenario, as the galaxies continued to self-enrich after the formation of the globular clusters., Comment: Published in ApJ Letters. It is difficult to visualize the uniformity of these strange clusters. Our best attempt is Fig. 2, where we compare them to globular clusters in Virgo galaxies

The ultra-diffuse galaxies DF2 and DF4 in the NGC1052 group share several unusual properties: they both have large sizes, rich populations of overluminous and large globular clusters, and very low velocity dispersions indicating little or no dark matter. It has been suggested that these galaxies were formed in the aftermath of high velocity encounters of gas rich galaxies, events that resemble the collision that created the bullet cluster but on much smaller scales. The gas separates from the dark matter in the collision and subsequent star formation leads to the formation of one or more dark matter-free galaxies. Here we show that the present-day line-of-sight distances and radial velocities of DF2 and DF4 are consistent with their joint formation in the aftermath of a single bullet-dwarf collision, around eight billion years ago. Moreover, we find that DF2 and DF4 are part of an apparent linear substructure of 7-11 large, low-luminosity objects. We propose that these all originated in the same event, forming a trail of dark matter-free galaxies that is more than 2 Mpc long and angled 7 +- 2 degrees from the line of sight. We also tentatively identify the highly dark matter-dominated remnants of the two progenitor galaxies that are expected at the leading edges of the trail., Comment: Nature, May 19 edition. See https://www.nature.com/articles/s41586-022-04665-6 . A pretty 40-orbit HST/ACS image of DF2 (used in Shen et al. 2021) can be found at https://www.pietervandokkum.com/publications

We explore how diffuse stellar light and globular clusters (GCs) can be used to trace the matter distribution of their host halo using an observational methodology. For this, we use 117 simulated dark matter (DM) haloes from the $(34.4~\rm cMpc)^3$ periodic volume of the E-MOSAICS project. For each halo, we compare the stellar surface brightness and GC projected number density maps to the surface densities of DM and total mass. We find that the dominant structures identified in the stellar light and in the GCs correspond closely with those from the DM and total mass. Our method is unaffected by the presence of satellites and its precision improves with fainter GC samples. We recover tight relations between the profiles of stellar surface brightness and GC number density to those of the DM, suggesting that the profile of DM can be accurately recovered from the stars and GCs ($\sigma\leq0.5~$dex). We quantify the projected morphology of DM, stars and GCs, and find that the stars and GCs are more flattened than the DM. Additionally, the semi-major axes of the distribution of stars and GCs are typically misaligned by $\sim 10~$degrees from that of DM. We demonstrate that deep imaging of diffuse stellar light and GCs can place constraints on the shape, profile and orientation of their host halo. These results extend down to haloes with central galaxies $M_{\star}\geq10^{10}~M_{\odot}$, and the analysis will be applicable to future data from the Euclid, Roman and the Rubin observatories., Comment: 17 pages, 15 figures; submitted to MNRAS (25/04/2022) - comments and suggestions welcome!

A key assumption in quasar absorption line studies of the circumgalactic medium (CGM) is that each absorption component maps to a spatially isolated "cloud" structure that has single valued properties (e.g. density, temperature, metallicity). We aim to assess and quantify the degree of accuracy underlying this assumption. We used adaptive mesh refinement hydrodynamic cosmological simulations of two $z=1$ dwarf galaxies and generated synthetic quasar absorption-line spectra of their CGM. For the SiII $\lambda 1260$ transition, and the CIV $\lambda\lambda1548, 1550$ and OVI $\lambda\lambda1031, 1037$ fine-structure doublets, we objectively determined which gas cells along a line-of-sight (LOS) contribute to detected absorption. We implemented a fast, efficient, and objective method to define individual absorption components in each absorption profile. For each absorption component, we quantified the spatial distribution of the absorbing gas. We studied a total of 1,302 absorption systems containing a total of 7,755 absorption components. 48% of SiII, 68% of CIV, and 72% of OVI absorption components arise from two or more spatially isolated "cloud" structures along the LOS. Spatially isolated "cloud" structures were most likely to have cloud-cloud LOS separations of 0.03$R_{vir}$, 0.11$R_{vir}$, and 0.13$R_{vir}$ for SiII, CIV, and OVI, respectively. There can be very little overlap between multi-phase gas structures giving rise to absorption components. If our results reflect the underlying reality of how absorption lines record CGM gas, they place tension on current observational analysis methods as they suggest that component-by-component absorption line formation is more complex than is assumed and applied for chemical-ionisation modelling., Comment: Submitted for publication in MNRAS on 02/24/2022

The formation and evolution of stellar clusters is intimately linked to that of their host galaxies. To study this connection, we present the EMP-Pathfinder suite of cosmological zoom-in Milky Way-mass simulations. These simulations contain a sub-grid description for stellar cluster formation and evolution, allowing us to study the simultaneous formation and evolution of stellar clusters alongside their host galaxies across cosmic time. As a key ingredient in these simulations, we include the physics of the multi-phase nature of the interstellar medium (ISM), which enables studies of how the presence of a cold, dense ISM affects cluster formation and evolution. We consider two different star formation prescriptions: a constant star formation efficiency per free-fall time, as well as an environmentally-dependent, turbulence-based prescription. We identify two key results drawn from these simulations. Firstly, we find that tidal shock-driven disruption caused by the graininess of the cold ISM produces old ($\tau>10~$Gyr) stellar cluster populations with properties that are in excellent agreement with the observed populations in the Milky Way and M31. Importantly, the addition of the cold ISM addresses the areas of disagreement found in previous simulations that lacked the cold gas phase. Secondly, the formation of stellar clusters is extremely sensitive to the baryonic physics that govern the properties of the cold, dense gas reservoir in the galaxy. This implies that the demographics of stellar cluster populations represent an important diagnostic tool for constraining baryonic physics models in upcoming galaxy formation simulations that also include a description of the cold ISM., Comment: 30 pages + appendices, 12 figures. MNRAS submitted, comments and suggestions welcome!

The mass function of globular cluster (GC) populations is a fundamental observable that encodes the physical conditions under which these massive stellar clusters formed and evolved. The high-mass end of star cluster mass functions are commonly described using a Schechter function, with an exponential truncation mass $M_{c,*}$. For the GC mass functions in the Virgo galaxy cluster, this truncation mass increases with galaxy mass ($M_{*}$). In this paper we fit Schechter mass functions to the GCs in the most massive galaxy group ($M_{\mathrm{200}} = 5.14 \times 10^{13} M_{\odot}$) in the E-MOSAICS simulations. The fiducial cluster formation model in E-MOSAICS reproduces the observed trend of $M_{c,*}$ with $M_{*}$ for the Virgo cluster. We therefore examine the origin of the relation by fitting $M_{c,*}$ as a function of galaxy mass, with and without accounting for mass loss by two-body relaxation, tidal shocks and/or dynamical friction. In the absence of these mass-loss mechanisms, the $M_{c,*}$-$M_{*}$ relation is flat above $M_* > 10^{10} M_{\odot}$. It is therefore the disruption of high-mass GCs in galaxies with $M_{*}\sim 10^{10} M_{\odot}$ that lowers the $M_{c,*}$ in these galaxies. High-mass GCs are able to survive in more massive galaxies, since there are more mergers to facilitate their redistribution to less-dense environments. The $M_{c,*}-M_*$ relation is therefore a consequence of both the formation conditions of massive star clusters and their environmentally-dependent disruption mechanisms., Comment: Accepted to MNRAS

It has been shown that ultra-diffuse galaxies (UDGs) have higher specific frequencies of globular clusters on average than other dwarf galaxies with similar luminosities. The UDG NGC5846-UDG1 is among the most extreme examples of globular cluster-rich galaxies found so far. Here we present new Hubble Space Telescope (HST) observations and analysis of this galaxy and its globular cluster system. We find that NGC5846-UDG1 hosts $54 \pm 9$ globular clusters, three to four times more than any previously known galaxy with a similar luminosity, and higher than reported in previous studies. With a galaxy luminosity of $L_{V,\mathrm{gal}} \approx 6 \times 10^7\,{\rm L}_{\odot}$ ($M_\star \approx 1.2 \times 10^8\,{\rm M}_\odot$) and a total globular cluster luminosity of $L_{V,\mathrm{GCs}} \approx 7.6 \times 10^6\,{\rm L}_{\odot}$, we find that the clusters currently comprise $\sim 13 \%$ of the total light. Taking into account the effects of mass loss from clusters during their formation and throughout their lifetime, we infer that most of the stars in the galaxy likely formed in globular clusters, and very little to no "normal" low-density star formation occurred. This result implies that the most extreme conditions during early galaxy formation promoted star formation in massive and dense clumps, in contrast to the dispersed star formation observed in galaxies today., Comment: Submitted to ApJL. Figure 1 shows the Hubble Space Telescope image of the galaxy, and Figure 5 shows the modeling results. Comments are very welcome!

We used adaptive mesh refinement hydrodynamic cosmological simulations of a $z=1$ Milky Way-type galaxy and a $z=0$ Dwarf galaxy and generated synthetic quasar absorption-line spectra of their circumgalactic medium (CGM). Our goal is to assess whether standard observational spectroscopic analysis methods accurately reproduce intrinsic column densities, metallicities [Si/H], and hydrogen densities $n_{H}$, in simulated absorption-line systems. Without knowledge of the intrinsic simulated properties (blind study), we analysed synthetic COS and HIRES spectra with fixed $S/N=30$ to determine the column densities, metallicity, and $n_{H}$, using Voigt profile fitting combined with Markov-Chain Monte-Carlo single-phase {\sc Cloudy} modelling techniques. To quantify the intrinsic simulated absorbing gas properties, we objectively determined which gas cells along a line of sight (LOS) contribute to detected absorption in the spectra and adopt the unweighted geometric mean of these gas cell properties. For this pilot study, we performed this experiment for five LOS in the Milky-Way galaxy and five LOS in the Dwarf galaxy. We found an average agreement between the "observed" and intrinsic metallicity overestimated within $0.8\sigma$ or $0.2$ dex for the "Milky-Way"' CGM and overestimated within $1.4\sigma$ or $0.2$ dex for the Dwarf galaxy CGM. We found that the spectroscopically-derived $n_{H}$ are underestimated within $0.8\sigma$ or $0.4$ dex of the intrinsic $n_{H}$ for the "Milky-Way" CGM and overestimated within $0.3\sigma$ or $0.3$ dex for the Dwarf galaxy CGM. The overall agreement suggests that, for single-phase ionisation modelling of systems where there is substantial spread in gas properties, global metallicity measurements from quasar absorption line studies are capturing the average metallicity and ionisation parameters in a given astrophysical environment., Comment: Submitted to MNRAS on August 6th, 2021. Accepted for publication on October 4th, 2021

Globular clusters (GCs) are bright objects that span a wide range of galactocentric distances, and are thus probes of the structure of dark matter (DM) haloes. In this work, we explore whether the projected radial profiles of GCs can be used to infer the structural properties of their host DM haloes. We use the simulated GC populations in a sample of 166 central galaxies from the $(34.4~\rm cMpc)^3$ periodic volume of the E-MOSAICS project. We find that more massive galaxies host stellar and GC populations with shallower density profiles that are more radially extended. In addition, the metal-poor GC subpopulations tend to have shallower and more extended profiles than the metal-rich subsamples, which we relate to the preferentially accreted origin of the metal-poor GCs. We find strong correlations between the slopes and effective radii of the radial profiles of the GC populations and the structural properties of the DM haloes, such as their power-law slopes, scale radii, and concentration parameters. Accounting for a dependence on the galaxy stellar mass decreases the scatter of the two-dimensional relations. This suggests that the projected number counts of GCs, combined with their galaxy mass, trace the density profile of the DM halo of their host galaxy. When applied to extragalactic GC systems, we recover the scale radii and the extent of the DM haloes of a sample of ETGs with uncertainties smaller than $0.2~\rm dex$. Thus, extragalactic GC systems provide a novel avenue to explore the structure of DM haloes beyond the Local Group., Comment: 19 pages, 11 figures and 2 tables; submitted to MNRAS, comments and/or suggestions are welcomed!

In addition to their low stellar densities, ultra-diffuse galaxies (UDGs) have a broad variety of dynamical mass-to-light ratios, ranging from dark matter (DM) dominated systems to objects nearly devoid of DM. To investigate the origin of this diversity, we develop a simple, semi-empirical model that predicts the structural evolution of galaxies, driven by feedback from massive star clusters, as a function of their departure from the mean SMHM relation. The model predicts that a galaxy located $\gtrsim 0.5$ dex above the mean relation at $M_{\rm halo}=10^{10}M_{\odot}$ will host a factor of $\sim 10-100$ larger globular cluster (GC) populations, and that feedback from these GCs drives a significant expansion of the stellar component and loss of DM compared to galaxies on the SMHM relation. This effect is stronger in haloes that collapse earlier and have enhanced star formation rates at $z\gtrsim2$, which leads to increased gas pressures, stellar clustering, and mean cluster masses, and significantly enhances the energy loading of galactic winds and its impact on the DM and stellar orbits. The impact on galaxy size and DM content can be large enough to explain observed galaxies that contain nearly the universal baryon fraction, as well as NGC1052-DF2 and DF4 and other isolated UDGs that contain almost no DM. The trend of increasing galaxy size with GC specific frequency observed in galaxy clusters also emerges naturally in the model. Our predictions can be tested with large and deep surveys of the stellar and GC populations in dwarfs and UDGs. Because stellar clustering drives the efficiency of galactic winds, it may be a dominant factor in the structural evolution of galaxies and should be included as an essential ingredient in galaxy formation models., Comment: 23 pages, 17 figures. Accepted for publication in MNRAS (November 5)

Dynamical models allow us to connect the motion of a set of tracers to the underlying gravitational potential, and thus to the total (luminous and dark) matter distribution. They are particularly useful for understanding the mass and spatial distribution of dark matter (DM) in a galaxy. Globular clusters (GCs) are an ideal tracer population in dynamical models, since they are bright and can be found far out into the halo of galaxies. We aim to test how well Jeans-Anisotropic-MGE (JAM) models using GCs (positions and line-of-sight velocities) as tracers can constrain the mass and radial distribution of DM halos. For this, we use the E-MOSAICS suite of 25 zoom-in simulations of L* galaxies. We find that the DM halo properties are reasonably well recovered by the JAM models. There is, however, a strong correlation between how well we recover the mass and the radial distribution of the DM and the number of GCs in the galaxy: the constraints get exponentially worse with fewer GCs, and at least 150 GCs are needed in order to guarantee that the JAM model will perform well. We find that while the data quality (uncertainty on the radial velocities) can be important, the number of GCs is the dominant factor in terms of the accuracy and precision of the measurements. This work shows promising results for these models to be used in extragalactic systems with a sample of more than 150 GCs., Comment: 18 pages, 13 figures. Accepted for publication in MNRAS

The ultra-diffuse galaxy (UDG) NGC1052-DF2 has a low dark matter content and hosts a very unusual globular cluster (GC) population, with a median luminosity $\sim4$ times higher than in most galaxies and containing about 5~per~cent of the galaxy's stars. We apply a theoretical model that predicts the initial cluster mass function as a function of the galactic environment to investigate the origin of DF2's peculiar GC system. Using the GC mass function, the model constrains the star-forming conditions in the galaxy during the formation of its GCs, $\sim9~\rm{Gyr}$ ago. We predict that the GCs formed in an environment with very high gas surface density, $\Sigma_{\rm ISM} \gtrsim 10^3 M_{\odot} \rm{pc}^{-2}$, and strong centrifugal support, $\Omega \gtrsim 0.7~\rm{Myr}^{-1}$, similar to nearby circum-nuclear starbursts and the central region of the Milky Way. These extreme conditions required to form the observed GC population imply a very high cluster formation efficiency of $\gtrsim 78$ per cent, and contrast strongly with the current diffuse nature of the galaxy. Since a nuclear starburst would lead to the rapid in-spiral of the GCs and is ruled out by the absence of a nuclear star cluster, we propose that the GCs plausibly formed during a major merger at $z\sim1.3$. The merger remnant must have undergone significant expansion of its stellar (and perhaps also its dark matter) component to reach its low present surface brightness, leading to the interesting possibility that it was the formation of DF2's extreme GC population that caused it to become a UDG. If true, this strong structural evolution would have important implications for understanding the origins of UDGs., Comment: 14 pages, 6 figures; Accepted for publication in MNRAS (June 17, 2021)

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 (Mgc) and the dark matter halo mass (Mhalo) 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.4$~cMpc 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 Mgc-Mhalo relation is linear for halo masses $>5\times10^{11}~Msun$, and is driven by the hierarchical assembly of galaxies. 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-halo mass relation of galaxies. Our fiducial model reproduces the observed Mgc-Mstar relation across the full mass range, which we argue is more physically relevant than the Mgc-Mhalo relation. We also explore the physical processes driving the observed constant value of $Mgc / Mhalo \sim 5\times10^{-5}$ and find that it is the result of a combination of cluster formation physics and cluster disruption., Comment: 13 pages, 9 figures, accepted for publication in MNRAS

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., Comment: 20 pages, 20 figures, resubmitted to MNRAS after accounting for referee's comments

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., Comment: 24 pages, 20 figures. Published in MNRAS

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., Comment: 20 pages, 9 figures, 2 tables; submitted to MNRAS (March 1, 2020); Figures 4, 5, and 9 show the main results of the paper

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., Comment: 13 pages, 7 figures. Accepted for publication in MNRAS (21 September 2020)

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., Comment: 24 pages, 14 figures. Accepted to MNRAS

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., Comment: 17 pages, 9 figures, MNRAS published

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., Comment: Matches published version

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), Comment: 17 pages, 22 figures; major revisions include clarification of the method, expanded comparison with simulations with a new figure, analysis of uncertainties in model as well as pressure support corrections, and a new table with nomenclature

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., Comment: APJ accepted; Updated Fig. 2 with more galaxy info

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., Comment: 18 pages, 7 figures, revised version accepted to ApJ Jan 2015

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., Comment: 23 pages, 19 figures, ApJ Accepted

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.

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., Comment: 25 pages plus full data table, 11 figures; Revised version accepted to The Astrophysical Journal (Oct 27, 2013)

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., Comment: 18 pages, 11 figures, accepted version, new figure added, results and conclusions unchanged

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., Comment: 17 pages, 7 figures, and 3 tables; accepted for publication in MNRAS

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., Comment: revised version includes additional galaxies and is the published version in Astrophysical Letters

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., Comment: 20 pages, 9 figures, submitted to ApJ

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., Comment: 40 pages, 18 figures, published in ApJ