Cosmological Zoom-in Simulations Research Articles

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

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

Hypernova signatures of the first stars in dwarf galaxies in the local group

ABSTRACT Observing the first generation of stars, Population III (Pop III), is still a challenge even with the JWST due to their faintness. Instead, searching for fossil records of Pop III stars in nearby dwarf galaxies provides an alternative method for studying their physical properties. It is intriguing that a star recently discovered in the Sculptor dwarf galaxy, named AS0039, is considered to show the unique signature of a Pop III star. The detailed abundance patterns of AS0039 are well matched with those predicted by nucleosynthesis models for Pop III exploding as an energetic hypernova (HN), confirming its potential to provide insight into the properties of the first stars. This study aims to explore the environmental conditions required for the formation of such a unique star using cosmological hydrodynamic zoom-in simulations on dwarf galaxies with a mass of $M_{\rm vir}\approx 10^8{\ \mathrm{ M}_\odot }$ at z = 0 while varying the fraction of Pop III stars that undergo HNe. Our simulations identify rapid gas inflow ($\dot{M}_{\rm gas}\sim 0.08{\ \mathrm{ M}_\odot }$ $\rm yr^{-1}$) as a possible factor in facilitating the formation of stars similar to AS0039. Alternatively, the delayed formation of subsequent Population II (Pop II) stars in the gas-enriched environment may lead to low-metallicity stars like AS0039. Additionally, using the a-sloth code, we investigate the probability of finding remnants of Pop II stars with HN signatures in nearby dwarf satellite galaxies. We suggest that the most likely dwarf galaxies to contain HN signatures are massive satellites with a probability of 40 per cent in the range of $M_{\rm peak}\approx 10^{10}\ \mathrm{ to} \ 10^{11}{\ \mathrm{ M}_\odot }$ and $M_{\ast }\approx 10^7 \ \mathrm{ to} \ 10^8{\ \mathrm{ M}_\odot }$, considering observational limitations.

Direct Collapse to Precursors of Supermassive Black Hole Seeds: Radiation-feedback-generated Outflows

We use high-resolution zoom-in cosmological simulations to model outflow triggered by radiation and thermal drivers around the central mass accumulation during direct collapse within the dark matter (DM) halo. The maximal resolution is 1.3 × 10−5 pc, and no restrictions are put on the geometry of the inflow/outflow. The central mass is considered prior to the formation of the supermassive black hole seed at a redshift of z ∼ 15.9 and can constitute either a supermassive star (SMS) of ∼105 M ⊙ surrounded by a growing accretion disk or a self-gravitating disk. The radiation transfer is modeled using the ray-tracing algorithm. Due to the high accretion rate of ∼1 M ⊙ yr−1 determined by the DM halo, accretion is mildly supercritical, resulting in mildly supercritical luminosity that has only a limited effect on the accretion rate, with a duty cycle of ∼0.9. We observe a fast development of hot cavities, which quickly extend into polar funnels and expand dense shells. Within the funnels, fast winds, ∼103 km s−1, are mass-loaded by the accreting gas. We follow the expanding shells to ∼1 pc, when the shell velocity remains substantially (∼5 times) above the escape speed. The ionization cones formed by the central UV/X-ray completely ionize the cavities. Extrapolating the outflow properties shows that the halo material outside the shell will have difficulty stopping it. We therefore conclude that the expanding wind-driven shell will break out of the central parsec and will reach the halo virial radius. Finally, the anisotropic accretion flow on subparsec scales will attenuate the UV/soft X-rays on the H2. Hence, the formation of funnels and powerful outflows around, e.g., SMSs can have interesting observational corollaries.

The inefficiency of stellar feedback in driving galactic outflows in massive galaxies at high redshift

ABSTRACT Recent observations indicate that galactic outflows are ubiquitous in high-redshift (high-z) galaxies, including normal star-forming galaxies, quasar hosts, and dusty star-forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star-formation histories and galactic outflow properties of galaxies in massive haloes ($10^{12}\, {\rm M}_{\odot }\ \lt\ M_{\rm vir}\ \lt\ 5\times 10^{12}\, {\rm M}_{\odot }$) at z ≳ 5.5 in three zoom-in cosmological simulations from the MassiveFIRE suite, as part of the Feedback In Realistic Environments (FIRE) project. The simulations were run with the FIRE-2 model, which does not include feedback from active galactic nuclei. The simulated galaxies resemble z > 4 DSFGs, with star-formation rates of $\sim\!{1000}\ {\rm M}_{\odot }\, \rm yr^{-1}$ and molecular gas masses of Mmol ∼ 1010 M⊙. However, the simulated galaxies are characterized by higher circular velocities than those observed in high-z DSFGs. The mass loading factors from stellar feedback are of the order of ∼0.1, implying that stellar feedback is inefficient in driving galactic outflows and gas is consumed by star formation on much shorter time-scales than it is expelled from the interstellar medium. We also find that stellar feedback is highly inefficient in self-regulating star formation in this regime, with an average integrated star formation efficiency (SFE) per dynamical time of 30 per cent. Finally, compared with FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high-z sample. We argue that both effects originate from the higher total and gas surface densities that characterize high-z massive systems.

The stellar halo in Local Group Hestia simulations

Stellar chemical abundances and kinematics provide key information for recovering the assembly history of galaxies. In this work we explore the chemo-chrono-kinematics of accreted and in situ stellar populations, by analyzing six M31/Milky Way (MW) analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We show that elemental abundances ([Fe/H], [Mg/Fe]) of merger debris in the stellar haloes are chemically distinct from the survived dwarf galaxies, in that they are [α/Fe]-enhanced and have lower metallicity in the same stellar mass range. Therefore, mergers debris have abundances expected for stars originating from dwarfs that had their star formation activity quenched at early times. Accreted stellar haloes, including individual debris, reveal [Fe/H] and [Mg/Fe] gradients in the E − Lz plane, with the most metal-rich, [α/Fe]-poor stars, which have formed in the inner parts of the disrupted systems before the merger, contributing mainly to the central regions of the host galaxies. This results in negative metallicity gradients in the accreted components of stellar haloes at z = 0, seen also for the individual merger debris. We suggest, therefore, that abundance measurements of halo stars in the inner MW will allow constraining better the parameters, such as the mass and merger time, of MW’s most massive merger Gaia-Sausage-Enceladus. The metallicity distribution functions (MDFs) of the individual debris show several peaks and the majority of debris have lower metallicity than the in situ stars in the prograde part of the E − Lz space. At the same time, non-rotating and retrograde accreted populations are very similar to the in situ stars in terms of [Fe/H] abundance. Prograde accreted stars show a prominent knee in the [Fe/H]–[Mg/Fe] plane, reaching up to solar [Mg/Fe], while retrograde stars typically contribute to the high-[Mg/Fe] sequence only. We find that the most metal-poor stars ([Fe/H] ≲ −1) of the HESTIA galaxies exhibit net rotation up to 80 km s−1, which is consistent with the Aurora population recently identified in the MW. At higher metallicities ([Fe/H] ≈ −0.5 ± 0.1) we detect a sharp transition (spin-up) from the turbulent phase to a regular disk-like rotation. Different merger debris appear similar in the [Fe/H]–[Mg/Fe] plane, thus making it difficult to identify individual events. However, combining a set of abundances, and especially stellar age, makes it possible to distinguish between different debris.

The stellar halo in Local Group Hestia simulations

Theory suggests that mergers play an important role in shaping galactic discs and stellar haloes, which was observationally confirmed in the Milky Way (MW) thanks to Gaia data. In this work, aiming to probe the contribution of mergers to the in situ stellar halo formation, we analyse six M 31 and MW analogues from the HESTIA suite of cosmological hydrodynamical zoom-in simulations of the Local Group. We found that all the HESTIA galaxies experience between one to four mergers with stellar mass ratios between 0.2 and 1 relative to the host at the time of the merger. These significant mergers, with a single exception, happened 7 − 11 Gyr ago. The overall impact of the most massive mergers in HESTIA is clearly seen as a sharp increase in the orbital eccentricity (and a corresponding decrease in the rotational velocity Vϕ) of pre-existing disc stars of the main progenitor, thus nicely reproducing the Splash-, Plume-like feature that was discovered in the MW. We do find a correlation between mergers and close pericentric passages of massive satellites and bursts of the star formation in the in situ component. Massive mergers sharply increase the disc velocity dispersion of the in situ stars; however, the latest significant merger often heats up the disc up to the numbers when the contribution of the previous ones is less prominent in the age-velocity dispersion relation. In HESTIA galaxies, the in situ halo is an important component of the inner stellar halo where its fraction is about 30 − 40%, while in the outer parts it typically does not exceed ≈5% beyond 15 kpc from the galactic centre. The simulations suggest that this component of the stellar haloes continues to grow well after mergers conclude; however, the most significant contribution comes from stars that formed recently before the merger. The orbital analysis of the HESTIA galaxies suggests that wedges in Rmax − Zmax (apocentre – maximum height from the mid-plane) space are mainly populated by the stars born in between significant mergers.

The stellar halo in Local Group Hestia simulations

Recent progress in understanding the assembly history of the Milky Way (MW) is driven by the tremendous amount of high-quality data delivered byGaia(ESA), revealing a number of substructures potentially linked to several ancient accretion events. In this work we aim to explore the phase-space structure of accreted stars by analysing six M31/MW analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We find that all HESTIA galaxies experience a few dozen mergers but only between one and four of those have stellar mass ratios > 0.2, relative to the host at the time of the merger. Depending on the halo definition, the most massive merger contributes from 20% to 70% of the total stellar halo mass. Individual merger remnants show diverse density distributions atz = 0, significantly overlapping with each other and with the in situ stars in theLz − E, (VR, Vϕ) and (R, vϕ) coordinates. Moreover, merger debris often shifts position in theLz − Espace with cosmic time due to the galactic mass growth and the non-axisymmetry of the potential. In agreement with previous works, we show that even individual merger debris exhibit a number of distinctLz − Efeatures. In the (VR, Vϕ) plane, all HESTIA galaxies reveal radially hot, non-rotating or weakly counter-rotating, Gaia-Sausage-like features, which are the remnants of the most recent significant mergers. We find an age gradient inLz − Espace for individual debris, where the youngest stars, formed in the inner regions of accreting systems, deposit to the innermost regions of the host galaxies. The bulk of these stars formed during the last stages of accretion, making it possible to use the stellar ages of the remnants to date the merger event. In action space (Jr, Jz, Jϕ), merger debris do not appear as isolated substructures, but are instead scattered over a large parameter area and overlap with the in situ stars. We suggest that accreted stars can be best identified usingJr > 0.2−0.3(104 kpc km s−1)0.5. We also introduce a new, purely kinematic space (Jz/Jr-orbital eccentricity), where different merger debris can be disentangled better from each other and from the in situ stars. Accreted stars have a broad distribution of eccentricities, peaking atϵ ≈ 0.6 − 0.9, and their mean eccentricity tends to be smaller for systems accreted more recently.

Fast Radio Bursts as Probes of Magnetic Fields in Galaxies at z < 0.5

We present a sample of nine fast radio bursts (FRBs) from which we derive magnetic field strengths of the host galaxies represented by normal, z < 0.5 star-forming galaxies with stellar masses M * ≈ 108–1010.5 M ⊙. We find no correlation between the FRB rotation measure (RM) and redshift, which indicates that the RM values are due mostly to the FRB host contribution. This assertion is further supported by a significant positive correlation (Spearman test probability P S < 0.05) found between the RM and the estimated host dispersion measure (DMhost; with Spearman rank correlation coefficient r S = +0.75). For these nine galaxies, we estimate their magnetic field strengths projected along the sight line ∣B ∥∣, finding a low median value of 0.5 μG. This implies the magnetic fields of our sample of hosts are weaker than those characteristic of the solar neighborhood (≈6 μG), but relatively consistent with a lower limit on the observed range of ≈2–10 μG for star-forming disk galaxies, especially as we consider reversals in the B-field, and that we are only probing B∥. We compare to RMs from simulated galaxies of the Auriga project—magneto-hydrodynamic cosmological zoom simulations—and find that the simulations predict the observed values to within a 95% confidence interval. Upcoming FRB surveys will provide hundreds of new FRBs with high-precision localizations, RMs, and imaging follow-up to support further investigation into the magnetic fields of a diverse population of z < 1 galaxies.

Quiescent Low-mass Galaxies Observed by JWST in the Epoch of Reionization

The surprising JWST discovery of a quiescent, low-mass (M ⋆ = 108.7 M ⊙) galaxy at redshift z = 7.3 (JADES-GS-z7-01-QU) represents a unique opportunity to study the imprint of feedback processes on early galaxy evolution. We build a sample of 130 low-mass (M ⋆ ≲ 109.5 M ⊙) galaxies from the serra cosmological zoom-in simulations, which show a feedback-regulated, bursty star formation history (SFH). The fraction of time spent in an active phase increases with the stellar mass from f duty ≈ 0.6 at M ⋆ ≈ 107.5 M ⊙ to ≈0.99 at M ⋆ ≥ 109 M ⊙, and it is in agreement with the value f duty ≈ 0.75 estimated for JADES-GS-z7-01-QU. On average, 30% of the galaxies are quiescent in the range 6 < z < 8.4; they become the dominant population at M ⋆ ≲ 108.3 M ⊙. However, none of these quiescent systems matches the spectral energy distribution of JADES-GS-z7-01-QU, unless their SFH is artificially truncated a few Myr after the main star formation peak. As supernova feedback can only act on a longer timescale (≳30 Myr), this implies that the observed abrupt quenching must be caused by a faster physical mechanism, such as radiation-driven winds from young massive stars and/or an active galactic nucleus.

Born this way: thin disc, thick disc, and isotropic spheroid formation in FIRE-2 Milky Way–mass galaxy simulations

ABSTRACT We investigate the formation of Milky Way–mass galaxies using FIRE-2 ΛCDM cosmological zoom-in simulations by studying the orbital evolution of stars formed in the main progenitor of the galaxy, from birth to the present day. We classify in situ stars as isotropic spheroid, thick-disc, and thin-disc according to their orbital circularities and show that these components are assembled in a time-ordered sequence from early to late times, respectively. All simulated galaxies experience an early phase of bursty star formation that transitions to a late-time steady phase. This transition coincides with the time that the inner CGM virializes. During the early bursty phase, galaxies have irregular morphologies and new stars are born on radial orbits; these stars evolve into an isotropic spheroidal population today. The bulk of thick-disc stars form at intermediate times, during a clumpy-disc ‘spin-up’ phase, slightly later than the peak of spheroid formation. At late times, once the CGM virializes and star formation ‘cools down,’ stars are born on circular orbits within a narrow plane. Those stars mostly inhabit thin discs today. Broadly speaking, stars with disc-like or spheroid-like orbits today were born that way. Mergers on to discs and secular processes do affect kinematics in our simulations, but play only secondary roles in populating thick-disc and in situ spheroid populations at z = 0. The age distributions of spheroid, thick disc, and thin disc populations scale self-similarly with the steady-phase transition time, which suggests that morphological age dating can be linked to the CGM virialization time in galaxies.

Spatially resolved mock observations of stellar kinematics: full radiative transfer treatment of simulated galaxies

ABSTRACT We present a framework to build realistic mock spectroscopic observations for state-of-the-art hydrodynamical simulations, using high spectral resolution stellar population models and full radiative transfer treatment with skirt. As a first application, we generate stellar continuum mock observations for the Auriga cosmological zoom simulations emulating integral-field observations from the Sydney–AOO Multi-object Integral Field Spectrograph (SAMI) Galaxy Survey. We perform spectral fitting on our synthetic cubes and compute the resulting rotation velocity (Vrot) and velocity dispersion within 1Re (σe) for a subset of the Auriga sample. We find that the kinematics produced by Auriga are in good agreement with the observations from the SAMI galaxy survey after taking into account the effects of dust and the systematics produced by the observation limitations. We also explore the effects of seeing convolution, inclination, and attenuation on the line-of-sight velocity distribution. For highly inclined galaxies, these effects can lead to an artificial decrease in the measured V/σ by nearly a factor of two (after inclination correction). We also demonstrate the utility of our method for high-redshift galaxies by emulating spatially resolved continuum spectra from the Large Early Galaxy Census (LEGA-C) survey and, looking forward, the European Extremely Large Telescope (E-ELT) High Angular Resolution Monolithic Optical and Near-infrared Integral field spectrograph (HARMONI). Our framework represents a crucial link between the ground truth for stellar populations and kinematics in simulations and the observed stellar continuum observations at low and high redshift.

Cosmic Sands: The Origin of Dusty, Star-forming Galaxies in the Epoch of Reionization

We present the Cosmic Sands suite of cosmological zoom-in simulations based on the simba galaxy formation model in order to study the buildup of the first massive and dusty galaxies in the early universe. Residing in the most massive halos, we find that the compact proto-massive galaxies undergo nearly continuous mergers with smaller subhalos, boosting star formation rates (SFRs) and the buildup of stellar mass. The galaxies are already appreciably chemically evolved by z = 7, with modeled dust masses comparable to those inferred from observations in the same epoch, except for the most extreme systems. We track gas accretion onto the galaxies to understand how extreme SFRs can be sustained by these early systems. We find that smooth gas accretion can maintain SFRs above 250 M ⊙ yr−1, but to achieve SFRs that boost galaxies well above the main sequence, a larger perturbation like a gas-rich major merger is necessary to trigger a starburst episode. Post-processing the Cosmic Sands simulations with dust RT, we find that, while the infrared luminosities of the most-dust-rich galaxies are comparable to local ULIRGs, they are substantially dimmer than classical z = 2 submillimeter galaxies. We end with a discussion on the possible reasons for this discrepancy at the highest masses and the future work we intend to carry out to study the chemical enrichment of the earliest dusty galaxies.

Strong Dark Matter Self-interactions Diversify Halo Populations within and surrounding the Milky Way

We perform a high-resolution cosmological zoom-in simulation of a Milky Way (MW)–like system, which includes a realistic Large Magellanic Cloud analog, using a large differential elastic dark matter self-interaction cross section that reaches ≈100 cm2 g−1 at relative velocities of ≈10 km s−1, motivated by the diverse and orbitally dependent central densities of dwarf galaxies within and surrounding the MW. We explore the effects of dark matter self-interactions on satellite, splashback, and isolated halos through their abundance, central densities, maximum circular velocities, orbital parameters, and correlations between these variables. We use an effective constant cross section model to analytically predict the stages of our simulated halos’ gravothermal evolution, demonstrating that deviations from the collisionless – relation can be used to select deeply core-collapsed halos, where is a halo’s maximum circular velocity, and is the radius at which it occurs. We predict that a sizable fraction (≈20%) of subhalos with masses down to ≈108 M ⊙ is deeply core collapsed in our SIDM model. Core-collapsed systems form ≈10% of the isolated halo population down to the same mass; these isolated, core-collapsed halos would host faint dwarf field galaxies with extremely steep central density profiles. Finally, most halos with masses above ≈109 M ⊙ are core-forming in our simulation. Our study thus demonstrates how self-interactions diversify halo populations in an environmentally dependent fashion within and surrounding MW-mass hosts, providing a compelling avenue to address the diverse dark matter distributions of observed dwarf galaxies.

Dark matter subhalo disturbance in the tidal field: Towards a new simulations accuracy

Dark matter (DM) identification is a pressing challenge for modern particle physics and cosmology. The success of this endeavor, however, is contingent on a number of parameters, some of which stem from the difficulty of predicting the DM sub-halo distribution in Milky Way-like galaxies and others from the effect of DM sub-halo signals on the gamma-ray instruments being utilized. This research compared the results of detailed simulations of individual DM sub-halos emerging in the external tidal fields with and without disc components and a baryonic bulge. The average DM particle is assumed to have a mass three times smaller than that assumed in cosmological zoom-in simulations of galaxy formation. The Via Lactea II simulation, which provides the basis for simulating the sub-haloes' tidal field in DM, was used to establish the baseline conditions for this investigation. A model of the Milky Way based on observations was used to simulate the tidal field; this model included a DM halo, a stellar bulge, and a stellar disc. When small (rp 24 kpc) percenters were included, the number of sub-haloes dropped to 39%, and the remaining sub-haloes were found to have a less dense orbital distribution. Since the total number of outer sub-haloes was only around 74% of the value anticipated by CDM, the addition of a stellar component to the tidal field had little effect on the sub-haloes with greater percenters. Previous research on the impacts of baryonic matter as mediated via the spurious disruption had reported ratios that were lower than the ones obtained in this study, suggesting that the latter's accuracy was to blame.

Rhapsody-Csimulations – anisotropic thermal conduction, black hole physics, and the robustness of massive galaxy cluster scaling relations

ABSTRACTWe present the Rhapsody-C simulations that extend the Rhapsody-G suite of massive galaxy clusters at the $M_{\rm vir}\sim 10^{15}\, {\rm M}_{\odot }$ scale with cosmological magneto-hydrodynamic zoom-in simulations that include anisotropic thermal conduction, modified supermassive black hole (SMBH) feedback, new SMBH seeding, and the SMBH orbital decay model. These modelling improvements have a dramatic effect on the SMBH growth, star formation, and gas depletion in the proto-clusters. We explore the parameter space of the models and report their effect on both star formation and the thermodynamics of the intracluster medium (ICM) as observed in X-ray and SZ observations. We report that the star formation in proto-clusters is strongly impacted by the choice of the SMBH seeding as well as the orbital decay of SMBHs. The feedback from AGNs is substantially boosted by the SMBH decay; its time evolution and impact range differ noticeably depending on the AGN energy injection scheme used. Compared to a mass-weighted injection whose energy remains confined close to the central SMBHs, a volume-weighted thermal energy deposition allows to heat the ICM out to large radii that severely quenches the star formation in proto-clusters. By flattening out temperature gradients in the ICM, anisotropic thermal conduction can reduce star formation early on but weakens and delays the AGN activity. Despite the dissimilarities found in the stellar and gaseous content of our haloes, the cluster scaling relations we report are surprisingly insensitive to the subresolution models used and are in good agreement with recent observational and numerical studies.

Public Data Release of the FIRE-2 Cosmological Zoom-in Simulations of Galaxy Formation

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

The connection between the escape of ionizing radiation and galaxy properties at z ∼ 3 in the Keck Lyman continuum spectroscopic survey

ABSTRACT The connection between the escape fraction of ionizing radiation (fesc) and the properties of galaxies, such as stellar mass ($\rm M_{\rm *}$), age, star-formation rate (SFR), and dust content, are key inputs for reionization models, but many of these relationships remain untested at high redshift. We present an analysis of a sample of 96 $z$ ∼ 3 galaxies from the Keck Lyman Continuum Spectroscopic Survey (KLCS). These galaxies have both sensitive Keck/LRIS spectroscopic measurements of the Lyman continuum (LyC) region, and multiband photometry that places constraints on stellar population parameters. We construct composite spectra from subsamples binned as a function of galaxy property and quantify the ionizing-photon escape for each composite. We find a significant anti-correlation between fesc and $\rm M_{\rm *}$, consistent with predictions from cosmological zoom-in simulations. We also find significant anti-correlation between fesc and E(B−V), encoding the underlying physics of LyC escape in our sample. We also find no significant correlation between fesc and either stellar age or specific SFR (= SFR/$\rm M_{\rm *}$), challenging interpretations that synchronize recent star formation and favorable conditions for ionizing escape. The galaxy properties now shown to correlate with fesc in the KLCS are Lyα equivalent width, UV Luminosity, $\rm M_{\rm *}$, SFR, and E(B−V), but not age or sSFR. This comprehensive analysis of galaxy properties and LyC escape at high redshift will be used to guide future models and observations of the reionization epoch.

Firefly: A Browser-based Interactive 3D Data Visualization Tool for Millions of Data Points

We present Firefly, a new browser-based interactive tool for visualizing 3D particle data sets. On a typical personal computer, Firefly can simultaneously render and enable real-time interactions with ≳10 million particles, and can interactively explore data sets with billions of particles using the included custom-built octree render engine. Once created, viewing a Firefly visualization requires no installation and is immediately usable in most modern internet browsers simply by visiting a URL. As a result, a Firefly visualization works out-of-the-box on most devices including smartphones and tablets. Firefly is primarily developed for researchers to explore their own data, but can also be useful to communicate results to researchers and/or collaborators and as an effective public outreach tool. Every element of the user interface can be customized and disabled, enabling easy adaptation of the same visualization for different audiences with little additional effort. Creating a new Firefly visualization is simple with the provided Python data preprocessor that translates input data to a Firefly-compatible format and provides helpful methods for hosting instances of Firefly both locally and on the internet. In addition to visualizing the positions of particles, users can visualize vector fields (e.g., velocities) and also filter and color points by scalar fields. We share three examples of Firefly applied to astronomical data sets: (1) the FIRE cosmological zoom-in simulations, (2) the SDSS galaxy catalog, and (3) Gaia Data Release 3. A gallery of additional interactive demos is available at alexbgurvi.ch/Firefly.

Chasing the impact of the Gaia-Sausage-Enceladus merger on the formation of the Milky Way thick disc

ABSTRACT We employ our Bayesian Machine Learning framework BINGO (Bayesian INference for Galactic archaeOlogy) to obtain high-quality stellar age estimates for 68 360 red giant and red clump stars present in the 17th data release of the Sloan Digital Sky Survey, the APOGEE-2 high-resolution spectroscopic survey. By examining the denoised age-metallicity relationship of the Galactic disc stars, we identify a drop in metallicity with an increase in [Mg/Fe] at an early epoch, followed by a chemical enrichment episode with increasing [Fe/H] and decreasing [Mg/Fe]. This result is congruent with the chemical evolution induced by an early-epoch gas-rich merger identified in the Milky Way-like zoom-in cosmological simulation Auriga. In the initial phase of the merger of Auriga 18 there is a drop in metallicity due to the merger diluting the metal content and an increase in the [Mg/Fe] of the primary galaxy. Our findings suggest that the last massive merger of our Galaxy, the Gaia-Sausage-Enceladus, was likely a significant gas-rich merger and induced a starburst, contributing to the chemical enrichment and building of the metal-rich part of the thick disc at an early epoch.

Spiral arms are metal freeways: azimuthal gas-phase metallicity variations in flocculent discs in the FIRE-2 cosmological zoom-in simulations

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