-
The OSIRIS Lens-Amplified Survey (OLAS) I: Dynamical Effects of Stellar Feedback in Low Mass Galaxies at z ~ 2
Authors:
Jessie Hirtenstein,
Tucker Jones,
Xin Wang,
Andrew Wetzel,
Kareem El-Badry,
Austin Hoag,
Tommaso Treu,
Marusa Bradac,
Takahiro Morishita
Abstract:
We introduce the OSIRIS Lens-Amplified Survey (OLAS), a kinematic survey of gravitationally lensed galaxies at cosmic noon taken with Keck adaptive optics. In this paper we present spatially resolved spectroscopy and nebular emission kinematic maps for 17 star forming galaxies with stellar masses 8 < log($M_*$/$M_{\odot}$) < 9.8 and redshifts 1.2 < z < 2.3. OLAS is designed to probe the stellar ma…
▽ More
We introduce the OSIRIS Lens-Amplified Survey (OLAS), a kinematic survey of gravitationally lensed galaxies at cosmic noon taken with Keck adaptive optics. In this paper we present spatially resolved spectroscopy and nebular emission kinematic maps for 17 star forming galaxies with stellar masses 8 < log($M_*$/$M_{\odot}$) < 9.8 and redshifts 1.2 < z < 2.3. OLAS is designed to probe the stellar mass ($M_*$) and specific star formation rate (sSFR) range where simulations suggest that stellar feedback is most effective at driving gaseous outflows that create galaxy-wide potential fluctuations which can generate dark matter cores. We compare our kinematic data with the trend between sSFR, $M_*$ and H$α$ velocity dispersion, $σ$, from the Feedback In Realistic Environments (FIRE) simulations. Our observations reveal a correlation between sSFR and sigma at fixed $M_*$ that is similar to the trend predicted by simulations: feedback from star formation drives star-forming gas and newly formed stars into more dispersion dominated orbits. The observed magnitude of this effect is in good agreement with the FIRE simulations, in which feedback alters the central density profiles of low mass galaxies, converting dark matter cusps into cores over time. Our data support the scenario that stellar feedback drives gaseous outflows and potential fluctuations, which in turn drive dark matter core formation in dwarf galaxies.
△ Less
Submitted 18 April, 2019; v1 submitted 28 November, 2018;
originally announced November 2018.
-
The Origins of the Circumgalactic Medium in the FIRE Simulations
Authors:
Z. Hafen,
C. -A. Faucher-Giguere,
D. Angles-Alcazar,
J. Stern,
D. Keres,
C. Hummels,
C. Esmerian,
S. Garrison-Kimmel,
K. El-Badry,
A. Wetzel,
T. K. Chan,
P. F. Hopkins,
N. Murray
Abstract:
We use a particle tracking analysis to study the origins of the circumgalactic medium (CGM), separating it into (1) accretion from the intergalactic medium (IGM), (2) wind from the central galaxy, and (3) gas ejected from other galaxies. Our sample consists of 21 FIRE-2 simulations, spanning the halo mass range log(Mh/Msun) ~ 10-12 , and we focus on z=0.25 and z=2. Owing to strong stellar feedback…
▽ More
We use a particle tracking analysis to study the origins of the circumgalactic medium (CGM), separating it into (1) accretion from the intergalactic medium (IGM), (2) wind from the central galaxy, and (3) gas ejected from other galaxies. Our sample consists of 21 FIRE-2 simulations, spanning the halo mass range log(Mh/Msun) ~ 10-12 , and we focus on z=0.25 and z=2. Owing to strong stellar feedback, only ~L* halos retain a baryon mass >~50% of their cosmic budget. Metals are more efficiently retained by halos, with a retention fraction >~50%. Across all masses and redshifts analyzed >~60% of the CGM mass originates as IGM accretion (some of which is associated with infalling halos). Overall, the second most important contribution is wind from the central galaxy, though gas ejected or stripped from satellites can contribute a comparable mass in ~L* halos. Gas can persist in the CGM for billions of years, resulting in well-mixed halo gas. Sight lines through the CGM are therefore likely to intersect gas of multiple origins. For low-redshift ~L* halos, cool gas (T<10^4.7 K) is distributed on average preferentially along the galaxy plane, however with strong halo-to-halo variability. The metallicity of IGM accretion is systematically lower than the metallicity of winds (typically by >~1 dex), although CGM and IGM metallicities depend significantly on the treatment of subgrid metal diffusion. Our results highlight the multiple physical mechanisms that contribute to the CGM and will inform observational efforts to develop a cohesive picture.
△ Less
Submitted 24 June, 2019; v1 submitted 28 November, 2018;
originally announced November 2018.
-
The geometric challenge of testing gravity with wide binaries
Authors:
Kareem El-Badry
Abstract:
Wide binaries provide promising laboratories for testing general relativity (GR) in the low-acceleration regime. Recent observational studies have found that the difference in the proper motions and/or radial velocities of the components of nearby wide binaries appear larger than predicted by Kepler's law's, indicating a potential breakdown of GR at low accelerations. These studies have not accoun…
▽ More
Wide binaries provide promising laboratories for testing general relativity (GR) in the low-acceleration regime. Recent observational studies have found that the difference in the proper motions and/or radial velocities of the components of nearby wide binaries appear larger than predicted by Kepler's law's, indicating a potential breakdown of GR at low accelerations. These studies have not accounted for projection effects owing to the different position of the two stars on the celestial sphere. I show that two stars in a wide binary with identical 3D space velocities often have significantly different proper motions and radial velocities purely due to projection effects. I construct a sample of simulated binaries that follow Kepler's laws and have similar phase-space distributions to the observed samples of nearby binaries. Beyond separations of $\sim$0.1 pc, direct comparison of the components' proper motions would suggest strong tensions with GR, even though the simulated binaries follow Kepler's laws by construction. The magnitude of the apparent disagreement is similar to that found observationally, suggesting that the apparent tension between observations and GR may largely be due to projection effects. I discuss prospects for constraining gravity at low accelerations with wide binaries. Robust tests of GR are possible with current data but require measurements of 3D velocities. Further work is also needed to model contamination from unbound moving groups and unrecognized hierarchical triples.
△ Less
Submitted 13 November, 2018; v1 submitted 31 October, 2018;
originally announced October 2018.
-
The wide binary fraction of solar-type stars: emergence of metallicity dependence at a < 200 AU
Authors:
Kareem El-Badry,
Hans-Walter Rix
Abstract:
We combine a catalog of wide binaries constructed from Gaia DR2 with [Fe/H] abundances from wide-field spectroscopic surveys to quantify how the binary fraction varies with metallicity over separations $50 \lesssim s/{\rm AU} \lesssim 50,000$. At a given distance, the completeness of the catalog is independent of metallicity, making it straightforward to constrain intrinsic variation with [Fe/H].…
▽ More
We combine a catalog of wide binaries constructed from Gaia DR2 with [Fe/H] abundances from wide-field spectroscopic surveys to quantify how the binary fraction varies with metallicity over separations $50 \lesssim s/{\rm AU} \lesssim 50,000$. At a given distance, the completeness of the catalog is independent of metallicity, making it straightforward to constrain intrinsic variation with [Fe/H]. The wide binary fraction is basically constant with [Fe/H] at large separations ($s \gtrsim 250$\,AU) but becomes quite rapidly anti-correlated with [Fe/H] at smaller separations: for $50 < s/{\rm AU} < 100$, the binary fraction at $\rm [Fe/H] = -1$ exceeds that at $\rm [Fe/H] = 0.5$ by a factor of 3, an anti-correlation almost as strong as that found for close binaries with $a < 10$ AU. Interpreted in terms of models where disk fragmentation is more efficient at low [Fe/H], our results suggest that $100 < a/{\rm AU} < 200$ is the separation below which a significant fraction of binaries formed via fragmentation of individual gravitationally unstable disks rather than through turbulent core fragmentation. We provide a public catalog of 8,407 binaries within 200 pc with spectroscopically-determined [Fe/H] for at least one component.
△ Less
Submitted 30 October, 2018; v1 submitted 18 September, 2018;
originally announced September 2018.
-
A Closer look at Bursty Star Formation with $L_{Hα}$ and $L_{UV}$ Distributions
Authors:
Najmeh Emami,
Brian Siana,
Daniel R. Weisz,
Benjamin D. Johnson,
Xiangcheng Ma,
Kareem El-Badry
Abstract:
We investigate the bursty star formation histories (SFHs) of dwarf galaxies using the distribution of log($L_{Hα}/L_{UV}$) of 185 local galaxies. We expand on the work of Weisz et al. 2012 to consider a wider range of SFHs and stellar metallicities, and show that there are large degeneracies in a periodic, top-hat burst model. We argue that all galaxies of a given mass have similar SFHs and we can…
▽ More
We investigate the bursty star formation histories (SFHs) of dwarf galaxies using the distribution of log($L_{Hα}/L_{UV}$) of 185 local galaxies. We expand on the work of Weisz et al. 2012 to consider a wider range of SFHs and stellar metallicities, and show that there are large degeneracies in a periodic, top-hat burst model. We argue that all galaxies of a given mass have similar SFHs and we can therefore include the $L_{Hα}$ distributions (subtracting the median trend with stellar mass, referred to as $Δ\text{log}(L_{Hα})$) in our analyses. $Δ\text{log}(L_{Hα})$ traces the amplitude of the bursts, and log($L_{Hα}/L_{UV}$) is a function of timescale, amplitude, and shape of the bursts. We examine the 2-dimensional distribution of these two indicators constrain the SFHs. We use exponentially rising/falling bursts to determine timescales ($e$-folding time, $τ$). We find that galaxies below $10^{7.5}$ M$_{\odot}$ undergo large (amplitudes of $\sim 100$) and rapid ($τ< 30$ Myr) bursts, while galaxies above $10^{8.5}$ M$_{\odot}$ experience smaller (maximum amplitudes $\sim 10$), slower ($τ\gtrsim 300$ Myr) bursts. We compare to the FIRE-2 hydrodynamical simulations and find that the burst amplitudes agree with observations, but they are too rapid in more massive galaxies ($M_* > 10^8$ M$_{\odot}$). Finally, we confirm that stochastic sampling of the stellar mass function can not reproduce the observed distributions unless the standard assumptions of cluster and stellar mass functions are changed. With the next generation of telescopes, measurements of $L_{UV}$ and $L_{Hα}$ will become available for dwarf galaxies at high-redshift, enabling similar analyses of galaxies in the early universe.
△ Less
Submitted 14 May, 2019; v1 submitted 17 September, 2018;
originally announced September 2018.
-
Imprints of white dwarf recoil in the separation distribution of Gaia wide binaries
Authors:
Kareem El-Badry,
Hans-Walter Rix
Abstract:
We construct from Gaia DR2 an extensive and very pure ($\lesssim 0.2\%$ contamination) catalog of wide binaries containing main-sequence (MS) and white dwarf (WD) components within 200 pc of the Sun. The public catalog contains, after removal of clusters and resolved higher-order multiples, $>$50,000 MS/MS, $>$3,000 WD/MS, and nearly 400 WD/WD binaries with projected separations of…
▽ More
We construct from Gaia DR2 an extensive and very pure ($\lesssim 0.2\%$ contamination) catalog of wide binaries containing main-sequence (MS) and white dwarf (WD) components within 200 pc of the Sun. The public catalog contains, after removal of clusters and resolved higher-order multiples, $>$50,000 MS/MS, $>$3,000 WD/MS, and nearly 400 WD/WD binaries with projected separations of $50 \lesssim s/{\rm AU} < 50,000$. Accounting for incompleteness and selection effects, we model the separation distribution of each class of binaries as a broken power-law, revealing marked differences between the three populations. The separation distribution of MS/MS systems is nearly consistent with a single power-law of slope $-1.6$ over at least $500 < s/{\rm AU} < 50,000$, with marginal steepening at $s > 10,000$ AU. In contrast, the separation distributions of WD/MS and WD/WD binaries show distinct breaks at $\sim$ 3,000 AU and $\sim$1,500 AU, respectively: they are flatter than the MS/MS distribution at small separations and steeper at large separations. Using binary population synthesis models, we show that these breaks are unlikely to be caused by external factors but can be explained if the WDs incur a kick of $\sim$ 0.75 km s$^{-1}$ during their formation, presumably due to asymmetric mass loss. The data rule out typical kick velocities above 2km s$^{-1}$. Our results imply that most wide binaries with separations exceeding a few thousand AU become unbound during post-MS evolution.
△ Less
Submitted 8 August, 2018; v1 submitted 16 July, 2018;
originally announced July 2018.
-
The Local Group on FIRE: Dwarf galaxy populations across a suite of hydrodynamic simulations
Authors:
Shea Garrison-Kimmel,
Philip F. Hopkins,
Andrew Wetzel,
James S. Bullock,
Michael Boylan-Kolchin,
Dusan Keres,
Claude-Andre Faucher-Giguere,
Kareem El-Badry,
Astrid Lamberts,
Eliot Quataert,
Robyn Sanderson
Abstract:
We present a new set of high-resolution hydrodynamic cosmological zoom-in simulations that apply the Feedback In Realistic Environments (FIRE) physics to both Local Group (LG)-like and isolated Milky Way (MW)-like volumes (ten host systems in total with baryonic particle mass $\simeq 3,500-7,000\,M_\odot$). We study the stellar mass functions, circular velocity or mass profiles, and velocity dispe…
▽ More
We present a new set of high-resolution hydrodynamic cosmological zoom-in simulations that apply the Feedback In Realistic Environments (FIRE) physics to both Local Group (LG)-like and isolated Milky Way (MW)-like volumes (ten host systems in total with baryonic particle mass $\simeq 3,500-7,000\,M_\odot$). We study the stellar mass functions, circular velocity or mass profiles, and velocity dispersions of the dwarf galaxy populations. The simulations reproduce the stellar mass function and central densities of MW satellite dwarfs for $M_\ast \geq10^{5.5}\,M_\odot$ and predict the existence of $\sim3$ unidentified galaxies with $M_\ast\sim10^5\,M_\odot$ within $300$ kpc of the MW. Overall, we find no evidence for the classical missing satellites or too-big-to-fail (TBTF) problems for satellite galaxies in our sample. Among the satellites, TBTF is resolved primarily by subhalo disruption and overall mass loss; central density profiles of subhalos are of secondary importance. For non-satellite galaxies, our LG-like simulations predict as many as $\sim10$ as-of-yet unseen galaxies at distances $0.3-1$ Mpc from both hosts, with $M_\ast\simeq10^{5-6}\,M_\odot$ (in halos with $V_\mathrm{max}\sim 20~\mathrm{km\,s}^{-1}$), albeit with large halo-to-halo variance. None of our simulations produces a compact, baryon-dominated, high-density dwarf elliptical-type galaxy (with $V_\mathrm{circ} \gtrsim 35~\mathrm{km\,s}^{-1}$ at $r<1\,$kpc), of which six may appear in the LG (but none in the MW). It may therefore remain a challenge to reproduce the full diversity of the dwarf population, including both the highest and lowest density systems.
△ Less
Submitted 11 June, 2018;
originally announced June 2018.
-
Formation, vertex deviation and age of the Milky Way's bulge: input from a cosmological simulation with a late-forming bar
Authors:
Victor P. Debattista,
Oscar A. Gonzalez,
Robyn E. Sanderson,
Kareem El-Badry,
Shea Garrison-Kimmel,
Andrew Wetzel,
Claude-André Faucher-Giguère,
Philip F. Hopkins
Abstract:
We present the late-time evolution of m12m, a cosmological simulation of a Milky Way-like galaxy from the FIRE project. The simulation forms a bar after redshift z = 0.2. We show that the evolution of the model exhibits behaviours typical of kinematic fractionation, with a bar weaker in older populations, an X-shape traced by the younger, metal-rich populations and a prominent X-shape in the edge-…
▽ More
We present the late-time evolution of m12m, a cosmological simulation of a Milky Way-like galaxy from the FIRE project. The simulation forms a bar after redshift z = 0.2. We show that the evolution of the model exhibits behaviours typical of kinematic fractionation, with a bar weaker in older populations, an X-shape traced by the younger, metal-rich populations and a prominent X-shape in the edge-on mean metallicity map. Because of the late formation of the bar in m12m, stars forming after 10 Gyr (z = 0.34) significantly contaminate the bulge, at a level higher than is observed at high latitudes in the Milky Way, implying that its bar cannot have formed as late as in m12m. We also study the model's vertex deviation of the velocity ellipsoid as a function of stellar metallicity and age in the equivalent of Baade's Window. The formation of the bar leads to a non-zero vertex deviation. We find that metal-rich stars have a large vertex deviation (~ 40 degrees), which becomes negligible for metal-poor stars, a trend also found in the Milky Way, despite not matching in detail. We demonstrate that the vertex deviation also varies with stellar age and is large for stars as old as 9 Gyr, while 13 Gyr old stars have negligible vertex deviation. When we exclude stars that have been accreted, the vertex deviation is not significantly changed, demonstrating that the observed variation of vertex deviation with metallicity is not necessarily due to an accreted population.
△ Less
Submitted 9 April, 2019; v1 submitted 30 May, 2018;
originally announced May 2018.
-
An Empirical Measurement of the Initial-Final Mass Relation with Gaia White Dwarfs
Authors:
Kareem El-Badry,
Hans-Walter Rix,
Daniel R. Weisz
Abstract:
We use data from Gaia DR2 to constrain the initial-final mass relation (IFMR) for field stars with initial masses $0.9 \lesssim m_{\rm in}/M_{\odot} \lesssim 8$. Precise parallaxes have revealed unprecedented substructure in the white dwarf (WD) cooling sequence on the color-magnitude diagram (CMD). Some of this substructure stems from the diversity of WD atmospheric compositions, but the CMD rema…
▽ More
We use data from Gaia DR2 to constrain the initial-final mass relation (IFMR) for field stars with initial masses $0.9 \lesssim m_{\rm in}/M_{\odot} \lesssim 8$. Precise parallaxes have revealed unprecedented substructure in the white dwarf (WD) cooling sequence on the color-magnitude diagram (CMD). Some of this substructure stems from the diversity of WD atmospheric compositions, but the CMD remains bimodal even when only spectroscopically-confirmed DA WDs are considered. We develop a generative model to predict the CMD for DA WDs as a function of the initial mass function (IMF), stellar age distribution, and a flexibly parameterized IFMR. We then fit the CMD of 1100 bright DA WDs within 100 pc, for which atmospheric composition and completeness are well-understood. The resulting best-fit IFMR flattens at $3.5 \lesssim m_{\rm in}/M_{\odot}\lesssim 5.5$, producing a secondary peak in the WD mass distribution at $m_{\rm WD} \sim 0.8 M_{\odot}$. Our IFMR is broadly consistent with weaker constraints obtained from binaries and star clusters in previous work but represents the clearest observational evidence obtained to date of theoretically-predicted non-linearity in the IFMR. A visibly bimodal CMD is only predicted for mixed-age stellar populations: in single-age clusters, more massive WDs reach the bottom of the cooling sequence before the first lower-mass WDs appear. This may explain why bimodal cooling sequences have thus far evaded detection in cluster CMDs.
△ Less
Submitted 5 June, 2018; v1 submitted 15 May, 2018;
originally announced May 2018.
-
The formation and hierarchical assembly of globular cluster populations
Authors:
Kareem El-Badry,
Eliot Quataert,
Daniel R. Weisz,
Nick Choksi,
Michael Boylan-Kolchin
Abstract:
We use a semi-analytic model for globular cluster (GC) formation built on dark matter merger trees to explore the relative role of formation physics and hierarchical assembly in determining the properties of GC populations. Many previous works have argued that the observed linear relation between total GC mass and halo mass points to a fundamental GC -- dark matter connection or indicates that GCs…
▽ More
We use a semi-analytic model for globular cluster (GC) formation built on dark matter merger trees to explore the relative role of formation physics and hierarchical assembly in determining the properties of GC populations. Many previous works have argued that the observed linear relation between total GC mass and halo mass points to a fundamental GC -- dark matter connection or indicates that GCs formed at very high redshift before feedback processes introduced nonlinearity in the baryon-to-dark matter mass relation. We demonstrate that at $M_{\rm vir}(z=0) \gtrsim 10^{11.5} M_{\odot}$, a constant ratio between halo mass and total GC mass is in fact an almost inevitable consequence of hierarchical assembly: by the central limit theorem, it is expected at $z=0$ independent of the GC-to-halo mass relation at the time of GC formation. The GC-to-halo mass relation at $M_{\rm vir}(z=0) < 10^{11.5} M_{\odot}$ is more sensitive to the details of the GC formation process. In our fiducial model, GC formation occurs in galaxies when the gas surface density exceeds a critical value. This model naturally predicts bimodal GC color distributions similar to those observed in nearby galaxies and reproduces the observed relation between GC system metallicity and halo mass. It predicts that the cosmic GC formation rate peaked at $z$ $\sim$ 4, too late for GCs to contribute significantly to the UV luminosity density during reionization.
△ Less
Submitted 2 November, 2018; v1 submitted 9 May, 2018;
originally announced May 2018.
-
Where are the most ancient stars in the Milky Way?
Authors:
Kareem El-Badry,
Joss Bland-Hawthorn,
Andrew Wetzel,
Eliot Quataert,
Daniel R. Weisz,
Michael Boylan-Kolchin,
Philip F. Hopkins,
Claude-André Faucher-Giguère,
Dušan Kereš,
Shea Garrison-Kimmel
Abstract:
The oldest stars in the Milky Way (MW) bear imprints of the Galaxy's early assembly history. We use FIRE cosmological zoom-in simulations of three MW-mass disk galaxies to study the spatial distribution, chemistry, and kinematics of the oldest surviving stars ($z_{\rm form} \gtrsim 5$) in MW-like galaxies. We predict the oldest stars to be less centrally concentrated at $z=0$ than stars formed at…
▽ More
The oldest stars in the Milky Way (MW) bear imprints of the Galaxy's early assembly history. We use FIRE cosmological zoom-in simulations of three MW-mass disk galaxies to study the spatial distribution, chemistry, and kinematics of the oldest surviving stars ($z_{\rm form} \gtrsim 5$) in MW-like galaxies. We predict the oldest stars to be less centrally concentrated at $z=0$ than stars formed at later times as a result of two processes. First, the majority of the oldest stars are not formed $\textit{in situ}$ but are accreted during hierarchical assembly. These $\textit{ex situ}$ stars are deposited on dispersion-supported, halo-like orbits but dominate over old stars formed $\textit{in situ}$ in the solar neighborhood, and in some simulations, even in the galactic center. Secondly, old stars formed $\textit{in situ}$ are driven outwards by bursty star formation and energetic feedback processes that create a time-varying gravitational potential at $z\gtrsim 2$, similar to the process that creates dark matter cores and expands stellar orbits in bursty dwarf galaxies. The total fraction of stars that are ancient is more than an order of magnitude higher for sight lines $\textit{away}$ from the bulge and inner halo than for inward-looking sight lines. Although the task of identifying specific stars as ancient remains challenging, we anticipate that million-star spectral surveys and photometric surveys targeting metal-poor stars already include hundreds of stars formed before $z=5$. We predict most of these targets to have higher metallicity ($-3 < \rm [Fe/H] < -2$) than the most extreme metal-poor stars.
△ Less
Submitted 9 July, 2018; v1 submitted 2 April, 2018;
originally announced April 2018.
-
No Assembly Required: Mergers are Mostly Irrelevant for the Growth of Low-mass Dwarf Galaxies
Authors:
Alex Fitts,
Michael Boylan-Kolchin,
James S. Bullock,
Daniel R. Weisz,
Kareem El-Badry,
Coral Wheeler,
Claude-André Faucher-Giguère,
Eliot Quataert,
Philip F. Hopkins,
Dušan Kereš,
Andrew Wetzel,
Chris Hayward
Abstract:
We investigate the merger histories of isolated dwarf galaxies based on a suite of 15 high-resolution cosmological zoom-in simulations, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_{\odot}$ (and M$_\star\sim10^5-10^7\,{\rm M}_{\odot}$) at $z=0$, from the Feedback in Realistic Environments (FIRE) project. The stellar populations of these dwarf galaxies at $z=0$ are formed essentially e…
▽ More
We investigate the merger histories of isolated dwarf galaxies based on a suite of 15 high-resolution cosmological zoom-in simulations, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_{\odot}$ (and M$_\star\sim10^5-10^7\,{\rm M}_{\odot}$) at $z=0$, from the Feedback in Realistic Environments (FIRE) project. The stellar populations of these dwarf galaxies at $z=0$ are formed essentially entirely "in situ": over 90$\%$ of the stellar mass is formed in the main progenitor in all but two cases, and all 15 of the galaxies have >70$\%$ of their stellar mass formed in situ. Virtually all galaxy mergers occur prior to $z\sim3$, meaning that accreted stellar populations are ancient. On average, our simulated dwarfs undergo 5 galaxy mergers in their lifetimes, with typical pre-merger galaxy mass ratios that are less than 1:10. This merger frequency is generally comparable to what has been found in dissipationless simulations when coupled with abundance matching. Two of the simulated dwarfs have a luminous satellite companion at $z=0$. These ultra-faint dwarfs lie at or below current detectability thresholds but are intriguing targets for next-generation facilities. The small contribution of accreted stars make it extremely difficult to discern the effects of mergers in the vast majority of dwarfs either photometrically or using resolved-star color-magnitude diagrams (CMDs). The important implication for near-field cosmology is that star formation histories of comparably massive galaxies derived from resolved CMDs should trace the build-up of stellar mass in one main system across cosmic time as opposed to reflecting the contributions of many individual star formation histories of merged dwarfs.
△ Less
Submitted 18 January, 2018;
originally announced January 2018.
-
Gas Kinematics in FIRE Simulated Galaxies Compared to Spatially Unresolved HI Observations
Authors:
Kareem El-Badry,
Jeremy Bradford,
Eliot Quataert,
Marla Geha,
Michael Boylan-Kolchin,
Daniel R. Weisz,
Andrew Wetzel,
Philip F. Hopkins,
T. K. Chan,
Alex Fitts,
Dušan Kereš,
Claude-André Faucher-Giguère
Abstract:
The shape of a galaxy's spatially unresolved, globally integrated 21-cm emission line depends on its internal gas kinematics: galaxies with rotation-supported gas disks produce double-horned profiles with steep wings, while galaxies with dispersion-supported gas produce Gaussian-like profiles with sloped wings. Using mock observations of simulated galaxies from the FIRE project, we show that one c…
▽ More
The shape of a galaxy's spatially unresolved, globally integrated 21-cm emission line depends on its internal gas kinematics: galaxies with rotation-supported gas disks produce double-horned profiles with steep wings, while galaxies with dispersion-supported gas produce Gaussian-like profiles with sloped wings. Using mock observations of simulated galaxies from the FIRE project, we show that one can therefore constrain a galaxy's gas kinematics from its unresolved 21-cm line profile. In particular, we find that the kurtosis of the 21-cm line increases with decreasing $V/σ$, and that this trend is robust across a wide range of masses, signal-to-noise ratios, and inclinations. We then quantify the shapes of 21-cm line profiles from a morphologically unbiased sample of $\sim$2000 low-redshift, HI-detected galaxies with $M_{\rm star} = 10^{7-11} M_{\odot}$ and compare to the simulated galaxies. At $M_{\rm star} \gtrsim 10^{10} M_{\odot}$, both the observed and simulated galaxies produce double-horned profiles with low kurtosis and steep wings, consistent with rotation-supported disks. Both the observed and simulated line profiles become more Gaussian-like (higher kurtosis and less-steep wings) at lower masses, indicating increased dispersion support. However, the simulated galaxies transition from rotation to dispersion support more strongly: at $M_{\rm star} = 10^{8-10}M_{\odot}$, most of the simulations produce more Gaussian-like profiles than typical observed galaxies with similar mass, indicating that gas in the low-mass simulated galaxies is, on average, overly dispersion-supported. Most of the lower-mass simulated galaxies also have somewhat lower gas fractions than the median of the observed population. The simulations nevertheless reproduce the observed line-width baryonic Tully-Fisher relation, which is insensitive to rotation vs. dispersion support.
△ Less
Submitted 14 March, 2018; v1 submitted 11 January, 2018;
originally announced January 2018.
-
The origin of the diverse morphologies and kinematics of Milky Way-mass galaxies in the FIRE-2 simulations
Authors:
Shea Garrison-Kimmel,
Philip F. Hopkins,
Andrew Wetzel,
Kareem El-Badry,
Robyn E. Sanderson,
James S. Bullock,
Xiangcheng Ma,
Freeke van de Voort,
Zachary Hafen,
Claude-André Faucher-Giguère,
Christopher C. Hayward,
Eliot Quataert,
Dusan Keres,
Michael Boylan-Kolchin
Abstract:
We use hydrodynamic cosmological zoom-in simulations from the FIRE project to explore the morphologies and kinematics of fifteen Milky Way (MW)-mass galaxies. Our sample ranges from compact, bulge-dominated systems with 90% of their stellar mass within 2.5 kpc to well-ordered disks that reach $\gtrsim15$ kpc. The gas in our galaxies always forms a thin, rotation-supported disk at $z=0$, with sizes…
▽ More
We use hydrodynamic cosmological zoom-in simulations from the FIRE project to explore the morphologies and kinematics of fifteen Milky Way (MW)-mass galaxies. Our sample ranges from compact, bulge-dominated systems with 90% of their stellar mass within 2.5 kpc to well-ordered disks that reach $\gtrsim15$ kpc. The gas in our galaxies always forms a thin, rotation-supported disk at $z=0$, with sizes primarily determined by the gas mass. For stars, we quantify kinematics and morphology both via the fraction of stars on disk-like orbits and with the radial extent of the stellar disk. In this mass range, stellar morphology and kinematics are poorly correlated with the properties of the halo available from dark matter-only simulations (halo merger history, spin, or formation time). They more strongly correlate with the gaseous histories of the galaxies: those that maintain a high gas mass in the disk after $z\sim1$ develop well-ordered stellar disks. The best predictor of morphology we identify is the spin of the gas in the halo at the time the galaxy formed 1/2 of its stars (i.e. the gas that builds the galaxy). High-$z$ mergers, before a hot halo emerges, produce some of the most massive bulges in the sample (from compact disks in gas-rich mergers), while later-forming bulges typically originate from internal processes, as satellites are stripped of gas before the galaxies merge. Moreover, most stars in $z=0$ MW-mass galaxies (even $z=0$ bulge stars) form in a disk: $\gtrsim$60-90% of stars begin their lives rotationally supported.
△ Less
Submitted 11 December, 2017;
originally announced December 2017.
-
Discovery and Characterization of 3000+ Main-Sequence Binaries from APOGEE Spectra
Authors:
Kareem El-Badry,
Yuan-Sen Ting,
Hans-Walter Rix,
Eliot Quataert,
Daniel R. Weisz,
Phillip Cargile,
Charlie Conroy,
David W. Hogg,
Maria Bergemann,
Chao Liu
Abstract:
We develop a data-driven spectral model for identifying and characterizing spatially unresolved multiple-star systems and apply it to APOGEE DR13 spectra of main-sequence stars. Binaries and triples are identified as targets whose spectra can be significantly better fit by a superposition of two or three model spectra, drawn from the same isochrone, than any single-star model. From an initial samp…
▽ More
We develop a data-driven spectral model for identifying and characterizing spatially unresolved multiple-star systems and apply it to APOGEE DR13 spectra of main-sequence stars. Binaries and triples are identified as targets whose spectra can be significantly better fit by a superposition of two or three model spectra, drawn from the same isochrone, than any single-star model. From an initial sample of $\sim$20,000 main-sequence targets, we identify $\sim$2,500 binaries in which both the primary and secondary star contribute detectably to the spectrum, simultaneously fitting for the velocities and stellar parameters of both components. We additionally identify and fit $\sim$200 triple systems, as well as $\sim$700 velocity-variable systems in which the secondary does not contribute detectably to the spectrum. Our model simplifies the process of simultaneously fitting single- or multi-epoch spectra with composite models and does not depend on a velocity offset between the two components of a binary, making it sensitive to traditionally undetectable systems with periods of hundreds or thousands of years. In agreement with conventional expectations, almost all the spectrally-identified binaries with measured parallaxes fall above the main sequence in the color-magnitude diagram. We find excellent agreement between spectrally and dynamically inferred mass ratios for the $\sim$600 binaries in which a dynamical mass ratio can be measured from multi-epoch radial velocities. We obtain full orbital solutions for 64 systems, including 14 close binaries within hierarchical triples. We make available catalogs of stellar parameters, abundances, mass ratios, and orbital parameters.
△ Less
Submitted 26 January, 2018; v1 submitted 23 November, 2017;
originally announced November 2017.
-
The origin of ultra diffuse galaxies: stellar feedback and quenching
Authors:
T. K. Chan,
D. Kereš,
A. Wetzel,
P. F. Hopkins,
C. -A. Faucher-Giguère,
K. El-Badry,
S. Garrison-Kimmel,
M. Boylan-Kolchin
Abstract:
We test if the cosmological zoom-in simulations of isolated galaxies from the FIRE project reproduce the properties of ultra diffuse galaxies. We show that stellar feedback-generated outflows that dynamically heat galactic stars, together with a passively aging stellar population after imposed quenching (from e.g. infall into a galaxy cluster), naturally reproduce the observed population of red UD…
▽ More
We test if the cosmological zoom-in simulations of isolated galaxies from the FIRE project reproduce the properties of ultra diffuse galaxies. We show that stellar feedback-generated outflows that dynamically heat galactic stars, together with a passively aging stellar population after imposed quenching (from e.g. infall into a galaxy cluster), naturally reproduce the observed population of red UDGs, without the need for high spin halos or dynamical influence from their host cluster. We reproduce the range of surface brightness, radius and absolute magnitude of the observed z=0 red UDGs by quenching simulated galaxies at a range of different times. They represent a mostly uniform population of dark matter-dominated galaxies with M_star ~1e8 Msun, low metallicity and a broad range of ages. The most massive simulated UDGs require earliest quenching and are therefore the oldest. Our simulations provide a good match to the central enclosed masses and the velocity dispersions of the observed UDGs (20-50 km/s). The enclosed masses of the simulated UDGs remain largely fixed across a broad range of quenching times because the central regions of their dark matter halos complete their growth early. A typical UDG forms in a dwarf halo mass range of Mh~4e10-1e11 Msun. The most massive red UDG in our sample requires quenching at z~3 when its halo reached Mh ~ 1e11 Msun. If it, instead, continues growing in the field, by z=0 its halo mass reaches > 5e11 Msun, comparable to the halo of an L* galaxy. If our simulated dwarfs are not quenched, they evolve into bluer low-surface brightness galaxies with mass-to-light ratios similar to observed field dwarfs. While our simulation sample covers a limited range of formation histories and halo masses, we predict that UDG is a common, and perhaps even dominant, galaxy type around Ms~1e8 Msun, both in the field and in clusters.
△ Less
Submitted 16 August, 2018; v1 submitted 13 November, 2017;
originally announced November 2017.
-
Signatures of unresolved binaries in stellar spectra: implications for spectral fitting
Authors:
Kareem El-Badry,
Hans-Walter Rix,
Yuan-Sen Ting,
Daniel R. Weisz,
Maria Bergemann,
Phillip Cargile,
Charlie Conroy,
Anna-Christina Eilers
Abstract:
The observable spectrum of an unresolved binary star system is a superposition of two single-star spectra. Even without a detectable velocity offset between the two stellar components, the combined spectrum of a binary system is in general different from that of either component, and fitting it with single-star models may yield inaccurate stellar parameters and abundances. We perform simple experi…
▽ More
The observable spectrum of an unresolved binary star system is a superposition of two single-star spectra. Even without a detectable velocity offset between the two stellar components, the combined spectrum of a binary system is in general different from that of either component, and fitting it with single-star models may yield inaccurate stellar parameters and abundances. We perform simple experiments with synthetic spectra to investigate the effect of unresolved main-sequence binaries on spectral fitting, modeling spectra similar to those collected by the APOGEE, GALAH, and LAMOST surveys. We find that fitting unresolved binaries with single-star models introduces systematic biases in the derived stellar parameters and abundances that are modest but certainly not negligible, with typical systematic errors of $300\,\rm K$ in $T_{\rm eff}$, 0.1 dex in $\log g$, and 0.1 dex in $[\rm Fe/H]$ for APOGEE-like spectra of solar-type stars. These biases are smaller for spectra at optical wavelengths than in the near-infrared. We show that biases can be corrected by fitting spectra with a binary model, which adds only two labels to the fit and includes single-star models as a special case. Our model provides a promising new method to constrain the Galactic binary population, including systems with single-epoch spectra and no detectable velocity offset between the two stars.
△ Less
Submitted 19 October, 2017; v1 submitted 12 September, 2017;
originally announced September 2017.
-
How To Model Supernovae in Simulations of Star and Galaxy Formation
Authors:
Philip F. Hopkins,
Andrew Wetzel,
Dusan Keres,
Claude-Andre Faucher-Giguere,
Eliot Quataert,
Michael Boylan-Kolchin,
Norman Murray,
Christopher C. Hayward,
Kareem El-Badry
Abstract:
We study the implementation of mechanical feedback from supernovae (SNe) and stellar mass loss in galaxy simulations, within the Feedback In Realistic Environments (FIRE) project. We present the FIRE-2 algorithm for coupling mechanical feedback, which can be applied to any hydrodynamics method (e.g. fixed-grid, moving-mesh, and mesh-less methods), and black hole as well as stellar feedback. This a…
▽ More
We study the implementation of mechanical feedback from supernovae (SNe) and stellar mass loss in galaxy simulations, within the Feedback In Realistic Environments (FIRE) project. We present the FIRE-2 algorithm for coupling mechanical feedback, which can be applied to any hydrodynamics method (e.g. fixed-grid, moving-mesh, and mesh-less methods), and black hole as well as stellar feedback. This algorithm ensures manifest conservation of mass, energy, and momentum, and avoids imprinting 'preferred directions' on the ejecta. We show that it is critical to incorporate both momentum and thermal energy of mechanical ejecta in a self-consistent manner, accounting for SNe cooling radii when they are not resolved. Using idealized simulations of single SN explosions, we show that the FIRE-2 algorithm, independent of resolution, reproduces converged solutions in both energy and momentum. In contrast, common 'fully-thermal' (energy-dump) or 'fully-kinetic' (particle-kicking) schemes in the literature depend strongly on resolution: when applied at mass resolution >100 solar masses, they diverge by orders-of-magnitude from the converged solution. In galaxy-formation simulations, this divergence leads to orders-of-magnitude differences in galaxy properties, unless those models are adjusted in a resolution-dependent way. We show that all models that individually time-resolve SNe converge to the FIRE-2 solution at sufficiently high resolution. However, in both idealized single-SN simulations and cosmological galaxy-formation simulations, the FIRE-2 algorithm converges much faster than other sub-grid models without re-tuning parameters.
△ Less
Submitted 11 November, 2018; v1 submitted 21 July, 2017;
originally announced July 2017.
-
Gas kinematics, morphology, and angular momentum in the FIRE simulations
Authors:
Kareem El-Badry,
Eliot Quataert,
Andrew Wetzel,
Philip F. Hopkins,
Daniel R. Weisz,
T. K. Chan,
Alex Fitts,
Michael Boylan-Kolchin,
Dušan Kereš,
Claude-André Faucher-Giguère,
Shea Garrison-Kimmel
Abstract:
We study the z=0 gas kinematics, morphology, and angular momentum content of isolated galaxies in a suite of cosmological zoom-in simulations from the FIRE project spanning $M_{\star}=10^{6-11}M_{\odot}$. Gas becomes increasingly rotationally supported with increasing galaxy mass. In the lowest-mass galaxies ($M_{\star}<10^{8}M_{\odot}$), gas fails to form a morphological disk and is primarily dis…
▽ More
We study the z=0 gas kinematics, morphology, and angular momentum content of isolated galaxies in a suite of cosmological zoom-in simulations from the FIRE project spanning $M_{\star}=10^{6-11}M_{\odot}$. Gas becomes increasingly rotationally supported with increasing galaxy mass. In the lowest-mass galaxies ($M_{\star}<10^{8}M_{\odot}$), gas fails to form a morphological disk and is primarily dispersion and pressure supported. At intermediate masses ($M_{\star}=10^{8-10}M_{\odot}$), galaxies display a wide range of gas kinematics and morphologies, from thin, rotating disks, to irregular spheroids with negligible net rotation. All the high-mass ($M_{\star}=10^{10-11}M_{\odot}$) galaxies form rotationally supported gas disks. Many of the halos whose galaxies fail to form disks harbor high angular momentum gas in their circumgalactic medium. The ratio of the specific angular momentum of gas in the central galaxy to that of the dark-matter halo increases significantly with galaxy mass, from $j_{\rm gas}/j_{\rm DM}\sim0.1$ at $M_{\star}=10^{6-7}M_{\odot}$ to $j_{\rm gas}/j_{\rm DM}\sim2$ at $M_{\star}=10^{10-11}M_{\odot}$. The reduced rotational support in the lowest-mass galaxies owes to (a) stellar feedback and the UV background suppressing the accretion of high-angular momentum gas at late times, and (b) stellar feedback driving large non-circular gas motions. We broadly reproduce the observed scaling relations between galaxy mass, gas rotation velocity, size, and angular momentum, but may somewhat underpredict the incidence of disky, high-angular momentum galaxies at the lowest observed masses ($M_{\star}=(10^{6}-2\times10^{7})M_{\odot}$). In our simulations, stars are uniformly less rotationally supported than gas. The common assumption that stars follow the same rotation curve as gas thus substantially overestimates galaxies' stellar angular momentum, particularly at low masses.
△ Less
Submitted 21 September, 2017; v1 submitted 29 May, 2017;
originally announced May 2017.
-
The statistical challenge of constraining the low-mass IMF in Local Group dwarf galaxies
Authors:
Kareem El-Badry,
Daniel R. Weisz,
Eliot Quataert
Abstract:
We use Monte Carlo simulations to explore the statistical challenges of constraining the characteristic mass ($m_c$) and width ($σ$) of a lognormal sub-solar initial mass function (IMF) in Local Group dwarf galaxies using direct star counts. For a typical Milky Way (MW) satellite ($M_{V} = -8$), jointly constraining $m_c$ and $σ$ to a precision of $\lesssim 20\%$ requires that observations be comp…
▽ More
We use Monte Carlo simulations to explore the statistical challenges of constraining the characteristic mass ($m_c$) and width ($σ$) of a lognormal sub-solar initial mass function (IMF) in Local Group dwarf galaxies using direct star counts. For a typical Milky Way (MW) satellite ($M_{V} = -8$), jointly constraining $m_c$ and $σ$ to a precision of $\lesssim 20\%$ requires that observations be complete to $\lesssim 0.2 M_{\odot}$, if the IMF is similar to the MW IMF. A similar statistical precision can be obtained if observations are only complete down to $0.4M_{\odot}$, but this requires measurement of nearly 100$\times$ more stars, and thus, a significantly more massive satellite ($M_{V} \sim -12$). In the absence of sufficiently deep data to constrain the low-mass turnover, it is common practice to fit a single-sloped power law to the low-mass IMF, or to fit $m_c$ for a lognormal while holding $σ$ fixed. We show that the former approximation leads to best-fit power law slopes that vary with the mass range observed and can largely explain existing claims of low-mass IMF variations in MW satellites, even if satellite galaxies have the same IMF as the MW. In addition, fixing $σ$ during fitting leads to substantially underestimated uncertainties in the recovered value of $m_c$ (by a factor of $\sim 4$ for typical observations). If the IMFs of nearby dwarf galaxies are lognormal and do vary, observations must reach down to $\sim m_c$ in order to robustly detect these variations. The high-sensitivity, near-infrared capabilities of JWST and WFIRST have the potential to dramatically improve constraints on the low-mass IMF. We present an efficient observational strategy for using these facilities to measure the IMFs of Local Group dwarf galaxies.
△ Less
Submitted 16 February, 2017; v1 submitted 9 January, 2017;
originally announced January 2017.
-
When the Jeans don't fit: How stellar feedback drives stellar kinematics and complicates dynamical modeling in low-mass galaxies
Authors:
Kareem El-Badry,
Andrew R. Wetzel,
Marla Geha,
Eliot Quataert,
Philip F. Hopkins,
Dusan Kereš,
T. K. Chan,
Claude-André Faucher-Giguère
Abstract:
In low-mass galaxies, stellar feedback can drive gas outflows that generate non-equilibrium fluctuations in the gravitational potential. Using cosmological zoom-in baryonic simulations from the Feedback in Realistic Environments (FIRE) project, we investigate how these fluctuations affect stellar kinematics and the reliability of Jeans dynamical modeling in low-mass galaxies. We find that stellar…
▽ More
In low-mass galaxies, stellar feedback can drive gas outflows that generate non-equilibrium fluctuations in the gravitational potential. Using cosmological zoom-in baryonic simulations from the Feedback in Realistic Environments (FIRE) project, we investigate how these fluctuations affect stellar kinematics and the reliability of Jeans dynamical modeling in low-mass galaxies. We find that stellar velocity dispersion and anisotropy profiles fluctuate significantly over the course of galaxies' starburst cycles. We therefore predict an observable correlation between star formation rate and stellar kinematics: dwarf galaxies with higher recent star formation rates should have systemically higher stellar velocity dispersions. This prediction provides an observational test of the role of stellar feedback in regulating both stellar and dark-matter densities in dwarf galaxies. We find that Jeans modeling, which treats galaxies as virialized systems in dynamical equilibrium, overestimates a galaxy's dynamical mass during periods of post-starburst gas outflow and underestimates it during periods of net inflow. Short-timescale potential fluctuations lead to typical errors of $\sim 20\%$ in dynamical mass estimates, even if full 3-dimensional stellar kinematics -- including the orbital anisotropy -- are known exactly. When orbital anisotropy is $\textit{not}$ known a priori, typical mass errors arising from non-equilibrium fluctuations in the potential are larger than those arising from the mass-anisotropy degeneracy. However, Jeans modeling alone $\textit{cannot}$ reliably constrain the orbital anisotropy, and problematically, it often favors anisotropy models that do not reflect the true profile. If galaxies completely lose their gas and cease forming stars, fluctuations in the potential subside, and Jeans modeling becomes much more reliable.
△ Less
Submitted 27 December, 2016; v1 submitted 13 October, 2016;
originally announced October 2016.
-
Breathing FIRE: How Stellar Feedback Drives Radial Migration, Rapid Size Fluctuations, and Population Gradients in Low-Mass Galaxies
Authors:
Kareem El-Badry,
Andrew R. Wetzel,
Marla Geha,
Philip F. Hopkins,
Dušan Kereš,
T. K. Chan,
Claude-André Faucher-Giguère
Abstract:
We examine the effects of stellar feedback and bursty star formation on low-mass galaxies ($M_{\rm star}=2\times10^6-5\times10^{10}{\rm M_{\odot}}$) using the FIRE (Feedback in Realistic Environments) simulations. While previous studies emphasized the impact of feedback on dark matter profiles, we investigate the impact on the stellar component: kinematics, radial migration, size evolution, and po…
▽ More
We examine the effects of stellar feedback and bursty star formation on low-mass galaxies ($M_{\rm star}=2\times10^6-5\times10^{10}{\rm M_{\odot}}$) using the FIRE (Feedback in Realistic Environments) simulations. While previous studies emphasized the impact of feedback on dark matter profiles, we investigate the impact on the stellar component: kinematics, radial migration, size evolution, and population gradients. Feedback-driven outflows/inflows drive significant radial stellar migration over both short and long timescales via two processes: (1) outflowing/infalling gas can remain star-forming, producing young stars that migrate $\sim1{\rm\,kpc}$ within their first $100 {\rm\,Myr}$, and (2) gas outflows/inflows drive strong fluctuations in the global potential, transferring energy to all stars. These processes produce several dramatic effects. First, galaxies' effective radii can fluctuate by factors of $>2$ over $\sim200 {\rm\,Myr}$, and these rapid size fluctuations can account for much of the observed scatter in radius at fixed $M_{\rm star}.$ Second, the cumulative effects of many outflow/infall episodes steadily heat stellar orbits, causing old stars to migrate outward most strongly. This age-dependent radial migration mixes---and even inverts---intrinsic age and metallicity gradients. Thus, the galactic-archaeology approach of calculating radial star-formation histories from stellar populations at $z=0$ can be severely biased. These effects are strongest at $M_{\rm star}\approx10^{7-9.6}{\rm M_{\odot}}$, the same regime where feedback most efficiently cores galaxies. Thus, detailed measurements of stellar kinematics in low-mass galaxies can strongly constrain feedback models and test baryonic solutions to small-scale problems in $Λ$CDM.
△ Less
Submitted 30 March, 2016; v1 submitted 3 December, 2015;
originally announced December 2015.
