Search

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

Milky Way-mass galaxies in the FIRE-2 simulations demonstrate two main modes of star formation. At high redshifts star formation occurs in a series of short and intense bursts, while at low redshifts star formation proceeds at a steady rate with a transition from one mode to another at times ranging from 3 to 7 Gyr ago for different galaxies. We analyse how the mode of star formation affects iron and alpha-element abundance. We find that the early bursty regime imprints a measurable pattern in stellar elemental abundances in the form of a "sideways chevron" shape on the [Fe/H] - [O/Fe] plane and the scatter in [O/Fe] at a given stellar age is higher than when a galaxy is in the steady regime. That suggests that the evolution of [O/Fe] scatter with age provides an estimate of the end of the bursty phase. We investigate the feasibility of observing of this effect by adding mock observational errors to a simulated stellar survey and find that the transition between the bursty and steady phase should be detectable in the Milky Way, although larger observational uncertainties make the transition shallower. We apply our method to observations of the Milky Way from the Second APOKASC Catalog and estimate that the transition to steady star formation in the Milky Way happened 7-8 Gyrs ago, earlier than transition times measured in the simulations., Comment: 15 pages; Accepted 17 January 2023 in MNRAS

Several recent simulations of galaxy formation predict two main phases of supermassive black hole (BH) accretion: an early, highly intermittent phase (during which BHs are under-massive relative to local scaling relations), followed by a phase of accelerated growth. We investigate physical factors that drive the transition in BH accretion in cosmological zoom-in simulations from the FIRE project, ranging from dwarf galaxies to galaxies sufficiently massive to host luminous quasars. The simulations model multi-channel stellar feedback, but neglect AGN feedback. We show that multiple physical properties, including halo mass, galaxy stellar mass, and depth of the central gravitational potential correlate with accelerated BH fueling: constant thresholds in these properties are typically crossed within ~0.1 Hubble time of accelerated BH fueling. Black hole masses increase sharply when the stellar surface density in the inner 1 kpc crosses a threshold Sigma1 ~ 10^9.5 Msun/kpc^2, a characteristic value above which gravity prevents stellar feedback from ejecting gas, and similar to the value above which galaxies are observed to quench. We further show that accelerated BH growth correlates with the emergence of long-lived, thin gas disks, as well as with virialization of the inner circumgalactic medium. The halo mass Mh ~ 10^12 Msun and stellar mass Mstar ~ 10^10.5 Msun at which BH growth accelerates correspond to ~L* galaxies. The fact that stellar feedback becomes inefficient at ejecting gas from the nucleus above this mass scale may play an important role in explaining why AGN feedback appears to be most important in galaxies above ~L*., Comment: 21 pages, 9 figures, submitted to MNRAS

We investigate the formation of Milky-Way-mass galaxies using FIRE-2 LCDM 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 onto 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., Comment: 16 pages, 10 figures; this version was accepted to MNRAS on June 9, 2023

Interstellar chemistry is important for galaxy formation, as it determines the rate at which gas can cool, and enables us to make predictions for observable spectroscopic lines from ions and molecules. We explore two central aspects of modelling the chemistry of the interstellar medium (ISM): (1) the effects of local stellar radiation, which ionises and heats the gas, and (2) the depletion of metals onto dust grains, which reduces the abundance of metals in the gas phase. We run high-resolution (400 M$_\odot$ per baryonic particle) simulations of isolated disc galaxies, from dwarfs to Milky Way-mass, using the FIRE galaxy formation models together with the CHIMES non-equilibrium chemistry and cooling module. In our fiducial model, we couple the chemistry to the stellar fluxes calculated from star particles using an approximate radiative transfer scheme, and we implement an empirical density-dependent prescription for metal depletion. For comparison, we also run simulations with a spatially uniform radiation field, and without metal depletion. Our fiducial model broadly reproduces observed trends in HI and H2 mass with stellar mass, and in line luminosity versus star formation rate for [CII] 158$\mu$m, [OI] 63$\mu$m, [OIII] 88$\mu$m, [NII] 122$\mu$m and H$\alpha$ 6563A. Our simulations with a uniform radiation field predict fainter luminosities, by up to an order of magnitude for [OIII] 88$\mu$m and H$\alpha$ 6563A, while ignoring metal depletion increases the luminosity of carbon and oxygen lines by a factor $\approx$2. However, the overall evolution of the galaxy is not strongly affected by local stellar fluxes or metal depletion, except in dwarf galaxies where the inclusion of local fluxes leads to weaker outflows and hence higher gas fractions., Comment: 29 pages, 24 figures (including appendices). Accepted for publication in MNRAS. Minor changes relative to v1

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 datasets 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/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 pre-processor (PDPP) 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 datasets: 1) the FIRE cosmological zoom-in simulations, 2) the SDSS galaxy catalog, and 3) Gaia DR3. A gallery of additional interactive demos is available at https://alexbgurvi.ch/Firefly., Comment: 26 pages, 9 figures

We study the sensitivity of the $z=0.1$ Lyman-$\alpha$ Forest observables, such as the column density distribution function (CDD), flux PDF, flux power spectrum, and line width distribution, to sub-grid models of active galactic nuclei (AGN) feedback using the Illustris and IllustrisTNG (TNG) cosmological simulations. The two simulations share an identical Ultraviolet Background (UVB) prescription and similar cosmological parameters, but TNG features an entirely reworked AGN feedback model. Due to changes in the AGN radio mode model, the original Illustris simulations have a factor of 2-3 fewer Lyman-$\alpha$ absorbers than TNG at column densities $N_{\rm HI}< 10^{15.5}$ cm$^{-2}$. We compare the simulated forest statistics to UV data from the Cosmic Origins Spectrograph (COS) and find that neither simulation can reproduce the slope of the absorber distribution. Both Illustris and TNG also produce significantly smaller line width distributions than observed in the COS data. We show that TNG is in much better agreement with the observed $z=0.1$ flux power spectrum than Illustris. We explore which statistics can disentangle the effects of AGN feedback from alternative UVB models by rescaling the UVB of Illustris to produce a CDD match to TNG. While this UVB rescaling is degenerate with the effect of AGN feedback on the CDD, the amplitude and shape of the flux PDF and 1D flux power spectrum change in a way distinct from a scaling of the UVB. Our study suggests that the $z=0.1$ Lyman-$\alpha$ forest observables can be used as a diagnostic of AGN feedback models., Comment: ApJL

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., Comment: 16 pages, 14 figures, submitted to MNRAS

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

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., Comment: 55 pages, 4 figures, 13 supplementary figures. 29 pages, 4 figures. Accepted for publication in Nature Astronomy. To appear on 14-February-2022. Published version: https://www.nature.com/articles/s41550-021-01598-4

Cosmic rays (CRs) are an important component in the interstellar medium (ISM), but their effect on the dynamics of the disk-halo interface (< 10 kpc from the disk) is still unclear. We study the influence of CRs on the gas above the disk with high-resolution FIRE-2 cosmological simulations of late-type Lstar galaxies at redshift around zero. We compare runs with and without CR feedback (with constant anisotropic diffusion around 3e29 cm^2/s and streaming). Our simulations capture the relevant disk halo interactions, including outflows, inflows, and galactic fountains. Extra-planar gas in all of the runs satisfies dynamical balance, where total pressure balances the weight of the overlying gas. While the kinetic pressure from non-uniform motion (>1-kpc scale) dominates in the midplane, thermal and bulk pressures (or CR pressure if included) take over at large heights. We find that with CR feedback, (1) the warm (1e4K) gas is slowly accelerated by CRs; (2) the hot (> 5e5K) gas scale height is suppressed; (3) the warm-hot (2e4-5e5K) medium becomes the most volume-filling phase in the disk-halo interface. We develop a novel conceptual model of the near-disk gas dynamics in low-redshift Lstar galaxies: With CRs, the disk-halo interface is filled with CR-driven warm winds and hot super-bubbles that are propagating into the CGM with a small fraction falling back to the disk. Without CRs, most outflows from hot superbubbles are trapped by the existing hot halo and gravity, so typically they form galactic fountains., Comment: 21 pages, 15 figures. Accepted for publication by MNRAS

We investigate thin and thick stellar disc formation in Milky-Way-mass galaxies using twelve FIRE-2 cosmological zoom-in simulations. All simulated galaxies experience an early period of bursty star formation that transitions to a late-time steady phase of near-constant star formation. Stars formed during the late-time steady phase have more circular orbits and thin-disc-like morphology at $z=0$, whilst stars born during the bursty phase have more radial orbits and thick-disc structure. The median age of thick-disc stars at $z=0$ correlates strongly with this transition time. We also find that galaxies with an earlier transition from bursty to steady star formation have a higher thin-disc fractions at $z=0$. Three of our systems have minor mergers with LMC-size satellites during the thin-disc phase. These mergers trigger short starbursts but do not destroy the thin disc nor alter broad trends between the star formation transition time and thin/thick disc properties. If our simulations are representative of the Universe, then stellar archaeological studies of the Milky Way (or M31) provide a window into past star-formation modes in the Galaxy. Current age estimates of the Galactic thick disc would suggest that the Milky Way transitioned from bursty to steady phase $\sim$6.5 Gyr ago; prior to that time the Milky Way likely lacked a recognisable thin disc., Comment: 15 pages, 11 figures, submitted to MNRAS

We carry out a comparative analysis of the metallicities from the stellar, neutral-gas, and ionized-gas components in the metal-rich spiral galaxy M83. We analyze spectroscopic observations taken with the Hubble Space Telescope (HST), the Large Binocular Telescope (LBT) and the Very Large Telescope (VLT). We detect a clear depletion of the HI gas, as observed from the HI column densities in the nuclear region of this spiral galaxy. We find column densities of log[$N$(HI) cm$^{-2}$] $<$ 20.0 at galactocentric distances of $<$ 0.18 kpc, in contrast to column densities of log[$N$(HI) cm$^{-2}$] $\sim$ 21.0 in the galactic disk, a trend observed in other nearby spiral galaxies. We measure a metallicity gradient of $-$0.03 $\pm$ 0.01 dex kpc$^{-1}$ for the ionized gas, comparable to the metallicity gradient of a local benchmark of 49 nearby star-forming galaxies of $-$0.026 $\pm$ 0.002 dex kpc$^{-1}$. Our co-spatial metallicity comparison of the multi-phase gas and stellar populations shows excellent agreement outside of the nucleus of the galaxy hinting at a scenario where the mixing of newly synthesized metals from the most massive stars in the star clusters takes longer than their lifetimes ($\sim$10 Myr). Finally, our work shows that caution must be taken when studying the metallicity gradient of the neutral-gas component in star-forming galaxies, since this can be strongly biased, as these environments can be dominated by molecular gas. In these regions the typical metallicity tracers can provide inaccurate abundances as they may trace both the neutral- and molecular-gas components., Comment: 29 pages, 14 tables, 9 figures, Accepted for Publication in ApJ

Understanding the rate at which stars form is central to studies of galaxy formation. Observationally, the star formation rates (SFRs) of galaxies are measured using the luminosity in different frequency bands, often under the assumption of a time-steady SFR in the recent past. We use star formation histories (SFHs) extracted from cosmological simulations of star-forming galaxies from the FIRE project to analyze the time-scales to which the H${\alpha}$ and far-ultraviolet (FUV) continuum SFR indicators are sensitive. In these simulations, the SFRs are highly time variable for all galaxies at high redshift, and continue to be bursty to z=0 in dwarf galaxies. When FIRE SFHs are partitioned into their bursty and time-steady phases, the best-fitting FUV time-scale fluctuates from its ~10 Myr value when the SFR is time-steady to >~100 Myr immediately following particularly extreme bursts of star formation during the bursty phase. On the other hand, the best-fitting averaging time-scale for H${\alpha}$ is generally insensitive to the SFR variability in the FIRE simulations and remains ~5 Myr at all times. These time-scales are shorter than the 100 Myr and 10 Myr time-scales sometimes assumed in the literature for FUV and H${\alpha}$, respectively, because while the FUV emission persists for stellar populations older than 100 Myr, the time-dependent luminosities are strongly dominated by younger stars. Our results confirm that the ratio of SFRs inferred using H${\alpha}$ vs. FUV can be used to probe the burstiness of star formation in galaxies., Comment: 14 pages, 10 figures, accepted to MNRAS

We use the FIRE-2 cosmological simulations to study the formation of a quasi-static, virial-temperature gas phase in the circumgalactic medium (CGM) at redshifts 0

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