Source author record

Sebastian Trujillo-Gomez

Sebastian Trujillo-Gomez appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2023arXiv

In-situ or accreted? Using deep learning to infer the origin of extragalactic globular clusters from observables

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.

preprint2022arXiv

A trail of dark matter-free galaxies from a bullet dwarf collision

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.

preprint2022arXiv

Constraining the formation of NGC1052-DF2 from its unusual globular cluster population

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, $Σ_{\rm ISM} \gtrsim 10^3 M_{\odot} \rm{pc}^{-2}$, and strong centrifugal support, $Ω\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.

preprint2022arXiv

Constraining the shape of dark matter haloes with globular clusters

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 ($σ\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.

preprint2022arXiv

Introducing EMP-Pathfinder: modelling the simultaneous formation and evolution of stellar clusters in their host galaxies

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 ($τ>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.

preprint2022arXiv

The emergence of dark matter-deficient ultra-diffuse galaxies driven by scatter in the stellar mass-halo mass relation and feedback from globular clusters

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.

preprint2021arXiv

Globular clusters as tracers of the dark matter halo: insights from the E-MOSAICS simulations

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.

preprint2021arXiv

NGC5846-UDG1: A galaxy formed mostly by star formation in massive, extremely dense clumps of gas

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.

preprint2021arXiv

The physics governing the upper truncation mass of the globular cluster mass function

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.

preprint2021arXiv

What to expect when using globular clusters as tracers of the total mass distribution in Milky Way-mass galaxies

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.

preprint2020arXiv

The globular cluster system mass-halo mass relation in the E-MOSAICS simulations

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.

preprint2020arXiv

Where did the globular clusters of the Milky Way form? Insights from the E-MOSAICS simulations

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.

preprint2015arXiv

An Extreme Metallicity, Large-Scale Outflow from a Star-Forming Galaxy at z ~ 0.4

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.

preprint2015arXiv

Direct Insights into Observational Absorption Line Analysis Methods of the Circumgalactic Medium Using Cosmological Simulations

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.

preprint2015arXiv

Low-mass galaxy assembly in simulations: regulation of early star formation by radiation from massive stars

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.

preprint2014arXiv

Halo Mass Dependence of HI and OVI Absorption: Evidence for Differential Kinematics

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.

preprint2014arXiv

Radiative feedback and the low efficiency of galaxy formation in low-mass haloes at high redshift

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.

preprint2013arXiv

Effects of baryon removal on the structure of dwarf spheroidal galaxies

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.

preprint2013arXiv

MAGIICAT III. Interpreting Self-Similarity of the Circumgalactic Medium with Virial Mass using MgII Absorption

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.

preprint2013arXiv

The Self-Similarity of the Circumgalactic Medium with Galaxy Virial Mass: Implications for Cold-Mode Accretion

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.

preprint2012arXiv

Quenched Cold Accretion of a Large Scale Metal-Poor Filament due to Virial Shocking in the Halo of a Massive z=0.7 Galaxy

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.

preprint2011arXiv

Galaxies in LCDM with Halo Abundance Matching: luminosity-velocity relation, baryonic mass-velocity relation, velocity function and clustering

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.