Showing 1–50 of 50 results for author: Loebman, S

The imprint of dark matter on the Galactic acceleration field

Authors: Arpit Arora, Robyn E. Sanderson, Sukanya Chakrabarti, Andrew Wetzel, Thomas Donlon II, Danny Horta, Sarah R. Loebman, Lina Necib, Micah Oeur

Abstract: Measurements of the accelerations of stars enabled by time-series extreme-precision spectroscopic observations, from pulsar timing, and from eclipsing binary stars in the Solar Neighborhood offer insights into the mass distribution of the Milky Way that do not rely on traditional equilibrium modeling. Given the measured accelerations, we can determine a total mass density, and from this, by accoun… ▽ More Measurements of the accelerations of stars enabled by time-series extreme-precision spectroscopic observations, from pulsar timing, and from eclipsing binary stars in the Solar Neighborhood offer insights into the mass distribution of the Milky Way that do not rely on traditional equilibrium modeling. Given the measured accelerations, we can determine a total mass density, and from this, by accounting for the mass in stars, gas, and dust, we can infer the amount of dark matter. Leveraging the FIRE-2 simulations of Milky Way-mass galaxies, we compare vertical acceleration profiles between cold dark matter (CDM) and self-interacting dark matter (SIDM) with constant cross-section of 1 cm$^2$ g$^{-1}$ across three halos with diverse assembly histories. Notably, significant asymmetries in vertical acceleration profiles near the midplane at fixed radii are observed in both CDM and SIDM, particularly in halos recently affected by mergers with satellites of Sagittarius/SMC-like masses or greater. These asymmetries offer a unique window into exploring the merger history of a galaxy. We show that SIDM halos consistently exhibit higher local stellar and dark matter densities and steeper vertical acceleration gradients, up to 30% steeper near the Solar Neighborhood. SIDM halos also manifest a more oblate halo shape in the Solar Neighborhood. Furthermore, enhanced precision in acceleration measurements and larger datasets promise to provide better constraints on the local dark matter density, complementing our understanding from kinematic analysis of their distribution within galaxies. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 14 pages, 5 figures, and 2 tables. Submitted to APJ

Spectacular nucleosynthesis from early massive stars

Authors: Alexander P. Ji, Sanjana Curtis, Nicholas Storm, Vedant Chandra, Kevin C. Schlaufman, Keivan G. Stassun, Alexander Heger, Marco Pignatari, Adrian M. Price-Whelan, Maria Bergemann, Guy S. Stringfellow, Carla Frohlich, Henrique Reggiani, Erika M. Holmbeck, Jamie Tayar, Shivani P. Shah, Emily J. Griffith, Chervin F. P. Laporte, Andrew R. Casey, Keith Hawkins, Danny Horta, William Cerny, Pierre Thibodeaux, Sam A. Usman, Joao A. S. Amarante , et al. (17 additional authors not shown)

Abstract: Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931… ▽ More Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 11 pages + 22 page appendix, accepted to ApJL

Disk settling and dynamical heating: histories of Milky Way-mass stellar disks across cosmic time in the FIRE simulations

Authors: Fiona McCluskey, Andrew Wetzel, Sarah R. Loebman, Jorge Moreno, Claude-Andre Faucher-Giguere, Philip F. Hopkins

Abstract: We study the kinematics of stars both at their formation and today within 14 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. We quantify the relative importance of cosmological disk settling and post-formation dynamical heating. We identify three eras: a Pre-Disk Era (typically >8 Gyr ago), when stars formed on dispersion-dominated orbits; an Early-Disk Era (~ 8 - 4… ▽ More We study the kinematics of stars both at their formation and today within 14 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. We quantify the relative importance of cosmological disk settling and post-formation dynamical heating. We identify three eras: a Pre-Disk Era (typically >8 Gyr ago), when stars formed on dispersion-dominated orbits; an Early-Disk Era (~ 8 - 4 Gyr ago), when stars started to form on rotation-dominated orbits but with high velocity dispersion, sigma_form; and a Late-Disk Era (< 4 Gyr ago), when stars formed with low sigma_form. sigma_form increased with time during the Pre-Disk Era, peaking ~ 8 Gyr ago, then decreased throughout the Early-Disk Era as the disk settled and remained low throughout the Late-Disk Era. By contrast, the velocity dispersion measured today, sigma_now, increases monotonically with age because of stronger post-formation heating for Pre-Disk stars. Importantly, most of sigma_now was in place at formation, not added post-formation, for stars younger than ~ 10 Gyr. We compare the evolution of the three velocity components: at all times, sigma_R,form > sigma_phi,form > sigma_Z,form. Post-formation heating primarily increased sigma_R at ages < 4 Gyr but acted nearly isotropically for older stars. The kinematics of young stars in FIRE-2 broadly agree with the range observed across the MW, M31, M33, and PHANGS-MUSE galaxies. The lookback time that the disk began to settle correlates with its dynamical state today: earlier-settling galaxies currently form colder disks. Including stellar cosmic-ray feedback does not significantly change disk rotational support at fixed stellar mass. △ Less

Submitted 21 November, 2023; v1 submitted 24 March, 2023; originally announced March 2023.

Comments: 23 pages, accepted for publication in MNRAS

What Causes The Formation of Disks and End of Bursty Star Formation?

Authors: Philip F. Hopkins, Alexander B. Gurvich, Xuejian Shen, Zachary Hafen, Michael Y. Grudic, Shalini Kurinchi-Vendhan, Christopher C. Hayward, Fangzhou Jiang, Matthew E. Orr, Andrew Wetzel, Dusan Keres, Jonathan Stern, Claude-Andre Faucher-Giguere, James Bullock, Coral Wheeler, Kareem El-Badry, Sarah R. Loebman, Jorge Moreno, Michael Boylan-Kolchin, Eliot Quataert

Abstract: As they grow, galaxies can transition from irregular/spheroidal with 'bursty' star formation histories (SFHs), to disky with smooth SFHs. But even in simulations, the direct physical cause of such transitions remains unclear. We therefore explore this in a large suite of numerical experiments re-running portions of cosmological simulations with widely varied physics, further validated with existin… ▽ More As they grow, galaxies can transition from irregular/spheroidal with 'bursty' star formation histories (SFHs), to disky with smooth SFHs. But even in simulations, the direct physical cause of such transitions remains unclear. We therefore explore this in a large suite of numerical experiments re-running portions of cosmological simulations with widely varied physics, further validated with existing FIRE simulations. We show that gas supply, cooling/thermodynamics, star formation model, Toomre scale, galaxy dynamical times, and feedback properties do not have a direct causal effect on these transitions. Rather, both the formation of disks and cessation of bursty star formation are driven by the gravitational potential, but in different ways. Disk formation is promoted when the mass profile becomes sufficiently centrally-concentrated in shape (relative to circularization radii): we show that this provides a well-defined dynamical center, ceases to support the global 'breathing modes' which can persist indefinitely in less-concentrated profiles and efficiently destroy disks, promotes orbit mixing to form a coherent angular momentum, and stabilizes the disk. Smooth SF is promoted by the potential or escape velocity (not circular velocity) becoming sufficiently large at the radii of star formation that cool, mass-loaded (momentum-conserving) outflows are trapped/confined near the galaxy, as opposed to escaping after bursts. We discuss the detailed physics, how these conditions arise in cosmological contexts, their relation to other correlated phenomena (e.g. inner halo virialization, vertical disk 'settling'), and observations. △ Less

Submitted 28 August, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: Updated to match MNRAS accepted version. 45 pages, 32 figures

Spiral Arms are Metal Freeways: Azimuthal Gas-Phase Metallicity Variations in Simulated Cosmological Zoom-in Flocculent Disks

Authors: Matthew E. Orr, Blakesley Burkhart, Andrew Wetzel, Philip F. Hopkins, Ivanna A. Escala, Allison L. Strom, Paul F. Goldsmith, Jorge L. Pineda, Christopher C. Hayward, Sarah R. Loebman

Abstract: We examine the azimuthal variations in gas-phase metallicity profiles in simulated Milky Way mass disk galaxies from the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite, which includes a sub-grid turbulent metal mixing model. We produce spatially resolved maps of the disks at $z \approx 0$ with pixel sizes ranging from 250 to 750~pc, analogous to modern integral f… ▽ More We examine the azimuthal variations in gas-phase metallicity profiles in simulated Milky Way mass disk galaxies from the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite, which includes a sub-grid turbulent metal mixing model. We produce spatially resolved maps of the disks at $z \approx 0$ with pixel sizes ranging from 250 to 750~pc, analogous to modern integral field unit (IFU) galaxy surveys, mapping the gas-phase metallicities in both the cold & dense gas and the ionized gas correlated with HII regions. We report that the spiral arms alternate in a pattern of metal rich and metal poor relative to the median metallicity on the order of $\lesssim 0.1$~dex, appearing generally in this sample of flocculent spirals. The pattern persists even in a simulation with different strengths of metal mixing, indicating that the pattern emerges from physics above the sub-grid scale. Local enrichment does not appear to be the dominant source of the azimuthal metallicity variations at $z \approx 0$: there is no correlation with local star formation on these spatial scales. Rather, the arms are moving inwards and outwards relative to each other, carrying their local metallicity gradients with them radially before mixing into the larger-scale interstellar medium. We propose that the arms act as freeways channeling relatively metal poor gas radially inwards, and relatively enriched gas radially outwards. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 16 Pages, 14 Figures, submitted to MNRAS

Magnetic Fields on FIRE: Comparing B-fields in the multiphase ISM and CGM of Simulated L$_*$ Galaxies to Observations

Authors: Sam B. Ponnada, Georgia V. Panopoulou, Iryna S. Butsky, Philip F. Hopkins, Sarah R. Loebman, Cameron Hummels, Suoqing Ji, Andrew Wetzel, Claude-André Faucher-Giguère, Christopher C. Hayward

Abstract: The physics of magnetic fields ($\textbf{B}$) and cosmic rays (CRs) have recently been included in simulations of galaxy formation. However, significant uncertainties remain in how these components affect galaxy evolution. To understand their common observational tracers, we analyze the magnetic fields in a set of high-resolution, magneto-hydrodynamic, cosmological simulations of Milky-Way-like ga… ▽ More The physics of magnetic fields ($\textbf{B}$) and cosmic rays (CRs) have recently been included in simulations of galaxy formation. However, significant uncertainties remain in how these components affect galaxy evolution. To understand their common observational tracers, we analyze the magnetic fields in a set of high-resolution, magneto-hydrodynamic, cosmological simulations of Milky-Way-like galaxies from the FIRE-2 project. We compare mock observables of magnetic field tracers for simulations with and without CRs to observations of Zeeman splitting and rotation/dispersion measures. We find reasonable agreement between simulations and observations in both the neutral and the ionised interstellar medium (ISM). We find that the simulated galaxies with CRs show weaker ISM $|\textbf{B}|$ fields on average compared to their magnetic-field-only counterparts. This is a manifestation of the effects of CRs in the diffuse, low density inner circum-galactic medium (CGM). We find that equipartition between magnetic and cosmic ray energy densities may be valid at large ($>$ 1 kpc) scales for typical ISM densities of Milky-Way-like galaxies, but not in their halos. Within the ISM, the magnetic fields in our simulated galaxies follow a power-law scaling with gas density. The scaling extends down to hydrogen number densities $<$ 300 cm$^{-3}$, in contrast to observationally-derived models, but consistent with the observational measurements. Finally, we generate synthetic rotation measure (RM) profiles for projections of the simulated galaxies and compare to observational constraints in the CGM. While consistent with upper limits, improved data are needed to detect the predicted CGM RMs at 10-200 kpc and better constrain theoretical predictions. △ Less

Submitted 28 August, 2022; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: 15 pages, 9 figures. Accepted in MNRAS 2022 August 26. Received 2022 August 25; in original form 2022 June 10

Disconnecting the Dots: Re-Examining the Nature of Stellar "Strings" in the Milky Way

Authors: Catherine Zucker, J. E. G. Peek, Sarah R. Loebman

Abstract: Recent analyses of Gaia data have resulted in the identification of new stellar structures, including a new class of extended stellar filaments called stellar "strings", first proposed by Kounkel and Covey 2019. We explore the spatial, kinematic, and chemical composition of strings to demonstrate that these newfound structures are largely inconsistent with being physical objects whose members shar… ▽ More Recent analyses of Gaia data have resulted in the identification of new stellar structures, including a new class of extended stellar filaments called stellar "strings", first proposed by Kounkel and Covey 2019. We explore the spatial, kinematic, and chemical composition of strings to demonstrate that these newfound structures are largely inconsistent with being physical objects whose members share a common origin. Examining the 3D spatial distribution of string members, we find that the spatial dispersion around the claimed string spine does not improve in the latest Gaia DR3 data release -- despite tangible gains in the signal-to-noise of the parallax measurements -- counter to expectations of a bona fide structure. Using the radial velocity dispersion of the strings (averaging $\rmσ_{V_r} = 16\;km\;s^{-1}$) to estimate their virial masses, we find that all strings are gravitationally unbound. Given the finding that the strings are dispersing, the reported stellar ages of the strings are typically $120\times$ larger than their measured dispersal times. Finally, we validate prior work that stellar strings are more chemically homogeneous than their local field stars, but show it is possible to obtain the same signatures of chemical homogeneity by drawing random samples of stars from spatially, temporally, and kinematically unrelated open clusters. Our results show that while some strings may be composed of real sub-structures, there is no consistent evidence for larger string-like connections over the sample. These results underline the need for caution in over-interpreting the significance of these strings and their role in understanding the star formation history of the Milky Way. △ Less

Submitted 12 August, 2022; v1 submitted 27 May, 2022; originally announced May 2022.

Comments: Accepted for publication in The Astrophysical Journal

Journal ref: ApJ 936 160 (2022)

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

FIRE-3: Updated Stellar Evolution Models, Yields, & Microphysics and Fitting Functions for Applications in Galaxy Simulations

Authors: Philip F. Hopkins, Andrew Wetzel, Coral Wheeler, Robyn Sanderson, Michael Y. Grudic, Omid Sameie, Michael Boylan-Kolchin, Matthew Orr, Xiangcheng Ma, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert, Kung-Yi Su, Jorge Moreno, Robert Feldmann, James S. Bullock, Sarah R. Loebman, Daniel Angles-Alcazar, Jonathan Stern, Lina Necib, Christopher C. Hayward

Abstract: Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numeric… ▽ More Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numerical methods, here we update the stellar evolution tracks used to determine stellar feedback inputs, e.g. stellar mass-loss (O/B and AGB), spectra (luminosities and ionization rates), and supernova rates (core-collapse and Ia), as well as detailed mass-dependent yields. We also update the low-temperature cooling and chemistry, to enable improved accuracy at $T \lesssim 10^{4}\,$K and densities $n\gg 1\,{\rm cm^{-3}}$, and the meta-galactic ionizing background. All of these synthesize newer empirical constraints on these quantities and updated stellar evolution and yield models from a number of groups, addressing different aspects of stellar evolution. To make the updated models as accessible as possible, we provide fitting functions for all of the relevant updated tracks, yields, etc, in a form specifically designed so they can be directly 'plugged in' to existing galaxy formation simulations. We also summarize the default FIRE-3 implementations of 'optional' physics, including spectrally-resolved cosmic rays and supermassive black hole growth and feedback. △ Less

Submitted 28 February, 2022; originally announced March 2022.

Comments: 26 pages, 12 figures. Submitted to MNRAS. Comments welcome

Public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation

Authors: Andrew Wetzel, Christopher C. Hayward, Robyn E. Sanderson, Xiangcheng Ma, Daniel Angles-Alcazar, Robert Feldmann, T. K Chan, Kareem El-Badry, Coral Wheeler, Shea Garrison-Kimmel, Farnik Nikakhtar, Nondh Panithanpaisal, Arpit Arora, Alexander B. Gurvich, Jenna Samuel, Omid Sameie, Viraj Pandya, Zachary Hafen, Cameron Hummels, Sarah Loebman, Michael Boylan-Kolchin, James S. Bullock, Claude-Andre Faucher-Giguere, Dusan Keres, Eliot Quataert , et al. (1 additional authors not shown)

Abstract: We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, includ… ▽ More We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from 3 suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way-mass galaxies, 5 SMC/LMC-mass galaxies, and 4 lower-mass galaxies including 1 ultra-faint; we release 39 snapshots across z = 0 - 10. The second comprises 4 massive galaxies, with 19 snapshots across z = 1 - 10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5 - 10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all Milky Way-mass galaxies, we release the formation coordinates and an "ex-situ" flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations. △ Less

Submitted 29 March, 2023; v1 submitted 14 February, 2022; originally announced February 2022.

Comments: 22 pages. Accepted for publication in ApJS, matches the published version, with updated references. Now includes up to 39 snapshots per simulation and formation coordinates for all star particles. Data available at http://flathub.flatironinstitute.org/fire

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

Reading the CARDs: the Imprint of Accretion History in the Chemical Abundances of the Milky Way's Stellar Halo

Authors: Emily C. Cunningham, Robyn E. Sanderson, Kathryn V. Johnston, Nondh Panithanpaisal, Melissa K. Ness, Andrew Wetzel, Sarah R. Loebman, Ivanna Escala, Danny Horta, Claude-André Faucher-Giguère

Abstract: In the era of large-scale spectroscopic surveys in the Local Group (LG), we can explore using chemical abundances of halo stars to study the star formation and chemical enrichment histories of the dwarf galaxy progenitors of the Milky Way (MW) and M31 stellar halos. In this paper, we investigate using the Chemical Abundance Ratio Distributions (CARDs) of seven stellar halos from the Latte suite of… ▽ More In the era of large-scale spectroscopic surveys in the Local Group (LG), we can explore using chemical abundances of halo stars to study the star formation and chemical enrichment histories of the dwarf galaxy progenitors of the Milky Way (MW) and M31 stellar halos. In this paper, we investigate using the Chemical Abundance Ratio Distributions (CARDs) of seven stellar halos from the Latte suite of FIRE-2 simulations. We attempt to infer galaxies' assembly histories by modelling the CARDs of the stellar halos of the Latte galaxies as a linear combination of template CARDs from disrupted dwarfs, with different stellar masses $M_{\star}$ and quenching times $t_{100}$. We present a method for constructing these templates using present-day dwarf galaxies. For four of the seven Latte halos studied in this work, we recover the mass spectrum of accreted dwarfs to a precision of $<10\%$. For the fraction of mass accreted as a function of $t_{100}$, we find residuals of $20-30\%$ for five of the seven simulations. We discuss the failure modes of this method, which arise from the diversity of star formation and chemical enrichment histories dwarf galaxies can take. These failure cases can be robustly identified by the high model residuals. Though the CARDs modeling method does not successfully infer the assembly histories in these cases, the CARDs of these disrupted dwarfs contain signatures of their unusual formation histories. Our results are promising for using CARDs to learn more about the histories of the progenitors of the MW and M31 stellar halos. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: 32 pages, 20 figures. Submitted to ApJ; comments welcome!

New families in our Solar neighborhood: applying Gaussian Mixture models for objective classification of structures in the Milky Way and in simulations

Authors: Farnik Nikakhtar, Robyn E. Sanderson, Andrew Wetzel, Sarah Loebman, Sanjib Sharma, Rachael Beaton, J. Ted Mackereth, Vijith Jacob Poovelil, Gail Zasowski, Ana Bonaca, Sarah Martell, Henrik Jonsson, Claude-Andre Faucher-Giguere

Abstract: The standard picture of galaxy formation motivates the decomposition of the Milky Way into 3--4 stellar populations with distinct kinematic and elemental abundance distributions: the thin disk, thick disk, bulge, and stellar halo. To test this idea, we construct a Gaussian mixture model (GMM) for both simulated and observed stars in the Solar neighborhood, using measured velocities and iron abunda… ▽ More The standard picture of galaxy formation motivates the decomposition of the Milky Way into 3--4 stellar populations with distinct kinematic and elemental abundance distributions: the thin disk, thick disk, bulge, and stellar halo. To test this idea, we construct a Gaussian mixture model (GMM) for both simulated and observed stars in the Solar neighborhood, using measured velocities and iron abundances (i.e., an augmented Toomre diagram) as the distributions to be decomposed. We compare results for the Gaia-APOGEE DR16 crossmatch catalog of the Solar neighborhood with those from a suite of synthetic Gaia-APOGEE crossmatches constructed from FIRE-2 cosmological simulations of Milky Way-mass galaxies. We find that in both the synthetic and real data, the best-fit GMM uses five independent components, some of whose properties resemble the standard populations predicted by galaxy formation theory. Two components can be identified unambiguously as the thin disk and another as the halo. However, instead of a single counterpart to the thick disk, there are three intermediate components with different age and alpha abundance distributions (although these data are not used to construct the model). We use decompositions of the synthetic data to show that the classified components indeed correspond to stars with different origins. By analogy with the simulated data, we show that our mixture model of the real Gaia-APOGEE crossmatch distinguishes the following components: (1) a classic thin disk of young stars on circular orbits (46%), (2) thin disk stars heated by interactions with satellites (22%), (3, 4) two components representing the velocity asymmetry of the alpha-enhanced thick disk (27%), and (5) a stellar halo consistent with early, massive accretion (4%). △ Less

Submitted 16 April, 2021; originally announced April 2021.

Comments: 20 pages, 14 figures, 1 table

Stellar Migration and Chemical Enrichment in the Milky Way Disc: A Hybrid Model

Authors: James W. Johnson, David H. Weinberg, Fiorenzo Vincenzo, Jonathan C. Bird, Sarah R. Loebman, Alyson M. Brooks, Thomas R. Quinn, Charlotte R. Christensen, Emily J. Griffith

Abstract: We develop a hybrid model of galactic chemical evolution that combines a multi-ring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradien… ▽ More We develop a hybrid model of galactic chemical evolution that combines a multi-ring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradient in outflow mass loading tuned to reproduce the observed metallicity gradient. With this approach, the model reproduces many of the striking qualitative features of the Milky Way disc's abundance structure: (i) the dependence of the [O/Fe]-[Fe/H] distribution on radius $R_\text{gal}$ and midplane distance $|z|$; (ii) the changing shapes of the [O/H] and [Fe/H] distributions with $R_\text{gal}$ and $|z|$; (iii) a broad distribution of [O/Fe] at sub-solar metallicity and changes in the [O/Fe] distribution with $R_\text{gal}$, $|z|$, and [Fe/H]; (iv) a tight correlation between [O/Fe] and stellar age for [O/Fe] $>$ 0.1; (v) a population of young and intermediate-age $α$-enhanced stars caused by migration-induced variability in the Type Ia supernova rate; (vi) non-monotonic age-[O/H] and age-[Fe/H] relations, with large scatter and a median age of $\sim$4 Gyr near solar metallicity. Observationally motivated models with an enhanced star formation rate $\sim$2 Gyr ago improve agreement with the observed age-[Fe/H] and age-[O/H] relations, but worsen agreement with the observed age-[O/Fe] relation. None of our models predict an [O/Fe] distribution with the distinct bimodality seen in the observations, suggesting that more dramatic evolutionary pathways are required. All code and tables used for our models are publicly available through the Versatile Integrator for Chemical Evolution (VICE; https://pypi.org/project/vice). △ Less

Submitted 17 March, 2021; originally announced March 2021.

Comments: 24 pages; 18 figures; to be submitted to MNRAS; comments welcome

3-D gas-phase elemental abundances across the formation histories of Milky Way-mass galaxies in the FIRE simulations: initial conditions for chemical tagging

Authors: Matthew A. Bellardini, Andrew Wetzel, Sarah R. Loebman, Claude-André Faucher-Giguère, Xiangcheng Ma, Robert Feldmann

Abstract: We use FIRE-2 simulations to examine 3-D variations of gas-phase elemental abundances of [O/H], [Fe/H], and [N/H] in 11 Milky Way (MW) and M31-mass galaxies across their formation histories at $z \leq 1.5$ ($t_{\rm lookback} \leq 9.4$ Gyr), motivated by characterizing the initial conditions of stars for chemical tagging. Gas within $1$ kpc of the disk midplane is vertically homogeneous to… ▽ More We use FIRE-2 simulations to examine 3-D variations of gas-phase elemental abundances of [O/H], [Fe/H], and [N/H] in 11 Milky Way (MW) and M31-mass galaxies across their formation histories at $z \leq 1.5$ ($t_{\rm lookback} \leq 9.4$ Gyr), motivated by characterizing the initial conditions of stars for chemical tagging. Gas within $1$ kpc of the disk midplane is vertically homogeneous to $\lesssim 0.008$ dex at all $z \leq 1.5$. We find negative radial gradients (metallicity decreases with galactocentric radius) at all times, which steepen over time from $\approx -0.01$ dex kpc$^{-1}$ at $z = 1$ ($t_{\rm lookback} = 7.8$ Gyr) to $\approx -0.03$ dex kpc$^{-1}$ at $z = 0$, and which broadly agree with observations of the MW, M31, and nearby MW/M31-mass galaxies. Azimuthal variations at fixed radius are typically $0.14$ dex at $z = 1$, reducing to $0.05$ dex at $z = 0$. Thus, over time radial gradients become steeper while azimuthal variations become weaker (more homogeneous). As a result, azimuthal variations were larger than radial variations at $z \gtrsim 0.8$ ($t_{\rm lookback} \gtrsim 6.9$ Gyr). Furthermore, elemental abundances are measurably homogeneous (to $\lesssim 0.05$ dex) across a radial range of $ΔR \approx 3.5$ kpc at $z \gtrsim 1$ and $ΔR \approx 1.7$ kpc at $z = 0$. We also measure full distributions of elemental abundances, finding typically negatively skewed normal distributions at $z \gtrsim 1$ that evolve to typically Gaussian distributions by $z = 0$. Our results on gas abundances inform the initial conditions for stars, including the spatial and temporal scales for applying chemical tagging to understand stellar birth in the MW. △ Less

Submitted 19 August, 2021; v1 submitted 11 February, 2021; originally announced February 2021.

Comments: 24 pages, 17 figures, Accepted for publication in MNRAS

Fiery Cores: Bursty and Smooth Star Formation Distributions across Galaxy Centers in Cosmological Zoom-in Simulations

Authors: Matthew E. Orr, H Perry Hatchfield, Cara Battersby, Christopher C. Hayward, Philip F. Hopkins, Andrew Wetzel, Samantha M. Benincasa, Sarah R. Loebman, Mattia C. Sormani, Ralf S. Klessen

Abstract: We present an analysis of the $R\lesssim 1.5$ kpc core regions of seven simulated Milky Way mass galaxies, from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, for a finely sampled period ($Δt = 2.2$ Myr) of 22 Myr at $z \approx 0$, and compare them with star formation rate (SFR) and gas surface density observations of the Milky Way's Central Molecular Zone (… ▽ More We present an analysis of the $R\lesssim 1.5$ kpc core regions of seven simulated Milky Way mass galaxies, from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, for a finely sampled period ($Δt = 2.2$ Myr) of 22 Myr at $z \approx 0$, and compare them with star formation rate (SFR) and gas surface density observations of the Milky Way's Central Molecular Zone (CMZ). Despite not being tuned to reproduce the detailed structure of the CMZ, we find that four of these galaxies are consistent with CMZ observations at some point during this 22 Myr period. The galaxies presented here are not homogeneous in their central structures, roughly dividing into two morphological classes; (a) several of the galaxies have very asymmetric gas and SFR distributions, with intense (compact) starbursts occurring over a period of roughly 10 Myr, and structures on highly eccentric orbits through the CMZ, whereas (b) others have smoother gas and SFR distributions, with only slowly varying SFRs over the period analyzed. In class (a) centers, the orbital motion of gas and star-forming complexes across small apertures ($R \lesssim 150$pc, analogously $|l|<1^\circ$ in the CMZ observations) contributes as much to tracers of star formation/dense gas appearing in those apertures, as the internal evolution of those structures does. These asymmetric/bursty galactic centers can simultaneously match CMZ gas and SFR observations, demonstrating that time-varying star formation can explain the CMZ's low star formation efficiency. △ Less

Submitted 26 January, 2021; originally announced January 2021.

Comments: 9 pages, 3 figures, accepted to ApJ Letters. One appended supplemental appendix and figure, which do not appear in the ApJ version

One-Two Quench: A Double Minor Merger Scenario

Authors: N. Nicole Sanchez, Michael Tremmel, Jessica K. Werk, Andrew Pontzen, Charlotte Christensen, Thomas Quinn, Sarah Loebman, Akaxia Cruz

Abstract: Using the N-body+Smoothed particle hydrodynamics code, ChaNGa, we identify two merger-driven processes\textemdash disk disruption and supermassive black hole (SMBH) feedback\textemdash which work together to quench L$^*$ galaxies for over 7 Gyr. Specifically, we examine the cessation of star formation in a simulated Milky Way (MW) analog, driven by an interaction with two minor satellites. Both in… ▽ More Using the N-body+Smoothed particle hydrodynamics code, ChaNGa, we identify two merger-driven processes\textemdash disk disruption and supermassive black hole (SMBH) feedback\textemdash which work together to quench L$^*$ galaxies for over 7 Gyr. Specifically, we examine the cessation of star formation in a simulated Milky Way (MW) analog, driven by an interaction with two minor satellites. Both interactions occur within $\sim$100 Myr of each other, and the satellites both have masses 5 to 20 times smaller than that of their MW-like host galaxy. Using the genetic modification process of \cite{Roth2016}, we generate a set of four zoom-in, MW-mass galaxies all of which exhibit unique star formation histories due to small changes to their assembly histories. In two of these four cases, the galaxy is quenched by $z = 1$. Because these are controlled modifications, we are able to isolate the effects of two closely-spaced minor merger events, the relative timing of which determines whether the MW-mass main galaxy quenches. This one-two punch works to: 1. fuel the primary halo's supermassive black hole (SMBH) at its peak accretion rate; and 2. disrupt the cold, gaseous disk of the host galaxy. The end result is that feedback from the SMBH thoroughly and abruptly ends the galaxy's star formation by $z\approx1$. We search for and find a similar quenching event in {\sc Romulus25}, a hydrodynamical $(25\,\mathrm{Mpc})^3$ volume simulation, demonstrating that the mechanism is common enough to occur even in a small sample of MW-mass quenched galaxies at $z=0$. △ Less

Submitted 25 May, 2021; v1 submitted 11 September, 2020; originally announced September 2020.

Comments: 12 pages, 9 figures, accepted to ApJ

Inside Out and Upside-Down: The Roles of Gas Cooling and Dynamical Heating in Shaping the Stellar Age-Velocity Relation

Authors: Jonathan C. Bird, Sarah R. Loebman, David H. Weinberg, Alyson Brooks, Thomas R. Quinn, Charlotte R. Christensen

Abstract: Kinematic studies of disk galaxies, using individual stars in the Milky Way or statistical studies of global disk kinematics over time, provide insight into how disks form and evolve. We use a high-resolution, cosmological zoom-simulation of a Milky Way-mass disk galaxy h277 to tie together local disk kinematics and the evolution of the disk over time. The present-day stellar age-velocity relation… ▽ More Kinematic studies of disk galaxies, using individual stars in the Milky Way or statistical studies of global disk kinematics over time, provide insight into how disks form and evolve. We use a high-resolution, cosmological zoom-simulation of a Milky Way-mass disk galaxy h277 to tie together local disk kinematics and the evolution of the disk over time. The present-day stellar age-velocity relationship (AVR) of h277 is nearly identical to that of the analogous solar-neighborhood measurement in the Milky Way. A crucial element of this success is the simulation's dynamically cold multi-phase ISM, which allows young stars to form with a low velocity dispersion ($σ_{\mathrm{birth}}$$\sim 6 - 8\ \mathrm{km}\, \mathrm{s}^{-1}$) at late times. Older stars are born kinematically hotter (i.e., the disk settles over time in an "upside-down" formation scenario), and are subsequently heated after birth. The disk also grows "inside-out", and many of the older stars in the solar neighborhood at $z=0$ are present because of radial mixing. We demonstrate that the evolution of $σ_{\mathrm{birth}}$ in h277 can be explained by the same model used to describe the general decrease in velocity dispersion observed in disk galaxies from $z\sim 2-3$ to the present-day, in which the disk evolves in quasi-stable equilibrium and the ISM velocity dispersion decreases over time due to a decreasing gas fraction. Thus, our results tie together local observations of the Milky Way's AVR with observed kinematics of high $z$ disk galaxies. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: Submitting to MNRAS. Comments welcomed. This simulation is publicly available. Contact Alyson Brooks for details

Pressure balance in the multiphase ISM of cosmologically simulated disk galaxies

Authors: Alexander B. Gurvich, Claude-André Faucher-Giguère, Alexander J. Richings, Philip F. Hopkins, Michael Y. Grudić, Zachary Hafen, Sarah Wellons, Jonathan Stern, Eliot Quataert, T. K. Chan, Matthew E. Orr, Dušan Kereš, Andrew Wetzel, Christopher C. Hayward, Sarah R. Loebman, Norman Murray

Abstract: Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disk galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. W… ▽ More Pressure balance plays a central role in models of the interstellar medium (ISM), but whether and how pressure balance is realized in a realistic multiphase ISM is not yet well understood. We address this question using a set of FIRE-2 cosmological zoom-in simulations of Milky Way-mass disk galaxies, in which a multiphase ISM is self-consistently shaped by gravity, cooling, and stellar feedback. We analyze how gravity determines the vertical pressure profile as well as how the total ISM pressure is partitioned between different phases and components (thermal, dispersion/turbulence, and bulk flows). We show that, on average and consistent with previous more idealized simulations, the total ISM pressure balances the weight of the overlying gas. Deviations from vertical pressure balance increase with increasing galactocentric radius and with decreasing averaging scale. The different phases are in rough total pressure equilibrium with one another, but with large deviations from thermal pressure equilibrium owing to kinetic support in the cold and warm phases, which dominate the total pressure near the midplane. Bulk flows (e.g., inflows and fountains) are important at a few disk scale heights, while thermal pressure from hot gas dominates at larger heights. Overall, the total midplane pressure is well-predicted by the weight of the disk gas, and we show that it also scales linearly with the star formation rate surface density (Sigma_SFR). These results support the notion that the Kennicutt-Schmidt relation arises because Sigma_SFR and the gas surface density (Sigma_g) are connected via the ISM midplane pressure. △ Less

Submitted 16 November, 2020; v1 submitted 26 May, 2020; originally announced May 2020.

Comments: 22 pages, 16 figures, accepted in MNRAS

Reproducing the CO-to-H$_2$ conversion factor in cosmological simulations of Milky Way-mass galaxies

Authors: Laura C. Keating, Alexander J. Richings, Norman Murray, Claude-Andre Faucher-Giguere, Philip F. Hopkins, Andrew Wetzel, Dusan Keres, Samantha Benincasa, Robert Feldmann, Sarah Loebman, Matthew E. Orr

Abstract: We present models of CO(1-0) emission from Milky Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results… ▽ More We present models of CO(1-0) emission from Milky Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which in our models sets the attenuation of the incident UV radiation field. At the resolution of these simulations, commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H$_2$ abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of $\sim 3$ pc in cold gas ($T < 300$ K) for both CO and H$_{2}$, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H$_{2}$ conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model which controls the amount of radiation that penetrates giant molecular clouds. △ Less

Submitted 28 September, 2020; v1 submitted 22 January, 2020; originally announced January 2020.

Comments: Accepted to MNRAS. 16 pages, 10 figures

Live Fast, Die Young: GMC lifetimes in the FIRE cosmological simulations of Milky Way-mass galaxies

Authors: Samantha M. Benincasa, Sarah R. Loebman, Andrew Wetzel, Philip F. Hopkins, Norman Murray, Matthew A. Bellardini, Claude-André Faucher-Giguère, Dávid Guszejnov, Matthew Orr

Abstract: We present the first measurement of the lifetimes of Giant Molecular Clouds (GMCs) in cosmological simulations at $z = 0$, using the Latte suite of FIRE-2 simulations of Milky Way-mass galaxies. We track GMCs with total gas mass $\gtrsim 10^5$ M$_\odot$ at high spatial ($\sim1$ pc), mass ($7100$ M$_{\odot}$), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution… ▽ More We present the first measurement of the lifetimes of Giant Molecular Clouds (GMCs) in cosmological simulations at $z = 0$, using the Latte suite of FIRE-2 simulations of Milky Way-mass galaxies. We track GMCs with total gas mass $\gtrsim 10^5$ M$_\odot$ at high spatial ($\sim1$ pc), mass ($7100$ M$_{\odot}$), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution of masses for massive GMCs in the Milky Way and nearby galaxies. We find GMC lifetimes of $5-7$ Myr, or 1-2 freefall times, on average, with less than 2$\%$ of clouds living longer than 20 Myr. We find decreasing GMC lifetimes with increasing virial parameter, and weakly increasing GMC lifetimes with galactocentric radius, implying that environment affects the evolutionary cycle of GMCs. However, our GMC lifetimes show no systematic dependence on GMC mass or amount of star formation. These results are broadly consistent with inferences from the literature and provide an initial investigation into ultimately understanding the physical processes that govern GMC lifetimes in a cosmological setting. △ Less

Submitted 2 December, 2020; v1 submitted 12 November, 2019; originally announced November 2019.

Comments: accepted to MNRAS

Evolution of giant molecular clouds across cosmic time

Authors: David Guszejnov, Michael Y. Grudić, Stella S. R. Offner, Michael Boylan-Kolchin, Claude-André Faucher-Giguère, Andrew Wetzel, Samantha M. Benincasa, Sarah Loebman

Abstract: Giant molecular clouds (GMCs) are well-studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of giant molecular clouds in a Milky Way-like galaxy and an LMC-like dwarf galaxy of the FIRE-2 (Feedback In Realis… ▽ More Giant molecular clouds (GMCs) are well-studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of giant molecular clouds in a Milky Way-like galaxy and an LMC-like dwarf galaxy of the FIRE-2 (Feedback In Realistic Environments) simulation suite, which have sufficient resolution to predict the bulk properties of GMCs in cosmological galaxy formation self-consistently. We show explicitly that the majority of star formation outside the galactic center occurs within self-gravitating gas structures that have properties consistent with observed bound GMCs. We find that the typical cloud bulk properties such as mass and surface density do not vary more than a factor of 2 in any systematic way after the first Gyr of cosmic evolution within a given galaxy from its progenitor. While the median properties are constant, the tails of the distributions can briefly undergo drastic changes, which can produce very massive and dense self-gravitating gas clouds. Once the galaxy forms, we identify only two systematic trends in bulk properties over cosmic time: a steady increase in metallicity produced by previous stellar populations and a weak decrease in bulk cloud temperatures. With the exception of metallicity we find no significant differences in cloud properties between the Milky Way-like and dwarf galaxies. These results have important implications for cosmological star and star cluster formation and put especially strong constraints on theories relating the stellar initial mass function to cloud properties. △ Less

Submitted 24 February, 2020; v1 submitted 2 October, 2019; originally announced October 2019.

Comments: 15 pages, 11 figures accepted by MNRAS

When Cold Radial Migration is Hot: Constraints from Resonant Overlap

Authors: Kathryne J. Daniel, David A. Schaffner, Fiona McCluskey, Codie Fiedler Kawaguchi, Sarah Loebman

Abstract: It is widely accepted that stars in a spiral disk, like the Milky Way's, can radially migrate on order a scale length over the disk's lifetime. With the exception of cold torquing, also known as "churning," processes that contribute to the radial migration of stars are necessarily associated with kinematic heating. Additionally, it is an open question whether or not an episode of cold torquing is… ▽ More It is widely accepted that stars in a spiral disk, like the Milky Way's, can radially migrate on order a scale length over the disk's lifetime. With the exception of cold torquing, also known as "churning," processes that contribute to the radial migration of stars are necessarily associated with kinematic heating. Additionally, it is an open question whether or not an episode of cold torquing is kinemically cold over long radial distances. This study uses a suite of analytically based simulations to investigate the dynamical response when stars are subject to cold torquing and are also resonant with an ultraharmonic. Model results demonstrate that these populations are kinematically heated and have RMS changes in orbital angular momentum around corotation that can exceed those of populations that do not experience resonant overlap. Thus, kinematic heating can occur during episodes of cold torquing. In a case study of a Milky Way-like disk with an exponential surface density profile and flat rotation curve, up to 40% of cold torqued stars in the solar cylinder experience resonant overlap. This fraction increases toward the galactic center. To first approximation, the maximum radial excursions from cold torquing depend only on the strength of the spiral pattern and the underlying rotation curve. This work places an upper limit to these excursions to be the distance between the ultraharmonics, otherwise radial migration near corotation can kinematically heat. The diffusion rate for kinematically cold radial migration is thus constrained by limiting the step size in the random walk approximation. △ Less

Submitted 23 July, 2019; originally announced July 2019.

Comments: 18 pages, 8 figures, 4 tables, Accepted to ApJ

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

The imprint of clump formation at high redshift. I. A disc alpha-abundance dichotomy

Authors: Adam J. Clarke, Victor P. Debattista, David L. Nidever, Sarah R. Loebman, Raymond C. Simons, Susan Kassin, Min Du, Melissa Ness, Deanne B. Fisher, Thomas R. Quinn, James Wadsley, Ken C. Freeman, Cristina C. Popescu

Abstract: The disc structure of the Milky Way is marked by a chemical dichotomy, with high-alpha and low-alpha abundance sequences, traditionally identified with the geometric thick and thin discs. This identification is aided by the old ages of the high-alpha stars, and lower average ages of the low-alpha ones. Recent large scale surveys such as APOGEE have provided a wealth of data on this chemical struct… ▽ More The disc structure of the Milky Way is marked by a chemical dichotomy, with high-alpha and low-alpha abundance sequences, traditionally identified with the geometric thick and thin discs. This identification is aided by the old ages of the high-alpha stars, and lower average ages of the low-alpha ones. Recent large scale surveys such as APOGEE have provided a wealth of data on this chemical structure, including showing that an identification of chemical and geometric thick discs is not exact, but the origin of the chemical dichotomy has remained unclear. Here we demonstrate that a dichotomy arises naturally if the early gas-rich disc fragments, leading to some fraction of the star formation occurring in clumps of the type observed in high-redshift galaxies. These clumps have high star formation rate density. They, therefore, enrich rapidly, moving from the low-alpha to the high-alpha sequence, while more distributed star formation produces the low-alpha sequence. We demonstrate that this model produces a chemically-defined thick disc that has many of the properties of the Milky Way's thick disc. Because clump formation is common in high redshift galaxies, we predict that chemical bimodalities are common in massive galaxies. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: Accepted to MNRAS

HALO7D II: The Halo Velocity Ellipsoid and Velocity Anisotropy with Distant Main Sequence Stars

Authors: Emily C. Cunningham, Alis J. Deason, Robyn E. Sanderson, Sangmo Tony Sohn, Jay Anderson, Puragra Guhathakurta, Constance M. Rockosi, Roeland P. van der Marel, Sarah R. Loebman, Andrew Wetzel

Abstract: The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way (MW) halo main sequence turnoff stars in the CANDELS fields. In this paper, the second in the HALO7D series, we present the proper motions for the HALO7D sample. We discuss our measurement methodology, which makes use… ▽ More The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way (MW) halo main sequence turnoff stars in the CANDELS fields. In this paper, the second in the HALO7D series, we present the proper motions for the HALO7D sample. We discuss our measurement methodology, which makes use of a Bayesian mixture modeling approach for creating the stationary reference frame of distant galaxies. Using the 3D kinematic HALO7D sample, we estimate the parameters of the halo velocity ellipsoid, $\langle v_φ \rangle, σ_r, σ_φ, σ_θ$, and the velocity anisotropy $β$. Using the full HALO7D sample, we find $β=0.63 \pm 0.05$ at $\langle r \rangle =24$ kpc. We also estimate the ellipsoid parameters for our sample split into three apparent magnitude bins; the posterior medians for these estimates of $β$, while consistent with one another, increase as a function of mean sample distance. Finally, we estimate $β$ in each of the individual HALO7D fields. We find that the velocity anisotropy $β$ can vary from field to field, which suggests that the halo is not phase mixed at $\langle r \rangle =24$ kpc. We explore the $β$ variation across the skies of two stellar halos from the \textit{Latte} suite of FIRE-2 simulations, finding that both simulated galaxies show $β$ variation over a similar range to the variation observed across the four HALO7D fields. The accretion histories of the two simulated galaxies result in different $β$ variation patterns; spatially mapping $β$ is thus a way forward in characterizing the accretion history of the Galaxy. △ Less

Submitted 29 October, 2018; originally announced October 2018.

Comments: 23 pages, 13 figures. Submitted to ApJ; comments welcome

Synthetic Gaia surveys from the FIRE cosmological simulations of Milky Way-mass galaxies

Authors: Robyn E. Sanderson, Andrew Wetzel, Sarah Loebman, Sanjib Sharma, Philip F. Hopkins, Shea Garrison-Kimmel, Claude-André Faucher-Giguère, Dušan Kereš, Eliot Quataert

Abstract: With Gaia Data Release 2, the astronomical community is entering a new era of multidimensional surveys of the Milky Way. This new phase-space view of our Galaxy demands new tools for comparing observations to simulations of Milky-Way-mass galaxies in a cosmological context, to test the physics of both dark matter and galaxy formation. We present ananke, a framework for generating synthetic phase-s… ▽ More With Gaia Data Release 2, the astronomical community is entering a new era of multidimensional surveys of the Milky Way. This new phase-space view of our Galaxy demands new tools for comparing observations to simulations of Milky-Way-mass galaxies in a cosmological context, to test the physics of both dark matter and galaxy formation. We present ananke, a framework for generating synthetic phase-space surveys from high-resolution baryonic simulations, and use it to generate a suite of synthetic surveys resembling Gaia DR2 in data structure, magnitude limits, and observational errors. We use three cosmological simulations of Milky-Way-mass galaxies from the Latte suite of the Feedback In Realistic Environments (FIRE) project, which feature self-consistent clustering of star formation in dense molecular clouds and thin stellar/gaseous disks in live cosmological halos with satellite dwarf galaxies and stellar halos. We select three solar viewpoints from each simulation to generate nine synthetic Gaia-like surveys. We sample synthetic stars by assuming each star particle (of mass 7070 $M_{\odot}$) represents a single stellar population. At each viewpoint, we compute dust extinction from the simulated gas metallicity distribution and apply a simple error model to produce a synthetic Gaia-like survey that includes both observational properties and a pointer to the generating star particle. We provide the complete simulation snapshot at $z = 0$ for each simulated galaxy. We describe data access points, the data model, and plans for future upgrades. These synthetic surveys provide a tool for the scientific community to test analysis methods and interpret Gaia data. △ Less

Submitted 19 March, 2020; v1 submitted 27 June, 2018; originally announced June 2018.

Comments: Matches accepted version. Data accompanying this paper are available at https://ananke.hub.yt and https://binder.flatironinstitute.org/~rsanderson/ananke. More info and updates at http://fire.northwestern.edu/ananke

The Impact of Stellar Migration on Disk Outskirts

Authors: Victor P. Debattista, Rok Roskar, Sarah R. Loebman

Abstract: Stellar migration, whether due to trapping by transient spirals (churning), or to scattering by non-axisymmetric perturbations, has been proposed to explain the presence of stars in outer disks. After a review of the basic theory, we present compelling, but not yet conclusive, evidence that churning has been important in the outer disks of galaxies with type II (down-bending) profiles, while scatt… ▽ More Stellar migration, whether due to trapping by transient spirals (churning), or to scattering by non-axisymmetric perturbations, has been proposed to explain the presence of stars in outer disks. After a review of the basic theory, we present compelling, but not yet conclusive, evidence that churning has been important in the outer disks of galaxies with type II (down-bending) profiles, while scattering has produced the outer disks of type III (up-bending) galaxies. In contrast, field galaxies with type I (pure exponential) profiles appear to not have experienced substantial migration. We conclude by suggesting work that would improve our understanding of the origin of outer disks. △ Less

Submitted 6 June, 2017; originally announced June 2017.

Comments: Invited review, Book chapter in "Outskirts of Galaxies", Eds. J. H. Knapen, J. C. Lee and A. Gil de Paz, Astrophysics and Space Science Library, Springer, in press 39 pages, 15 figures

Reliability of the measured velocity anisotropy of the Milky Way stellar halo

Authors: Kohei Hattori, Monica Valluri, Sarah R. Loebman, Eric F. Bell

Abstract: Determining the velocity distribution of halo stars is essential for estimating the mass of the Milky Way and for inferring its formation history. Since the stellar halo is a dynamically hot system, the velocity distribution of halo stars is well described by the 3-dimensional velocity dispersions $(σ_r, σ_θ, σ_φ)$, or by the velocity anisotropy parameter $β=1-(σ_θ^2+σ_φ^2)/(2σ_r^2)$. Direct measu… ▽ More Determining the velocity distribution of halo stars is essential for estimating the mass of the Milky Way and for inferring its formation history. Since the stellar halo is a dynamically hot system, the velocity distribution of halo stars is well described by the 3-dimensional velocity dispersions $(σ_r, σ_θ, σ_φ)$, or by the velocity anisotropy parameter $β=1-(σ_θ^2+σ_φ^2)/(2σ_r^2)$. Direct measurements of $(σ_r, σ_θ, σ_φ)$ consistently suggest $β=0.5$-$0.7$ for nearby halo stars. In contrast, the value of $β$ at large Galactocentric radius $r$ is still controversial, since reliable proper motion data are available for only a handful of stars. In the last decade, several authors have tried to estimate $β$ for distant halo stars by fitting the observed line-of-sight velocities at each radius with simple velocity distribution models (local fitting methods). Some results of local fitting methods imply $β<0$ at $r \gtrsim 20 \;\rm{kpc}$, which is inconsistent with recent predictions from cosmological simulations. Here we perform mock-catalogue analyses to show that the estimates of $β$ based on local fitting methods are reliable only at $r \leq 15 \;\rm{kpc}$ with the current sample size ($\sim10^3$ stars at a given radius). As $r$ increases, the line-of-sight velocity (corrected for the Solar reflex motion) becomes increasingly closer to the Galactocentric radial velocity, so that it becomes increasingly more difficult to estimate tangential velocity dispersion $(σ_θ, σ_φ)$ from line-of-sight velocity distribution. Our results suggest that the forthcoming Gaia data will be crucial for understanding the velocity distribution of halo stars at $r \geq 20\;\rm{kpc}$. △ Less

Submitted 30 April, 2017; v1 submitted 20 April, 2017; originally announced April 2017.

Comments: [Deleted non-relevant figures] ApJ submitted. Comments welcome. 20 pages (14 pages + 6 pages for Appendix). 13 figures. Main result: Fig 7. Schematic diagram: Fig 9. Companion paper to Loebman et al. (2017)

Beta Dips in the Gaia Era: Simulation Predictions of the Galactic Velocity Anisotropy Parameter for Stellar Halos

Authors: Sarah R. Loebman, Monica Valluri, Kohei Hattori, Victor P. Debattista, Eric F. Bell, Greg Stinson, Charlotte R. Christensen, Alyson Brooks, Thomas R. Quinn, Fabio Governato

Abstract: The velocity anisotropy parameter, beta, is a measure of the kinematic state of orbits in the stellar halo which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for beta as a function of radius from three suites of simulations, including accretion only and cosmological hydrodynamic simulations. We find that both types of simulations are consisten… ▽ More The velocity anisotropy parameter, beta, is a measure of the kinematic state of orbits in the stellar halo which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for beta as a function of radius from three suites of simulations, including accretion only and cosmological hydrodynamic simulations. We find that both types of simulations are consistent and predict strong radial anisotropy (<beta>~0.7) for Galactocentric radii greater than 10 kpc. Previous observations of beta for the MW's stellar halo claim a detection of an isotropic or tangential "dip" at r~20 kpc. Using the N-body+SPH simulations, we investigate the temporal persistence, population origin, and severity of "dips" in beta. We find dips in the in situ stellar halo are long-lived, while dips in the accreted stellar halo are short-lived and tied to the recent accretion of satellite material. We also find that a major merger as early as z~1 can result in a present day low (isotropic to tangential) value of beta over a wide range of radii and angular expanse. While all of these mechanisms are plausible drivers for the beta dip observed in the MW, in the simulations, each mechanism has a unique metallicity signature associated with it, implying that future spectroscopic surveys could distinguish between them. Since an accurate knowledge of beta(r) is required for measuring the mass of the MW halo, we note significant transient dips in beta could cause an overestimate of the halo's mass when using spherical Jeans equation modeling. △ Less

Submitted 3 February, 2018; v1 submitted 20 April, 2017; originally announced April 2017.

Comments: 14 pages, 7 figures, accepted to ApJ, companion paper to Hattori et al. (2017), appendix with Gaia specific predictions

The DESI Experiment Part II: Instrument Design

Authors: DESI Collaboration, Amir Aghamousa, Jessica Aguilar, Steve Ahlen, Shadab Alam, Lori E. Allen, Carlos Allende Prieto, James Annis, Stephen Bailey, Christophe Balland, Otger Ballester, Charles Baltay, Lucas Beaufore, Chris Bebek, Timothy C. Beers, Eric F. Bell, José Luis Bernal, Robert Besuner, Florian Beutler, Chris Blake, Hannes Bleuler, Michael Blomqvist, Robert Blum, Adam S. Bolton, Cesar Briceno , et al. (268 additional authors not shown)

Abstract: DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from… ▽ More DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use. △ Less

Submitted 13 December, 2016; v1 submitted 31 October, 2016; originally announced November 2016.

The DESI Experiment Part I: Science,Targeting, and Survey Design

Authors: DESI Collaboration, Amir Aghamousa, Jessica Aguilar, Steve Ahlen, Shadab Alam, Lori E. Allen, Carlos Allende Prieto, James Annis, Stephen Bailey, Christophe Balland, Otger Ballester, Charles Baltay, Lucas Beaufore, Chris Bebek, Timothy C. Beers, Eric F. Bell, José Luis Bernal, Robert Besuner, Florian Beutler, Chris Blake, Hannes Bleuler, Michael Blomqvist, Robert Blum, Adam S. Bolton, Cesar Briceno , et al. (268 additional authors not shown)

Abstract: DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure… ▽ More DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions. △ Less

Submitted 13 December, 2016; v1 submitted 31 October, 2016; originally announced November 2016.

Detailed Chemical Abundances in the r-Process-Rich Ultra-Faint Dwarf Galaxy Reticulum 2

Authors: Ian U. Roederer, Mario Mateo, John I. Bailey III, Yingyi Song, Eric F. Bell, Jeffrey D. Crane, Sarah Loebman, David L. Nidever, Edward W. Olszewski, Stephen A. Shectman, Ian B. Thompson, Monica Valluri, Matthew G. Walker

Abstract: The ultra-faint dwarf galaxy Reticulum 2 (Ret 2) was recently discovered in images obtained by the Dark Energy Survey. We have observed the four brightest red giants in Ret 2 at high spectral resolution using the Michigan/Magellan Fiber System. We present detailed abundances for as many as 20 elements per star, including 12 elements heavier than the Fe group. We confirm previous detection of high… ▽ More The ultra-faint dwarf galaxy Reticulum 2 (Ret 2) was recently discovered in images obtained by the Dark Energy Survey. We have observed the four brightest red giants in Ret 2 at high spectral resolution using the Michigan/Magellan Fiber System. We present detailed abundances for as many as 20 elements per star, including 12 elements heavier than the Fe group. We confirm previous detection of high levels of r-process material in Ret 2 (mean [Eu/Fe]=+1.69+/-0.05) found in three of these stars (mean [Fe/H]=-2.88+/-0.10). The abundances closely match the r-process pattern found in the well-studied metal-poor halo star CS22892-052. Such r-process-enhanced stars have not been found in any other ultra-faint dwarf galaxy, though their existence has been predicted by at least one model. The fourth star in Ret 2 ([Fe/H]=-3.42+/-0.20) contains only trace amounts of Sr ([Sr/Fe]=-1.73+/-0.43) and no detectable heavier elements. One r-process enhanced star is also enhanced in C (natal [C/Fe]=+1.1). This is only the third such star known, which suggests that the nucleosynthesis sites leading to C and r-process enhancements are decoupled. The r-process-deficient star is enhanced in Mg ([Mg/Fe]=+0.81+/-0.14), and the other three stars show normal levels of alpha-enhancement (mean [Mg/Fe]=+0.34+/-0.03). The abundances of other alpha and Fe-group elements closely resemble those in ultra-faint dwarf galaxies and metal-poor halo stars, suggesting that the nucleosynthesis that led to the large r-process enhancements either produced no light elements or produced light-element abundance signatures indistinguishable from normal supernovae. △ Less

Submitted 15 January, 2016; originally announced January 2016.

Comments: Accepted for publication in the Astronomical Journal. 12 pages, 6 figures, 8 tables

Imprints of radial migration on the Milky Way's metallicity distribution functions

Authors: Sarah R. Loebman, Victor P. Debattista, David L. Nidever, Michael R. Hayden, Jon A. Holtzman, Adam J. Clarke, Rok Roskar, Monica Valluri

Abstract: Recent analysis of the SDSS-III/APOGEE Data Release 12 stellar catalogue has revealed that the Milky Way's metallicity distribution function (MDF) changes shape as a function of radius, transitioning from being negatively skewed at small Galactocentric radii to positively skewed at large Galactocentric radii. Using a high resolution, N-body+SPH simulation, we show that the changing skewness arises… ▽ More Recent analysis of the SDSS-III/APOGEE Data Release 12 stellar catalogue has revealed that the Milky Way's metallicity distribution function (MDF) changes shape as a function of radius, transitioning from being negatively skewed at small Galactocentric radii to positively skewed at large Galactocentric radii. Using a high resolution, N-body+SPH simulation, we show that the changing skewness arises from radial migration - metal-rich stars form in the inner disk and subsequently migrate to the metal-poorer outer disk. These migrated stars represent a large fraction (> 50%) of the stars in the outer disk; they populate the high metallicity tail of the MDFs and are, in general, more metal-rich than the surrounding outer disk gas. The simulation also reproduces another surprising APOGEE result: the spatially invariant high-[alpha/Fe] MDFs. This arises in the simulation from the migration of a population formed within a narrow range of radii (3.2+/-1.2 kpc) and time (8.8+/-0.6 Gyr ago), rather than from spatially extended star formation in a homogeneous medium at early times. These results point toward the crucial role radial migration has played in shaping our Milky Way. △ Less

Submitted 8 February, 2016; v1 submitted 19 November, 2015; originally announced November 2015.

Comments: 5 pages, 5 figures, published in ApJ Letters

Magellan/M2FS Spectroscopy of Tucana 2 and Grus 1

Authors: Matthew G. Walker, Mario Mateo, Edward W. Olszewski, Sergey E. Koposov, Vasily Belokurov, Prashin Jethwa, David L. Nidever, Vincent Bonnivard, John I. Bailey III, Eric F. Bell, Sarah R. Loebman

Abstract: We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of $147$ stellar targets along the line of sight to the newly-discovered `ultrafaint' stellar systems Tucana 2 (Tuc 2) and Grus 1 (Gru 1). Based on simultaneous estimates of line-of-sight velocity and stellar-atmospheric parameters, we identify 8 and 7 stars as probable members of Tuc 2 and and Gru 1,… ▽ More We present results from spectroscopic observations with the Michigan/Magellan Fiber System (M2FS) of $147$ stellar targets along the line of sight to the newly-discovered `ultrafaint' stellar systems Tucana 2 (Tuc 2) and Grus 1 (Gru 1). Based on simultaneous estimates of line-of-sight velocity and stellar-atmospheric parameters, we identify 8 and 7 stars as probable members of Tuc 2 and and Gru 1, respectively. Our sample for Tuc 2 is sufficient to resolve an internal velocity dispersion of $8.6_{-2.7}^{+4.4}$ km s$^{-1}$ about a mean of $-129.1_{-3.5}^{+3.5}$ km s$^{-1}$ (solar rest frame), and to estimate a mean metallicity of [Fe/H]= $-2.23_{-0.12}^{+0.18}$. These results place Tuc 2 on chemodynamical scaling relations followed by dwarf galaxies, suggesting a dominant dark matter component with dynamical mass $2.7_{-1.3}^{+3.1}\times 10^6$ $\mathrm{M}_{\odot}$ enclosed within the central $\sim 160$ pc, and dynamical mass-to-light ratio $1900_{-900}^{+2200}$ $\mathrm{M}_{\odot}/L_{V,\odot}$. For Gru 1 we estimate a mean velocity of $-140.5_{-1.6}^{+2.4}$ km s$^{-1}$ and a mean metallicity of [Fe/H]=$-1.42_{-0.42}^{+0.55}$, but our sample does not resolve Gru 1's velocity dispersion. The radial coordinates of Tuc 2 and Gru 1 in Galactic phase space suggest that their orbits are among the most energetic within distance $\leq 300$ kpc. Moreover, their proximity to each other in this space arises naturally if both objects are trailing the Large Magellanic Cloud. △ Less

Submitted 14 December, 2015; v1 submitted 19 November, 2015; originally announced November 2015.

Comments: replaced with ApJ-accepted version, all spectra and data products (including samples from posterior PDFs) are available at http://www.andrew.cmu.edu/user/mgwalker/tuc2gru1_dataproducts.tar.gz

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

The continued optical to mid-IR evolution of V838 Monocerotis

Authors: S. R. Loebman, J. P. Wisniewski, S. J. Schmidt, A. F. Kowalski, R. K. Barry, K. S. Bjorkman, H. B. Hammel, S. L. Hawley, L. Hebb, M. M. Kasliwal, D. K. Lynch, R. W. Russell, M. L. Sitko, P. Szkody

Abstract: The eruptive variable V838 Monocerotis gained notoriety in 2002 when it brightened nine magnitudes in a series of three outbursts and then rapidly evolved into an extremely cool supergiant. We present optical, near-IR, and mid-IR spectroscopic and photometric observations of V838 Monocerotis obtained between 2008 and 2012 at the Apache Point Observatory 3.5m, NASA IRTF 3m, and Gemini South 8m tele… ▽ More The eruptive variable V838 Monocerotis gained notoriety in 2002 when it brightened nine magnitudes in a series of three outbursts and then rapidly evolved into an extremely cool supergiant. We present optical, near-IR, and mid-IR spectroscopic and photometric observations of V838 Monocerotis obtained between 2008 and 2012 at the Apache Point Observatory 3.5m, NASA IRTF 3m, and Gemini South 8m telescopes. We contemporaneously analyze the optical & IR spectroscopic properties of V838 Monocerotis to arrive at a revised spectral type L3 supergiant and effective temperature Teff~2000--2200 K. Because there are no existing optical observational data for L supergiants in the optical, we speculate that V838 Monocerotis may represent the prototype for L supergiants in this wavelength regime. We find a low level of Halpha emission present in the system, consistent with interaction between V838 Monocerotis and its B3V binary; however, we cannot rule out a stellar collision as the genesis event, which could result in the observed Halpha activity. Based upon a two-component blackbody fit to all wavelengths of our data, we conclude that, as of 2009, a shell of ejecta surrounded V838 Monocerotis at a radius of R=263+/-10 AU with a temperature of T=285+/-2 K. This result is consistent with IR interferometric observations from the same era and predictions from the Lynch et al. model of the expanding system, which provides a simple framework for understanding this complicated system. △ Less

Submitted 8 September, 2014; originally announced September 2014.

Comments: 6 pages, 2 tables, 6 figures; accepted to AJ

The Milky Way Tomography with SDSS. V. Mapping the Dark Matter Halo

Authors: Sarah R. Loebman, Zeljko Ivezic, Thomas R. Quinn, Jo Bovy, Charlotte R. Christensen, Mario Juric, Rok Roskar, Alyson M. Brooks, Fabio Governato

Abstract: We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding 10 kpc and galactocentric distances exceeding 20 kpc. Our analysis utilizes Jeans equations to genera… ▽ More We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding 10 kpc and galactocentric distances exceeding 20 kpc. Our analysis utilizes Jeans equations to generate two-dimensional acceleration maps throughout the volume; this approach is thoroughly tested on a cosmologically derived N-body+SPH simulation of a MW-like galaxy. We show that the known accelerations (gradients of the gravitational potential) can be successfully recovered in such a realistic system. Leveraging the baryonic gravitational potential derived by Bovy & Rix (2013), we show that the gravitational potential implied by the SDSS observations cannot be explained, assuming Newtonian gravity, by visible matter alone: the gravitational force experienced by stars at galactocentric distances of 20 kpc is as much as three times stronger than what can be attributed to purely visible matter. We also show that the SDSS data provide a strong constraint on the shape of the dark matter halo potential. Within galactocentric distances of 20 kpc, the dark matter halo potential is well described as an oblate halo with axis ratio qDM=0.7+/-0.1; this corresponds to an axis ratio qDM=0.4+/-0.1 for the dark matter density distribution. Because of our precise two-dimensional measurements of the acceleration of the halo stars, we can reject several MOND models as an explanation of the observed behavior. △ Less

Submitted 24 November, 2014; v1 submitted 22 August, 2014; originally announced August 2014.

Comments: 15 pages, 21 figures, accepted to ApJ. If viewing using Preview, some figures may appear blurry. If possible, please view using Adobe PDF Viewer

Faint dwarfs as a test of DM models: WDM vs. CDM

Authors: Fabio Governato, Daniel Weisz, Andrew Pontzen, Sarah Loebman, Darren Reed, Alyson M. Brooks, Peter Behroozi, Charlotte Christensen, Piero Madau, Lucio Mayer, Sijing Shen, Matthew Walker, Thomas Quinn, Benjamin W. Keller, James Wadsley

Abstract: We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $Λ$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. In… ▽ More We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $Λ$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra. △ Less

Submitted 23 November, 2014; v1 submitted 30 June, 2014; originally announced July 2014.

Comments: 13 pages including Appendix on Color Magnitude Diagrams. Accepted by MNRAS. Added one plot and details on ChaNGa implementation. Reduced number of citations after editorial request

Journal ref: MNRAS, 2015, 448, 792

Effects of Inclination on Measuring Velocity Dispersion and Implications for Black Holes

Authors: Jillian Bellovary, Kelly Holley-Bockelmann, Kayhan Gültekin, Charlotte Christensen, Fabio Governato, Alyson Brooks, Sarah Loebman, Ferah Munshi

Abstract: The relation of central black hole mass and stellar spheroid velocity dispersion (the M-$σ$ relation) is one of the best-known and tightest correlations linking black holes and their host galaxies. There has been much scrutiny concerning the difficulty of obtaining accurate black hole measurements, and rightly so; however, it has been taken for granted that measurements of velocity dispersion are… ▽ More The relation of central black hole mass and stellar spheroid velocity dispersion (the M-$σ$ relation) is one of the best-known and tightest correlations linking black holes and their host galaxies. There has been much scrutiny concerning the difficulty of obtaining accurate black hole measurements, and rightly so; however, it has been taken for granted that measurements of velocity dispersion are essentially straightforward. We examine five disk galaxies from cosmological SPH simulations and find that line-of-sight effects due to galaxy orientation can affect the measured $σ$ by 30%, and consequently black hole mass predictions by up to 1.0 dex. Face-on orientations correspond to systematically lower velocity dispersion measurements, while more edge-on orientations give higher velocity dispersions, due to contamination by disk stars when measuring line of sight quantities. We caution observers that the uncertainty of velocity dispersion measurements is at least 20 km/s, and can be much larger for moderate inclinations. This effect may account for some of the scatter in the locally measured M-$σ$ relation, particularly at the low-mass end. We provide a method for correcting observed $σ_{\rm los}$ values for inclination effects based on observable quantities. △ Less

Submitted 22 September, 2014; v1 submitted 1 May, 2014; originally announced May 2014.

Comments: 11 pages, 7 figures, replaced with accepted version

A Survey of Satellite Galaxies around NGC 4258

Authors: Meghin Spencer, Sarah Loebman, Peter Yoachim

Abstract: We conduct a survey of satellite galaxies around the nearby spiral NGC 4258 by combining spectroscopic observations from the Apache Point Observatory 3.5-meter telescope with SDSS spectra. New spectroscopy is obtained for 15 galaxies. Of the 47 observed objects, we categorize 8 of them as probable satellites, 8 as possible satellites, and 17 as unlikely to be satellites. We do not speculate on the… ▽ More We conduct a survey of satellite galaxies around the nearby spiral NGC 4258 by combining spectroscopic observations from the Apache Point Observatory 3.5-meter telescope with SDSS spectra. New spectroscopy is obtained for 15 galaxies. Of the 47 observed objects, we categorize 8 of them as probable satellites, 8 as possible satellites, and 17 as unlikely to be satellites. We do not speculate on the membership of the remaining 14 galaxies due to a lack of velocity and distance information. Radially integrating our best fit NFW profile for NGC 4258 yields a total mass of 1.8e12 Msun within 200 kpc. We find that the angular distribution of the satellites appears to be random, and not preferentially aligned with the disk of NGC 4258. In addition, many of the probable satellite galaxies have blue u-r colors and appear to be star-forming irregulars in SDSS images; this stands in contrast to the low number of blue satellites in the Milky Way and M31 systems at comparable distances. △ Less

Submitted 28 April, 2014; originally announced April 2014.

Comments: accepted to ApJ

The Pressure of the Star Forming ISM in Cosmological Simulations

Authors: Ferah Munshi, Charlotte Christensen, Thomas R. Quinn, Fabio Governato, James Wadsley, Sarah Loebman, Sijing Shen

Abstract: We examine the pressure of the star-forming interstellar medium (ISM) of Milky-Way sized disk galaxies using fully cosmological SPH+N-body, high resolution simulations. These simulations include explicit treatment of metal-line cooling in addition to dust and self-shielding, $\mathrm{H_{2}}$ based star formation. The 4 simulated halos have masses ranging from a few times $10^{10}$ to nearly… ▽ More We examine the pressure of the star-forming interstellar medium (ISM) of Milky-Way sized disk galaxies using fully cosmological SPH+N-body, high resolution simulations. These simulations include explicit treatment of metal-line cooling in addition to dust and self-shielding, $\mathrm{H_{2}}$ based star formation. The 4 simulated halos have masses ranging from a few times $10^{10}$ to nearly $10^{12}$ solar masses. Using a kinematic decomposition of these galaxies into present-day bulge and disk components, we find that the typical pressure of the star-forming ISM in the present-day bulge is higher than that in the present-day disk by an order of magnitude. We also find that pressure of the star-forming ISM at high redshift is on average, higher than ISM pressures at low redshift. This explains the why the bulge forms at higher pressures: the disk assembles at lower redshift, when the ISM is lower pressure and the bulge forms at high redshift, when the ISM is at higher pressure. If ISM pressure and IMF variation are tied together as suggested in studies like \cite{Conroy2012}, these results could indicate a time-dependent IMF in Milky-Way like systems, as well as a different IMF in the bulge and the disk. △ Less

Submitted 5 December, 2013; originally announced December 2013.

Comments: Accepted 12/3/2013 to ApJL

Exploring the Variable Sky with LINEAR. III. Classification of Periodic Light Curves

Authors: Lovro Palaversa, Željko Ivezić, Laurent Eyer, Domagoj Ruždjak, Davor Sudar, Mario Galin, Andrea Kroflin, Martina Mesarić, Petra Munk, Dijana Vrbanec, Hrvoje Božić, Sarah Loebman, Branimir Sesar, Lorenzo Rimoldini, Nicholas Hunt-Walker, Jacob VanderPlas, David Westman, J. Scott Stuart, Andrew C. Becker, Gregor Srdoč, Przemyslaw Wozniak, Hakeem Oluseyi

Abstract: We describe the construction of a highly reliable sample of approximately 7,000 optically faint periodic variable stars with light curves obtained by the asteroid survey LINEAR across 10,000 sq.deg of northern sky. Majority of these variables have not been cataloged yet. The sample flux limit is several magnitudes fainter than for most other wide-angle surveys; the photometric errors range from ~0… ▽ More We describe the construction of a highly reliable sample of approximately 7,000 optically faint periodic variable stars with light curves obtained by the asteroid survey LINEAR across 10,000 sq.deg of northern sky. Majority of these variables have not been cataloged yet. The sample flux limit is several magnitudes fainter than for most other wide-angle surveys; the photometric errors range from ~0.03 mag at $r=15$ to ~0.20 mag at r=18. Light curves include on average 250 data points, collected over about a decade. Using SDSS-based photometric recalibration of the LINEAR data for about 25 million objects, we selected ~200,000 most probable candidate variables and visually confirmed and classified approximately 7,000 periodic variables using phased light curves. The reliability and uniformity of visual classification across eight human classifiers was calibrated and tested using a SDSS Stripe 82 region variable star catalog, and verified using an unsupervised machine learning approach. The resulting sample of periodic LINEAR variables is dominated by 3,900 RR Lyrae stars and 2,700 eclipsing binary stars of all subtypes, and includes small fractions of relatively rare populations such as asymptotic giant branch stars and SX Phoenicis stars. We discuss the distribution of these mostly uncataloged variables in various diagrams constructed with optical-to-infrared SDSS, 2MASS and WISE photometry, and with LINEAR light curve features. An interesting side result is a robust and precise quantitative description of a strong correlation between the light-curve period and color/spectral type for close and contact eclipsing binary stars. These large samples of robustly classified variable stars will enable detailed statistical studies of the Galactic structure and physics of binary and other stars, and we make them publicly available. △ Less

Submitted 1 August, 2013; originally announced August 2013.

Comments: 71 pages, 21 figures. For additional machine readable tables and light curves, see http://www.astro.washington.edu/users/ivezic/linear/PaperIII/PLV.html . Accepted by AJ

The Effects of Radial Migration on the Vertical Structure of Galactic Discs

Authors: Rok Roškar, Victor P. Debattista, Sarah R. Loebman

Abstract: We present evidence that isolated growing discs, subject to internal spiral perturbations, thicken due to both heating \emph{and} radial migration. We show this by demonstrating that the thickness and vertical velocity dispersions of coeval stars depend on their age as well as the change in their radii. While the disc thickens due to internal processes, we find that this induces only a minor amoun… ▽ More We present evidence that isolated growing discs, subject to internal spiral perturbations, thicken due to both heating \emph{and} radial migration. We show this by demonstrating that the thickness and vertical velocity dispersions of coeval stars depend on their age as well as the change in their radii. While the disc thickens due to internal processes, we find that this induces only a minor amount of flaring. We further demonstrate the consequences of such thickening on the structural properties of stellar populations and find that they qualitatively agree with recent studies of the Milky Way disc. △ Less

Submitted 9 May, 2013; v1 submitted 8 November, 2012; originally announced November 2012.

Comments: revised version, accepted for publication

Constraints on the Shape of the Milky Way Dark Matter Halo from Jeans Equations Applied to SDSS Data

Authors: Sarah R. Loebman, Zeljko Ivezic, Thomas R. Quinn, Fabio Governato, Alyson M. Brooks, Charlotte R. Christensen, Mario Juric

Abstract: We search for evidence of dark matter in the Milky Way by utilizing the stellar number density distribution and kinematics measured by the Sloan Digital Sky Survey (SDSS) to heliocentric distances exceeding ~10 kpc. We employ the cylindrically symmetric form of Jeans equations and focus on the morphology of the resulting acceleration maps, rather than the normalization of the total mass as done in… ▽ More We search for evidence of dark matter in the Milky Way by utilizing the stellar number density distribution and kinematics measured by the Sloan Digital Sky Survey (SDSS) to heliocentric distances exceeding ~10 kpc. We employ the cylindrically symmetric form of Jeans equations and focus on the morphology of the resulting acceleration maps, rather than the normalization of the total mass as done in previous, mostly local, studies. Jeans equations are first applied to a mock catalog based on a cosmologically derived N-body + SPH simulation, and the known acceleration (gradient of gravitational potential) is successfully recovered. The same simulation is also used to quantify the impact of dark matter on the total acceleration. We use Galfast, a code designed to quantitatively reproduce SDSS measurements and selection effects, to generate a synthetic stellar catalog. We apply Jeans equations to this catalog and produce two-dimensional maps of stellar acceleration. These maps reveal that in a Newtonian framework, the implied gravitational potential cannot be explained by visible matter alone. The acceleration experienced by stars at galactocentric distances of ~20 kpc is three times larger than what can be explained by purely visible matter. The application of an analytic method for estimating the dark matter halo axis ratio to SDSS data implies an oblate halo with q_DM = 0.47 +/- 0.14 within the same distance range. These techniques can be used to map the dark matter halo to much larger distances from the Galactic center using upcoming deep optical surveys, such as LSST. △ Less

Submitted 12 September, 2012; originally announced September 2012.

Comments: accepted to ApJ Letters, six pages, five figures

Reproducing the Stellar Mass/Halo Mass Relation in Simulated LCDM Galaxies: Theory vs Observational Estimates

Authors: Ferah Munshi, F. Governato, A. M. Brooks, C. Christensen, S. Shen, S. Loebman, B. Moster, T. Quinn, J. Wadsley

Abstract: We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for a new sample of simulated field galaxies using fully cosmological, LCDM, high resolution SPH + N-Body simulations.These simulations include an explicit treatment of metal line cooling, dust and self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos have masses… ▽ More We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for a new sample of simulated field galaxies using fully cosmological, LCDM, high resolution SPH + N-Body simulations.These simulations include an explicit treatment of metal line cooling, dust and self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos have masses ranging from a few times 10^8 to nearly 10^12 solar masses. At z=0 our simulated galaxies have a baryon content and morphology typical of field galaxies. Over a stellar mass range of 2.2 x 10^3 to 4.5 x 10^10 solar masses, we find extremely good agreement between the SHM ratio in simulations and the present-day predictions from the statistical Abundance Matching Technique presented in Moster et al. (2012). This improvement over past simulations is due to a number systematic factors, each decreasing the SHM ratios: 1) gas outflows that reduce the overall SF efficiency but allow for the formation of a cold gas component 2) estimating the stellar masses of simulated galaxies using artificial observations and photometric techniques similar to those used in observations and 3) accounting for a systematic, up to 30 percent overestimate in total halo masses in DM-only simulations, due to the neglect of baryon loss over cosmic times. Our analysis suggests that stellar mass estimates based on photometric magnitudes can underestimate the contribution of old stellar populations to the total stellar mass, leading to stellar mass errors of up to 50 percent for individual galaxies. These results highlight the importance of using proper techniques to compare simulations with observations and reduce the perceived tension between the star formation efficiency in galaxy formation models and in real galaxies. △ Less

Submitted 6 September, 2012; originally announced September 2012.

Comments: Submitted to ApJ 9 pages, 5 figures

Milky Way Tomography IV: Dissecting Dust

Authors: Michael Berry, Željko Ivezić, Branimir Sesar, Mario Jurić, Edward F. Schlafly, Jillian Bellovary, Douglas Finkbeiner, Dijana Vrbanec, Timothy C. Beers, Keira J. Brooks, Donald P. Schneider, Robert R. Gibson, Amy Kimball, Lynne Jones, Peter Yoachim, Simon Krughoff, Andrew J. Connolly, Sarah Loebman, Nicholas A. Bond, David Schlegel, Julianne Dalcanton, Brian Yanny, Steven R. Majewski, Gillian R. Knapp, James E. Gunn , et al. (6 additional authors not shown)

Abstract: We use SDSS photometry of 73 million stars to simultaneously obtain best-fit main-sequence stellar energy distribution (SED) and amount of dust extinction along the line of sight towards each star. Using a subsample of 23 million stars with 2MASS photometry, whose addition enables more robust results, we show that SDSS photometry alone is sufficient to break degeneracies between intrinsic stellar… ▽ More We use SDSS photometry of 73 million stars to simultaneously obtain best-fit main-sequence stellar energy distribution (SED) and amount of dust extinction along the line of sight towards each star. Using a subsample of 23 million stars with 2MASS photometry, whose addition enables more robust results, we show that SDSS photometry alone is sufficient to break degeneracies between intrinsic stellar color and dust amount when the shape of extinction curve is fixed. When using both SDSS and 2MASS photometry, the ratio of the total to selective absorption, $R_V$, can be determined with an uncertainty of about 0.1 for most stars in high-extinction regions. These fits enable detailed studies of the dust properties and its spatial distribution, and of the stellar spatial distribution at low Galactic latitudes. Our results are in good agreement with the extinction normalization given by the Schlegel et al. (1998, SFD) dust maps at high northern Galactic latitudes, but indicate that the SFD extinction map appears to be consistently overestimated by about 20% in the southern sky, in agreement with Schlafly et al. (2010). The constraints on the shape of the dust extinction curve across the SDSS and 2MASS bandpasses support the models by Fitzpatrick (1999) and Cardelli et al. (1989). For the latter, we find an $R_V=3.0\pm0.1$(random) $\pm0.1$(systematic) over most of the high-latitude sky. At low Galactic latitudes (|b|<5), we demonstrate that the SFD map cannot be reliably used to correct for extinction as most stars are embedded in dust, rather than behind it. We introduce a method for efficient selection of candidate red giant stars in the disk, dubbed "dusty parallax relation", which utilizes a correlation between distance and the extinction along the line of sight. We make these best-fit parameters, as well as all the input SDSS and 2MASS data, publicly available in a user-friendly format. △ Less

Submitted 21 November, 2011; originally announced November 2011.

Comments: Submitted to ApJ, 55 pages, 37 figures

Implications of Radial Migration for Stellar Population Studies

Authors: R. Roškar, V. P. Debattista, S. R. Loebman, Ž. Ivezić, T. R. Quinn

Abstract: Recent theoretical work suggests that it may be common for stars in the disks of spiral galaxies to migrate radially across significant distances in the disk. Such migrations are a result of resonant scattering with spiral arms and move the guiding centers of the stars while preserving the circularities of their orbits. Migration can therefore efficiently mix stars in all parts of the Galactic dis… ▽ More Recent theoretical work suggests that it may be common for stars in the disks of spiral galaxies to migrate radially across significant distances in the disk. Such migrations are a result of resonant scattering with spiral arms and move the guiding centers of the stars while preserving the circularities of their orbits. Migration can therefore efficiently mix stars in all parts of the Galactic disk. We are rapidly approaching an important confluence of theory and observation, where we may soon be able to uncover signatures of such processes in our own Milky Way. The resolution and robustness of the physical modeling in simulations has improved drastically, while observational datasets are increasing in depth and astrometric accuracy. Here, we discuss the results from our idealized N-body/SPH simulations of disk formation and evolution, emphasizing specifically the observational consequences of stellar migration on the solar neighborhood and the vertical structure of the disk. We demonstrate that radial mixing of stars is a crucial dynamical process that we must try to understand if we are to draw significant conclusions about our Galactic history from the properties of stars in our vicinity. △ Less

Submitted 6 January, 2011; originally announced January 2011.

Comments: Contributed talk proceedings from Cool Stars 16 Meeting

The Genesis of the Milky Way's Thick Disk via Stellar Migration

Authors: Sarah R. Loebman, Rok Roskar, Victor P. Debattista, Zeljko Ivezic, Thomas R. Quinn, James Wadsley

Abstract: The separation of the Milky Way disk into a thin and thick component is supported by differences in the spatial, kinematic and metallicity distributions of their stars. These differences have led to the view that the thick disk formed early via a cataclysmic event and constitutes fossil evidence of the hierarchical growth of the Milky Way. We show here, using N-body simulations, how a double-expon… ▽ More The separation of the Milky Way disk into a thin and thick component is supported by differences in the spatial, kinematic and metallicity distributions of their stars. These differences have led to the view that the thick disk formed early via a cataclysmic event and constitutes fossil evidence of the hierarchical growth of the Milky Way. We show here, using N-body simulations, how a double-exponential vertical structure, with stellar populations displaying similar dichotomies can arise purely through internal evolution. In this picture, stars migrate radially, while retaining nearly circular orbits, as described by Sellwood & Binney (2002). As stars move outwards they populate a thickened component. Such stars found at the present time in the solar neighborhood formed early in the disk's history at smaller radii where stars are more metal-poor and alpha-enhanced, leading to the properties observed for thick disk stars. Classifying stars as members of the thin or thick disk by either velocity or metallicity leads to an apparent separation in the other property as observed. This scenario is supported by the SDSS observation that stars in the transition region do not show any correlation between rotational velocity and metallicity. The good qualitative agreement between our simulation and observations in the Milky Way hint that the thick disk may be a ubiquitous galaxy feature generated by stellar migration. Nonetheless, we cannot exclude that some fraction of the thick disk is a fossil of a past more violent history, nor can this scenario explain thick disks in all galaxies. △ Less

Submitted 29 September, 2010; originally announced September 2010.

SDSS Observations of the Milky Way vs. N-body Models: A Comparison of Stellar Distributions in the Position-Velocity-Metallicity Space

Authors: S. Loebman, R. Roskar, Z. Ivezic

Abstract: The data obtained by the recent modern sky surveys enable detailed studies of the stellar distribution in the multi-dimensional space spanned by spatial coordinates, velocity and metallicity, from the solar neighborhood all the way out to the outer Milky Way halo. While these results represent exciting observational breakthroughs, their interpretation is not simple. For example, traditional deco… ▽ More The data obtained by the recent modern sky surveys enable detailed studies of the stellar distribution in the multi-dimensional space spanned by spatial coordinates, velocity and metallicity, from the solar neighborhood all the way out to the outer Milky Way halo. While these results represent exciting observational breakthroughs, their interpretation is not simple. For example, traditional decomposition of the thin and thick disks predicts a strong correlation in metallicity and kinematics at $\sim$1 kpc from the Galactic plane; however, recent SDSS--based work has demonstrated an absence of this correlation for disk stars. Instead, the variation of the metallicity and rotational velocity distributions can be modeled using non--Gaussian functions that retain their shapes and only shift as the distance from the mid--plane increases. To fully contextualize these recent observational results, a detailed comparison with sophisticated numerical models is necessary. Modern simulations have sufficient resolution and physical detail to study the formation of stellar disks and spheroids over a large baseline of masses and cosmic ages. We discuss preliminary comparisons of various observed maps and N--body model predictions and find them encouraging. In particular, the N--body disk models of Roškar et al. \cite{Roskar 2008} reproduce a change of disk scale height reminiscent of thin/thick disk decomposition, as well as metallicity and rotational velocity gradients, while not inducing a correlation of the latter two quantities, in qualitative agreement with SDSS observations. △ Less

Submitted 30 October, 2008; v1 submitted 28 October, 2008; originally announced October 2008.

Comments: Presented at the "Classification and Discovery in Large Astronomical Surveys" meeting, Ringberg Castle, 14-17 October, 2008