Showing 1–50 of 82 results for author: Bailin, J

First resolved stellar halo kinematics of a MW-mass galaxy outside the Local Group: A flat counter-rotating halo in NGC 4945

Authors: Camila Beltrand, Antonela Monachesi, Richard D'Souza, Eric F. Bell, Roelof S. de Jong, Facundo A. Gomez, Jeremy Bailin, In Sung Jang, Adam Smercina

Abstract: Stellar halos of galaxies, primarily formed through the accretion of smaller objects, are important to understand the hierarchical mass assembly of galaxies. However, the inner regions of stellar halos in disk galaxies are predicted to have an in-situ component that is expected to be prominent along the major axis. Kinematic information is crucial to disentangle the contribution of the in-situ com… ▽ More Stellar halos of galaxies, primarily formed through the accretion of smaller objects, are important to understand the hierarchical mass assembly of galaxies. However, the inner regions of stellar halos in disk galaxies are predicted to have an in-situ component that is expected to be prominent along the major axis. Kinematic information is crucial to disentangle the contribution of the in-situ component from the accreted stellar halos. The low surface brightness of stellar halos makes it inaccessible with traditional integrated light spectroscopy. In this work, using a novel technique, we study the kinematics of the stellar halo of the edge-on galaxy NGC 4945. We couple new deep Multi Unit Spectroscopic Explorer spectroscopic observations with existing Hubble Space Telescope imaging data to spectroscopically measure the line-of-sight (LOS) heliocentric velocity and velocity dispersion in two fields at a galactocentric distance of 12.2 kpc (outer disk field) and 34.6 kpc (stellar halo field) along NGC 4945 major axis, by stacking individual spectra of red giant branch and asymptotic giant branch stars. We obtain a LOS velocity and dispersion of 673+/-11 km/s and 73+/-14 km/s, respectively, for the outer disk field. This is consistent with the mean HI velocity of the disk at that distance. For the halo field we obtain a LOS velocity and dispersion of 519+/-12 km/s and 42+/-22 km/s. The halo fields' velocity measurement is within ~40 km/s from the systemic LOS velocity of NGC 4945, which is 563 km/s, suggesting that its stellar halo at 34.6 kpc along the major axis is counter-rotating and is of likely accretion origin. This provides the first ever kinematic measurement of the stellar halo of a Milky Way-mass galaxy outside the Local Group from its resolved stellar population, and establishes a powerful technique for measuring the velocity field of the stellar halos of nearby galaxies. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: Accepted to A&A. 23 pages, 12 figures, 2 appendices

Spatial Variations of Stellar Elemental Abundances in FIRE Simulations of Milky Way-Mass Galaxies: Patterns Today Mostly Reflect Those at Formation

Authors: Russell L. Graf, Andrew Wetzel, Matthew A. Bellardini, Jeremy Bailin

Abstract: Spatial patterns of stellar elemental abundances encode rich information about a galaxy's formation history. We analyze the radial, vertical, and azimuthal variations of metals in stars, both today and at formation, in the FIRE-2 cosmological simulations of Milky Way (MW)-mass galaxies, and we compare with the MW. The radial gradient today is steeper (more negative) for younger stars, which agrees… ▽ More Spatial patterns of stellar elemental abundances encode rich information about a galaxy's formation history. We analyze the radial, vertical, and azimuthal variations of metals in stars, both today and at formation, in the FIRE-2 cosmological simulations of Milky Way (MW)-mass galaxies, and we compare with the MW. The radial gradient today is steeper (more negative) for younger stars, which agrees with the MW, although radial gradients are shallower in FIRE-2. Importantly, this age dependence was present already at birth: radial gradients today are only modestly ($\lesssim$ 0.01 dex kpc$^{-1}$) shallower than at birth. Disk vertical settling gives rise to negative vertical gradients across all stars, but vertical gradients of mono-age stellar populations are weak. Similar to the MW, vertical gradients in FIRE-2 are shallower at larger radii, but they are overall shallower in FIRE-2. This vertical dependence was present already at birth: vertical gradients today are only modestly ($\lesssim$ 0.1 dex kpc$^{-1}$) shallower than at birth. Azimuthal scatter is nearly constant with radius, and it is nearly constant with age $\lesssim$ 8 Gyr ago, but increases for older stars. Azimuthal scatter is slightly larger ($\lesssim$ 0.04 dex) today than at formation. Galaxies with larger azimuthal scatter have a stronger radial gradient, implying that azimuthal scatter today arises primarily from radial redistribution of gas and stars. Overall, spatial variations of stellar metallicities show only modest differences between formation and today; spatial variations today primarily reflect the conditions of stars at birth, with spatial redistribution of stars after birth contributing secondarily. △ Less

Submitted 27 February, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: 24 pages, submitted to ApJ

Origins of the Evil Eye: M64's Stellar Halo Reveals the Recent Accretion of an SMC-mass Satellite

Authors: Adam Smercina, Eric F. Bell, Paul A. Price, Jeremy Bailin, Julianne J. Dalcanton, Roelof S. de Jong, Richard D'Souza, Katya Gozman, In Sung Jang, Antonela Monachesi, David Nidever, Colin T. Slater

Abstract: M64, often called the "Evil Eye" galaxy, is unique among local galaxies. Beyond its dramatic, dusty nucleus, it also hosts an outer gas disk that counter-rotates relative to its stars. The mass of this outer disk is comparable to the gas content of the Small Magellanic Cloud (SMC), prompting the idea that it was likely accreted in a recent minor merger. Yet, detailed follow-up studies of M64's out… ▽ More M64, often called the "Evil Eye" galaxy, is unique among local galaxies. Beyond its dramatic, dusty nucleus, it also hosts an outer gas disk that counter-rotates relative to its stars. The mass of this outer disk is comparable to the gas content of the Small Magellanic Cloud (SMC), prompting the idea that it was likely accreted in a recent minor merger. Yet, detailed follow-up studies of M64's outer disk have shown no evidence of such an event, leading to other interpretations, such as a "flyby" interaction with the distant diffuse satellite Coma P. We present Subaru Hyper Suprime-Cam observations of M64's stellar halo, which resolve its stellar populations and reveal a spectacular radial shell feature, oriented $\sim$30$^{\circ}$ relative to the major axis and along the rotation axis of the outer gas disk. The shell is $\sim$45 kpc southeast of M64, while a similar but more diffuse plume to the northwest extends to $>$100 kpc. We estimate a stellar mass and metallicity for the southern shell of $M_{\star} {=} 1.80~{\pm}~0.54{\times}10^8~M_{\odot}$ and [M/H] $=$ $-$1.0, respectively, and a similar mass of $1.42~{\pm}~0.71{\times}10^8 M_{\odot}$ for the northern plume. Taking into account the accreted material in M64's inner disk, we estimate a total stellar mass for the progenitor satellite of $M_{\rm \star,prog}~{\simeq}~5{\times}10^8~M_{\odot}$. These results suggest that M64 is in the final stages of a minor merger with a gas-rich satellite strikingly similar to the SMC, in which M64's accreted counter-rotating gas originated, and which is responsible for the formation of its dusty inner star-forming disk. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: 12 pages, 5 figures, accepted for publication in ApJ Letters

Saying Hallo to M94's Stellar Halo: Investigating the Accretion History of the Largest Pseudobulge Host in the Local Universe

Authors: Katya Gozman, Eric F. Bell, Adam Smercina, Paul Price, Jeremy Bailin, Roelof S. de Jong, Richard D'Souza, In Sung Jang, Antonela Monachesi, Colin Slater

Abstract: It is not yet settled how the combination of secular processes and merging gives rise to the bulges and pseudobulges of galaxies. The nearby ($D\sim$ 4.2 Mpc) disk galaxy M94 (NGC 4736) has the largest pseudobulge in the local universe, and offers a unique opportunity for investigating the role of merging in the formation of its pseudobulge. We present a first ever look at M94's stellar halo, whic… ▽ More It is not yet settled how the combination of secular processes and merging gives rise to the bulges and pseudobulges of galaxies. The nearby ($D\sim$ 4.2 Mpc) disk galaxy M94 (NGC 4736) has the largest pseudobulge in the local universe, and offers a unique opportunity for investigating the role of merging in the formation of its pseudobulge. We present a first ever look at M94's stellar halo, which we expect to contain a fossil record of M94's past mergers. Using Subaru's Hyper Suprime-Cam, we resolve and identify red giant branch (RGB) stars in M94's halo, finding two distinct populations. After correcting for completeness through artificial star tests, we can measure the radial profile of each RGB population. The metal-rich RGB stars show an unbroken exponential profile to a radius of 30 kpc that is a clear continuation of M94's outer disk. M94's metal poor stellar halo is detectable over a wider area and clearly separates from its metal-rich disk. By integrating the halo density profile, we infer a total accreted stellar mass of $\sim 2.8 \times 10^8 M_\odot$, with a median metallicity of [M/H] $=-$1.4. This indicates that M94's most-massive past merger was with a galaxy similar to, or less massive than, the Small Magellanic Cloud. Few nearby galaxies have had such a low-mass dominant merger; therefore we suggest that M94's pseudobulge was not significantly impacted by merging. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 19 pages, 12 figures, 1 table; Published in ApJ on April 14, 2023

The imprint of bursty star formation on alpha-element abundance patterns in Milky Way-like galaxies

Authors: Hanna Parul, Jeremy Bailin, Andrew Wetzel, Alexander B. Gurvich, Claude-André Faucher-Giguère, Zachary Hafen, Jonathan Stern, Owain Snaith

Abstract: Milky Way-mass galaxies in the FIRE-2 simulations demonstrate two main modes of star formation. At high redshifts star formation occurs in a series of short and intense bursts, while at low redshifts star formation proceeds at a steady rate with a transition from one mode to another at times ranging from 3 to 7 Gyr ago for different galaxies. We analyse how the mode of star formation affects iron… ▽ More Milky Way-mass galaxies in the FIRE-2 simulations demonstrate two main modes of star formation. At high redshifts star formation occurs in a series of short and intense bursts, while at low redshifts star formation proceeds at a steady rate with a transition from one mode to another at times ranging from 3 to 7 Gyr ago for different galaxies. We analyse how the mode of star formation affects iron and alpha-element abundance. We find that the early bursty regime imprints a measurable pattern in stellar elemental abundances in the form of a "sideways chevron" shape on the [Fe/H] - [O/Fe] plane and the scatter in [O/Fe] at a given stellar age is higher than when a galaxy is in the steady regime. That suggests that the evolution of [O/Fe] scatter with age provides an estimate of the end of the bursty phase. We investigate the feasibility of observing of this effect by adding mock observational errors to a simulated stellar survey and find that the transition between the bursty and steady phase should be detectable in the Milky Way, although larger observational uncertainties make the transition shallower. We apply our method to observations of the Milky Way from the Second APOKASC Catalog and estimate that the transition to steady star formation in the Milky Way happened 7-8 Gyrs ago, earlier than transition times measured in the simulations. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 15 pages; Accepted 17 January 2023 in MNRAS

Orientations of DM Halos in FIRE-2 Milky Way-mass Galaxies

Authors: Jay Baptista, Robyn Sanderson, Dan Huber, Andrew Wetzel, Omid Sameie, Michael Boylan-Kolchin, Jeremy Bailin, Philip F. Hopkins, Claude-Andre Faucher-Giguere, Sukanya Chakrabarti, Drona Vargya, Nondh Panithanpaisal, Arpit Arora, Emily Cunningham

Abstract: The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way's DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FI… ▽ More The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way's DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FIRE-2 Milky Way-mass galaxies to explore the evolution of the DM halo's orientation with radius and time, with or without a major merger with a Large Magellanic Cloud (LMC) analog, and when varying the DM model. In three of the four CDM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20 degrees beyond about 30 galactocentric kpc, reaching a maximum of 30--90 degrees depending on the individual halo's formation history. In identical simulations using a model of self-interacting DM with $σ= 1 \, \mathrm{cm}^2 \, \mathrm{g}^{-1}$, the halo remains aligned with the stellar disk out to $\sim$200--400 kpc. Interactions with massive satellites ($M \gtrsim 4 \times 10^{10} \, \rm{M_\odot}$ at pericenter; $M \gtrsim 3.3 \times 10^{10} \, \rm{M_\odot}$ at infall) affect the orientation of the halo significantly, aligning the halo's major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of SIDM if the effects of massive satellites can be accounted for. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Submitted to ApJ, 19 pages, 12 figures, 3 tables

The observable properties of galaxy accretion events in Milky Way-like galaxies in the FIRE-2 cosmological simulations

Authors: Danny Horta, Emily C. Cunningham, Robyn E. Sanderson, Kathryn V. Johnston, Nondh Panithanpaisal, Arpit Arora, Lina Necib, Andrew Wetzel, Jeremy Bailin, Claude-André Faucher-Giguère

Abstract: In the $Λ$-Cold Dark Matter model of the Universe, galaxies form in part through accreting satellite systems. Previous work have built an understanding of the signatures of these processes contained within galactic stellar halos. This work revisits that picture using seven Milky Way-like galaxies in the \textit{Latte} suite of FIRE-2 cosmological simulations. The resolution of these simulations al… ▽ More In the $Λ$-Cold Dark Matter model of the Universe, galaxies form in part through accreting satellite systems. Previous work have built an understanding of the signatures of these processes contained within galactic stellar halos. This work revisits that picture using seven Milky Way-like galaxies in the \textit{Latte} suite of FIRE-2 cosmological simulations. The resolution of these simulations allows a comparison of contributions from satellites above M$_{*}$$\gtrsim$10$\times$$^{7}$M$_{\odot}$, enabling the analysis of observable properties for disrupted satellites in a fully self-consistent and cosmological context. Our results show that, the time of accretion and the stellar mass of an accreted satellite are fundamental parameters that in partnership dictate the resulting spatial distribution, orbital energy, and [$α$/Fe]-[Fe/H] compositions of the stellar debris of such mergers $at$ $present$ $day$. These parameters also govern the resulting dynamical state of an accreted galaxy at $z=0$, leading to the expectation that the inner regions of the stellar halo (R$_{\mathrm{GC}}$ $\lesssim$30 kpc) should contain fully phase-mixed debris from both lower and higher mass satellites. In addition, we find that a significant fraction of the lower mass satellites accreted at early times deposit debris in the outer halo (R$_{\mathrm{GC}}$ $>$50 kpc) that are $not$ fully phased-mixed, indicating that they could be identified in kinematic surveys. Our results suggest that, as future surveys become increasingly able to map the outer halo of our Galaxy, they may reveal the remnants of long-dead dwarf galaxies whose counterparts are too faint to be seen $in$ $situ$ in higher redshift surveys. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: Submitted for publication in ApJ

Coupling Semi-Analytic and N-body Galaxies (CoSANG) for cosmological stellar halo simulations I- Methods and the structure of dark matter halos

Authors: Shahram Talei, Krista Mccord, Jeremy Bailin, Darren J. Croton, Alexandra Mannings, Michael Sitarz, Annelia Anderson, Brooke Bailey, Alexis Rollins

Abstract: We present CoSANG (Coupling Semi-Analytic and N-body Galaxies), a new hybrid model for cosmological dark matter and stellar halo simulations. In this approach a collisionless model (Gadget3) is used for gravitational interactions while a coupled semi-analytic model (SAGE) calculates baryonic effects at each time-step. This live self-consistent interaction at each time-step is the key difference be… ▽ More We present CoSANG (Coupling Semi-Analytic and N-body Galaxies), a new hybrid model for cosmological dark matter and stellar halo simulations. In this approach a collisionless model (Gadget3) is used for gravitational interactions while a coupled semi-analytic model (SAGE) calculates baryonic effects at each time-step. This live self-consistent interaction at each time-step is the key difference between CoSANG and traditional semi-analytic models that are mainly used for post-processing. By accounting for the gravitational effect of the baryons, CoSANG can overcome some of the deficiencies of pure N-body simulations, while being less computationally expensive than hydrodynamic simulations. Moreover CoSANG can produce stellar halo populations via tagging tracer dark matter particles. We demonstrate the performance and dynamical accuracy of this approach using both controlled test simulations and a set of three cosmological zoom-in simulations of Milky Way (MW) mass halos. We simulate each target halo both without the coupling (dark matter only, hereafter DMO) and with the coupling (CoSANG). We compare the internal structure (subhalo distribution, shape and orientation) of the halos. The following changes are observed in the CoSANG model compared to the DMO model: 1) the total number of subhaloes close to the center of the halo is reduced, 2) the $V_{\mathrm{max}}$ distribution peaks at a lower value and lies below the DMO model, 3) the axis ratio is smaller. The difference between DMO and CoSANG simulations is more significant in early forming halos. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 24 pages, 15 figures, submitted to MNRAS, comments are welcome

Ultrafaint Dwarf Galaxy Candidates in the M81 Group: Signatures of Group Accretion

Authors: Eric F. Bell, Adam Smercina, Paul A. Price, Richard D'Souza, Jeremy Bailin, Roelof S. de Jong, Katya Gozman, In Sung Jang, Antonela Monachesi, Oleg Y. Gnedin, Colin T. Slater

Abstract: The faint and ultrafaint dwarf galaxies in the Local Group form the observational bedrock upon which our understanding of small-scale cosmology rests. In order to understand whether this insight generalizes, it is imperative to use resolved-star techniques to discover similarly faint satellites in nearby galaxy groups. We describe our search for ultrafaint galaxies in the M81 group using deep grou… ▽ More The faint and ultrafaint dwarf galaxies in the Local Group form the observational bedrock upon which our understanding of small-scale cosmology rests. In order to understand whether this insight generalizes, it is imperative to use resolved-star techniques to discover similarly faint satellites in nearby galaxy groups. We describe our search for ultrafaint galaxies in the M81 group using deep ground-based resolved-star data sets from Subaru's Hyper Suprime-Cam. We present one new ultrafaint dwarf galaxy in the M81 group and identify five additional extremely low surface brightness candidate ultrafaint dwarfs that reach deep into the ultrafaint regime to $M_V \sim -6$ (similar to current limits for Andromeda satellites). These candidates' luminosities and sizes are similar to known Local Group dwarf galaxies Tucana B, Canes Venatici I, Hercules, and Boötes I. Most of these candidates are likely to be real, based on tests of our techniques on blank fields. Intriguingly, all of these candidates are spatially clustered around NGC 3077, which is itself an M81 group satellite in an advanced state of tidal disruption. This is somewhat surprising, as M81 itself and its largest satellite M82 are both substantially more massive than NGC 3077 and by virtue of their greater masses, would have been expected to host as many or more ultrafaint candidates. These results lend considerable support to the idea that satellites of satellites are an important contribution to the growth of satellite populations around Milky Way-mass galaxies. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: The Astrophysical Journal Letters; in press. 11 pages, 4 figures, 1 table

New Globular Cluster Candidates in the M81 group

Authors: Jiaming Pan, Eric F. Bell, Adam Smercina, Paul Price, Colin T. Slater, Jeremy Bailin, Roelof S. de Jong, Richard D'Souza, In Sung Jang, Antonela Monachesi

Abstract: 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 prev… ▽ More 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. △ Less

Submitted 14 June, 2022; originally announced June 2022.

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

Authors: Jenna Samuel, Andrew Wetzel, Isaiah Santistevan, Erik Tollerud, Jorge Moreno, Michael Boylan-Kolchin, Jeremy Bailin, Bhavya Pardasani

Abstract: 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 s… ▽ More 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. △ Less

Submitted 17 June, 2022; v1 submitted 14 March, 2022; originally announced March 2022.

Comments: 17 pages, 13 figures, + appendices. Accepted by MNRAS

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

Authors: Matthew A. Bellardini, Andrew Wetzel, Sarah R. Loebman, Jeremy Bailin

Abstract: 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 co… ▽ More 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. △ Less

Submitted 9 June, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: 22 pages, 16 figures, Accepted for publication in MNRAS

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

Authors: Jeremy Bailin, Ryker von Klar

Abstract: 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 p… ▽ More 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. △ Less

Submitted 27 October, 2021; originally announced October 2021.

Comments: Accepted for publication in ApJ. Machine readable full version of Table 1 will be uploaded to CDS once full volume/page number is known; until then it is available at http://simulated-galaxies.ua.edu/gcsigpred_table1.txt

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

Authors: Preet B. Patel, Sarah R. Loebman, Andrew Wetzel, Claude-André Faucher-Giguère, Kareem El-Badry, Jeremy Bailin

Abstract: 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 explor… ▽ More 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. △ Less

Submitted 23 March, 2022; v1 submitted 15 October, 2021; originally announced October 2021.

Comments: v1: 14 pages, 8 figures. Submitted to MNRAS on Oct. 18th, 2021 | v2: 15 pages, 8 figures. Accepted by MNRAS on Mar. 17th, 2022

Predicting the spectrum of UGC 2885, Rubin's Galaxy with machine learning

Authors: Benne W. Holwerda, John F. Wu, William C. Keel, Jason Young, Ren Mullins, Joannah Hinz, K. E. Saavik Ford, Pauline Barmby, Rupali Chandar, Jeremy Bailin, Josh Peek, Tim Pickering, Torsten Böker

Abstract: Wu & Peek (2020) predict SDSS-quality spectra based on Pan-STARRS broad-band \textit{grizy} images using machine learning (ML). In this letter, we test their prediction for a unique object, UGC 2885 ("Rubin's galaxy"), the largest and most massive, isolated disk galaxy in the local Universe ($D<100$ Mpc). After obtaining the ML predicted spectrum, we compare it to all existing spectroscopic inform… ▽ More Wu & Peek (2020) predict SDSS-quality spectra based on Pan-STARRS broad-band \textit{grizy} images using machine learning (ML). In this letter, we test their prediction for a unique object, UGC 2885 ("Rubin's galaxy"), the largest and most massive, isolated disk galaxy in the local Universe ($D<100$ Mpc). After obtaining the ML predicted spectrum, we compare it to all existing spectroscopic information that is comparable to an SDSS spectrum of the central region: two archival spectra, one extracted from the VIRUS-P observations of this galaxy, and a new, targeted MMT/Binospec observation. Agreement is qualitatively good, though the ML prediction prefers line ratios slightly more towards those of an active galactic nucleus (AGN), compared to archival and VIRUS-P observed values. The MMT/Binospec nuclear spectrum unequivocally shows strong emission lines except H$β$, the ratios of which are consistent with AGN activity. The ML approach to galaxy spectra may be a viable way to identify AGN supplementing NIR colors. How such a massive disk galaxy ($M^* = 10^{11}$ M$_\odot$), which uncharacteristically shows no sign of interaction or mergers, manages to fuel its central AGN remains to be investigated. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: 9 pages, 5 figures, submitted to ApJL

The Galaxy Progenitors of Stellar Streams around Milky Way-mass Galaxies in the FIRE Cosmological Simulations

Authors: Nondh Panithanpaisal, Robyn E. Sanderson, Andrew Wetzel, Emily C. Cunningham, Jeremy Bailin, Claude-André Faucher-Giguère

Abstract: Stellar streams record the accretion history of their host galaxy. We present a set of simulated streams from disrupted dwarf galaxies in 13 cosmological simulations of Milky Way (MW)-mass galaxies from the FIRE-2 suite at $z=0$, including 7 isolated Milky Way-mass systems and 6 hosts resembling the MW-M31 pair (full dataset at: https://flathub.flatironinstitute.org/sapfire). In total, we identify… ▽ More Stellar streams record the accretion history of their host galaxy. We present a set of simulated streams from disrupted dwarf galaxies in 13 cosmological simulations of Milky Way (MW)-mass galaxies from the FIRE-2 suite at $z=0$, including 7 isolated Milky Way-mass systems and 6 hosts resembling the MW-M31 pair (full dataset at: https://flathub.flatironinstitute.org/sapfire). In total, we identify 106 simulated stellar streams, with no significant differences in the number of streams and masses of their progenitors between the isolated and paired environments. We resolve simulated streams with stellar masses ranging from $\sim 5\times10^5$ up to $\sim 10^{9} M_\odot$, similar to the mass range between the Orphan and Sagittarius streams in the MW. We confirm that present-day simulated satellite galaxies are good proxies for stellar stream progenitors, with similar properties including their stellar mass function, velocity dispersion, [Fe/H] and [$α$/H] evolution tracks, and orbital distribution with respect to the galactic disk plane. Each progenitor's lifetime is marked by several important timescales: its infall, star-formation quenching, and stream-formation times. We show that the ordering of these timescales is different between progenitors with stellar masses higher and lower than $\sim 2\times10^6 M_\odot$. Finally, we show that the main factor controlling the rate of phase-mixing, and therefore fading, of tidal streams from satellite galaxies in MW-mass hosts is non-adiabatic evolution of the host potential. Other factors commonly used to predict phase-mixing timescales, such as progenitor mass and orbital circularity, show virtually no correlation with the number of dynamical times required for a stream to become phase-mixed. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Comments: 26 pages, 19 figures, public dataset available at https://flathub.flatironinstitute.org/sapfire

Planes of satellites around Milky Way/M31-mass galaxies in the FIRE simulations and comparisons with the Local Group

Authors: Jenna Samuel, Andrew Wetzel, Sierra Chapman, Erik Tollerud, Philip F. Hopkins, Michael Boylan-Kolchin, Jeremy Bailin, Claude-André Faucher-Giguère

Abstract: We examine the prevalence, longevity, and causes of planes of satellite dwarf galaxies, as observed in the Local Group. We use 14 Milky Way/Andromeda-(MW/M31) mass host galaxies from the FIRE-2 simulations. We select the 14 most massive satellites by stellar mass within 300 kpc of each host and correct for incompleteness from the foreground galactic disc when comparing to the MW. We find that MW-l… ▽ More We examine the prevalence, longevity, and causes of planes of satellite dwarf galaxies, as observed in the Local Group. We use 14 Milky Way/Andromeda-(MW/M31) mass host galaxies from the FIRE-2 simulations. We select the 14 most massive satellites by stellar mass within 300 kpc of each host and correct for incompleteness from the foreground galactic disc when comparing to the MW. We find that MW-like planes as spatially thin and/or kinematically coherent as observed are uncommon, but they do exist in our simulations. Spatially thin planes occur in 1-2 per cent of snapshots during $z=0-0.2$, and kinematically coherent planes occur in 5 per cent of snapshots. These planes are generally transient, surviving for less than 500 Myr. However, if we select hosts with an LMC-like satellite near first pericentre, the fraction of snapshots with MW-like planes increases dramatically to 7-16 per cent, with lifetimes of 0.7-1 Gyr, likely because of group accretion of satellites. We find that M31's satellite distribution is much more common: M31's satellites lie within about 1 sigma of the simulation median for every plane metric we consider. We find no significant difference in average satellite planarity for isolated hosts versus hosts in LG-like pairs. Baryonic and dark matter-only simulations exhibit similar levels of planarity, even though baryonic subhaloes are less centrally concentrated within their host haloes. We conclude that planes of satellites are not a strong challenge to LCDM cosmology. △ Less

Submitted 20 July, 2021; v1 submitted 16 October, 2020; originally announced October 2020.

Comments: 15 pages, 8 figures. Main results in Figures 2 and 7. Accepted by MNRAS

Is NGC 300 a pure exponential disk galaxy?

Authors: In Sung Jang, Roelof S. de Jong, I. Minchev, Eric F. Bell, Antonela Monachesi, Benne W. Holwerda, Jeremy Bailin, Adam Smercina, Richard D'Souza

Abstract: 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… ▽ More 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. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 6 pages, 5 figures, 1 appendix, accepted for publication in A&A Letters

Tracing the anemic stellar halo of M101

Authors: In Sung Jang, Roelof S. de Jong, Benne W. Holwerda, Antonela Monachesi, Eric F. Bell, Jeremy Bailin

Abstract: 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.… ▽ More 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. △ Less

Submitted 31 January, 2020; originally announced January 2020.

Comments: 20 pages, 17 figures, accepted for publication in A&A

Journal ref: A&A 637, A8 (2020)

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

Authors: Isaiah B. Santistevan, Andrew Wetzel, Kareem El-Badry, Joss Bland-Hawthorn, Michael Boylan-Kolchin, Jeremy Bailin, Claude-Andre Faucher-Giguere, Samantha Benincasa

Abstract: 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 in… ▽ More 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. △ Less

Submitted 28 July, 2020; v1 submitted 9 January, 2020; originally announced January 2020.

Comments: 20 pages, 8 figures, accepted for publication in MNRAS

The Saga of M81: Global View of a Massive Stellar Halo in Formation

Authors: Adam Smercina, Eric F. Bell, Paul A. Price, Colin T. Slater, Richard D'Souza, Jeremy Bailin, Roelof S. de Jong, In Sung Jang, Antonela Monachesi, David Nidever

Abstract: Recent work has shown that Milky Way-mass galaxies display an incredible range of stellar halo properties, yet the origin of this diversity is unclear. The nearby galaxy M81 $-$ currently interacting with M82 and NGC 3077 $-$ sheds unique light on this problem. We present a Subaru Hyper Suprime-Cam survey of the resolved stellar populations around M81, revealing M81's stellar halo in never-before-… ▽ More Recent work has shown that Milky Way-mass galaxies display an incredible range of stellar halo properties, yet the origin of this diversity is unclear. The nearby galaxy M81 $-$ currently interacting with M82 and NGC 3077 $-$ sheds unique light on this problem. We present a Subaru Hyper Suprime-Cam survey of the resolved stellar populations around M81, revealing M81's stellar halo in never-before-seen detail. We resolve the halo to unprecedented $V$-band equivalent surface brightnesses of 33 mag arcsec$^{-2}$, and produce the first-ever global stellar mass density map for a Milky Way-mass stellar halo outside of the Local Group. Using the minor axis, we confirm M81's halo as one of the lowest mass and metal-poorest known ($M_{\star} \simeq 1.16{\times}10^9 M_{\odot}$, [Fe/H] $\simeq {-}1.2$) $-$ indicating a relatively quiet prior accretion history. Yet, our global halo census finds that tidally unbound material from M82 and NGC 3077 provides a substantial infusion of metal-rich material ($M_{\star} \simeq 5.4{\times}10^8$ $M_{\odot}$, [Fe/H] $\simeq {-}$0.9). We further show that, following the accretion of its massive satellite M82 (and the LMC-like NGC 3077), M81 will host one of the most massive and metal-rich stellar halos in the nearby universe. Thus, the saga of M81: following a passive history, M81's merger with M82 will completely transform its halo from a low-mass, anemic halo rivaling the MW, to a metal-rich behemoth rivaled only by systems such as M31. This dramatic transformation indicates that the observed diversity in stellar halo properties is primarily driven by diversity in the largest mergers these galaxies have experienced. △ Less

Submitted 26 October, 2020; v1 submitted 31 October, 2019; originally announced October 2019.

Comments: 25 pages, 15 figures, 3 tables, accepted to ApJ

Globular Cluster Intrinsic Iron Abundance Spreads: I. Catalog

Authors: Jeremy Bailin

Abstract: 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 resol… ▽ More 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. △ Less

Submitted 25 September, 2019; originally announced September 2019.

Comments: ApJS in press

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

Authors: Jenna Samuel, Andrew Wetzel, Erik Tollerud, Shea Garrison-Kimmel, Sarah Loebman, Kareem El-Badry, Philip F. Hopkins, Michael Boylan-Kolchin, Claude-Andre Faucher-Giguere, James S. Bullock, Samantha Benincasa, Jeremy Bailin

Abstract: 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… ▽ More 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. △ Less

Submitted 17 March, 2020; v1 submitted 25 April, 2019; originally announced April 2019.

Comments: 18 pages, 10 figures, plus appendices. Main results in figures 2, 3, and 4. Accepted version

Search for a metallicity spread in the multiple population Large Magellanic Cloud cluster NGC 1978

Authors: Andrés E. Piatti, Jeremy Bailin

Abstract: We report on the spread of [Fe/H] values in the massive Large Magellanic Cloud cluster NGC1978, recently confirmed to harbor multiple populations of nearly the same age. We used accurate Stromgren photometry of carefully selected cluster red giant branch stars along with a high-dispersion spectroscopy-based calibration of the metallicity-sensitive index m1. Once accounting for the photometry quali… ▽ More We report on the spread of [Fe/H] values in the massive Large Magellanic Cloud cluster NGC1978, recently confirmed to harbor multiple populations of nearly the same age. We used accurate Stromgren photometry of carefully selected cluster red giant branch stars along with a high-dispersion spectroscopy-based calibration of the metallicity-sensitive index m1. Once accounting for the photometry quality, assessed from extensive artificial stars tests to trace the photometric uncertainties as a function of the position to the cluster's center as well as the stellar brightness, and those from the metallicity calibration, we found that NGC1978 exhibits a small metallicity spread of 0.035 dex (+- 0.019-0.023), depending on whether stars with individual sigma[Fe/H] < 0.15 dex or those located in the cluster's outer areas are considered. Such a spread in [Fe/H] is consistent with a cluster formation model with self-enrichment, if mass loss higher than 90 per cent due to stellar evolutionary and galactic tidal effects is assumed. Nevertheless, scenarios where the apparent [Fe/H] variation reflects CN abundance anomalies or less extreme mass loss models with environmentally-dependent self enrichment should not be ruled out. △ Less

Submitted 3 December, 2018; originally announced December 2018.

Comments: 11 pages, 4 figures. Accepted for publication in The Astronomical Journal

A Model for Clumpy Self-Enrichment in Globular Clusters

Authors: Jeremy Bailin

Abstract: Detailed observations of globular clusters (GCs) have revealed evidence of self-enrichment: some of the heavy elements that we see in stars today were produced by cluster stars themselves. Moreover, GCs have internal subpopulations with different elemental abundances, including, in some cases, in elements such as iron that are produced by supernovae. This paper presents a theoretical model for GC… ▽ More Detailed observations of globular clusters (GCs) have revealed evidence of self-enrichment: some of the heavy elements that we see in stars today were produced by cluster stars themselves. Moreover, GCs have internal subpopulations with different elemental abundances, including, in some cases, in elements such as iron that are produced by supernovae. This paper presents a theoretical model for GC formation motivated by observations of Milky Way star forming regions and simulations of star formation, where giant molecular clouds fragment into multiple clumps which undergo star formation at slightly different times. Core collapse supernovae from earlier-forming clumps can enrich later-forming clumps to the degree that the ejecta can be retained within the gravitational potential well, resulting in subpopulations with different total metallicities once the clumps merge to form the final cluster. The model matches the mass-metallicity relation seen in GC populations around massive elliptical galaxies, and predicts metallicity spreads within clusters in excellent agreement with those seen in Milky Way GCs, even for those whose internal abundance spreads are so large that their entire identity as a GC is in question. The internal metallicity spread serves as an excellent measurement of how much self-enrichment has occurred in a cluster, a result that is very robust to variation in the model parameters. △ Less

Submitted 4 July, 2018; originally announced July 2018.

Comments: ApJ, in press. 18 pages. Code instantiating the model is at doi:10.5281/zenodo.1253347

Haloes at the ragged edge: The importance of the splashback radius

Authors: O. N. Snaith, J. Bailin, A. Knebe, G. Stinson, J. Wadsley, H. Couchman

Abstract: We have explored the outskirts of dark matter haloes out to 2.5 times the virial radius using a large sample of halos drawn from Illustris, along with a set of zoom simulations (MUGS). Using these, we make a systematic exploration of the shape profile beyond R$_{vir}$. In the mean sphericity profile of Illustris halos we identify a dip close to the virial radius, which is robust across a broad ran… ▽ More We have explored the outskirts of dark matter haloes out to 2.5 times the virial radius using a large sample of halos drawn from Illustris, along with a set of zoom simulations (MUGS). Using these, we make a systematic exploration of the shape profile beyond R$_{vir}$. In the mean sphericity profile of Illustris halos we identify a dip close to the virial radius, which is robust across a broad range of masses and infall rates. The inner edge of this feature may be related to the virial radius and the outer edge with the splashback radius. Due to the high halo-to-halo variation this result is visible only on average. However, in four individual halos in the MUGS sample, a decrease in the sphericity and a subsequent recovery is evident close to the splashback radius. We find that this feature persists for several Gyr, growing with the halo. This feature appears at the interface between the spherical halo density distribution and the filamentary structure in the environment. The shape feature is strongest when there is a high rate of infall, implying that the effect is due to the mixing of accreting and virializing material. The filamentary velocity field becomes rapidly mixed in the halo region inside the virial radius, with the area between this and the splashback radius serving as the transition region. We also identify a long-lasting and smoothly evolving splashback region in the radial density gradient in many of the MUGS halos. △ Less

Submitted 21 August, 2017; originally announced August 2017.

Comments: 22 pages, 15 figures, Accepted for publication in MNRAS

d1005+68: A New Faint Dwarf Galaxy in the M81 Group

Authors: Adam Smercina, Eric F. Bell, Colin T. Slater, Paul A. Price, Jeremy Bailin, Antonela Monachesi

Abstract: We present the discovery of d1005+68, a new faint dwarf galaxy in the M81 Group, using observations taken with the Subaru Hyper Suprime-Cam. d1005+68's color-magnitude diagram is consistent with a distance of $3.98_{-0.43}^{+0.39}$ Mpc, establishing group membership. We derive an absolute $V$-band magnitude, from stellar isochrone fitting, of $M_{V} = -7.94_{-0.50}^{+0.38}$, with a half-light radi… ▽ More We present the discovery of d1005+68, a new faint dwarf galaxy in the M81 Group, using observations taken with the Subaru Hyper Suprime-Cam. d1005+68's color-magnitude diagram is consistent with a distance of $3.98_{-0.43}^{+0.39}$ Mpc, establishing group membership. We derive an absolute $V$-band magnitude, from stellar isochrone fitting, of $M_{V} = -7.94_{-0.50}^{+0.38}$, with a half-light radius of $r_{h} = 188_{-41}^{+39}$ pc. These place d1005+68 within the radius-luminosity locus of Local Group and M81 satellites and among the faintest confirmed satellites outside the Local Group. Assuming an age of 12 Gyr, d1005+68's red giant branch is best fit by an isochrone of [Fe/H] $= -1.90 \pm 0.24$. It has a projected separation from nearby M81 satellite BK5N of only 5 kpc. As this is well within BK5N's virial radius, we speculate that d1005+68 may be a satellite of BK5N. If confirmed, this would make d1005+68 one of the first detected satellites-of-a-satellite. △ Less

Submitted 22 June, 2017; v1 submitted 21 June, 2017; originally announced June 2017.

Comments: 7 pages, 4 figures, 1 table, additional affiliations included

Galaxies grow their bulges and black holes in diverse ways

Authors: Eric F. Bell, Antonela Monachesi, Benjamin Harmsen, Roelof S. de Jong, Jeremy Bailin, David J. Radburn-Smith, Richard D'Souza, Benne W. Holwerda

Abstract: Galaxies with Milky Way-like stellar masses have a wide range of bulge and black hole masses; in turn, these correlate with other properties such as star formation history. While many processes may drive bulge formation, major and minor mergers are expected to play a crucial role. Stellar halos offer a novel and robust measurement of galactic merger history; cosmologically-motivated models predict… ▽ More Galaxies with Milky Way-like stellar masses have a wide range of bulge and black hole masses; in turn, these correlate with other properties such as star formation history. While many processes may drive bulge formation, major and minor mergers are expected to play a crucial role. Stellar halos offer a novel and robust measurement of galactic merger history; cosmologically-motivated models predict that mergers with larger satellites produce more massive, higher metallicity stellar halos, reproducing the recently-observed stellar halo metallicity-mass relation. We quantify the relationship between stellar halo mass and bulge or black hole prominence using a sample of eighteen Milky Way-mass galaxies with newly-available measurements of (or limits on) stellar halo properties. There is an order of magnitude range in bulge mass, and two orders of magnitude in black hole mass, at a given stellar halo mass (or, equivalently, merger history). Galaxies with low mass bulges show a wide range of quiet merger histories, implying formation mechanisms that do not require intense merging activity. Galaxies with massive 'classical' bulges and central black holes also show a wide range of merger histories. While three of these galaxies have massive stellar halos consistent with a merger origin, two do not - merging appears to have had little impact in making these two massive 'classical' bulges. Such galaxies may be ideal laboratories to study massive bulge formation through pathways such as early gas-rich accretion, violent disk instabilities or misaligned infall of gas throughout cosmic time. △ Less

Submitted 20 February, 2017; originally announced February 2017.

Comments: Astrophysical Journal Letters, in press. 6 pages, 3 figures, 1 table

Globular Cluster Systems in Brightest Cluster Galaxies. III: Beyond Bimodality

Authors: W. E. Harris, S. M. Ciccone, G. M. Eadie, O. Y. Gnedin, D. Geisler, B. Rothberg, J. Bailin

Abstract: We present new deep photometry of the rich globular cluster (GC) systems around the Brightest Cluster Galaxies UGC 9799 (Abell 2052) and UGC 10143 (Abell 2147), obtained with the HST ACS and WFC3 cameras. For comparison, we also present new reductions of similar HST/ACS data for the Coma supergiants NGC 4874 and 4889. All four of these galaxies have huge cluster populations (to the radial limits o… ▽ More We present new deep photometry of the rich globular cluster (GC) systems around the Brightest Cluster Galaxies UGC 9799 (Abell 2052) and UGC 10143 (Abell 2147), obtained with the HST ACS and WFC3 cameras. For comparison, we also present new reductions of similar HST/ACS data for the Coma supergiants NGC 4874 and 4889. All four of these galaxies have huge cluster populations (to the radial limits of our data, comprising from 12000 to 23000 clusters per galaxy). The metallicity distribution functions (MDFs) of the GCs can still be matched by a bimodal-Gaussian form where the metal-rich and metal-poor modes are separated by ~0.8 dex, but the internal dispersions of each mode are so large that the total MDF becomes very broad and nearly continuous from [Fe/H] = -2.4 to Solar. There are, however, significant differences between galaxies in the relative numbers of \emph{metal-rich} clusters, suggesting that they underwent significantly different histories of mergers with massive, gas-rich halos. Lastly, the proportion of metal-poor GCs rises especially rapidly outside projected radii R > 4 R_eff, suggesting the importance of accreted dwarf satellites in the outer halo. Comprehensive models for the formation of GCs as part of the hierarchical formation of their parent galaxies will be needed to trace the systematic change in structure of the MDF with galaxy mass, from the distinctly bimodal form in smaller galaxies up to the broad continuum that we see in the very largest systems. △ Less

Submitted 21 December, 2016; originally announced December 2016.

Comments: In press for Astrophysical Journal

Characterizing dw1335-29, a recently discovered dwarf satellite of M83

Authors: Andreia Carrillo, Eric Bell, Jeremy Bailin, Antonela Monachesi, Roelof S. de Jong, Benjamin Harmsen, Colin Slater

Abstract: The number, distribution, and properties of dwarf satellites are crucial probes of the physics of galaxy formation at low masses and the response of satellite galaxies to the tidal and gas dynamical effects of their more massive parent.To make progress, it is necessary to augment and solidify the census of dwarf satellites of galaxies outside the Local Group. Müller et al. (2015) presented 16 dwar… ▽ More The number, distribution, and properties of dwarf satellites are crucial probes of the physics of galaxy formation at low masses and the response of satellite galaxies to the tidal and gas dynamical effects of their more massive parent.To make progress, it is necessary to augment and solidify the census of dwarf satellites of galaxies outside the Local Group. Müller et al. (2015) presented 16 dwarf galaxy candidates near M83, but lacking reliable distances, it is unclear which candidates are M83 satellites. Using red giant branch stars from the HST/GHOSTS survey in conjunction with ground-based images from VLT/VIMOS, we confirm that one of the candidates, dw1335-29-- with a projected distance of 26 kpc from M83 and a distance modulus of $(m - M)_0 = 28.5^{+0.3}_{-0.1}$ -- is a satellite of M83. We estimate an absolute magnitude $M_V = -10.1 \pm{0.4}$, an ellipticity of $0.40^{+0.14}_{-0.22}$, a half light radius of $656^{+121}_{-170 }$ pc, and [Fe/H] = $-1.3^{+0.3}_{-0.4}$. Owing to dw1335-29's somewhat irregular shape and possible young stars, we classify this galaxy as a dwarf irregular or transition dwarf. This is curious, as with a projected distance of 26 kpc from M83, dw1335-29 is expected to lack recent star formation. Further study of M83's dwarf population will reveal if star formation in its satellites is commonplace (suggesting a lack of a hot gas envelope for M83 that would quench star formation) or rare (suggesting that dw1335-29 has a larger M83-centric distance, and is fortuitously projected to small radii). △ Less

Submitted 30 November, 2016; originally announced November 2016.

Comments: 7 pages, 5 figures, accepted for publication in MNRAS

Diverse Stellar Haloes in Nearby Milky Way-Mass Disc Galaxies

Authors: Benjamin Harmsen, Antonela Monachesi, Eric F. Bell, Roelof S. de Jong, Jeremy Bailin, David J. Radburn-Smith, Benne W. Holwerda

Abstract: We have examined the resolved stellar populations at large galactocentric distances along the minor axis (from 10 kpc up to between 40 and 75 kpc), with limited major axis coverage, of six nearby highly-inclined Milky Way-mass disc galaxies using HST data from the GHOSTS survey. We select red giant branch stars to derive stellar halo density profiles. The projected minor axis density profiles can… ▽ More We have examined the resolved stellar populations at large galactocentric distances along the minor axis (from 10 kpc up to between 40 and 75 kpc), with limited major axis coverage, of six nearby highly-inclined Milky Way-mass disc galaxies using HST data from the GHOSTS survey. We select red giant branch stars to derive stellar halo density profiles. The projected minor axis density profiles can be approximated by power laws with projected slopes of between $-2$ and $-3.7$ and a diversity of stellar halo masses of $1-6\times 10^{9}M_{\odot}$, or $2-14\%$ of the total galaxy stellar masses. The typical intrinsic scatter around a smooth power law fit is $0.05-0.1$ dex owing to substructure. By comparing the minor and major axis profiles, we infer projected axis ratios $c/a$ at $\sim 25$ kpc between $0.4-0.75$. The GHOSTS stellar haloes are diverse, lying between the extremes charted out by the (rather atypical) haloes of the Milky Way and M31. We find a strong correlation between the stellar halo metallicities and the stellar halo masses. We compare our results with cosmological models, finding good agreement between our observations and accretion-only models where the stellar haloes are formed by the disruption of dwarf satellites. In particular, the strong observed correlation between stellar halo metallicity and mass is naturally reproduced. Low-resolution hydrodynamical models have unrealistically high stellar halo masses. Current high-resolution hydrodynamical models appear to predict stellar halo masses somewhat higher than observed but with reasonable metallicities, metallicity gradients and density profiles. △ Less

Submitted 16 November, 2016; originally announced November 2016.

Comments: 26 pages, 17 figures. Accepted for publication in MNRAS

HI4PI: A full-sky HI survey based on EBHIS and GASS

Authors: HI4PI Collaboration, N. Ben Bekhti, L. Flöer, R. Keller, J. Kerp, D. Lenz, B. Winkel, J. Bailin, M. R. Calabretta, L. Dedes, H. A. Ford, B. K. Gibson, U. Haud, S. Janowiecki, P. M. W. Kalberla, F. J. Lockman, N. M. McClure-Griffiths, T. Murphy, H. Nakanishi, D. J. Pisano, L. Staveley-Smith

Abstract: Measurement of the Galactic neutral atomic hydrogen (HI) column density, NHI, and brightness temperatures, Tb, is of high scientific value for a broad range of astrophysical disciplines. In the past two decades, one of the most-used legacy HI datasets has been the Leiden/Argentine/Bonn Survey (LAB). We release the HI 4$π$ survey (HI4PI), an all-sky database of Galactic HI, which supersedes the LAB… ▽ More Measurement of the Galactic neutral atomic hydrogen (HI) column density, NHI, and brightness temperatures, Tb, is of high scientific value for a broad range of astrophysical disciplines. In the past two decades, one of the most-used legacy HI datasets has been the Leiden/Argentine/Bonn Survey (LAB). We release the HI 4$π$ survey (HI4PI), an all-sky database of Galactic HI, which supersedes the LAB survey. The HI4PI survey is based on data from the recently completed first coverage of the Effelsberg-Bonn HI Survey (EBHIS) and from the third revision of the Galactic All-Sky Survey (GASS). EBHIS and GASS share similar angular resolution and match well in sensitivity. Combined, they are ideally suited to be a successor to LAB. The new HI4PI survey outperforms the LAB in angular resolution (16.2', FWHM) and sensitivity (RMS: 43 mK). Moreover, it has full spatial sampling and thus overcomes a major drawback of LAB, which severely undersamples the sky. We publish all-sky column density maps of the neutral atomic hydrogen in the Milky Way, along with full spectroscopic data, in several map projections including HEALPix. △ Less

Submitted 19 October, 2016; originally announced October 2016.

Comments: 14 pages, 8 figures, A&A in press, corresponding authors: B. Winkel and P. M. W. Kalberla

The impact of baryons on the direct detection of dark matter

Authors: Chris Kelso, Christopher Savage, Monica Valluri, Katherine Freese, Gregory S. Stinson, Jeremy Bailin

Abstract: 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 complic… ▽ More 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. △ Less

Submitted 7 November, 2016; v1 submitted 18 January, 2016; originally announced January 2016.

Comments: Minor changes to match JCAP version. 21 pages, 9 figures

Report number: NORDITA-2015-15, CETUP2015-030

Journal ref: JCAP 1608 (2016) 071

Origin of the Metallicity Distribution in the Thick Disc

Authors: Maider S. Miranda, Kate Pilkington, Brad K. Gibson, Christopher B. Brook, Patricia Sánchez-Blázquez, Ivan Minchev, Christopher Gareth Few, Rory Smith, Rosa Domínguez-Tenreiro, Aura Obreja, Jeremy Bailin, Greg S. Stinson

Abstract: 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 simulat… ▽ More 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. △ Less

Submitted 14 December, 2015; originally announced December 2015.

Comments: 15 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 587, A10 (2016)

The history of stellar metallicity in a simulated disc galaxy

Authors: O. N. Snaith, J. Bailin, B. K. Gibson, E. F. Bell, G. Stinson, M. Valluri, J. Wadsley, H. Couchman

Abstract: We explore the chemical distribution of stars in a simulated galaxy. Using simulations of the same initial conditions but with two different feedback schemes (MUGS and MaGICC), we examine the features of the age-metallicity relation (AMR), and the three-dimensional age-metallicity-[O/Fe] distribution, both for the galaxy as a whole and decomposed into disc, bulge, halo, and satellites. The MUGS si… ▽ More We explore the chemical distribution of stars in a simulated galaxy. Using simulations of the same initial conditions but with two different feedback schemes (MUGS and MaGICC), we examine the features of the age-metallicity relation (AMR), and the three-dimensional age-metallicity-[O/Fe] distribution, both for the galaxy as a whole and decomposed into disc, bulge, halo, and satellites. The MUGS simulation, which uses traditional supernova feedback, is replete with chemical substructure. This sub- structure is absent from the MaGICC simulation, which includes early feedback from stellar winds, a modified IMF and more efficient feedback. The reduced amount of substructure is due to the almost complete lack of satellites in MaGICC. We identify a significant separation between the bulge and disc AMRs, where the bulge is considerably more metal-rich with a smaller spread in metallicity at any given time than the disc. Our results suggest, however, that identifying the substructure in observations will require exquisite age resolution, on the order of 0.25 Gyr. Certain satellites show exotic features in the AMR, even forming a 'sawtooth' shape of increasing metallicity followed by sharp declines which correspond to pericentric passages. This fact, along with the large spread in stellar age at a given metallicity, compromises the use of metallicity as an age indicator, although alpha abundance provides a more robust clock at early times. This may also impact algorithms that are used to reconstruct star formation histories from resolved stellar populations, which frequently assume a monotonically-increasing AMR. △ Less

Submitted 8 December, 2015; originally announced December 2015.

Comments: 25 pages, 14 Figures, 2 Tables. Accepted for publication in MNRAS

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

Authors: Monica Valluri, Sarah R. Loebman, Jeremy Bailin, Adam Clarke, Victor P. Debattista, Greg Stinson

Abstract: 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 sim… ▽ More 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. △ Less

Submitted 18 October, 2015; originally announced October 2015.

Comments: 2 pages, 1 figure, to appear in: The General Assembly of Stellar Halos: Structure, Origin and Evolution Proceedings IAU Symposium No. 317, 2015, M. Arnaboldi, A. Bragaglia, M. Rejkuba & D. Romano eds

Extragalactic archeology with the GHOSTS Survey I. - Age-resolved disk structure of nearby low-mass galaxies

Authors: David Streich, Roelof S. de Jong, Jeremy Bailin, Eric F. Bell, Benne W. Holwerda, Ivan Minchev, Antonela Monachesi, David J. Radburn-Smith

Abstract: We study the individual evolution histories of three nearby low-mass edge-on galaxies (IC 5052, NGC4244, and NGC5023). Using resolved stellar populations, we constructed star count density maps for populations of different ages and analyzed the change of structural parameters with stellar age within each galaxy. We do not detect a separate thick disk in any of the three galaxies, even though our… ▽ More We study the individual evolution histories of three nearby low-mass edge-on galaxies (IC 5052, NGC4244, and NGC5023). Using resolved stellar populations, we constructed star count density maps for populations of different ages and analyzed the change of structural parameters with stellar age within each galaxy. We do not detect a separate thick disk in any of the three galaxies, even though our observations cover a wider range in equivalent surface brightness than any integrated light study. While scale heights increase with age, each population can be well described by a single disk. Two of the galaxies contain a very weak additional component, which we identify as the faint halo. The mass of these faint halos is lower than 1% of the mass of the disk. The three galaxies show low vertical heating rates, which are much lower than the heating rate of the Milky Way. This indicates that heating agents, such as giant molecular clouds and spiral structure, are weak in low-mass galaxies. All populations in the three galaxies exhibit no or only little flaring. While this finding is consistent with previous integrated light studies, it poses strong constraints on galaxy simulations, where strong flaring is often found as a result of interactions or radial migration. △ Less

Submitted 22 September, 2015; originally announced September 2015.

Comments: Accepted for publication in A&A

Journal ref: A&A 585, A97 (2016)

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

Authors: Antonela Monachesi, Eric F. Bell, David Radburn-Smith, Jeremy Bailin, Roelof S. de Jong, Benne Holwerda, David Streich, Grace Silverstein

Abstract: 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… ▽ More 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. △ Less

Submitted 19 January, 2016; v1 submitted 23 July, 2015; originally announced July 2015.

Comments: Accepted by MNRAS, 31 pages, 18 figures. Appendix added and some editions to match accepted version. Conclusions unchanged

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

Authors: William E. Harris, Warren Morningstar, Oleg Y. Gnedin, Heather O'Halloran, John P. Blakeslee, Bradley C. Whitmore, Patrick Cote, Douglas Geisler, Eric W. Peng, Jeremy Bailin, Barry Rothberg, Robert Cockcroft, Regina Barber DeGraaff

Abstract: 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, w… ▽ More 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. △ Less

Submitted 23 October, 2014; originally announced October 2014.

Comments: To appear in Astrophysical Journal, December 2014

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

Authors: John M. Cannon, Megan Johnson, Kristen B. W. McQuinn, Erik D. Alfvin, Jeremy Bailin, H. Alyson Ford, Léo Girardi, Alec S. Hirschauer, Steven Janowiecki, John J. Salzer, Angela Van Sistine, Andrew Dolphin, E. C. Elson, Baerbel Koribalski, Paola Marigo, Jessica L. Rosenberg, Philip Rosenfield, Evan D. Skillman, Aparna Venkatesan, Steven R. Warren

Abstract: 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 relat… ▽ More 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. △ Less

Submitted 14 April, 2014; originally announced April 2014.

Comments: Accepted for publication in the Astrophysical Journal Letters

The Role of Feedback in Shaping the Structure of the Interstellar Medium

Authors: A. P. Walker, B. K. Gibson, K. Pilkington, C. B. Brook, P. Dutta, S. Stanimirovic, G. S. Stinson, J. Bailin

Abstract: We present an analysis of the role of feedback in shaping the neutral hydrogen (HI) content of simulated disc galaxies. For our analysis, we have used two realisations of two separate Milky Way-like (~L*) discs - one employing a conservative feedback scheme (MUGS), the other significantly more energetic (MaGICC). To quantify the impact of these schemes, we generate zeroth moment (surface density)… ▽ More We present an analysis of the role of feedback in shaping the neutral hydrogen (HI) content of simulated disc galaxies. For our analysis, we have used two realisations of two separate Milky Way-like (~L*) discs - one employing a conservative feedback scheme (MUGS), the other significantly more energetic (MaGICC). To quantify the impact of these schemes, we generate zeroth moment (surface density) maps of the inferred HI distribution; construct power spectra associated with the underlying structure of the simulated cold ISM, in addition to their radial surface density and velocity dispersion profiles. Our results are compared with a parallel, self-consistent, analysis of empirical data from THINGS (The HI Nearby Galaxy Survey). Single power-law fits (P~k^gamma) to the power spectra of the stronger-feedback (MaGICC) runs (over spatial scales corresponding to 0.5 kpc to 20 kpc) result in slopes consistent with those seen in the THINGS sample (gamma = -2.5). The weaker-feedback (MUGS) runs exhibit shallower power law slopes (gamma = -1.2). The power spectra of the MaGICC simulations are more consistent though with a two-component fit, with a flatter distribution of power on larger scales (i.e., gamma = -1.4 for scales in excess of 2 kpc) and a steeper slope on scales below 1 kpc (gamma = -5), qualitatively consistent with empirical claims, as well as our earlier work on dwarf discs. The radial HI surface density profiles of the MaGICC discs show a clear exponential behaviour, while those of the MUGS suite are essentially flat; both behaviours are encountered in nature, although the THINGS sample is more consistent with our stronger (MaGICC) feedback runs. △ Less

Submitted 3 March, 2014; originally announced March 2014.

Comments: Accepted for publication in MNRAS; 9 pages; 7 figures

Systematic Problems With Using Dark Matter Simulations to Model Stellar Halos

Authors: Jeremy Bailin, Eric F. Bell, Monica Valluri, Greg S. Stinson, Victor P. Debattista, H. M. P. Couchman, James Wadsley

Abstract: 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 assu… ▽ More 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. △ Less

Submitted 21 January, 2014; originally announced January 2014.

Comments: ApJ, in press

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

Authors: David Streich, Roelof S. de Jong, Jeremy Bailin, Paul Goudfrooij, David Radburn-Smith, Marija Vlajic

Abstract: 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… ▽ More 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. △ Less

Submitted 13 January, 2014; originally announced January 2014.

Comments: accepted for publication in A&A

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

Authors: David J. Radburn-Smith, Roelof S. de Jong, David Streich, Eric F. Bell, Julianne J. Dalcanton, Andrew E. Dolphin, Adrienne M. Stilp, Antonela Monachesi, Benne W. Holwerda, Jeremy Bailin

Abstract: 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^… ▽ More 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. △ Less

Submitted 6 December, 2013; originally announced December 2013.

Comments: 12 pages, 9 figures. Accepted for publication in ApJ

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

Authors: Antonela Monachesi, Eric F. Bell, David J. Radburn-Smith, Roelof S. de Jong, Jeremy Bailin, Julianne J. Dalcanton, Benne W. Holwerda, H. Alyson Ford, David Streich, Marija Vlajic, Daniel B. Zucker

Abstract: 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… ▽ More 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. △ Less

Submitted 2 December, 2013; originally announced December 2013.

Comments: 8 pages, 7 figures. Accepted for publication in ApJ

Constraining Sub-Grid Physics with High-Redshift Spatially-Resolved Metallicity Distributions

Authors: B. K. Gibson, K. Pilkington, C. B. Brook, G. S. Stinson, J. Bailin

Abstract: We examine the role of energy feedback in shaping the distribution of metals within cosmological hydrodynamical simulations of L* disc galaxies. While negative abundance gradients today provide a boundary condition for galaxy evolution models, in support of inside-out disc growth, empirical evidence as to whether abundance gradients steepen or flatten with time remains highly contradictory. We mad… ▽ More We examine the role of energy feedback in shaping the distribution of metals within cosmological hydrodynamical simulations of L* disc galaxies. While negative abundance gradients today provide a boundary condition for galaxy evolution models, in support of inside-out disc growth, empirical evidence as to whether abundance gradients steepen or flatten with time remains highly contradictory. We made use of a suite of L* discs, realised with and without `enhanced' feedback. All the simulations were produced using the smoothed particle hydrodynamics code Gasoline, and their in situ gas-phase metallicity gradients traced from redshift z~2 to the present-day. Present-day age-metallicity relations and metallicity distribution functions were derived for each system. The `enhanced' feedback models, which have been shown to be in agreement with a broad range of empirical scaling relations, distribute energy and re-cycled ISM material over large scales and predict the existence of relatively `flat' and temporally invariant abundance gradients. Enhanced feedback schemes reduce significantly the scatter in the local stellar age-metallicity relation and, especially, the [O/Fe]-[Fe/H] relation. The local [O/Fe] distribution functions for our L* discs show clear bimodality, with peaks at [O/Fe]=-0.05 and +0.05 (for stars with [Fe/H]>-1), consistent with our earlier work on dwarf discs. Our results with `enhanced' feedback are inconsistent with our earlier generation of simulations realised with `conservative' feedback. We conclude that spatially-resolved metallicity distributions, particularly at high-redshift, offer a unique and under-utilised constraint on the uncertain nature of stellar feedback processes. △ Less

Submitted 14 May, 2013; v1 submitted 10 April, 2013; originally announced April 2013.

Comments: Accepted for publication in A&A; 5 pages; 3 figures; inset/legend to figure fixed (misaligned symbol fixed!)

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

Authors: Antonela Monachesi, Eric F. Bell, David Radburn-Smith, Marija Vlajic, Roelof S. de Jong, Jeremy Bailin, Julianne J. Dalcanton, Benne W. Holwerda, David Streich

Abstract: 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 ma… ▽ More 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. △ Less

Submitted 11 February, 2013; originally announced February 2013.

Comments: 16 pages, 7 figures. Accepted to ApJ

Halo orbits in cosmological disk galaxies: tracers of formation history

Authors: Monica Valluri, Victor P. Debattista, Gregory S. Stinson, Jeremy Bailin, Thomas R. Quinn, H. M. P. Couchman, James Wadsley

Abstract: 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 clos… ▽ More 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. △ Less

Submitted 28 February, 2013; v1 submitted 18 January, 2013; originally announced January 2013.

Comments: Revised version accepted by ApJ, Feb 28th 2013. Minor changes, not affecting conclusions

The luminosity function of diverse satellite galaxy systems

Authors: S. Nickerson, G. Stinson, H. M. P. Couchman, J. Bailin, J. Wadsley

Abstract: The high-resolution, SPH galaxies of the McMaster Unbiased Galaxy Survey (MUGS) are used to examine the satellite systems of sixteen model host galaxies. Each galaxy has a different mass, angular momentum and merger history that yield a rich set of satellite luminosity functions. With new observations of distant satellite systems, we can compare these luminosity functions to satellite systems beyo… ▽ More The high-resolution, SPH galaxies of the McMaster Unbiased Galaxy Survey (MUGS) are used to examine the satellite systems of sixteen model host galaxies. Each galaxy has a different mass, angular momentum and merger history that yield a rich set of satellite luminosity functions. With new observations of distant satellite systems, we can compare these luminosity functions to satellite systems beyond the Local Group. We find that the luminosity functions of our simulations compare well to observations when the luminosity functions are scaled according to host mass. We use the recently-found relationship between dwarf satellites and host mass in distant satellite systems (Trentham & Tully 2009) to normalize a theoretical, complete luminosity function for the Milky Way (Koposov et al. 2008). The luminosity function of satellites, expressed as a function of the host mass, is given by dN/dM_V= 3.5 x M_host^(0.91) x 10^(0.1 x M_V-10.2), where host mass is given in M_sun. The mass of a host galaxy can be used to predict the number of dwarf satellites and even when considering spiral and elliptical hosts separately this relation holds. △ Less

Submitted 8 November, 2012; originally announced November 2012.

Comments: 7 pages, 5 figures, 1 table, accepted by MNRAS

The Distribution of Metals in Cosmological Hydrodynamical Simulations of Dwarf Disk Galaxies

Authors: K. Pilkington, B. K. Gibson, C. B. Brook, F. Calura, G. S. Stinson, R. J. Thacker, L. Michel-Dansac, J. Bailin, H. M. P. Couchman, J. Wadsley, T. R. Quinn, A. Maccio

Abstract: We examine the chemical properties of 5 cosmological hydrodynamical simulations of an M33-like disc galaxy which have been shown to be consistent with the morphological characteristics and bulk scaling relations expected of late-type spirals. These simulations are part of the Making Galaxies In a Cosmological Context (MaGICC) Project, in which stellar feedback is tuned to match the stellar mass --… ▽ More We examine the chemical properties of 5 cosmological hydrodynamical simulations of an M33-like disc galaxy which have been shown to be consistent with the morphological characteristics and bulk scaling relations expected of late-type spirals. These simulations are part of the Making Galaxies In a Cosmological Context (MaGICC) Project, in which stellar feedback is tuned to match the stellar mass -- halo mass relationship. Each realisation employed identical initial conditions and assembly histories, but differed from one another in their underlying baryonic physics prescriptions, including (a) the efficiency with which each supernova energy couples to the ISM, (b) the impact of feedback associated with massive star radiation pressure, (c) the role of the minimum shut-off time for radiative cooling of Type II SNe remnants, (d) the treatment of metal diffusion, and (e) varying the IMF. Our analysis focusses on the resulting stellar metallicity distribution functions (MDFs) in each simulated (analogous) `solar neighbourhood' and central `bulge' region. We compare the simulated MDFs' skewness, kurtosis, and dispersion (inter-quartile, inter-decile, inter-centile, and inter-tenth-percentile regions) with that of the empirical solar neighbourhood MDF and Local Group dwarfs. We find that the MDFs of the simulated discs are more negatively skewed, with higher kurtosis, than those observed locally. We can trace this difference to the simulations' tight and correlated age-metallicity relations (compared with that of the Milky Way), suggesting that these relations within `dwarf' discs might be steeper than in L* discs and/or the degree of stellar orbital re-distribution and migration inferred locally has not been captured in their entirety, at the resolution of our simulations. The important role of metal diffusion in ameliorating the over-production of extremely metal-poor stars is highlighted. △ Less

Submitted 21 May, 2012; originally announced May 2012.

Comments: Accepted for publication in MNRAS; 12 pages; 7 figures