No Arabic abstract
We use FIRE-2 simulations to examine 3-D variations of gas-phase elemental abundances of [O/H], [Fe/H], and [N/H] in 11 Milky Way (MW) and M31-mass galaxies across their formation histories at $z leq 1.5$ ($t_{rm lookback} leq 9.4$ Gyr), motivated by characterizing the initial conditions of stars for chemical tagging. Gas within $1$ kpc of the disk midplane is vertically homogeneous to $lesssim 0.008$ dex at all $z leq 1.5$. We find negative radial gradients (metallicity decreases with galactocentric radius) at all times, which steepen over time from $approx -0.01$ dex kpc$^{-1}$ at $z = 1$ ($t_{rm lookback} = 7.8$ Gyr) to $approx -0.03$ dex kpc$^{-1}$ at $z = 0$, and which broadly agree with observations of the MW, M31, and nearby MW/M31-mass galaxies. Azimuthal variations at fixed radius are typically $0.14$ dex at $z = 1$, reducing to $0.05$ dex at $z = 0$. Thus, over time radial gradients become steeper while azimuthal variations become weaker (more homogeneous). As a result, azimuthal variations were larger than radial variations at $z gtrsim 0.8$ ($t_{rm lookback} gtrsim 6.9$ Gyr). Furthermore, elemental abundances are measurably homogeneous (to $lesssim 0.05$ dex) across a radial range of $Delta R approx 3.5$ kpc at $z gtrsim 1$ and $Delta R approx 1.7$ kpc at $z = 0$. We also measure full distributions of elemental abundances, finding typically negatively skewed normal distributions at $z gtrsim 1$ that evolve to typically Gaussian distributions by $z = 0$. Our results on gas abundances inform the initial conditions for stars, including the spatial and temporal scales for applying chemical tagging to understand stellar birth in the MW.
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.
Stellar streams record the accretion history of their host galaxy. We present a set of simulated streams from disrupted dwarf galaxies in 13 cosmological simulations of Milky Way (MW)-mass galaxies from the FIRE-2 suite at $z=0$, including 7 isolated Milky Way-mass systems and 6 hosts resembling the MW-M31 pair (full dataset at: https://flathub.flatironinstitute.org/sapfire). In total, we identify 106 simulated stellar streams, with no significant differences in the number of streams and masses of their progenitors between the isolated and paired environments. We resolve simulated streams with stellar masses ranging from $sim 5times10^5$ up to $sim 10^{9} M_odot$, similar to the mass range between the Orphan and Sagittarius streams in the MW. We confirm that present-day simulated satellite galaxies are good proxies for stellar stream progenitors, with similar properties including their stellar mass function, velocity dispersion, [Fe/H] and [$alpha$/H] evolution tracks, and orbital distribution with respect to the galactic disk plane. Each progenitors lifetime is marked by several important timescales: its infall, star-formation quenching, and stream-formation times. We show that the ordering of these timescales is different between progenitors with stellar masses higher and lower than $sim 2times10^6 M_odot$. Finally, we show that the main factor controlling the rate of phase-mixing, and therefore fading, of tidal streams from satellite galaxies in MW-mass hosts is non-adiabatic evolution of the host potential. Other factors commonly used to predict phase-mixing timescales, such as progenitor mass and orbital circularity, show virtually no correlation with the number of dynamical times required for a stream to become phase-mixed.
We present the first measurement of the lifetimes of Giant Molecular Clouds (GMCs) in cosmological simulations at $z = 0$, using the Latte suite of FIRE-2 simulations of Milky Way-mass galaxies. We track GMCs with total gas mass $gtrsim 10^5$ M$_odot$ at high spatial ($sim1$ pc), mass ($7100$ M$_{odot}$), and temporal (1 Myr) resolution. Our simulated GMCs are consistent with the distribution of masses for massive GMCs in the Milky Way and nearby galaxies. We find GMC lifetimes of $5-7$ Myr, or 1-2 freefall times, on average, with less than 2$%$ of clouds living longer than 20 Myr. We find decreasing GMC lifetimes with increasing virial parameter, and weakly increasing GMC lifetimes with galactocentric radius, implying that environment affects the evolutionary cycle of GMCs. However, our GMC lifetimes show no systematic dependence on GMC mass or amount of star formation. These results are broadly consistent with inferences from the literature and provide an initial investigation into ultimately understanding the physical processes that govern GMC lifetimes in a cosmological setting.
We study star formation histories (SFHs) of $simeq500$ dwarf galaxies (stellar mass $M_ast = 10^5 - 10^9,M_odot$) from FIRE-2 cosmological zoom-in simulations. We compare dwarfs around individual Milky Way (MW)-mass galaxies, dwarfs in Local Group (LG)-like environments, and true field (i.e. isolated) dwarf galaxies. We reproduce observed trends wherein higher-mass dwarfs quench later (if at all), regardless of environment. We also identify differences between the environments, both in terms of satellite vs. central and LG vs. individual MWvs. isolated dwarf central. Around the individual MW-mass hosts, we recover the result expected from environmental quenching: central galaxies in the near field have more extended SFHs than their satellite counterparts, with the former more closely resemble isolated (true field) dwarfs (though near-field centrals are still somewhat earlier forming). However, this difference is muted in the LG-like environments, where both near-field centrals and satellites have similar SFHs, which resemble satellites of single MW-mass hosts. This distinction is strongest for $M_ast = 10^6 - 10^7,M_odot$ but exists at other masses. Our results suggest that the paired halo nature of the LG may regulate star formation in dwarf galaxies even beyond the virial radii of the MW and Andromeda. Caution is needed when comparing zoom-in simulations targeting isolated dwarf galaxies against observed dwarf galaxies in the LG.
We examine the prevalence, longevity, and causes of planes of satellite dwarf galaxies, as observed in the Local Group. We use 14 Milky Way/Andromeda-(MW/M31) mass host galaxies from the FIRE-2 simulations. We select the 14 most massive satellites by stellar mass within 300 kpc of each host and correct for incompleteness from the foreground galactic disc when comparing to the MW. We find that MW-like planes as spatially thin and/or kinematically coherent as observed are uncommon, but they do exist in our simulations. Spatially thin planes occur in 1-2 per cent of snapshots during $z=0-0.2$, and kinematically coherent planes occur in 5 per cent of snapshots. These planes are generally transient, surviving for less than 500 Myr. However, if we select hosts with an LMC-like satellite near first pericentre, the fraction of snapshots with MW-like planes increases dramatically to 7-16 per cent, with lifetimes of 0.7-1 Gyr, likely because of group accretion of satellites. We find that M31s satellite distribution is much more common: M31s satellites lie within about 1 sigma of the simulation median for every plane metric we consider. We find no significant difference in average satellite planarity for isolated hosts versus hosts in LG-like pairs. Baryonic and dark matter-only simulations exhibit similar levels of planarity, even though baryonic subhaloes are less centrally concentrated within their host haloes. We conclude that planes of satellites are not a strong challenge to LCDM cosmology.