Source author record

Jeremy Bailin

Jeremy Bailin 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

27works
4topics
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

27 published item(s)

preprint2026arXiv

The origin of strong $α$-element bimodalities in FIRE simulations of Milky Way-mass galaxies

One of the Milky Way's characteristic features is a strongly bimodal distribution of $α$-process elements, such as Mg, at fixed [Fe/H] in stellar abundances. We examine patterns in [Mg/Fe] versus [Fe/H] in FIRE-2 simulations of Milky Way-mass galaxies. Out of 16 galaxies, 4 are capable of producing a strongly bimodal distribution. In all four galaxies, the high-$α$ population corresponds to an older, radially-compact, thick disk, and the low-$α$ population corresponds to a younger, radially-extended, thin disk, similar to the MW.The transition from high- to low-$α$ took $0.3-1.2\Gyr$ and began $5.5-6.5\Gyr$ ago. [Mg/Fe] decreased at relatively fixed [Fe/H], both in the galaxy overall and at fixed radii: Fe enrichment nearly balanced gas accretion (and therefore dilution), but Mg enrichment was weaker. Importantly, this transition occurred during a period of relatively low gas fraction ($5-15\%$), immediately after a rapid decline in star formation (halving within a few hundred Myr), which caused an increase in Fe-producing white-dwarf supernovae relative to Mg-producing core-collapse supernovae. Only one case coincided with a major merger coalescence. We find similar trends in measuring stars by their current radius and by their birth radius, therefore, radial redistribution did not play a dominant role in the formation of a bimodality or its spatial dependence today. Overall, in FIRE-2, strong $α$-element bimodalities are relatively uncommon ($\sim25\%$), often not associated with a major merger, and arise primarily from a rapid decline in star formation during relatively low gas fraction.

preprint2022arXiv

3D elemental abundances of stars at formation across the histories of Milky Way-mass galaxies in the FIRE simulations

We characterize the 3-D spatial variations of [Fe/H], [Mg/H], and [Mg/Fe] in stars at the time of their formation, across 11 simulated Milky Way (MW)- and M31-mass galaxies in the FIRE-2 simulations, to inform initial conditions for chemical tagging. The overall scatter in [Fe/H] within a galaxy decreased with time until $\approx 7$ Gyr ago, after which it increased to today: this arises from a competition between a reduction of azimuthal scatter and a steepening of the radial gradient in abundance over time. The radial gradient is generally negative, and it steepened over time from an initially flat gradient $\gtrsim 12$ Gyr ago. The strength of the present-day abundance gradient does not correlate with when the disk `settled'; instead, it best correlates with the radial velocity dispersion within the galaxy. The strength of azimuthal variation is nearly independent of radius, and the 360 degree scatter decreased over time, from $\lesssim 0.17$ dex at $t_{\rm lb} = 11.6$ Gyr to $\sim 0.04$ dex at present day. Consequently, stars at $t_{\rm lb} \gtrsim 8$ Gyr formed in a disk with primarily azimuthal scatter in abundances. All stars formed in a vertically homogeneous disk, $Δ$[Fe/H] $\leq 0.02$ dex within $1$ kpc of the galactic midplane, with the exception of the young stars in the inner $\approx 4$ kpc at $z \sim 0$. These results generally agree with our previous analysis of gas-phase elemental abundances, which reinforces the importance of cosmological disk evolution and azimuthal scatter in the context of stellar chemical tagging. We provide analytic fits to our results for use in chemical-tagging analyses.

preprint2022arXiv

Extinguishing the FIRE: environmental quenching of satellite galaxies around Milky Way-mass hosts in simulations

The star formation and gas content of satellite galaxies around the Milky Way (MW) and Andromeda (M31) are depleted relative to more isolated galaxies in the Local Group (LG) at fixed stellar mass. We explore the environmental regulation of gas content and quenching of star formation in $z=0$ galaxies at $M*=10^{5-10}\rm{M}_{\odot}$ around 14 MW-mass hosts from the FIRE-2 simulations. Lower-mass satellites ($M*\lesssim10^7\rm{M}_{\odot}$) are mostly quiescent and higher-mass satellites ($M*\gtrsim10^8\rm{M}_{\odot}$) are mostly star-forming, with intermediate-mass satellites ($M*\approx10^{7-8}\rm{M}_{\odot}$) split roughly equally between quiescent and star-forming. Hosts with more gas in their circumgalactic medium have a higher quiescent fraction of massive satellites ($M*=10^{8-9}\rm{M}_{\odot}$). We find no significant dependence on isolated versus paired (LG-like) host environments, and the quiescent fractions of satellites around MW-mass and LMC-mass hosts from the FIRE-2 simulations are remarkably similar. Environmental effects that lead to quenching can also occur as preprocessing in low-mass groups prior to MW infall. Lower-mass satellites typically quenched before MW infall as central galaxies or rapidly during infall into a low-mass group or a MW-mass galaxy. Most intermediate- to high-mass quiescent satellites have experienced $\geq1-2$ pericentre passages ($\approx2.5-5$ Gyr) within a MW-mass halo. Most galaxies with $M*\gtrsim10^{6.5}\rm{M}_{\odot}$ did not quench before falling into a host, indicating a possible upper mass limit for isolated quenching. The simulations reproduce the average trend in the LG quiescent fraction across the full range of satellite stellar masses. Though the simulations are consistent with the SAGA survey's quiescent fraction at $M*\gtrsim10^8\rm{M}_{\odot}$, they do not generally reproduce SAGA's turnover at lower masses.

preprint2022arXiv

New Globular Cluster Candidates in the M81 group

The study of outer halo globular cluster (GC) populations can give insight into galaxy merging, globular cluster accretion and the origin of GCs. We use archival Subaru Hyper Suprime-Cam (HSC) data in concert with space-based GALEX, IRAC and Gaia EDR3 data to select candidate Globular clusters (GCs) in the outer halo of the M81 group for confirmation and future study. We use a small sample of previously-discovered GCs to tune our selection criteria, finding that bright already-known GCs in the M81 group have sizes that are typically slightly larger than the Subaru PSF in our fields. In the optical bands, GCs appear to have colours that are only slightly different from stars. The inclusion of archival IRAC data yields dramatic improvements in colour separation, as the long wavelength baseline aids somewhat in the separation from stars and clearly separates GCs from many compact background galaxies. We show that some previously spectroscopically-identified GCs in the M81 group are instead foreground stars or background galaxies. GCs close to M82 have radial velocities suggesting that they fell into the M81 group along with M82. The overall M81 GC luminosity function is similar to the Milky Way and M31. M81's outer halo GCs are similar to the Milky Way in their metallicities and numbers, and much less numerous than M31's more metal-rich outer halo GC population. These properties reflect differences in the three galaxies' merger histories, highlighting the possibility of using outer halo GCs to trace merger history in larger samples of galaxies.

preprint2022arXiv

Predictions for Complex Distributions of Stellar Elemental Abundances in Low-Mass Galaxies

We investigate stellar elemental abundance patterns at z = 0 in 8 low-mass (M_* = 10^6 - 10^9 M_sun) galaxies in the Feedback in Realistic Environments (FIRE-2) cosmological simulations. Using magnesium (Mg) as a representative alpha-element, we explore stellar abundance patterns in [Mg/Fe] versus [Fe/H], which follow an overall monotonic trend that evolved slowly over time. Beyond this, we explore 3 notable secondary features in enrichment (found in three different case-study galaxies) that arise from a galaxy merger or bursty star formation. First, we observe a secondary track with a lower [Mg/Fe] than the main trend. At z = 0, stars from this track are predominantly found within 2-6 kpc of the center; they were accreted in a 1:3 total-mass-ratio merger ~ 0.4 Gyr ago. Second, we find a distinct elemental bi-modality that forms following a strong burst in star formation in a galaxy at t_lookback ~ 10 Gyr. This burst quenched star formation for ~ 0.66 Gyr, allowing Ia supernovae to enrich the system with iron before star formation resumed. Third, we examine stripes in enrichment that run roughly orthogonal to the dominant [Mg/Fe] versus [Fe/H] trend; these stripes correspond to short bursts of star formation during which core-collapse supernovae enrich the surrounding medium with Mg (and Fe) on short timescales. If observed, these features would substantiate the utility of elemental abundances in revealing the assembly and star formation histories of dwarf galaxies. We explore the observability of these features for upcoming spectroscopic studies. Our results show that precise measurements of elemental abundance patterns can reveal critical events in the formation histories of low-mass galaxies.

preprint2021arXiv

Globular Cluster Intrinsic Iron Abundance Spreads: II. Protocluster Metallicities and the Age-Metallicity Relations of Milky Way Progenitors

Intrinsic iron abundance spreads in globular clusters, although usually small, are very common, and are signatures of self enrichment: some stars within the cluster have been enriched by supernova ejecta from other stars within the same cluster. We use the Bailin (2018) self enrichment model to predict the relationship between properties of the protocluster -- its mass and the metallicity of the protocluster gas cloud -- and the final observable properties today -- its current metallicity and the internal iron abundance spread. We apply this model to an updated catalog of Milky Way globular clusters where the initial mass and/or the iron abundance spread is known to reconstruct their initial metallicities. We find that with the exception of the known anomalous bulge cluster Terzan 5 and three clusters strongly suspected to be nuclear star clusters from stripped dwarf galaxies, the model provides a good lens for understanding their iron spreads and initial metallicities. We then use these initial metallicities to construct age-metallicity relations for kinematically-identified major accretion events in the Milky Way's history. We find that using the initial metallicity instead of the current metallicity does not alter the overall picture of the Milky Way's history, since the difference is usually small, but does provide information that can help distinguish which accretion event some individual globular clusters with ambiguous kinematics should be associated with, and points to potential complexity within the accretion events themselves.

preprint2020arXiv

A profile in FIRE: resolving the radial distributions of satellite galaxies in the Local Group with simulations

While many tensions between Local Group (LG) satellite galaxies and LCDM cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with Mstar > 10^5 Msun around 8 isolated Milky Way- (MW) mass host galaxies and 4 hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with Mstar >~ 10^5 Msun. The simulations broadly agree with LG observations, spanning the radial profiles around the MW and M31. This agreement does not depend strongly on satellite mass, even at distances <~ 100 kpc. Host-to-host variation dominates the scatter in satellite counts within 300 kpc of the hosts, while time variation dominates scatter within 50 kpc. More massive host galaxies within our sample have fewer satellites at small distances, likely because of enhanced tidal destruction of satellites via the baryonic disks of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhalos in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at Mstar >~ 10^5 Msun: we predict 2-10 such satellites to be discovered around the MW and possibly 6-9 around M31. To provide cosmological context, we compare our results with the radial profiles of satellites around MW analogs in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems.

preprint2020arXiv

Is NGC 300 a pure exponential disk galaxy?

NGC 300 is a low-mass disk galaxy in the Sculptor group. In the literature, it has been identified as a pure exponential disk galaxy, as its luminosity profile could be well fitted with a single exponential law over many disk scale lengths (Type I). We investigate the stellar luminosity distribution of NGC 300 using $Hubble$ $Space$ $Telescope$ (HST) archive data, reaching farther and deeper than any other previous studies. Color magnitude diagrams show a significant population of old red giant branch (RGB) stars in all fields out to $R\sim19$ kpc ($32'$), as well as younger populations in the inner regions. We construct the density profiles of the young, intermediate-aged, and old stellar populations. We find two clear breaks in the density profiles of the old RGB and intermediate-aged stars: one down-bending (Type II) at $R\sim5.9$ kpc, and another up-bending (Type III) at $R\sim8.3$ kpc. Moreover, the old RGB stars exhibit a negative radial color gradient with an up-bending at $R\sim8$~kpc, beyond which the stellar populations are uniformly old ($>$7~Gyr) and metal-poor ($\rm[Fe/H] = -1.6^{+0.2}_{-0.4}$ dex). The outer stellar component at $R\gtrapprox8$ kpc is, therefore, well separated from the inner disk in terms of the stellar density and stellar populations. While our results cast doubt on the currently established wisdom that NGC\,300 is a pure exponential disk galaxy, a more detailed survey should be carried out to identify the outskirts as either a disk or a stellar halo.

preprint2020arXiv

The Formation Times and Building Blocks of Milky Way-mass Galaxies in the FIRE Simulations

Surveys of the Milky Way (MW) and M31 enable detailed studies of stellar populations across ages and metallicities, with the goal of reconstructing formation histories across cosmic time. These surveys motivate key questions for galactic archaeology in a cosmological context: when did the main progenitor of a MW/M31-mass galaxy form, and what were the galactic building blocks that formed it? We investigate the formation times and progenitor galaxies of MW/M31-mass galaxies using the FIRE-2 cosmological simulations, including 6 isolated MW/M31-mass galaxies and 6 galaxies in Local Group (LG)-like pairs at z = 0. We examine main progenitor "formation" based on two metrics: (1) transition from primarily ex-situ to in-situ stellar mass growth and (2) mass dominance compared to other progenitors. We find that the main progenitor of a MW/M31-mass galaxy emerged typically at z ~ 3-4 (11.6-12.2 Gyr ago), while stars in the bulge region (inner 2 kpc) at z = 0 formed primarily in a single main progenitor at z < 5 (< 12.6 Gyr ago). Compared with isolated hosts, the main progenitors of LG-like paired hosts emerged significantly earlier (Δz ~ 2, Δt ~ 1.6 Gyr), with ~ 4x higher stellar mass at all z > 4 (> 12.2 Gyr ago). This highlights the importance of environment in MW/M31-mass galaxy formation, especially at early times. Overall, about 100 galaxies with M_star > 10^5 M_sun formed a typical MW/M31-mass system. Thus, surviving satellites represent a highly incomplete census (by ~ 5x) of the progenitor population.

preprint2020arXiv

Tracing the anemic stellar halo of M101

Models of galaxy formation in a cosmological context predict that massive disk galaxies should have structured extended stellar halos. Recent studies in integrated light, however, report a few galaxies, including the nearby disk galaxy M101, that have no measurable stellar halos to the detection limit. We aim to quantify the stellar content and structure of M101's outskirts by resolving its stars. We present the photometry of its stars based on deep F606W and F814W images taken with Hubble Space Telescope as part of the GHOSTS survey. The constructed CMDs of stars reach down to two magnitudes below the tip of the red giant branch. We derived radial number density profiles of the bright red giant branch (RGB) stars. The mean color of the RGB stars at $R \sim$ 40 -- 60 kpc is similar to those of metal-poor globular clusters in the Milky Way. We also derived radial surface brightness profiles using the public image data provided by the Dragonfly team. Both the radial number density and surface brightness profiles were converted to radial mass density profiles and combined. We find that the mass density profiles show a weak upturn at the very outer region, where surface brightness is as faint as $μ_g\approx 34$ mag arcsec$^{-1}$. An exponential disk + power-law halo model on the mass density profiles finds the total stellar halo mass of $M_{halo}=8.2_{-2.2}^{+3.5}\times 10^7M_\odot$. The total stellar halo mass does not exceed $M_{halo} = 3.2 \times 10^8$ $M_{\odot}$ when strongly truncated disk models are considered. Combining the halo mass with the total stellar mass of M101, we obtain the stellar halo mass fraction of $M_{halo}/M_{gal} = 0.20_{-0.08}^{+0.10}\%$ with an upper limit of 0.78\%. We compare the halo properties of M101 with those of six GHOSTS survey galaxies as well as the Milky Way and M31 and find that M101 has an anemic stellar halo.

preprint2019arXiv

Globular Cluster Intrinsic Iron Abundance Spreads: I. Catalog

We present an up-to-date catalog of intrinsic iron abundance spreads in the 55 Milky Way globular clusters for which sufficiently precise spectroscopic measurements are available. Our method combines multiple datasets when possible to improve the statistics, taking into account the fact that different methods and instruments can lead to systematically offset metallicities. Only high spectral resolution (R>14,000) studies that measure the equivalent widths of individual iron lines are found to have uncertainties on the metallicities of the individual stars that can be calibrated sufficiently well for the intrinsic dispersion to be separated cleanly from random measurement error. The median intrinsic iron spread is found to be 0.045 dex, which is small but unambiguously measured to be non-zero in most cases. There is large variation between clusters, but more luminous globular clusters, above 10^5 L_sun, have increasingly large iron spreads on average; no trend between iron spread and metallicity is found.

preprint2016arXiv

The GHOSTS survey. II. The diversity of Halo Color and Metallicity Profiles of Massive Disk Galaxies

We study the stellar halo color properties of six nearby massive highly inclined disk galaxies using Hubble Space Telescope Advanced Camera for Surveys and Wide Field Camera 3 observations in both F606W and F814W filters from the GHOSTS survey. The observed fields, placed both along the minor and major axis of each galaxy, probe the stellar outskirts out to projected distances of ~ 50-70 kpc from their galactic centre along the minor axis. The 50% completeness levels of the color magnitude diagrams are typically at two mag below the tip of the red giant branch. We find that all galaxies have extended stellar halos out to ~ 50 kpc and two out to ~ 70 kpc. We determined the halo color distribution and color profile for each galaxy using the median colors of stars in the RGB. Within each galaxy we find variations in the median colors as a function of radius which likely indicates population variations, reflecting that their outskirts were built from several small accreted objects. We find that half of the galaxies (NGC 0891, NGC 4565, and NGC 7814) present a clear negative color gradient, reflecting a declining metallicity in their halos; the other have no significant color or population gradient. In addition, notwithstanding the modest sample size of galaxies, there is no strong correlation between their halo color/metallicity or gradient with galaxy's properties such as rotational velocity or stellar mass. The diversity in halo color profiles observed in the GHOSTS galaxies qualitatively supports the predicted galaxy-to-galaxy scatter in halo stellar properties; a consequence of the stochasticity inherent in the assembling history of galaxies.

preprint2016arXiv

The impact of baryons on the direct detection of dark matter

The spatial and velocity distributions of dark matter particles in the Milky Way Halo affect the signals expected to be observed in searches for dark matter. Results from direct detection experiments are often analyzed assuming a simple isothermal distribution of dark matter, the Standard Halo Model (SHM). Yet there has been skepticism regarding the validity of this simple model due to the complicated gravitational collapse and merger history of actual galaxies. In this paper we compare the SHM to the results of cosmological hydrodynamical simulations of galaxy formation to investigate whether or not the SHM is a good representation of the true WIMP distribution in the analysis of direct detection data. We examine two Milky Way-like galaxies from the MaGICC cosmological simulations (a) with dark matter only and (b) with baryonic physics included. The inclusion of baryons drives the shape of the DM halo to become more spherical and makes the velocity distribution of dark matter particles less anisotropic especially at large heliocentric velocities, thereby making the SHM a better fit. We also note that we do not find a significant disk-like rotating dark matter component in either of the two galaxy halos with baryons that we examine, suggesting that dark disks are not a generic prediction of cosmological hydrodynamical simulations. We conclude that in the Solar neighborhood, the SHM is in fact a good approximation to the true dark matter distribution in these cosmological simulations (with baryons) which are reasonable representations of the Milky Way, and hence can also be used for the purpose of dark matter direct detection calculations.

preprint2015arXiv

Origin of the Metallicity Distribution in the Thick Disc

Aims. Using a suite of cosmological chemodynamical disc galaxy simulations, we assess how (a) radial metallicity gradients evolve with scaleheight; (b) the vertical metallicity gradients change through the thick disc; and (c) the vertical gradient of the stellar rotation velocity varies through the disc. We compare with the Milky Way to search for analogous trends. Methods. We analyse five simulated spiral galaxies with masses comparable to the Milky Way. The simulations span a range of star formation and energy feedback strengths and prescriptions, particle- and grid-based hydrodynamical implementations, as well as initial conditions/assembly history. Results. Consistently, we find that the steeper, negative, radial metallicity gradients seen in the mid-plane flatten with increasing height away from the plane. In simulations with stronger (and/or more spatially-extended) feedback, the negative radial gradients invert, becoming positive for heights in excess of 1 kpc. Such behaviour is consistent with that inferred from recent observations. Our measurements of the vertical metallicity gradients show no clear correlation with galactocentric radius, and are in good agreement with those observed in the Milky Way's thick disc (locally). Conclusions. Simulations employing stronger/more extended feedback prescriptions possess radial and vertical metallicity and kinematic gradients more in line with recent observations. The inverted, positive, radial metallicity gradients seen in the simulated thick stellar discs originate from a population of younger, more metal-rich, stars formed in-situ, superimposed upon a background population of older migrators from the inner disc; the contrast provided by the former increases radially, due to the inside-out growth of the disc. A similar behaviour may be responsible for the same flattening seen in the radial gradients with scaleheight in the Milky Way.

preprint2015arXiv

Stellar orbital properties as diagnostics of the origin of the stellar halo

We examine metallicities, ages and orbital properties of halo stars in a Milky-Way like disk galaxy formed in the cosmological hydrodynamical MaGICC simulations. Halo stars were either accreted from satellites or they formed in situ in the disk or bulge of the galaxy and were then kicked up into the halo ("in situ/ kicked-up" stars). Regardless of where they formed both types show surprisingly similar orbital properties: the majority of stars of both types are on short-axis tubes with the same sense of rotation as the disk -- implying that a large fraction of satellites are accreted onto the halo with the same sense of angular momentum as the disk.

preprint2014arXiv

Discovery of a Gas-Rich Companion to the Extremely Metal-Poor Galaxy DDO 68

We present HI spectral-line imaging of the extremely metal-poor galaxy DDO 68. This system has a nebular oxygen abundance of only 3% Z$_{\odot}$, making it one of the most metal-deficient galaxies known in the local volume. Surprisingly, DDO 68 is a relatively massive and luminous galaxy for its metal content, making it a significant outlier in the mass-metallicity and luminosity-metallicity relationships. The origin of such a low oxygen abundance in DDO 68 presents a challenge for models of the chemical evolution of galaxies. One possible solution to this problem is the infall of pristine neutral gas, potentially initiated during a gravitational interaction. Using archival HI spectral-line imaging obtained with the Karl G. Jansky Very Large Array, we have discovered a previously unknown companion of DDO 68. This low-mass (M$_{\rm HI}$ $=$ 2.8$\times$10$^{7}$ M$_{\odot}$), recently star-forming (SFR$_{\rm FUV}$ $=$ 1.4$\times$10$^{-3}$ M$_{\odot}$ yr$^{-1}$, SFR$_{\rm Hα}$ $<$ 7$\times$10$^{-5}$ M$_{\odot}$ yr$^{-1}$) companion has the same systemic velocity as DDO 68 (V$_{\rm sys}$ $=$ 506 km s$^{-1}$; D $=$ 12.74$\pm$0.27 Mpc) and is located at a projected distance of 42 kpc. New HI maps obtained with the 100m Robert C. Byrd Green Bank Telescope provide evidence that DDO 68 and this companion are gravitationally interacting at the present time. Low surface brightness HI gas forms a bridge between these objects.

preprint2014arXiv

Globular Cluster Systems in Brightest Cluster Galaxies: A Near-Universal Luminosity Function?

We present the first results from our HST Brightest Cluster Galaxy (BCG) survey of seven central supergiant cluster galaxies and their globular cluster (GC) systems. We measure a total of 48000 GCs in all seven galaxies, representing the largest single GC database. We find that a log-normal shape accurately matches the observed luminosity function (LF) of the GCs down to the GCLF turnover point, which is near our photometric limit. In addition, the LF has a virtually identical shape in all seven galaxies. Our data underscore the similarity in the formation mechanism of massive star clusters in diverse galactic environments. At the highest luminosities (log L > 10^7 L_Sun) we find small numbers of "superluminous" objects in five of the galaxies; their luminosity and color ranges are at least partly consistent with those of UCDs (Ultra-Compact Dwarfs). Lastly, we find preliminary evidence that in the outer halo (R > 20 kpc), the LF turnover point shows a weak dependence on projected distance, scaling as L_0 ~ R^-0.2, while the LF dispersion remains nearly constant.

preprint2014arXiv

On the relation between metallicity and RGB color in HST/ACS data

The determination of stellar metallicity and its gradient in external galaxies is a difficult task, but crucial for the understanding of galaxy formation and evolution. The color of the Red Giant Branch (RGB) can be used to determine metallicities of stellar populations that have only shallow photometry. We will quantify the relation between metallicity and color in the widely used HST ACS filters F606W and F814W. We use a sample of globular clusters from the ACS Globular Cluster Survey and measure their RGB color at given absolute magnitudes to derive the color-metallicity relation. We especially investigate the scatter and the uncertainties in this relation and show its limitations. There is a clear relation between metallicity and RGB color. A comparison with isochrones shows reasonably good agreement with BaSTI models, a small offset to Dartmouth models, and a larger offset to Padua models. Even for the best globular cluster data available, the metallicity of a simple stellar population can be determined from the RGB alone only with an accuracy of 0.3dex for [M/H]< -1, and 0.15dex for [M/H]> -1. For mixed populations, as they are observed in external galaxies, the uncertainties will be even larger due to uncertainties in extinction, age, etc. Therefore caution is necessary when interpreting photometric metallicities.

preprint2014arXiv

Systematic Problems With Using Dark Matter Simulations to Model Stellar Halos

The limits of available computing power have forced models for the structure of stellar halos to adopt one or both of the following simplifying assumptions: (1) stellar mass can be "painted" onto dark matter particles in progenitor satellites; (2) pure dark matter simulations that do not form a luminous galaxy can be used. We estimate the magnitude of the systematic errors introduced by these assumptions using a controlled set of stellar halo models where we independently vary whether we look at star particles or painted dark matter particles, and whether we use a simulation in which a baryonic disk galaxy forms or a matching pure dark matter simulation that does not form a baryonic disk. We find that the "painting" simplification reduces the halo concentration and internal structure, predominantly because painted dark matter particles have different kinematics than star particles even when both are buried deep in the potential well of the satellite. The simplification of using pure dark matter simulations reduces the concentration further, but increases the internal structure, and results in a more prolate stellar halo. These differences can be a factor of 1.5-7 in concentration (as measured by the half-mass radius) and 2-7 in internal density structure. Given this level of systematic uncertainty, one should be wary of overinterpreting differences between observations and the current generation of stellar halo models based on dark matter only simulations when such differences are less than an order of magnitude.

preprint2013arXiv

Constraining the age of the NGC 4565 HI Disk Warp: Determining the Origin of Gas Warps

We have mapped the distribution of young and old stars in the gaseous HI warp of NGC 4565. We find a clear correlation of young stars (<600 Myr) with the warp, but no coincident old stars (>1 Gyr), which places an upper limit on the age of the structure. The formation rate of the young stars, which increased ~300 Myr ago relative to the surrounding regions, is (6.3 +2.5/-1.5) x 10^-5 M_sol/yr/kpc^2. This implies a ~60+/-20 Gyr depletion time of the HI warp, similar to the timescales calculated for the outer HI disks of nearby spiral galaxies. While some stars associated with the warp fall into the asymptotic giant branch (AGB) region of the color magnitude diagram, where stars could be as old as 1 Gyr, further investigation suggests that they may be interlopers rather than real AGB stars. We discuss the implications of these age constraints for the formation of HI warps, and the gas fueling of disk galaxies.

preprint2013arXiv

GHOSTS I: A New Faint very Isolated Dwarf Galaxy at D = 12 +/- 2 Mpc

We report the discovery of a new faint dwarf galaxy, GHOSTS I, using HST/ACS data from one of our GHOSTS (Galaxy Halos, Outer disks, Substructure, Thick disk, and Star clusters) fields. Its detected individual stars populate an approximately one magnitude range of its luminosity function (LF). Using synthetic color-magnitude diagrams (CMDs) to compare with the galaxy's CMD, we find that the colors and magnitudes of GHOSTS I's individual stars are most consistent with being young helium-burning and asymptotic giant branch stars at a distance of 12 +/- 2 Mpc. Morphologically, GHOSTS I appears to be actively forming stars, so we tentatively classify it as a dwarf irregular (dIrr) galaxy, although future HST observations deep enough to resolve a larger magnitude range in its LF are required to make a more secure classification. GHOSTS I's absolute magnitude is $M_V = -9.85^{+ 0.40}_{- 0.33}$, making it one of the least luminous dIrr galaxies known, and its metallicity is lower than [Fe/H] =-1.5 dex. The half-light radius of GHOSTS I is 226 +/- 38 pc and its ellipticity is 0.47 +/- 0.07, similar to Milky Way and M31 dwarf satellites at comparable luminosity. There are no luminous massive galaxies or galaxy clusters within ~ 4 Mpc from GHOSTS I that could be considered as its host, making it a very isolated dwarf galaxy in the Local Universe.

preprint2013arXiv

Halo orbits in cosmological disk galaxies: tracers of formation history

We analyze the orbits of stars and dark matter particles in the halo of a disk galaxy formed in a cosmological hydrodynamical simulation. The halo is oblate within the inner ~20 kpc and triaxial beyond this radius. About 43% of orbits are short axis tubes - the rest belong to orbit families that characterize triaxial potentials (boxes, long-axis tubes and chaotic orbits), but their shapes are close to axisymmetric. We find no evidence that the self-consistent distribution function of the nearly oblate inner halo is comprised primarily of axisymmetric short-axis tube orbits. Orbits of all families, and both types of particles are highly eccentric with mean eccentricity >0.6. We find that randomly selected samples of halo stars show no substructure in "integrals of motion" space. However individual accretion events can be clearly identified in plots of metallicity versus formation time. Dynamically young tidal debris is found primarily on a single type of orbit. However, stars associated with older satellites become chaotically mixed during the formation process (possibly due to scattering by the central bulge and disk, and baryonic processes), and appear on all four types of orbits. We find that the tidal debris in cosmological hydrodynamical simulations experiences significantly more chaotic evolution than in collisionless simulations, making it much harder to identify individual progenitors using phase space coordinates alone. However by combining information on stellar ages and chemical abundances with the orbital properties of halo stars in the underlying self-consistent potential, the identification of progenitors is likely to be possible.

preprint2013arXiv

Testing galaxy formation models with the GHOSTS survey: The color profile of M81's stellar halo

We study the properties of the stellar populations in M81's outermost part, which hereafter we will term the stellar halo, using HST ACS/WFC observations of 19 fields from the GHOSTS survey. The observed fields probe the stellar halo out to a projected distance of ~ 50 kpc from the galactic center. Each field was observed in both F606W and F814W filters. The 50% completeness levels of the color magnitude diagrams (CMDs) are typically at 2 mag below the tip of the red giant branch (TRGB). Fields at distances closer than 15 kpc show evidence of disk-dominated populations whereas fields at larger distances are mostly populated by halo stars. The RGB of the M81's halo CMDs is well matched with isochrones of ~ 10 Gyr and metallicities [Fe/H] ~ -1.2 dex, suggesting that the dominant stellar population of M81's halo has a similar age and metallicity. The halo of M81 is characterized by a color distribution of width ~ 0.4 mag and an approximately constant median value of (F606W - F814W) ~ 1 mag measured using stars within the magnitude range 23.7 < F814W < 25.5. When considering only fields located at galactocentric radius R > 15 kpc, we detect no color gradient in the stellar halo of M81. We place a limit of 0.03+/-0.11 mag difference between the median color of RGB M81 halo stars at ~ 15 and at 50 kpc, corresponding to a metallicity difference of 0.08+/-0.35 dex over that radial range for an assumed constant age of 10 Gyr. We compare these results with model predictions for the colors of stellar halos formed purely via accretion of satellite galaxies. When we analyze the cosmologically motivated models in the same way as the HST data, we find that they predict no color gradient for the stellar halos, in good agreement with the observations.

preprint2011arXiv

The Resolved Stellar Halo of NGC 253

We have obtained Magellan/IMACS and HST/ACS imaging data that resolve red giant branch stars in the stellar halo of the starburst galaxy NGC 253. The HST data cover a small area, and allow us to accurately interpret the ground-based data, which cover 30% of the halo to a distance of 30 kpc, allowing us to make detailed quantitative measurements of the global properties and structure of a stellar halo outside of the Local Group. The geometry of the halo is significantly flattened in the same sense as the disk, with a projected axis ratio of b/a ~ 0.35 +/- 0.1. The total stellar mass of the halo is estimated to be M_halo ~ 2.5 +/- 1.5 x 10^9 M_sun, or 6% of the total stellar mass of the galaxy, and has a projected radial dependence that follows a power law of index -2.8 +/- 0.6, corresponding to a three-dimensional power law index of ~ -4. The total luminosity and profile shape that we measure for NGC 253 are somewhat larger and steeper than the equivalent values for the Milky Way and M31, but are well within the scatter of model predictions for the properties of stellar halos built up in a cosmological context. Structure within the halo is seen at a variety of scales: there is small kpc-scale density variation and a large shelf-like feature near the middle of the field. The techniques that have been developed will be essential for quantitatively comparing our upcoming larger sample of observed stellar halos to models of halo formation.

preprint2010arXiv

Consequences of cosmic microwave background-regulated star formation

It has been hypothesized that the cosmic microwave background (CMB) provides a temperature floor for collapsing protostars that can regulate the process of star formation and result in a top-heavy initial mass function at high metallicity and high redshift. We examine whether this hypothesis has any testable observational consequences. First we determine, using a set of hydrodynamic galaxy formation simulations, that the CMB temperature floor would have influenced the majority of stars formed at redshifts between z=3 and 6, and probably even to higher redshift. Five signatures of CMB-regulated star formation are: (1) a higher supernova rate than currently predicted at high redshift; (2) a systematic discrepancy between direct and indirect measurements of the high redshift star formation rate; (3) a lack of surviving globular clusters that formed at high metallicity and high redshift; (4) a more rapid rise in the metallicity of cosmic gas than is predicted by current simulations; and (5) an enhancement in the abundances of alpha elements such as O and Mg at metallicities -2 < [Fe/H] < -0.5. Observations are not presently able to either confirm or rule out the presence of these signatures. However, if correct, the top-heavy IMF of high-redshift high-metallicity globular clusters could provide an explanation for the observed bimodality of their metallicity distribution.

preprint2009arXiv

A Blue Tilt in the Globular Cluster System of the Milky Way-like Galaxy NGC 5170

Here we present HST/ACS imaging, in the B and I bands, of the edge-on Sb/Sc galaxy NGC 5170. Excluding the central disk region region, we detect a 142 objects with colours and sizes typical of globular clusters (GCs). Our main result is the discovery of a `blue tilt' (a mass-metallicity relation), at the 3sigma level, in the metal-poor GC subpopulation of this Milky Way like galaxy. The tilt is consistent with that seen in massive elliptical galaxies and with the self enrichment model of Bailin & Harris. For a linear mass-metallicity relation, the tilt has the form Z ~ L^{0.42 +/- 0.13}. We derive a total GC system population of 600 +/- 100, making it much richer than the Milky Way. However when this number is normalised by the host galaxy luminosity or stellar mass it is similar to that of M31. Finally, we report the presence of a potential Ultra Compact Dwarf of size ~ 6 pc and luminosity M_I ~ -12.5, assuming it is physically associated with NGC 5170.

preprint2009arXiv

Stochastic self-enrichment, pre-enrichment, and the formation of globular clusters

We develop a model for stochastic pre-enrichment and self-enrichment in globular clusters (GCs) during their formation process. GCs beginning their formation have an initial metallicity determined by the pre-enrichment of their surrounding protocloud, but can also undergo internal self-enrichment during formation. Stochastic variations in metallicity arise because of the finite numbers of supernova. We construct an analytic formulation of the combined effects of pre-enrichment and self-enrichment and use Monte Carlo models to verify that the model accurately encapsulates the mean metallicity and metallicity spread among real GCs. The predicted metallicity spread due to self-enrichment alone, a robust prediction of the model, is much smaller than the observed spread among real GCs. This result rules out self-enrichment as a significant contributor to the metal content in most GCs, leaving pre-enrichment as the viable alternative. Self-enrichment can, however, be important for clusters with masses well above 10^6 Msun, which are massive enough to hold in a significant fraction of their SN ejecta even without any external pressure confinement. This transition point corresponds well to the mass at which a mass-metallicity relationship ("blue tilt") appears in the metal-poor cluster sequence in many large galaxies. We therefore suggest that self-enrichment is the primary driver for the mass-metallicity relation. Other predictions from our model are that the cluster-to-cluster metallicity spread decreases amongst the highest mass clusters; and that the red GC sequence should also display a more modest mass-metallicity trend if it can be traced to similarly high mass.