Showing 1–21 of 21 results for author: Orr, M E

Any Way the Wind Blows: Quantifying Superbubbles and their Outflows in Simulated Galaxies across $z \approx 0-3$

Authors: Lori E. Porter, Matthew E. Orr, Blakesley Burkhart, Andrew Wetzel, Dušan Kereš, Claude-André Faucher-Giguère, Philip F. Hopkins

Abstract: We present an investigation of clustered stellar feedback in the form of superbubbles identified within eleven galaxies from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, at both cosmic noon (1 < z < 3) and in the local Universe. We study the spatially-resolved multiphase outflows that these supernovae drive, comparing our findings with recent theory and ob… ▽ More We present an investigation of clustered stellar feedback in the form of superbubbles identified within eleven galaxies from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, at both cosmic noon (1 < z < 3) and in the local Universe. We study the spatially-resolved multiphase outflows that these supernovae drive, comparing our findings with recent theory and observations. These simulations consist of five LMC-mass galaxies and six Milky Way-mass progenitors (with a minimum baryonic particle mass of $m_{b.min} = 7100 M_{\odot}$), for which we calculate the local mass and energy loading factors on 750~pc scales from the identified outflows. We also characterize the multiphase morphology and properties of the identified superbubbles, including the `shell' of cool ($T<10^5$ K) gas and break out of energetic hot ($T>10^5$ K) gas when the shell bursts. For all galaxies, the outflow mass, momentum, and energy fluxes appear to reach their peak during the identified superbubbles, and we investigate the effects on the interstellar medium (ISM), circumgalactic medium (CGM), and subsequent star formation rates. We find that these simulations, regardless of redshift, have mass-loading factors and momentum fluxes in the cool gas that largely agree with recent observations. Lastly, we also investigate how methodological choices in measuring outflows can affect loading factors for galactic winds. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Objects May Be Closer Than They Appear: Significant Host Galaxy Dispersion Measures of Fast Radio Bursts in Zoom-in Simulations

Authors: Matthew E. Orr, Blakesley Burkhart, Wenbin Lu, Sam B. Ponnada, Cameron B. Hummels

Abstract: We investigate the contribution of host galaxies to the overall Dispersion Measures (DMs) for Fast Radio Bursts (FRBs) using the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite. We calculate DMs from every star particle in the simulated L* galaxies by ray-tracing through their multi-phase interstellar medium (ISM), summing the line-of-sight free thermal electron c… ▽ More We investigate the contribution of host galaxies to the overall Dispersion Measures (DMs) for Fast Radio Bursts (FRBs) using the Feedback in Realistic Environments (FIRE-2) cosmological zoom-in simulation suite. We calculate DMs from every star particle in the simulated L* galaxies by ray-tracing through their multi-phase interstellar medium (ISM), summing the line-of-sight free thermal electron column for all gas elements within $\pm$20 kpc of the galactic mid-plane. At $z=0$, we find average (median) host-galaxy DMs of 74 (43) and 210 (94) pc cm$^{-3}$ for older ($\gtrsim$10 Myr) and younger ($\lesssim$10 Myr) stellar populations, respectively. Inclination raises the median DM measured for older populations ($\gtrsim$10 Myr) in the simulations by a factor of $\sim$2, but generally does not affect the younger stars deeply embedded in H{\small II} regions except in extreme edge-on cases (inclination $\gtrsim 85^\circ$). In kinematically disturbed snapshots ($z = 1$ in FIRE), the average (median) host-galaxy DMs are higher: 80 (107) and 266 (795) pc cm$^{-3}$ for older ($\gtrsim$10 Myr) and younger ($\lesssim$10 Myr) stellar populations, respectively. FIRE galaxies tend to have higher DM values than cosmological simulations such as IllustrisTNG. As a result, FRB host galaxies may be closer (lower redshift) than previously inferred. Furthermore, constraining host-galaxy DM distributions may help significantly constrain FRB progenitor models. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 17 pages, 12 figures, submitted to ApJ Letters

What Causes The Formation of Disks and End of Bursty Star Formation?

Authors: Philip F. Hopkins, Alexander B. Gurvich, Xuejian Shen, Zachary Hafen, Michael Y. Grudic, Shalini Kurinchi-Vendhan, Christopher C. Hayward, Fangzhou Jiang, Matthew E. Orr, Andrew Wetzel, Dusan Keres, Jonathan Stern, Claude-Andre Faucher-Giguere, James Bullock, Coral Wheeler, Kareem El-Badry, Sarah R. Loebman, Jorge Moreno, Michael Boylan-Kolchin, Eliot Quataert

Abstract: 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 existin… ▽ More 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. △ Less

Submitted 28 August, 2023; v1 submitted 19 January, 2023; originally announced January 2023.

Comments: Updated to match MNRAS accepted version. 45 pages, 32 figures

Spiral Arms are Metal Freeways: Azimuthal Gas-Phase Metallicity Variations in Simulated Cosmological Zoom-in Flocculent Disks

Authors: Matthew E. Orr, Blakesley Burkhart, Andrew Wetzel, Philip F. Hopkins, Ivanna A. Escala, Allison L. Strom, Paul F. Goldsmith, Jorge L. Pineda, Christopher C. Hayward, Sarah R. Loebman

Abstract: We examine the azimuthal variations in gas-phase metallicity profiles in simulated Milky Way mass disk 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 disks at $z \approx 0$ with pixel sizes ranging from 250 to 750~pc, analogous to modern integral f… ▽ More We examine the azimuthal variations in gas-phase metallicity profiles in simulated Milky Way mass disk 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 disks at $z \approx 0$ with pixel sizes ranging from 250 to 750~pc, analogous to modern integral field unit (IFU) galaxy surveys, mapping the gas-phase metallicities in both the cold & dense gas and the ionized gas correlated with HII regions. We report that the spiral arms alternate in a pattern of metal rich and metal poor relative to the median metallicity on the order of $\lesssim 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 \approx 0$: there is no correlation with local star formation on these spatial scales. Rather, the arms are moving 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. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 16 Pages, 14 Figures, submitted to MNRAS

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

Authors: Lori E. Porter, Matthew E. Orr, Blakesley Burkhart, Andrew Wetzel, Xiangcheng Ma, Philip F. Hopkins, Andrew Emerick

Abstract: 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,… ▽ More 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. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 14 pages, 10 figures, submitted to MNRAS

Probing Hot Gas Components of Circumgalactic Medium in Cosmological Simulations with the Thermal Sunyaev-Zel'dovich Effect

Authors: Junhan Kim, Sunil Golwala, James G. Bartlett, Stefania Amodeo, Nicholas Battaglia, Andrew J. Benson, J. Colin Hill, Philip F. Hopkins, Cameron B. Hummels, Emily Moser, Matthew E. Orr

Abstract: The thermal Sunyaev-Zel'dovich (tSZ) effect is a powerful tool with the potential for constraining directly the properties of the hot gas that dominates dark matter halos because it measures pressure and thus thermal energy density. Studying this hot component of the circumgalactic medium (CGM) is important because it is strongly impacted by star-formation and active galactic nuclei (AGN) activity… ▽ More The thermal Sunyaev-Zel'dovich (tSZ) effect is a powerful tool with the potential for constraining directly the properties of the hot gas that dominates dark matter halos because it measures pressure and thus thermal energy density. Studying this hot component of the circumgalactic medium (CGM) is important because it is strongly impacted by star-formation and active galactic nuclei (AGN) activity in galaxies, participating in the feedback loop that regulates star and black hole mass growth in galaxies. We study the tSZ effect across a wide halo mass range using three cosmological hydrodynamical simulations: Illustris-TNG, EAGLE, and FIRE-2. Specifically, we present the scaling relation between tSZ signal and halo mass and radial profiles of gas density, temperature, and pressure for all three simulations. The analysis includes comparisons to Planck tSZ observations and to the thermal pressure profile inferred from the Atacama Cosmology Telescope (ACT) measurements. We compare these tSZ data to the simulations to interpret the measurements in terms of feedback and accretion processes in the CGM. We also identify as-yet unobserved potential signatures of these processes that may be visible in future measurements, which will have the capability of measuring tSZ signals to even lower masses. We also perform internal comparisons between runs with different physical assumptions. We conclude: (1) there is strong evidence for the impact of feedback at $R_{500}$ but that this impact decreases by $5R_{500}$, and (2) the thermodynamic profiles of the CGM are highly dependent on the implemented model, such as cosmic-ray or AGN feedback prescriptions. △ Less

Submitted 28 October, 2021; originally announced October 2021.

Comments: 21 pages, 10 figures, submitted to ApJ, comments welcome

Bursting Bubbles: Feedback from Clustered SNe and the Trade-off Between Turbulence and Outflows

Authors: Matthew E. Orr, Drummond B. Fielding, Christopher C. Hayward, Blakesley Burkhart

Abstract: We present an analytic model for clustered supernovae (SNe) feedback in galaxy disks, incorporating the dynamical evolution of superbubbles formed from spatially overlapping SNe remnants. We propose two realistic outcomes for the evolution of superbubbles in galactic disks: (1) the expansion velocity of the shock front falls below the turbulent velocity dispersion of the ISM in the galaxy disk, wh… ▽ More We present an analytic model for clustered supernovae (SNe) feedback in galaxy disks, incorporating the dynamical evolution of superbubbles formed from spatially overlapping SNe remnants. We propose two realistic outcomes for the evolution of superbubbles in galactic disks: (1) the expansion velocity of the shock front falls below the turbulent velocity dispersion of the ISM in the galaxy disk, whereupon the superbubble stalls and fragments, depositing its momentum entirely within the galaxy disk, or (2) the superbubble grows in size to reach the gas scale height, breaking out of the galaxy disk and driving galactic outflows/fountains. In either case, we find that superbubble breakup/breakout almost always occurs before the last Type-II SN ($\lesssim$40 Myr) in the recently formed star cluster, assuming a standard high-end IMF slope, and scalings between stellar lifetimes and masses. The threshold between these two cases implies a break in the effective strength of feedback in driving turbulence within galaxies, and a resulting change in the scalings of, e.g., star formation rates with gas surface density (the Kennicutt-Schmidt relation) and the star formation efficiency in galaxy disks. △ Less

Submitted 2 May, 2022; v1 submitted 29 September, 2021; originally announced September 2021.

Comments: 24 pages, 11 figures, accepted to ApJ

Bursting Bubbles: Clustered Supernova Feedback in Local and High-redshift Galaxies

Authors: Matthew E. Orr, Drummond B. Fielding, Christopher C. Hayward, Blakesley Burkhart

Abstract: We compare an analytic model for the evolution of supernova-driven superbubbles with observations of local and high-redshift galaxies, and the properties of intact HI shells in local star-forming galaxies. Our model correctly predicts the presence of superwinds in local star-forming galaxies (e.g., NGC 253) and the ubiquity of outflows near $z \sim 2$. We find that high-redshift galaxies may `capt… ▽ More We compare an analytic model for the evolution of supernova-driven superbubbles with observations of local and high-redshift galaxies, and the properties of intact HI shells in local star-forming galaxies. Our model correctly predicts the presence of superwinds in local star-forming galaxies (e.g., NGC 253) and the ubiquity of outflows near $z \sim 2$. We find that high-redshift galaxies may `capture' 20-50\% of their feedback momentum in the dense ISM (with the remainder escaping into the nearby CGM), whereas local galaxies may contain $\lesssim$10\% of their feedback momentum from the central starburst. Using azimuthally averaged galaxy properties, we predict that most superbubbles stall and fragment \emph{within} the ISM, and that this occurs at, or near, the gas scale height. We find a consistent interpretation in the observed HI bubble radii and velocities, and predict that most will fragment within the ISM, and that those able to break-out originate from short dynamical time regions (where the dynamical time is shorter than feedback timescales). Additionally, we demonstrate that models with constant star cluster formation efficiency per Toomre mass are inconsistent with the occurrence of outflows from high-$z$ starbursts and local circumnuclear regions. △ Less

Submitted 21 January, 2022; v1 submitted 29 September, 2021; originally announced September 2021.

Comments: 10 pages, 4 figures, accepted to ApJL

Journal ref: ApJL 924 L28 (2022)

Which AGN Jets Quench Star Formation in Massive Galaxies?

Authors: Kung-Yi Su, Philip F. Hopkins, Greg L. Bryan, Rachel S. Somerville, Christopher C. Hayward, Daniel Anglés-Alcázar, Claude-André Faucher-Giguère, Sarah Wellons, Jonathan Stern, Bryan A. Terrazas, T. K. Chan, Matthew E. Orr, Cameron Hummels, Robert Feldmann, Dušan Kereš

Abstract: Without additional heating, radiative cooling of gas in the halos of massive galaxies (Milky Way and above) produces cold gas or stars in excess of that observed. Previous work suggested that AGN jets are likely required, but the form of jet energy required to quench remains unclear. This is particularly challenging for galaxy simulations, in which the resolution is orders of magnitude coarser tha… ▽ More Without additional heating, radiative cooling of gas in the halos of massive galaxies (Milky Way and above) produces cold gas or stars in excess of that observed. Previous work suggested that AGN jets are likely required, but the form of jet energy required to quench remains unclear. This is particularly challenging for galaxy simulations, in which the resolution is orders of magnitude coarser than necessary to form and evolve the jet. On such scales, the uncertain parameters include: jet energy form (kinetic, thermal, and cosmic ray (CR) energy), energy, momentum, and mass flux, magnetic field strength and geometry, jet precession angle and period, opening-angle, and duty cycle. We investigate all of these parameters in a $10^{14}\,{\rm M}_{\odot}$ halo using high-resolution non-cosmological MHD simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We explore which scenarios match observational constraints and show that CR-dominated jets can most efficiently quench the central galaxy through a combination of CR pressure support and a modification of the thermal instability. Jets with most energy in mildly relativistic ($\sim$ MeV or $\sim10^{10}$ K) thermal plasma work, but require a factor $\sim 10$ larger energy input. For a fixed energy flux, jets with higher specific energy (longer cooling times) quench more effectively. For this halo size, kinetic jets are less efficient in quenching unless they have wide opening or precession angles. Magnetic fields play a minor role except when the magnetic flux reaches $\gtrsim 10^{44}$ erg s$^{-1}$ in a kinetic jet model, which causes the jet cocoon to significantly widen, and the quenching to become explosive. We conclude that the criteria for a successful jet model are an optimal energy flux and a sufficiently wide jet cocoon with long enough cooling time at the cooling radius. △ Less

Submitted 3 February, 2021; originally announced February 2021.

Comments: 29 pages, 20 figures

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

Authors: Matthew E. Orr, H Perry Hatchfield, Cara Battersby, Christopher C. Hayward, Philip F. Hopkins, Andrew Wetzel, Samantha M. Benincasa, Sarah R. Loebman, Mattia C. Sormani, Ralf S. Klessen

Abstract: 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's Central Molecular Zone (… ▽ More 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'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's low star formation efficiency. △ Less

Submitted 26 January, 2021; originally announced January 2021.

Comments: 9 pages, 3 figures, accepted to ApJ Letters. One appended supplemental appendix and figure, which do not appear in the ApJ version

Pressure balance in the multiphase ISM of cosmologically simulated disk galaxies

Authors: Alexander B. Gurvich, Claude-André Faucher-Giguère, Alexander J. Richings, Philip F. Hopkins, Michael Y. Grudić, Zachary Hafen, Sarah Wellons, Jonathan Stern, Eliot Quataert, T. K. Chan, Matthew E. Orr, Dušan Kereš, Andrew Wetzel, Christopher C. Hayward, Sarah R. Loebman, Norman Murray

Abstract: 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. W… ▽ More 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. △ Less

Submitted 16 November, 2020; v1 submitted 26 May, 2020; originally announced May 2020.

Comments: 22 pages, 16 figures, accepted in MNRAS

Reproducing the CO-to-H$_2$ conversion factor in cosmological simulations of Milky Way-mass galaxies

Authors: Laura C. Keating, Alexander J. Richings, Norman Murray, Claude-Andre Faucher-Giguere, Philip F. Hopkins, Andrew Wetzel, Dusan Keres, Samantha Benincasa, Robert Feldmann, Sarah Loebman, Matthew E. Orr

Abstract: We present models of CO(1-0) emission from Milky Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results… ▽ More We present models of CO(1-0) emission from Milky Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1-0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which in our models sets the attenuation of the incident UV radiation field. At the resolution of these simulations, commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H$_2$ abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of $\sim 3$ pc in cold gas ($T < 300$ K) for both CO and H$_{2}$, is able to reproduce both the observed CO(1-0) luminosity and inferred CO-to-H$_{2}$ conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model which controls the amount of radiation that penetrates giant molecular clouds. △ Less

Submitted 28 September, 2020; v1 submitted 22 January, 2020; originally announced January 2020.

Comments: Accepted to MNRAS. 16 pages, 10 figures

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

Authors: Maan H. Hani, Christopher C. Hayward, Matthew E. Orr, Sara L. Ellison, Paul Torrey, Norm Murray, Andrew Wetzel, Claude-André Faucher-Giguère

Abstract: The correlation between galaxies' integrated stellar masses and star formation rates (the `star formation main sequence'; 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'; rSFMS). In this work, we demonstrate that the rSFMS emerges naturally in… ▽ More The correlation between galaxies' integrated stellar masses and star formation rates (the `star formation main sequence'; 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'; 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' 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. △ Less

Submitted 4 February, 2020; v1 submitted 9 December, 2019; originally announced December 2019.

Comments: Accepted for publication in MNRAS Letters

Swirls of FIRE: Spatially Resolved Gas Velocity Dispersions and Star Formation Rates in FIRE-2 Disk Environments

Authors: Matthew E. Orr, Christopher C. Hayward, Anne M. Medling, Philip F. Hopkins, Norman Murray, Jorge L. Pineda, Claude-André Faucher-Giguère, Dušan Kereš, Kung-Yi Su

Abstract: We study the spatially resolved (sub-kpc) gas velocity dispersion ($σ$)--star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way mass disk galaxies at late times. In agreement with observations, we find a relatively flat relationship, with $σ\approx 15-30$ km/s in neutral gas across 3 dex in SFRs. We show th… ▽ More We study the spatially resolved (sub-kpc) gas velocity dispersion ($σ$)--star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way mass disk galaxies at late times. In agreement with observations, we find a relatively flat relationship, with $σ\approx 15-30$ km/s in neutral gas across 3 dex in SFRs. We show that higher dense gas fractions (ratios of dense gas to neutral gas) and SFRs are correlated at constant $σ$. Similarly, lower gas fractions (ratios of gas to stellar mass) are correlated with higher $σ$ at constant SFR. The limits of the $σ$-$Σ_{\rm SFR}$ relation correspond to the onset of strong outflows. We see evidence of "on-off" cycles of star formation in the simulations, corresponding to feedback injection timescales of 10-100 Myr, where SFRs oscillate about equilibrium SFR predictions. Finally, SFRs and velocity dispersions in the simulations agree well with feedback-regulated and marginally stable gas disk (Toomre's $Q =1$) model predictions, and the data effectively rule out models assuming that gas turns into stars at (low) constant efficiency (i.e., ${\rm 1\%}$ per free-fall time). And although the simulation data do not entirely exclude gas accretion/gravitationally powered turbulence as a driver of $σ$, it appears to be strongly subdominant to stellar feedback in the simulated galaxy disks. △ Less

Submitted 31 October, 2019; originally announced November 2019.

Comments: 18 pages, 12 figures, submitted to MNRAS

Cosmic Rays or Turbulence can Suppress Cooling Flows (Where Thermal Heating or Momentum Injection Fail)

Authors: Kung-Yi Su, Philip F. Hopkins, Christopher C. Hayward, Claude-André Faucher-Giguère, Dušan Kereš, Xiangcheng Ma, Matthew E. Orr, T. K. Chan, Victor H. Robles

Abstract: The quenching `maintenance' and `cooling flow' problems are important from the Milky Way through massive cluster elliptical galaxies. Previous work has shown that some source of energy beyond that from stars and pure magnetohydrodynamic processes is required, perhaps from AGN, but even the qualitative form of this energetic input remains uncertain. Different scenarios include thermal `heating,' di… ▽ More The quenching `maintenance' and `cooling flow' problems are important from the Milky Way through massive cluster elliptical galaxies. Previous work has shown that some source of energy beyond that from stars and pure magnetohydrodynamic processes is required, perhaps from AGN, but even the qualitative form of this energetic input remains uncertain. Different scenarios include thermal `heating,' direct wind or momentum injection, cosmic ray heating or pressure support, or turbulent `stirring' of the intra-cluster medium (ICM). We investigate these in $10^{12}-10^{14}\,{\rm M}_{\odot}$ halos using high-resolution non-cosmological simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, including simplified toy energy-injection models, where we arbitrarily vary the strength, injection scale, and physical form of the energy. We explore which scenarios can quench without violating observational constraints on energetics or ICM gas. We show that turbulent stirring in the central $\sim100\,$kpc, or cosmic-ray injection, can both maintain a stable low-SFR halo for $>$Gyr timescales with modest energy input, by providing a non-thermal pressure which stably lowers the core density and cooling rates. In both cases, associated thermal-heating processes are negligible. Turbulent stirring preserves cool-core features while mixing condensed core gas into the hotter halo and is by far the most energy efficient model. Pure thermal heating or nuclear isotropic momentum injection require vastly larger energy, are less efficient in lower-mass halos, easily over-heat cores, and require fine-tuning to avoid driving unphysical temperature gradients or gas expulsion from the halo center. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: 22 pages, 16 figures

A Simple Non-equilibrium Feedback Model for Galaxy-Scale Star Formation: Delayed Feedback and SFR Scatter

Authors: Matthew E. Orr, Christopher C. Hayward, Philip F. Hopkins

Abstract: We explore a class of simple non-equilibrium star formation models within the framework of a feedback-regulated model of the ISM, applicable to kiloparsec-scale resolved star formation relations (e.g. Kennicutt-Schmidt). Combining a Toomre-Q-dependent local star formation efficiency per free-fall time with a model for delayed feedback, we are able to match the normalization and scatter of resolved… ▽ More We explore a class of simple non-equilibrium star formation models within the framework of a feedback-regulated model of the ISM, applicable to kiloparsec-scale resolved star formation relations (e.g. Kennicutt-Schmidt). Combining a Toomre-Q-dependent local star formation efficiency per free-fall time with a model for delayed feedback, we are able to match the normalization and scatter of resolved star formation scaling relations. In particular, this simple model suggests that large ($\sim$dex) variations in star formation rates (SFRs) on kiloparsec scales may be due to the fact that supernova feedback is not instantaneous following star formation. The scatter in SFRs at constant gas surface density in a galaxy then depends on the properties of feedback and when we observe its star-forming regions at various points throughout their collapse/star formation "cycles". This has the following important observational consequences: (1) the scatter and normalization of the Kennicutt-Schmidt relation are relatively insensitive to the local (small-scale) star formation efficiency, (2) but gas depletion times and velocity dispersions are; (3) the scatter in and normalization of the Kennicutt-Schmidt relation is a sensitive probe of the feedback timescale and strength; (4) even in a model where $\tilde Q_{\rm gas}$ deterministically dictates star formation locally, time evolution, variation in local conditions (e.g., gas fractions and dynamical times), and variations between galaxies can destroy much of the observable correlation between SFR and $\tilde Q_{\rm gas}$ in resolved galaxy surveys. Additionally, this model exhibits large scatter in SFRs at low gas surface densities, in agreement with observations of flat outer HI disk velocity dispersion profiles. △ Less

Submitted 30 April, 2019; v1 submitted 22 October, 2018; originally announced October 2018.

Comments: 15 pages, 6 figures, accepted by MNRAS (04/25/2019)

The failure of stellar feedback, magnetic fields, conduction, and morphological quenching in maintaining red galaxies

Authors: Kung-Yi Su, Philip F. Hopkins, Christopher C. Hayward, Xiangcheng Ma, Claude-André Faucher-Giguère, Dušan Kereš, Matthew E. Orr, Victor H. Robles

Abstract: The quenching "maintenance'" and related "cooling flow" problems are important in galaxies from Milky Way mass through clusters. We investigate this in halos with masses $\sim 10^{12}-10^{14}\,{\rm M}_{\odot}$, using non-cosmological high-resolution hydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model. We specifically focus on physics present without… ▽ More The quenching "maintenance'" and related "cooling flow" problems are important in galaxies from Milky Way mass through clusters. We investigate this in halos with masses $\sim 10^{12}-10^{14}\,{\rm M}_{\odot}$, using non-cosmological high-resolution hydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model. We specifically focus on physics present without AGN, and show that various proposed "non-AGN" solution mechanisms in the literature, including Type Ia supernovae, shocked AGB winds, other forms of stellar feedback (e.g. cosmic rays), magnetic fields, Spitzer-Braginskii conduction, or "morphological quenching" do not halt or substantially reduce cooling flows nor maintain "quenched" galaxies in this mass range. We show that stellar feedback (including cosmic rays from SNe) alters the balance of cold/warm gas and the rate at which the cooled gas within the galaxy turns into stars, but not the net baryonic inflow. If anything, outflowing metals and dense gas promote additional cooling. Conduction is important only in the most massive halos, as expected, but even at $\sim 10^{14}\,{\rm M}_{\odot}$ reduces inflow only by a factor $\sim 2$ (owing to saturation effects and anisotropic suppression). Changing the morphology of the galaxies only slightly alters their Toomre-$Q$ parameter, and has no effect on cooling (as expected), so has essentially no effect on cooling flows or maintaining quenching. This all supports the idea that additional physics, e.g., AGN feedback, must be important in massive galaxies. △ Less

Submitted 24 September, 2018; originally announced September 2018.

Comments: 16 pages, 12 figures

Discrete Effects in Stellar Feedback: Individual Supernovae, Hypernovae, and IMF Sampling in Dwarf Galaxies

Authors: Kung-Yi Su, Philip F. Hopkins, Christopher C. Hayward, Xiangcheng Ma, Michael Boylan-Kolchin, Daniel Kasen, Dušan Kereš, Claude-Andrè Faucher-Giguère, Matthew E. Orr

Abstract: Using high-resolution simulations from the FIRE-2 (Feedback In Realistic Environments) project, we study the effects of discreteness in stellar feedback processes on the evolution of galaxies and the properties of the interstellar medium (ISM). We specifically consider the discretization of supernovae (SNe), including hypernovae (HNe), and sampling the initial mass function (IMF). We study these p… ▽ More Using high-resolution simulations from the FIRE-2 (Feedback In Realistic Environments) project, we study the effects of discreteness in stellar feedback processes on the evolution of galaxies and the properties of the interstellar medium (ISM). We specifically consider the discretization of supernovae (SNe), including hypernovae (HNe), and sampling the initial mass function (IMF). We study these processes in cosmological simulations of dwarf galaxies with $z=0$ stellar masses $M_{\ast}\sim 10^{4}-3\times10^{6}\,M_\odot$ (halo masses $\sim 10^{9}-10^{10}\,M_\odot$). We show that the discrete nature of individual SNe (as opposed to a model in which their energy/momentum deposition is continuous over time, similar to stellar winds) is crucial in generating a reasonable ISM structure and galactic winds and in regulating dwarf stellar masses. However, once SNe are discretized, accounting for the effects of IMF sampling on continuous mechanisms such as radiative feedback and stellar mass-loss (as opposed to adopting IMF-averaged rates) has weak effects on galaxy-scale properties. We also consider the effects of rare HNe events with energies $\sim 10^{53}\,{\rm erg}$. The effects of HNe are similar to the effects of clustered explosions of SNe -- which are already captured in our default simulation setup -- and do not quench star formation (provided that the HNe do not dominate the total SNe energy budget), which suggests that HNe yield products should be observable in ultra-faint dwarfs today. △ Less

Submitted 7 December, 2017; originally announced December 2017.

Comments: 9 pages, 4 figures

SDSS-IV MaNGA: Constraints on the Conditions for Star Formation in Galaxy Discs

Authors: David V. Stark, Kevin A. Bundy, Matthew E. Orr, Philip F. Hopkins, Kyle Westfall, Matthew Bershady, Cheng Li, Dmitry Bizyaev, Karen L. Masters, Anne-Marie Weijmans, Ivan Lacerna, Daniel Thomas, Niv Drory, Renbin Yan, Kai Zhang

Abstract: Regions of disc galaxies with widespread star formation tend to be both gravitationally unstable and self-shielded against ionizing radiation, whereas extended outer discs with little or no star formation tend to be stable and unshielded on average. We explore what drives the transition between these two regimes, specifically whether discs first meet the conditions for self-shielding (parameterize… ▽ More Regions of disc galaxies with widespread star formation tend to be both gravitationally unstable and self-shielded against ionizing radiation, whereas extended outer discs with little or no star formation tend to be stable and unshielded on average. We explore what drives the transition between these two regimes, specifically whether discs first meet the conditions for self-shielding (parameterized by dust optical depth, $τ$) or gravitational instability (parameterized by a modified version of Toomre's instability parameters, $Q_{\rm thermal}$, which quantifies the stability of a gas disc that is thermally supported at $T=10^4$ K). We first introduce a new metric formed by the product of these quantities, $Q_{\rm thermal}τ$, which indicates whether the conditions for disk instability or self-shielding are easier to meet in a given region of a galaxy, and we discuss how $Q_{\rm thermal}τ$ can be constrained even in the absence of direct gas information. We then analyse a sample of 13 galaxies with resolved gas measurements and find that on average galaxies will reach the threshold for disk instabilities ($Q_{\rm thermal}<1$) before reaching the threshold for self-shielding ($τ>1$). Using integral field spectroscopic observations of a sample of 236 galaxies from the MaNGA survey, we find that the value of $Q_{\rm thermal}τ$ in star-forming discs is consistent with similar behavior. These results support a scenario where disc fragmentation and collapse occurs before self-shielding, suggesting that gravitational instabilities are the primary condition for widespread star formation in galaxy discs. Our results support similar conclusions based on recent galaxy simulations. △ Less

Submitted 19 November, 2017; v1 submitted 31 October, 2017; originally announced November 2017.

Comments: 12 pages, 6 figures, accepted for publication in MNRAS

Stacked star formation rate profiles of bursty galaxies exhibit 'coherent' star formation

Authors: Matthew E. Orr, Christopher C. Hayward, Erica J. Nelson, Philip F. Hopkins, Claude-André Faucher-Giguère, Dušan Kereš, T. K. Chan, Denise M. Schmitz, Tim B. Miller

Abstract: In a recent work based on 3200 stacked H$α$ maps of galaxies at $z \sim 1$, Nelson et al.~find evidence for `coherent star formation': the stacked SFR profiles of galaxies above (below) the 'star formation main sequence' (MS) are above (below) that of galaxies on the MS at all radii. One might interpret this result as inconsistent with highly bursty star formation and evidence that galaxies evolve… ▽ More In a recent work based on 3200 stacked H$α$ maps of galaxies at $z \sim 1$, Nelson et al.~find evidence for `coherent star formation': the stacked SFR profiles of galaxies above (below) the 'star formation main sequence' (MS) are above (below) that of galaxies on the MS at all radii. One might interpret this result as inconsistent with highly bursty star formation and evidence that galaxies evolve smoothly along the MS rather than crossing it many times. We analyze six simulated galaxies at $z\sim1$ from the Feedback in Realistic Environments (FIRE) project in a manner analogous to the observations to test whether the above interpretations are correct. The trends in stacked SFR profiles are qualitatively consistent with those observed. However, SFR profiles of individual galaxies are much more complex than the stacked profiles: the former can be flat or even peak at large radii because of the highly clustered nature of star formation in the simulations. Moreover, the SFR profiles of individual galaxies above (below) the MS are not systematically above (below) those of MS galaxies at all radii. We conclude that the time-averaged coherent star formation evident stacks of observed galaxies is consistent with highly bursty, clumpy star formation of individual galaxies and is not evidence that galaxies evolve smoothly along the MS. △ Less

Submitted 28 September, 2017; originally announced September 2017.

Comments: 7 pages, 4 figures, accepted for publication in ApJL

What FIREs Up Star Formation: the Emergence of the Kennicutt-Schmidt Law from Feedback

Authors: Matthew E. Orr, Christopher C. Hayward, Philip Hopkins, T. K. Chan, Claude-André Faucher-Giguère, Robert Feldmann, Dušan Kereš, Norman Murray, Eliot Quataert

Abstract: We present an analysis of the global and spatially-resolved Kennicutt-Schmidt (KS) star formation relation in the FIRE (Feedback In Realistic Environments) suite of cosmological simulations, including halos with $z = 0$ masses ranging from $10^{10}$ -- $10^{13}$ M$_{\odot}$. We show that the KS relation emerges and is robustly maintained due to the effects of feedback on local scales regulating st… ▽ More We present an analysis of the global and spatially-resolved Kennicutt-Schmidt (KS) star formation relation in the FIRE (Feedback In Realistic Environments) suite of cosmological simulations, including halos with $z = 0$ masses ranging from $10^{10}$ -- $10^{13}$ M$_{\odot}$. We show that the KS relation emerges and is robustly maintained due to the effects of feedback on local scales regulating star-forming gas, independent of the particular small-scale star formation prescriptions employed. We demonstrate that the time-averaged KS relation is relatively independent of redshift and spatial averaging scale, and that the star formation rate surface density is weakly dependent on metallicity and inversely dependent on orbital dynamical time. At constant star formation rate surface density, the `Cold \& Dense' gas surface density (gas with $T < 300$~K and $n > 10$~cm$^{-3}$, used as a proxy for the molecular gas surface density) of the simulated galaxies is $\sim$0.5~dex less than observed at $\sim$kpc scales. This discrepancy may arise from underestimates of the local column density at the particle-scale for the purposes of shielding in the simulations. Finally, we show that on scales larger than individual giant molecular clouds, the primary condition that determines whether star formation occurs is whether a patch of the galactic disk is thermally Toomre-unstable (not whether it is self-shielding): once a patch can no longer be thermally stabilized against fragmentation, it collapses, becomes self-shielding, cools, and forms stars, regardless of epoch or environment. △ Less

Submitted 30 May, 2018; v1 submitted 6 January, 2017; originally announced January 2017.

Comments: 23 pages, 16 figures, accepted to MNRAS, fourth revision