BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Andrew Wetzel

Andrew Wetzel contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.COastro-ph.SRhep-phastro-ph.IMMachine Learning

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
35works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

35 published item(s)

preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

The JWST Resolved Stellar Populations Early Release Science Program II. Survey Overview

We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

We characterize the 3-D spatial variations of [Fe/H], [Mg/H], and [Mg/Fe] in stars at the time of their formation, across 11 simulated Milky Way (MW)- and M31-mass galaxies in the FIRE-2 simulations, to inform initial conditions for chemical tagging. The overall scatter in [Fe/H] within a galaxy decreased with time until $\approx 7$ Gyr ago, after which it increased to today: this arises from a competition between a reduction of azimuthal scatter and a steepening of the radial gradient in abundance over time. The radial gradient is generally negative, and it steepened over time from an initially flat gradient $\gtrsim 12$ Gyr ago. The strength of the present-day abundance gradient does not correlate with when the disk `settled&#39;; 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Amplified J-factors in the Galactic Center for velocity-dependent darkmatter annihilation in FIRE simulations

We use FIRE-2 zoom cosmological simulations of Milky Way size galaxy halos to calculate astrophysical J-factors for dark matter annihilation and indirect detection studies. In addition to velocity-independent (s-wave) annihilation cross sections $σ_v$, we also calculate effective J-factors for velocity-dependent models, where the annihilation cross section is either either p-wave ($\propto v^2/c^2$) or d-wave ($\propto v^4/c^4$). We use 12 pairs of simulations, each run with dark-matter-only (DMO) physics and FIRE-2 physics. We observe FIRE runs produce central dark matter velocity dispersions that are systematically larger than in DMO runs by factors of $\sim 2.5-4$. They also have a larger range of central ($\sim 400$ pc) dark matter densities than the DMO runs ($ρ_{\rm FIRE}/ρ_{\rm DMO} \simeq 0.5 - 3$) owing to the competing effects of baryonic contraction and feedback. At 3 degrees from the Galactic Center, FIRE J-factors are $5-50$ (p-wave) and $15-500$ (d-wave) times higher than in the DMO runs. The change in s-wave signal at 3 degrees is more modest and can be higher or lower ($\sim 0.3-6$), though the shape of the emission profile is flatter (less peaked towards the Galactic Center) and more circular on the sky in FIRE runs. Our results for s-wave are broadly consistent with the range of assumptions in most indirect detection studies. We observe p-wave J-factors that are significantly enhanced compared to most past estimates. We find that thermal models with p-wave annihilation may be within range of detection in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dissipative Dark Matter on FIRE: II. Observational signatures and constraints from local dwarf galaxies

We analyze the first set of cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modelled as a constant fractional energy dissipation ($f_{\rm diss}=0.5$). In this paper, we examine the properties of dwarf galaxies with $M_{\ast} \sim 10^{5}\operatorname{-}10^{9}\,{\rm M}_{\odot}$ in both isolation and within Milky Way-mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of stellar/neutral gas disk formation in dSIDM with $(σ/m)\leq 1\,{\rm cm^2\,g^{-1}}$ but they are still consistent with observed galaxy sizes and masses. In addition, as a result of the steeper central density profiles developed in dSIDM, the sub-kpc circular velocities of isolated dwarfs in models with $(σ/m)\geq 0.1\,{\rm cm^2\,g^{-1}}$ are enhanced by about a factor of two, which are still consistent with the measured stellar velocity dispersions of Local Group dwarfs but in tension with the HI rotation curves of more massive field dwarfs. Meanwhile, for satellites of the simulated Milky Way-mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help address the missing compact dwarf satellites in CDM. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with constant cross-section $(σ/m) \gtrsim 0.1\,{\rm cm^2\,g^{-1}}$ (assuming $f_{\rm diss}=0.5$) are effectively ruled out in bright dwarfs ($M_{\rm halo}\sim 10^{11}\,{\rm M}_{\odot}$) by circular velocity constraints. However, models with lower effective cross-sections (at this halo mass/velocity scale) are still viable and can give rise to non-trivial observable signatures.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Evidence for a Vast Prograde Stellar Stream in the Solar Vicinity

Massive dwarf galaxies that merge with the Milky Way on prograde orbits can be dragged into the disk plane before being completely disrupted. Such mergers can contribute to an accreted stellar disk and a dark matter disk. We present evidence for Nyx, a vast new stellar stream in the vicinity of the Sun, that may provide the first indication that such an event occurred in the Milky Way. We identify about 500 stars that have coherent radial and prograde motion in this stream using a catalog of accreted stars built by applying deep learning methods to the second Gaia data release. Nyx is concentrated within $\pm 2$ kpc of the Galactic midplane and spans the full radial range studied (6.5-9.5 kpc). The kinematics of Nyx stars are distinct from those of both the thin and thick disk. In particular, its rotational speed lags the disk by $\sim 80$ km/s and its stars follow more eccentric orbits. A small number of Nyx stars have chemical abundances or inferred ages; from these, we deduce that Nyx stars have a peak metallicity of [Fe/H] $\sim -0.5$ and ages $\sim $10-13 Gyr. Taken together with the kinematic observations, these results strongly favor the interpretation that Nyx is the remnant of a disrupted dwarf galaxy. To further justify this interpretation, we explicitly demonstrate that metal-rich, prograde streams like Nyx can be found in the disk plane of Milky Way-like galaxies using the FIRE hydrodynamic simulations. Future spectroscopic studies will be able to validate whether Nyx stars originate from a single progenitor.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Exploring metallicity-dependent rates of Type Ia supernovae and their impact on galaxy formation

Type Ia supernovae play a critical role in stellar feedback and elemental enrichment in galaxies. Recent transient surveys like the All-Sky Automated Survey for Supernova (ASAS-SN) and the Dark Energy Survey (DES) find that the specific Ia rate at z ~ 0 may be ~ 15-50 times higher in lower-mass galaxies than at Milky Way-mass. Independently, Milky Way observations show that the close-binary fraction of solar-type stars is higher at lower metallicity. Motivated by these observations, we use the FIRE-2 cosmological zoom-in simulations to explore the impact of varying Ia rate models, including metallicity dependence, on galaxies across a range of stellar masses: 10^7 Msun - 10^{11} Msun. First, we benchmark our simulated star-formation histories (SFHs) against observations. We show that assumed SFHs and stellar mass functions play a major role in determining the degree of tension between observations and metallicity-independent Ia rate models, and potentially cause ASAS-SN and DES observations to be much more consistent with each other than might naively appear. Models in which the Ia rate increases with decreasing metallicity (as ~ Z^{-0.5} to Z^{-1}) provide significantly better agreement with observations. Encouragingly, these increases in Ia rate (> 10 times in low-mass galaxies) do not significantly impact galaxy stellar masses and morphologies: effective radii, axis ratios, and v/sigma remain largely unaffected except for our most extreme rate models. We explore implications for both [Fe/H] and [alpha/Fe] enrichment: metallicity-dependent Ia rate models can improve agreement with observed stellar mass-metallicity relations in low-mass galaxies. Our results demonstrate that a wide range of metallicity-dependent Ia models are viable for galaxy formation and motivate future work in this area.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark matter and modified gravity. Apprehension grew after several cosmological simulations using the conventional model failed to form adequate numerical analogues with comparable internal characteristics (stellar masses, sizes, velocity dispersions and morphologies). Here we show that the standard paradigm naturally produces galaxies lacking dark matter with internal characteristics in agreement with observations. Using a state-of-the-art cosmological simulation and a meticulous galaxy-identification technique, we find that extreme close encounters with massive neighbours can be responsible for this. We predict that approximately 30 percent of massive central galaxies (with at least 1e11 solar masses in stars) harbour at least one dark-matter-deficient satellite (with 1e8 - 1e9 solar masses in stars). This distinctive class of galaxies provides an additional layer in our understanding of the role of interactions in shaping galactic properties. Future observations surveying galaxies in the aforementioned regime will provide a crucial test of this scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gas infall and radial transport in cosmological simulations of Milky Way-mass disks

Observations indicate that a continuous supply of gas is needed to maintain observed star formation rates in large, disky galaxies. To fuel star formation, gas must reach the inner regions of such galaxies. Despite its crucial importance for galaxy evolution, how and where gas joins galaxies is poorly constrained observationally and is rarely explored in fully cosmological simulations. To investigate gas accretion in the vicinity of galaxies, we analyze the FIRE-2 cosmological zoom-in simulations for 4 Milky Way mass galaxies (M_halo ~ 10E12 solar masses), focusing on simulations with cosmic ray physics. We find that at z~0, gas approaches the disk with angular momentum similar to the gaseous disk edge and low radial velocities, piling-up near the edge and settling into full rotational support. Accreting gas moves predominantly parallel to the disk with small but nonzero vertical velocity components, and joins the disk largely in the outskirts as opposed to &#34;raining&#34; down onto the disk. Once in the disk, gas trajectories are complex, being dominated by spiral arm induced oscillations and feedback. However, time and azimuthal averages show clear but slow net radial infall with transport speeds of 1-3 km/s and net mass fluxes through the disk on the order of one solar mass per year, comparable to the star formation rates of the galaxies and decreasing towards galactic center as gas is sunk into star formation. These rates are slightly higher in simulations without cosmic rays (1-7 km/s, ~4-5 solar masses per year). We find overall consistency of our results with observational constraints and discuss prospects of future observations of gas flows in and around galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Great Balls of FIRE I: The formation of star clusters across cosmic time in a Milky Way-mass galaxy

The properties of young star clusters formed within a galaxy are thought to vary in different interstellar medium (ISM) conditions, but the details of this mapping from galactic to cluster scales are poorly understood due to the large dynamic range involved in galaxy and star cluster formation. We introduce a new method for modeling cluster formation in galaxy simulations: mapping giant molecular clouds (GMCs) formed self-consistently in a FIRE-2 MHD galaxy simulation onto a cluster population according to a GMC-scale cluster formation model calibrated to higher-resolution simulations, obtaining detailed properties of the galaxy&#39;s star clusters in mass, metallicity, space, and time. We find $\sim 10\%$ of all stars formed in the galaxy originate in gravitationally-bound clusters overall, and this fraction increases in regions with elevated $Σ_{\rm gas}$ and $Σ_{\rm SFR}$, because such regions host denser GMCs with higher star formation efficiency. These quantities vary systematically over the history of the galaxy, driving variations in cluster formation. The mass function of bound clusters varies -- no single Schechter-like or power-law distribution applies at all times. In the most extreme episodes, clusters as massive as $7\times 10^6 M_\odot$ form in massive, dense clouds with high star formation efficiency. The initial mass-radius relation of young star clusters is consistent with an environmentally-dependent 3D density that increases with $Σ_{\rm gas}$ and $Σ_{\rm SFR}$. The model does not reproduce the age and metallicity statistics of old ($>11\rm Gyr$) globular clusters found in the Milky Way, possibly because it forms stars more slowly at $z>3$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Hot-mode accretion and the physics of thin-disk galaxy formation

We use FIRE simulations to study disk formation in z~0, Milky Way-mass galaxies, and conclude that a key ingredient for the formation of thin stellar disks is the ability for accreting gas to develop an aligned angular momentum distribution via internal cancellation *prior* to joining the galaxy. Among galaxies with a high fraction (>70%) of their young stars in a thin disk (h/R~0.1) we find that: (i) hot, virial-temperature gas dominates the inflowing gas mass on halo scales (>~20 kpc), with radiative losses offset by compression heating; (ii) this hot accretion proceeds until angular momentum support slows inward motion, at which point the gas cools to T~10^4 K or less; (iii) prior to cooling, the accreting gas develops an angular momentum distribution that is aligned with the galaxy disk, and while cooling transitions from a quasi-spherical spatial configuration to a more flattened, disk-like configuration. We show that the existence of this &#34;rotating cooling flow&#34; accretion mode is strongly correlated with the fraction of stars forming in a thin disk among a sample of 17 z~0 galaxies spanning a halo mass range of 10^10.5 solar masses to 10^12 solar masses, or a stellar mass range 10^8 solar masses to 10^11 solar masses. Notably, galaxies with a thick disk or irregular morphology do not undergo significant angular momentum alignment of gas prior to accretion and show no correspondence between halo gas cooling and flattening. Our results suggest that rotating cooling flows (or, more generally, rotating subsonic flows) that become coherent and angular momentum-supported prior to accretion onto the galaxy are likely a necessary condition for the formation of thin, star-forming disk galaxies in a LambdaCDM universe.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Rapid disc settling and the transition from bursty to steady star formation in Milky Way-mass galaxies

Recent observations and simulations indicate substantial evolution in the properties of galaxies with time, wherein rotationally-supported and steady thin discs (like those frequently observed in the local universe) emerge from galaxies that are clumpy, irregular, and have bursty star formation rates (SFRs). To better understand the progenitors of local disc galaxies we carry out an analysis of three FIRE-2 simulated galaxies with a mass similar to the Milky Way at redshift z=0. We show that all three galaxies transition from bursty to steady SFRs at a redshift between z=0.5 and z=0.8, and that this transition coincides with a rapid (< ~1 Gyr) emergence of a rotationally-supported interstellar medium (ISM).In the late phase with steady SFR, the rotational energy comprises > ~90% of the total kinetic + thermal energy in the ISM, and is roughly half the gravitational energy. By contrast, during the early phase with bursty star formation, the ISM has a quasi-spheroidal morphology and its energy budget is dominated by quasi-isotropic flows including turbulence and coherent inflows/outflows. This result, that rotational support is subdominant at early times, challenges the common application of equilibrium disc models to the high-redshift progenitors of Milky Way-like galaxies. We further find that the formation of a rotation-supported ISM coincides with the formation of a thermal energy-supported inner circumgalactic medium (CGM). Before this transition, the inner CGM is also supported by turbulence and coherent flows, indicating that at early times there is no clear boundary between the ISM and inner CGM.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Shapes of Milky-Way-Mass Galaxies with Self-Interacting Dark Matter

Self-interacting dark matter (SIDM) models offer one way to reconcile inconsistencies between observations and predictions from collisionless cold dark matter (CDM) models on dwarf-galaxy scales. In order to incorporate the effects of both baryonic and SIDM interactions, we study a suite of cosmological-baryonic simulations of Milky-Way (MW)-mass galaxies from the Feedback in Realistic Environments (FIRE-2) project where we vary the SIDM self-interaction cross-section $σ/m$. We compare the shape of the main dark matter (DM) halo at redshift $z=0$ predicted by SIDM simulations (at $σ/m=0.1$, $1$, and $10$ cm$^2$ g$^{-1}$) with CDM simulations using the same initial conditions. In the presence of baryonic feedback effects, we find that SIDM models do not produce the large differences in the inner structure of MW-mass galaxies predicted by SIDM-only models. However, we do find that the radius where the shape of the total mass distribution begins to differ from that of the stellar mass distribution is dependent on $σ/m$. This transition could potentially be used to set limits on the SIDM cross-section in the MW.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spatially Resolved Gas-phase Metallicity in FIRE-2 Dwarfs: Late-Time Evolution of Metallicity Relations in Simulations with Feedback and Mergers

We present an analysis of spatially resolved gas-phase metallicity relations in five dwarf galaxies ($M_{halo} \approx 10^{11} M_\odot$, $M_\star \approx 10^{8.8}-10^{9.6} M_\odot$) from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, which include an explicit model for sub-grid turbulent mixing of metals in gas, near $z\approx 0$, over a period of 1.4 Gyrs, and compare our findings with observations. While these dwarf galaxies represent a diverse sample, we find that all simulated galaxies match the observed mass-metallicity (MZR) and mass-metallicity gradient (MZGR) relations. We note that in all five galaxies, the metallicities are effectively identical between phases of the interstellar medium (ISM), with 95$\%$ being within $\pm$0.1 dex between various ISM phases, including the cold and dense gas ($T < 500$ K and $n_{\rm H} > 1$ cm$^{-3}$), ionized gas (near the H$α$ $T \approx 10^4$ K ridge-line), and nebular regions (ionized gas where the 10 Myr-averaged star formation rate is non-zero). We find that most of the scatter in relative metallicity between cold and dense gas and ionized gas/nebular regions can be attributed to either local starburst events or metal-poor inflows. We also note the presence of a major merger in one of our galaxies, m11e, with a substantial impact on the metallicity distribution in the spatially resolved map, showing two strong metallicity peaks and triggering a starburst in the main galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The In-situ Origins of Dwarf Stellar Outskirts in FIRE-2

Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies ($10^{8.5}<M_\star/M_\odot<10^{9.6}$) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structure, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in-situ, largely via the outward migration of disk stars. However, there also exists a large population of &#34;non-disky&#34; dwarfs in FIRE-2 that lack a well-defined disk/halo and do not resemble the observed dwarf population. These non-disky dwarfs tend to be either more gas poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy&#39;s intrinsic shape -- which is a direct quantification of the distribution of the galaxy&#39;s stellar content -- to interrogate the feedback implementation in simulated galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The relationship between age, metallicity, and abundances for disk stars in a simulated Milky Way galaxy

Observations of the Milky Way&#39;s low-$α$ disk show that at fixed metallicity, [Fe/H], several element abundance, [X/Fe], correlate with age, with unique slopes and small scatters around the age-[X/Fe] relations. In this study, we turn to simulations to explore the age-[X/Fe] relations for the elements C, N, O, Mg, Si, S, and Ca that are traced in a FIRE-2 cosmological zoom-in simulation of a Milky Way-like galaxy, m12i, and understand what physical conditions give rise to the observed age-[X/Fe] trends. We first explore the distributions of mono-age populations in their birth and current locations, [Fe/H], and [X/Fe], and find evidence for inside-out radial growth for stars with ages < 7 Gyr. We then examine the age-[X/Fe] relations across m12i&#39;s disk and find that the direction of the trends agree with observations, apart from C, O, and Ca, with remarkably small intrinsic scatters, $σ_{int}$ (0.01-0.04 dex). This $σ_{int}$ measured in the simulations is also metallicity-dependent, with $σ_{int}$ $\approx$ 0.025 dex at [Fe/H]=-0.25 dex versus $σ_{int}$ $\approx$ 0.015 dex at [Fe/H]=0 dex, and a similar metallicity dependence is seen in the GALAH survey for the elements in common. Additionally, we find that $σ_{int}$ is higher in the inner galaxy, where stars are older and formed in less chemically-homogeneous environments. The age-[X/Fe] relations and the small scatter around them indicate that simulations capture similar chemical enrichment variance as observed in the Milky Way, arising from stars sharing similar element abundances at a given birth place and time.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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&#39;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&#39;s low star formation efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Galactic Mass Estimates using Dwarf Galaxies as Kinematic Tracers

New mass estimates and cumulative mass profiles with Bayesian credible regions (c.r.) for the Milky Way (MW) are found using the Galactic Mass Estimator (GME) code and dwarf galaxy (DG) kinematic data from multiple sources. GME takes a hierarchical Bayesian approach to simultaneously estimate the true positions and velocities of the DGs, their velocity anisotropy, and the model parameters for the Galaxy&#39;s total gravitational potential. In this study, we incorporate meaningful prior information from past studies and simulations. The prior distributions for the physical model are informed by the results of Eadie & Juric (2019), which used globular clusters instead of DGs, as well as by the subhalo distributions of the Ananke Gaia-like surveys from Feedback In Realistic Environments-2 (Fire-2) cosmological simulations (see Sanderson et al. 2020). Using DGs beyond 45 kpc, we report median and 95% c.r estimates for $r_{200}$ = 212.8 (191.12,238.44) kpc, and for the total enclosed mass $M_{200}$ = 1.19 (0.87,1.68)$\times10^{12}M_{\odot}$ (adopting $Δ_c=200$). Median mass estimates at specific radii are also reported (e.g., $M(<50\text{ kpc})=0.52\times10^{12}M_{\odot}$ and $M(100\text{ kpc})=0.78\times10^{12}M_{\odot}$). Estimates are comparable to other recent studies using GAIA DR2 and DGs, but notably different from the estimates of Eadie & Juric (2019). We perform a sensitivity analysis to investigate whether individual DGs and/or a more massive Large Magellanic Cloud (LMC) on the order of $10^{11}M_{\odot}$ may be affecting our mass estimates. We find possible supporting evidence for the idea that some DGs are affected by a massive LMC and are not in equilibrium with the MW.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Reading the CARDs: the Imprint of Accretion History in the Chemical Abundances of the Milky Way&#39;s Stellar Halo

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&#39; 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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The effects of LMC-mass environments on their dwarf satellite galaxies in the FIRE simulations

Characterizing the predicted environments of dwarf galaxies like the Large Magellanic Cloud (LMC) is becoming increasingly important as next generation surveys push sensitivity limits into this low-mass regime at cosmological distances. We study the environmental effects of LMC-mass halos ($M_{200m} \sim 10^{11}$ M$_\odot$) on their populations of satellites ($M_\star \geq 10^4$ M$_\odot$) using a suite of zoom-in simulations from the Feedback In Realistic Environments (FIRE) project. Our simulations predict significant hot coronas with $T\sim10^5$ K and $M_\text{gas}\sim10^{9.5}$ M$_\odot$. We identify signatures of environmental quenching in dwarf satellite galaxies, particularly for satellites with intermediate mass ($M_\star = 10^{6-7}$ M$_\odot$). The gas content of such objects indicates ram-pressure as the likely quenching mechanism, sometimes aided by star formation feedback. Satellites of LMC-mass hosts replicate the stellar mass dependence of the quiescent fraction found in satellites of MW mass hosts (i.e. that the quiescent fraction increases as stellar mass decreases). Satellites of LMC-mass hosts have a wider variety of quenching times when compared to the strongly bi-modal distribution of quenching times of nearby centrals. Finally, we identify significant tidal stellar structures around four of our six LMC-analogs, suggesting that stellar streams may be common. These tidal features originated from satellites on close orbits, extend to $\sim$80 kpc from the central galaxy, and contain $\sim10^{6-7}$ M$_\odot$~of stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A dark matter profile to model diverse feedback-induced core sizes of $Λ$CDM haloes

We analyze the cold dark matter density profiles of 54 galaxy halos simulated with FIRE-2 galaxy formation physics, each resolved within $0.5\%$ of the halo virial radius. These halos contain galaxies with masses that range from ultra-faint dwarfs ($M_\star \simeq 10^{4.5} M_{\odot}$) to the largest spirals ($M_\star \simeq 10^{11} M_{\odot}$) and have density profiles that are both cored and cuspy. We characterize our results using a new analytic density profile that extends the standard Einasto form to allow for a pronounced constant-density core in the resolved innermost radius. With one additional core-radius parameter, $r_{c}$, this &#34;core-Einasto&#34; profile is able to characterize the shape and normalization of our feedback-impacted dark matter halos. In order to enable comparisons with observations, we provide fitting functions for $r_{c}$ and other profile parameters as a function of both $M_\star$ and $M_{\star}/M_{\rm halo}$. In agreement with similar studies done in the literature, we find that dark matter core formation is most efficient at the characteristic stellar-mass to halo-mass ratio $M_\star/M_{\rm halo} \simeq 5 \times 10^{-3}$, or $M_{\star} \sim 10^9 \, M_{\odot}$, with cores that are roughly the size of the galaxy half-light radius, $r_{c} \simeq 1-5$ kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with $M_{\star}/M_{\rm halo} < 5\times 10^{-4}$ or $M_\star \lesssim 10^6 \, M_{\odot}$. For Milky Way-size galaxies, baryonic contraction often makes halos significantly more concentrated and dense at the stellar half-light radius than dark matter only runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of $\simeq 0.5-2$ kpc in size. Recent evidence for a ${\sim} 2$ kpc core in the Milky Way&#39;s dark matter halo is consistent with this expectation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Relationship Between Stellar Metallicity Gradients and Galaxy Age in Dwarf Galaxies

We explore the origin of stellar metallicity gradients in simulated and observed dwarf galaxies. We use FIRE-2 cosmological baryonic zoom-in simulations of 26 isolated galaxies as well as existing observational data for 10 Local Group dwarf galaxies. Our simulated galaxies have stellar masses between $10^{5.5}$ and $10^{8.6} \msun$. Whilst gas-phase metallicty gradients are generally weak in our simulated galaxies, we find that stellar metallicity gradients are common, with central regions tending to be more metal-rich than the outer parts. The strength of the gradient is correlated with galaxy-wide median stellar age, such that galaxies with younger stellar populations have flatter gradients. Stellar metallicty gradients are set by two competing processes: (1) the steady &#34;puffing&#34; of old, metal-poor stars by feedback-driven potential fluctuations, and (2) the accretion of extended, metal-rich gas at late times, which fuels late-time metal-rich star formation. If recent star formation dominates, then extended, metal-rich star formation washes out pre-existing gradients from the &#34;puffing&#34; process. We use published results from ten Local Group dwarf galaxies to show that a similar relationship between age and stellar metallicity-gradient strength exists among real dwarfs. This suggests that observed stellar metallicity gradients may be driven largely by the baryon/feedback cycle rather than by external environmental effects.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cataloging Accreted Stars within Gaia DR2 using Deep Learning

The goal of this study is to present the development of a machine learning based approach that utilizes phase space alone to separate the Gaia DR2 stars into two categories: those accreted onto the Milky Way from those that are in situ. Traditional selection methods that have been used to identify accreted stars typically rely on full 3D velocity, metallicity information, or both, which significantly reduces the number of classifiable stars. The approach advocated here is applicable to a much larger portion of Gaia DR2. A method known as &#34;transfer learning&#34; is shown to be effective through extensive testing on a set of mock Gaia catalogs that are based on the FIRE cosmological zoom-in hydrodynamic simulations of Milky Way-mass galaxies. The machine is first trained on simulated data using only 5D kinematics as inputs and is then further trained on a cross-matched Gaia/RAVE data set, which improves sensitivity to properties of the real Milky Way. The result is a catalog that identifies around 767,000 accreted stars within Gaia DR2. This catalog can yield empirical insights into the merger history of the Milky Way and could be used to infer properties of the dark matter distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamical heating of the X-ray emitting intracluster medium: the roles of merger shocks and turbulence dissipation

The diffuse plasma inside clusters of galaxies has X-ray emitting temperatures of a few keV. The physical mechanisms that heat this intracluster medium (ICM) to such temperatures include the accretion shock at the periphery of a galaxy cluster, the shocks driven by merger events, as well as a somewhat overlooked mechanism -- the dissipation of intracluster turbulent motions. We study the relative role of these heating mechanisms using galaxy clusters in Lagrangian tracer particle re-simulations of the Omega500 cosmological simulation. We adopt a novel analysis method of decomposing the temperature increase at each time step into the contribution from dissipative heating and that from adiabatic heating. In the high-resolution spatial-temporal map of these heating rates, merger tracks are clearly visible, demonstrating the dominant role of merger events in heating the ICM. The dissipative heating contributed by each merger event is extended in time and also occurs in the rarefaction regions, suggesting the importance of heating by the dissipation of merger-induced turbulence. Quantitative analysis shows that turbulence heating, rather than direct heating at merger shocks, dominates the temperature increase of the ICM, especially at inner radii $r < r_{\rm 500c}$. In addition, we find that many merger shocks can propagate with almost constant velocity to very large radii $r \gg r_{\rm 500c}$, some even reach and join with the accretion shock and becoming the outer boundary of the ICM. Altogether, these results suggest that the ICM is heated more in an `inside-out&#39; fashion rather than `outside-in&#39; as depicted in the classical smooth accretion picture.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effects of Stellar Feedback on Stellar and Gas Kinematics of Star-Forming Galaxies at 0.6<z<1.0

Recent zoom-in cosmological simulations have shown that stellar feedback can flatten the inner density profile of the dark matter halo in low-mass galaxies. A correlation between the stellar/gas velocity dispersion ($σ_{star}$, $σ_{gas}$) and the specific star formation rate (sSFR) is predicted as an observational test of the role of stellar feedback in re-shaping the dark matter density profile. In this work we test the validity of this prediction by studying a sample of star-forming galaxies at $0.6<z<1.0$ from the LEGA-C survey, which provides high signal-to-noise measurements of stellar and gas kinematics. We find that a weak but significant correlation between $σ_{star}$ (and $σ_{gas}$) and sSFR indeed exists for galaxies in the lowest mass bin (M$_\ast\sim10^{10}\,$M$_\odot$). This correlation, albeit with a $\sim$35% scatter, holds for different tracers of star formation, and becomes stronger with redshift. This result generally agrees with the picture that at higher redshifts star formation rate was generally higher, and galaxies at M$_\ast\lesssim10^{10}\,$M$_\odot$ have not yet settled into a disk. As a consequence, they have shallower gravitational potentials more easily perturbed by stellar feedback. The observed correlation between $σ_{star}$ (and $σ_{gas}$) and sSFR supports the scenario predicted by cosmological simulations, in which feedback-driven outflows cause fluctuations in the gravitation potential that flatten the density profiles of low-mass galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Evolution of giant molecular clouds across cosmic time

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

No missing photons for reionization: moderate ionizing photon escape fractions from the FIRE-2 simulations

We present the escape fraction of hydrogen ionizing photons (f_esc) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z>5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 halos in M_vir~10^8--10^{12} M_sun (M_star~10^4--10^{10} M_sun) at z=5--12. We find the sample average <f_esc> increases with halo mass for M_vir~10^8--10^{9.5} M_sun, becomes nearly constant for M_vir~10^{9.5}--10^{11} M_sun, and decreases at M_vir>10^{11} M_sun. Equivalently, <f_esc> increases with stellar mass up to M_star~10^8 M_sun and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate <f_esc>~0.2 for galaxies at M_star~10^8 M_sun. Nearly half of the escaped ionizing photons come from stars 1--3 Myr old and the rest from stars 3--10 Myr old. Binaries only have a modest effect, boosting <f_esc> by ~25--35% and the number of escaped photons by 60--80%. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of <f_esc> at the high-mass end is due to dust attenuation, while at the low-mass end, <f_esc> decreases owing to inefficient star formation (and hence feedback). At fixed mass, <f_esc> tends to increase with redshift. Our simulations produce sufficient ionizing photons for cosmic reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Self-consistent proto-globular cluster formation in cosmological simulations of high-redshift galaxies

We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z>5 galaxies from the FIRE project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over 10^4 Msun pc^-2, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at birth. These high-pressure clouds are compressed by feedback-driven winds and/or collisions of smaller clouds/gas streams in highly gas-rich, turbulent environments. The newly formed clusters follow a power-law mass function of dN/dM~M^-2. The cluster formation efficiency is similar across galaxies with stellar masses of ~10^7-10^10 Msun at z>5. The age spread of cluster stars is typically a few Myrs and increases with cluster mass. The metallicity dispersion of cluster members is ~0.08 dex in [Z/H] and does not depend on cluster mass significantly. Our findings support the scenario that present-day old globular clusters (GCs) were formed during relatively normal star formation in high-redshift galaxies. Simulations with a stricter/looser star formation model form a factor of a few more/fewer bound clusters per stellar mass formed, while the shape of the mass function is unchanged. Simulations with a lower local star formation efficiency form more stars in bound clusters. The simulated clusters are larger than observed GCs due to finite resolution. Our simulations are among the first cosmological simulations that form bound clusters self-consistently in a wide range of high-redshift galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Variations in the slope of the resolved star-forming main sequence: a tool for constraining the mass of star-forming regions

The correlation between galaxies&#39; integrated stellar masses and star formation rates (the `star formation main sequence&#39;; SFMS) is a well-established scaling relation. Recently, surveys have found a relationship between the star formation rate and stellar mass surface densities on kpc and sub-kpc scales (the `resolved SFMS&#39;; rSFMS). In this work, we demonstrate that the rSFMS emerges naturally in FIRE-2 zoom-in simulations of Milky Way-mass galaxies. We make SFR and stellar mass maps of the simulated galaxies at a variety of spatial resolutions and star formation averaging time-scales and fit the rSFMS using multiple methods from the literature. While the absolute value of the SFMS slope depends on the fitting method, the slope is steeper for longer star formation time-scales and lower spatial resolutions regardless of the fitting method employed. We present a toy model that quantitatively captures the dependence of the simulated galaxies&#39; rSFMS slope on spatial resolution and use it to illustrate how this dependence can be used to constrain the characteristic mass of star-forming clumps.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Stars made in outflows may populate the stellar halo of the Milky Way

We study stellar-halo formation using six Milky Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that $5-40\%$ of the outer ($50-300$ kpc) stellar halo in each system consists of $\textit{in-situ}$ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by super-bubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at $z=0$.In the outer halo, outflow stars are more spatially homogeneous, metal rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to $\sim 10 \%$ of our kinematically-identified halo stars were born in outflows; the fraction rises to as high as $\sim 40\%$ for the most metal-rich local halo stars ([Fe/H] $> -0.5$). We conclude that the Milky Way stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell $\&$ Sandange, who suggested that halo stars formed in rapidly $\textit{infalling}$ gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme velocity stars in the Milky Way.

0 citations0 reviewsNo code linkedOpen access