Showing 1–50 of 258 results for author: Faucher-Giguère, C

Elevated UV luminosity density at Cosmic Dawn explained by non-evolving, weakly-mass dependent star formation efficiency

Authors: Robert Feldmann, Michael Boylan-Kolchin, James S. Bullock, Onur Çatmabacak, Claude-André Faucher-Giguère, Christopher C. Hayward, Dušan Kereš, Alexandres Lazar, Lichen Liang, Jorge Moreno, Pascal A. Oesch, Eliot Quataert, Xuejian Shen, Guochao Sun

Abstract: Recent observations with the James Webb Space Telescope (JWST) have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the Big Bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox-HR… ▽ More Recent observations with the James Webb Space Telescope (JWST) have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the Big Bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox-HR, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox-HR re-simulates the cosmic volume (L = 22.1 cMpc) of the original FIREbox run with eight times higher mass resolution (m_b ~ 7800 M_sun), but with identical physics, down to z ~ 6. FIREbox-HR predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at z ~ 6 - 14, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox-HR, the SFE - halo mass relation for intermediate mass halos (M_halo ~ 10^9 - 10^11 M_sun) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE - halo mass relation lead to a larger contribution from lower mass halos at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE - halo mass relation inferred from FIREbox-HR allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at z > 12 will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 26 pages, 14 figures, 5 tables, submitted to MNRAS, comments welcome

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.

Angular momentum transfer in cosmological simulations of Milky Way-mass discs

Authors: Cameron W. Trapp, Dušan Kereš, Philip F. Hopkins, Claude-André Faucher-Giguère

Abstract: Fueling star formation in large, disky galaxies, requires a continuous supply of gas accreting into star-forming regions. Previously, we characterized this accretion in 4 Milky Way mass galaxies ($M_{\rm halo}\sim10^{12}M_{\odot}$) in the FIRE-2 cosmological zoom-in simulations, focusing on runs with cosmic ray physics. At $z\sim0$, gas within the inner circumgalactic medium (iCGM) approaches the… ▽ More Fueling star formation in large, disky galaxies, requires a continuous supply of gas accreting into star-forming regions. Previously, we characterized this accretion in 4 Milky Way mass galaxies ($M_{\rm halo}\sim10^{12}M_{\odot}$) in the FIRE-2 cosmological zoom-in simulations, focusing on runs with cosmic ray physics. At $z\sim0$, gas within the inner circumgalactic medium (iCGM) approaches the disk with comparable angular momentum (AM) to the disk edge, joining in the outer half of the gaseous disk. Within the disk, gas moves inward at velocities of $\sim$1-5~km~s$^{-1}$ while fully rotationally supported. In this study, we analyze the torques that drive these flows. In all cases, we find that the torques in disks enable gas accreted near the disk edge to transport inwards and fuel star formation in the central few kpc. The primary sources of torque come from gravity, hydrodynamical forces, and the sub-grid $P dV$ work done by supernova (SNe) remnants interacting with gas on $\lesssim$10 pc scales. These SNe remnant interactions provide a major source of torque on the gas, inducing negative torques within the inner disk and positive torques in the outer disk. The gas-gas gravitational, hydro, and "feedback" torques all transfer AM outward to where accreting gas is joining the disk, playing an important role in driving inflows and regulating disk structure. Gravitational torques from stars and dark matter provide an AM sink within the innermost regions of the disk and iCGM, respectively. Torques from viscous shearing, magnetic forces, stellar winds, and radiative transfer are largely insignificant. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 19 pages, 17 figures

Confronting the Diversity Problem: The Limits of Galaxy Rotation Curves as a tool to Understand Dark Matter Profiles

Authors: Isabel S. Sands, Philip F. Hopkins, Xuejian Shen, Michael Boylan-Kolchin, James Bullock, Claude-Andre Faucher-Giguere, Francisco J. Mercado, Jorge Moreno, Lina Necib, Xiaowei Ou, Sarah Wellons, Andrew Wetzel

Abstract: While galaxy rotation curves provide one of the most powerful methods for measuring dark matter profiles in the inner regions of rotation-supported galaxies, at the dwarf scale there are factors that can complicate this analysis. Given the expectation of a universal profile in dark matter-only simulations, the diversity of observed rotation curves has become an often-discussed issue in Lambda Cold… ▽ More While galaxy rotation curves provide one of the most powerful methods for measuring dark matter profiles in the inner regions of rotation-supported galaxies, at the dwarf scale there are factors that can complicate this analysis. Given the expectation of a universal profile in dark matter-only simulations, the diversity of observed rotation curves has become an often-discussed issue in Lambda Cold Dark Matter cosmology on galactic scales. We analyze a suite of Feedback in Realistic Environments (FIRE) simulations of $10^{10}-10^{12}$ $M_\odot$ halos with standard cold dark matter, and compare the true circular velocity to rotation curve reconstructions. We find that, for galaxies with well-ordered gaseous disks, the measured rotation curve may deviate from true circular velocity by at most 10% within the radius of the disk. However, non-equilibrium behavior, non-circular motions, and non-thermal and non-kinetic stresses may cause much larger discrepancies of 50% or more. Most rotation curve reconstructions underestimate the true circular velocity, while some reconstructions transiently over-estimate it in the central few kiloparsecs due to dynamical phenomena. We further demonstrate that the features that contribute to these failures are not always visibly obvious in HI observations. If such dwarf galaxies are included in galaxy catalogs, they may give rise to the appearance of "artificial" rotation curve diversity that does not reflect the true variation in underlying dark matter profiles. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 20 pages, 10 figures

Predictions for CO emission and the CO-to-H$_2$ conversion factor in galaxy simulations with non-equilibrium chemistry

Authors: Oliver A. Thompson, Alexander J. Richings, Brad K. Gibson, Claude-André Faucher-Giguère, Robert Feldmann, Christopher C. Hayward

Abstract: Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H$_2$ gas mass of the region using the CO-to-H$_2$ conversion factor (Xco). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulat… ▽ More Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H$_2$ gas mass of the region using the CO-to-H$_2$ conversion factor (Xco). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulations of isolated disc galaxies, ranging from dwarfs to Milky Way-mass galaxies, that were run using the FIRE-2 subgrid models coupled to the CHIMES non-equilibrium chemistry solver. We use the non-equilibrium abundances from the simulations, and we also compare to results using abundances assuming equilibrium, which we calculate from the simulation in post-processing. Our non-equilibrium simulations are able to reproduce the relation between CO and H$_2$ column densities, and the relation between Xco and metallicity, seen within observations of the Milky Way. We also compare to the xCOLD GASS survey, and find agreement with their data to our predicted CO luminosities at fixed star formation rate. We also find the multivariate function used by xCOLD GASS overpredicts the H$_2$ mass for our simulations, motivating us to suggest an alternative multivariate function of our fitting, though we caution that this fitting is uncertain due to the limited range of galaxy conditions covered by our simulations. We also find that the non-equilibrium chemistry has little effect on the conversion factor (<5\%) for our high-mass galaxies, though still affects the H$_2$ mass and Lco by $\approx$25\%. △ Less

Submitted 12 June, 2024; v1 submitted 21 March, 2024; originally announced March 2024.

Comments: 19 pages, 16 figures (including appendices). Accepted for publication in MNRAS

The High-Redshift Gas-Phase Mass-Metallicity Relation in FIRE-2

Authors: Andrew Marszewski, Guochao Sun, Claude-André Faucher-Giguère, Christopher C. Hayward, Robert Feldmann

Abstract: The unprecedented infrared spectroscopic capabilities of JWST have provided high-quality interstellar medium (ISM) metallicity measurements and enabled characterization of the gas-phase mass-metallicity relation (MZR) for galaxies at $z \gtrsim 5$ for the first time. We analyze the gas-phase MZR and its evolution in a high-redshift suite of FIRE-2 cosmological zoom-in simulations at $z=5-12$ and f… ▽ More The unprecedented infrared spectroscopic capabilities of JWST have provided high-quality interstellar medium (ISM) metallicity measurements and enabled characterization of the gas-phase mass-metallicity relation (MZR) for galaxies at $z \gtrsim 5$ for the first time. We analyze the gas-phase MZR and its evolution in a high-redshift suite of FIRE-2 cosmological zoom-in simulations at $z=5-12$ and for stellar masses $M_* \sim 10^6-10^{10} \rm{M}_\odot$. These simulations implement a multi-channel stellar feedback model and produce broadly realistic galaxy properties, including when evolved to $z=0$. The simulations predict very weak redshift evolution of the MZR over the redshift range studied, with the normalization of the MZR increasing by less than $0.01$ dex as redshift decreases from $z = 12$ to $z=5$. The median MZR in the simulations is well-approximated as a constant power-law relation across this redshift range given by $\log(Z/Z_\odot) = 0.37\log(M_*/\rm{M}_\odot) - 4.3$. We find good agreement between our best-fit model and recent observations made by JWST at high redshift. The weak evolution of the MZR at $z > 5$ contrasts with the evolution at $z \lesssim 3$, where increasing normalization of the MZR with decreasing redshift is observed and predicted by most models. The FIRE-2 simulations predict increasing scatter in the gas-phase MZR with decreasing stellar mass, in qualitative agreement with some observations. △ Less

Submitted 23 May, 2024; v1 submitted 13 March, 2024; originally announced March 2024.

Comments: 13 pages, 6 figures

The Cosmic Ultraviolet Baryon Survey (CUBS) VII: on the warm-hot circumgalactic medium probed by O VI and Ne VIII at 0.4 $\lesssim$ z $\lesssim$ 0.7

Authors: Zhijie Qu, Hsiao-Wen Chen, Sean D. Johnson, Gwen C. Rudie, Fakhri S. Zahedy, David DePalma, Joop Schaye, Erin T. Boettcher, Sebastiano Cantalupo, Mandy C. Chen, Claude-André Faucher-Giguère, Jennifer I-Hsiu Li, John S. Mulchaey, Patrick Petitjean, Marc Rafelski

Abstract: This paper presents a newly established sample of 103 unique galaxies or galaxy groups at $0.4\lesssim z\lesssim 0.7$ from the Cosmic Ultraviolet Baryon Survey (CUBS) for studying the warm-hot circumgalactic medium (CGM) probed by both O VI and Ne VIII absorption. The galaxies and associated neighbors are identified at $< 1$ physical Mpc from the sightlines toward 15 CUBS QSOs at… ▽ More This paper presents a newly established sample of 103 unique galaxies or galaxy groups at $0.4\lesssim z\lesssim 0.7$ from the Cosmic Ultraviolet Baryon Survey (CUBS) for studying the warm-hot circumgalactic medium (CGM) probed by both O VI and Ne VIII absorption. The galaxies and associated neighbors are identified at $< 1$ physical Mpc from the sightlines toward 15 CUBS QSOs at $z_{\rm QSO}\gtrsim 0.8$. A total of 30 galaxies or galaxy groups exhibit associated O VI $λλ$ 1031, 1037 doublet absorption within a line-of-sight velocity interval of $\pm250$ km/s, while the rest show no trace of O VI to a detection limit of $\log N_{\rm OVI}/{\rm cm^{-2}}\approx13.7$. Meanwhile, only five galaxies or galaxy groups exhibit the Ne VIII $λλ$ 770,780 doublet absorption, down to a limiting column density of $\log N_{\rm NeVIII}/{\rm cm^{-2}}\approx14.0$. These O VI- and Ne VIII-bearing halos reside in different galaxy environments with stellar masses ranging from $\log M_{\rm star}/M_\odot \approx 8$ to $\approx11.5$. The warm-hot CGM around galaxies of different stellar masses and star formation rates exhibits different spatial profiles and kinematics. In particular, star-forming galaxies with $\log M_{\rm star}/M_\odot\approx9-11$ show a significant concentration of metal-enriched warm-hot CGM within the virial radius, while massive quiescent galaxies exhibit flatter radial profiles of both column densities and covering fractions. In addition, the velocity dispersion of O VI absorption is broad with $σ_v > 40$ km/s for galaxies of $\log M_{\rm star}/M_\odot>9$ within the virial radius, suggesting a more dynamic warm-hot halo around these galaxies. Finally, the warm-hot CGM probed by O VI and Ne VIII is suggested to be the dominant phase in sub-$L^*$ galaxies with $\log M_{\rm star}/M_\odot\approx9-10$ based on their high ionization fractions in the CGM. △ Less

Submitted 14 February, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: Submitted to ApJ after addressing the referee's comments; 28 pages, 16 figures

Inflow and outflow properties, not total gas fractions, drive the evolution of the mass-metallicity relation

Authors: Luigi Bassini, Robert Feldmann, Jindra Gensior, Claude-André Faucher-Giguère, Elia Cenci, Jorge Moreno, Mauro Bernardini, Lichen Liang

Abstract: Observations show a tight correlation between the stellar mass of galaxies and their gas-phase metallicity (MZR). This relation evolves with redshift, with higher-redshift galaxies being characterized by lower metallicities. Understanding the physical origin of the slope and redshift evolution of the MZR may provide important insight into the physical processes underpinning it: star formation, fee… ▽ More Observations show a tight correlation between the stellar mass of galaxies and their gas-phase metallicity (MZR). This relation evolves with redshift, with higher-redshift galaxies being characterized by lower metallicities. Understanding the physical origin of the slope and redshift evolution of the MZR may provide important insight into the physical processes underpinning it: star formation, feedback, and cosmological inflows. While theoretical models ascribe the shape of the MZR to the lower efficiency of galactic outflows in more massive galaxies, what drives its evolution remains an open question. In this letter, we analyze how the MZR evolves over $z=0-3$, combining results from the FIREbox cosmological volume simulation with analytical models. Contrary to a frequent assertion in the literature, we find that the evolution of the gas fraction does not contribute significantly to the redshift evolution of the MZR. Instead, we show that the latter is driven by the redshift-dependence of the inflow metallicity, outflow metallicity, and mass loading factor, whose relative importance depends on stellar mass. These findings also suggest that the evolution of the MZR is not explained by galaxies moving along a fixed surface in the space spanned by stellar mass, gas phase metallicity, and star formation rate. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: 6 pages, 4 figures. Submitted for publication on MNRAS

Ne VIII in the warm-hot circumgalactic medium of FIRE simulations and in observations

Authors: Nastasha A. Wijers, Claude-André Faucher-Giguère, Jonathan Stern, Lindsey Byrne, Imran Sultan

Abstract: The properties of warm-hot gas around $\sim L_{*}$ galaxies can be studied with absorption lines from highly ionized metals. We predict Ne VIII column densities from cosmological zoom-in simulations of halos with masses in $\sim 10^{12}$ and $\sim 10^{13}\,\mathrm{M}_{\odot}$ from the FIRE project. Ne VIII traces the volume-filling, virial-temperature gas in $\sim 10^{12}\,\mathrm{M}_{\odot}$ halo… ▽ More The properties of warm-hot gas around $\sim L_{*}$ galaxies can be studied with absorption lines from highly ionized metals. We predict Ne VIII column densities from cosmological zoom-in simulations of halos with masses in $\sim 10^{12}$ and $\sim 10^{13}\,\mathrm{M}_{\odot}$ from the FIRE project. Ne VIII traces the volume-filling, virial-temperature gas in $\sim 10^{12}\,\mathrm{M}_{\odot}$ halos. In $\sim 10^{13}\,\mathrm{M}_{\odot}$ halos the Ne VIII gas is clumpier, and biased towards the cooler part of the warm-hot phase. We compare the simulations to observations by the CASBaH and CUBS surveys. We show that when inferring halo masses from stellar masses to compare simulated and observed halos, it is important to account for the scatter in the stellar-mass-halo-mass relation, especially at $\,\mathrm{M}_{\star} \gtrsim 10^{10.5} \,\mathrm{M}_{\odot}$. Median Ne VIII columns in the fiducial FIRE-2 model are about as high as observed upper limits allow, while the simulations analyzed do not reproduce the highest observed columns. This suggests that the median Ne VIII profiles predicted by the simulations are consistent with observations, but that the simulations may underpredict the scatter. We find similar agreement with analytical models that assume a product of the halo gas fraction and metallicity (relative to solar) $\sim 0.1$-$0.3$, indicating that observations are consistent with plausible CGM temperatures, metallicities, and gas masses. Variants of the FIRE simulations with a modified supernova feedback model and/or AGN feedback included (as well as some other cosmological simulations from the literature) more systematically underpredict Ne VIII columns. The circumgalactic Ne VIII observations therefore provide valuable constraints on simulations that otherwise predict realistic galaxy properties. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 19 pages + references and appendices, submitted to ApJ

A Dusty Locale: Evolution of Galactic Dust Populations from Milky Way to Dwarf-Mass Galaxies

Authors: Caleb R. Choban, Dušan Kereš, Karin M. Sandstrom, Philip F. Hopkins, Christopher C. Hayward, Claude-André Faucher-Giguère

Abstract: Observations indicate dust populations vary between galaxies and within them, suggesting a complex life cycle and evolutionary history. Here we investigate the evolution of galactic dust populations across cosmic time using a suite of cosmological zoom-in simulations from the Feedback in Realistic Environments (FIRE) project, spanning $M_{\rm vir}=10^{9-12}M_{\odot};\,M_{*}=10^{6-11}\,M_{\odot}$.… ▽ More Observations indicate dust populations vary between galaxies and within them, suggesting a complex life cycle and evolutionary history. Here we investigate the evolution of galactic dust populations across cosmic time using a suite of cosmological zoom-in simulations from the Feedback in Realistic Environments (FIRE) project, spanning $M_{\rm vir}=10^{9-12}M_{\odot};\,M_{*}=10^{6-11}\,M_{\odot}$. Our simulations incorporate a dust evolution model that accounts for the dominant sources of dust production, growth, and destruction and follows the evolution of specific dust species. All galactic dust populations in our suite exhibit similar evolutionary histories, with gas-dust accretion being the dominant producer of dust mass for all but the most metal-poor galaxies. Similar to previous works, we find the onset of efficient gas-dust accretion occurs above a `critical' metallicity threshold ($Z_{\rm crit}$). Due to this threshold, our simulations reproduce observed trends between galactic D/Z and metallicity and element depletion trends in the ISM. However, we find $Z_{\rm crit}$ varies between dust species due to differences in key element abundances, dust physical properties, and life cycle processes resulting in $Z_{\rm crit}\sim0.05Z_{\odot},\,0.2Z_{\odot},\,0.5Z_{\odot}$ for metallic iron, silicates, and carbonaceous dust, respectively. These variations could explain the lack of small carbonaceous grains observed in the Magellanic Clouds. We also find a delay between the onset of gas-dust accretion and when a dust population reaches equilibrium, which we call the equilibrium timescale ($τ_{\rm eq}$). The relation between $τ_{\rm eq}$ and the metal enrichment timescale of a galaxy, determined by its recent evolutionary history, can contribute to the scatter in the observed relation between galactic D/Z and metallicity. △ Less

Submitted 19 March, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: 24 pages, 18 figures, accepted for publication in MNRAS. Minor changes

The HI covering fraction of Lyman Limit Systems in FIRE haloes

Authors: Lucas Tortora, Robert Feldmann, Mauro Bernardini, Claude-André Faucher-Giguère

Abstract: Atomic hydrogen (HI) serves a crucial role in connecting galactic-scale properties such as star formation with the large-scale structure of the Universe. While recent numerical simulations have successfully matched the observed covering fraction of HI near Lyman Break Galaxies (LBGs) and in the foreground of luminous quasars at redshifts $z \lesssim 3$, the low-mass end remains as-of-yet unexplore… ▽ More Atomic hydrogen (HI) serves a crucial role in connecting galactic-scale properties such as star formation with the large-scale structure of the Universe. While recent numerical simulations have successfully matched the observed covering fraction of HI near Lyman Break Galaxies (LBGs) and in the foreground of luminous quasars at redshifts $z \lesssim 3$, the low-mass end remains as-of-yet unexplored in observational and computational surveys. We employ a cosmological, hydrodynamical simulation (FIREbox) supplemented with zoom-in simulations (MassiveFIRE) from the Feedback In Realistic Environments (FIRE) project to investigate the HI covering fraction of Lyman Limit Systems ($N_{\mathrm{HI}} \gtrsim 10^{17.2}$ cm$^{-2}$) across a wide range of redshifts ($z=0-6$) and halo masses ($10^8-10^{13} M_{\odot}$ at $z=0$, $10^8-10^{11} M_{\odot}$ at $z=6$) in the absence of feedback from active galactic nuclei (AGN). We find that the covering fraction inside haloes exhibits a strong increase with redshift, with only a weak dependence on halo mass for higher-mass haloes. For massive haloes ($M_{\mathrm{vir}} \sim 10^{11}-10^{12} M_{\odot}$), the radial profiles showcase scale-invariance and remain independent of mass. The radial dependence is well-captured by a fitting function. The covering fractions in our simulations are in good agreement with measurements of the covering fraction in LBGs. Our comprehensive analysis unveils a complex dependence with redshift and halo mass for haloes with $M_{\mathrm{vir}} \lesssim 10^{10} M_{\odot}$ that future observations aim to constrain, providing key insights into the physics of structure formation and gas assembly. △ Less

Submitted 29 November, 2023; originally announced November 2023.

TESS photometry of the nova eruption in V606 Vul: asymmetric photosphere and multiple ejections?

Authors: Kirill V. Sokolovsky, Elias Aydi, Konstantin Malanchev, Colin J. Burke, Koji Mukai, J. L. Sokoloski, Brian D. Metzger, Kirill E. Atapin, Aleksandre A. Belinski, Yu-Ching Chen, Laura Chomiuk, Pavol A. Dubovsky, Claude-Andre Faucher-Giguere, Rebekah A. Hounsell, Natalia P. Ikonnikova, Vsevolod Yu. Lander, Junyao Li, Justin D. Linford, Amy J. Mioduszewski, Isabella Molina, Ulisse Munari, Sergey A. Potanin, Robert M. Quimby, Michael P. Rupen, Simone Scaringi , et al. (46 additional authors not shown)

Abstract: Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The ligh… ▽ More Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The lightcurve covers the first of two major peaks of V606 Vul that was reached 19 days after the start of the eruption. The nova reached its brightest visual magnitude V=9.9 in its second peak 64 days after the eruption onset, following the completion of Sector 41 observations. To increase the confidence level of the extracted lightcurve, we performed the analysis using four different codes implementing the aperture photometry (Lightkurve, VaST) and image subtraction (TESSreduce, tequila_shots) and find good agreement between them. We performed ground-based photometric and spectroscopic monitoring to complement the TESS data. The TESS lightcurve reveals two features: periodic variations (0.12771d, 0.01mag average peak-to-peak amplitude) that disappeared when the source was within 1mag of peak optical brightness and a series of isolated mini-flares (with peak-to-peak amplitudes of up to 0.5mag) appearing at seemingly random times. We interpret the periodic variations as the result of azimuthal asymmetry of the photosphere engulfing the nova-hosting binary that was distorted by and rotating with the binary. Whereas we use spectra to associate the two major peaks in the nova lightcurve with distinct episodes of mass ejection, the origin of mini-flares remains elusive. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 31 pages, 10 figures, submitted to ApJ; comments welcome

Dense stellar clump formation driven by strong quasar winds in the FIRE cosmological hydrodynamic simulations

Authors: Jonathan Mercedes-Feliz, Daniel Anglés-Alcázar, Boon Kiat Oh, Christopher C. Hayward, Rachel K. Cochrane, Alexander J. Richings, Claude-André Faucher-Giguère, Sarah Wellons, Bryan A. Terrazas, Jorge Moreno, Kung Yi Su, Philip F. Hopkins

Abstract: We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star forming galaxy at $z \sim 2$ under the presence of strong quasar winds. Our simulations include multi-phase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disk winds. We show that powerfu… ▽ More We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star forming galaxy at $z \sim 2$ under the presence of strong quasar winds. Our simulations include multi-phase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disk winds. We show that powerful quasar winds can have a global negative impact on galaxy growth while in the strongest cases triggering the formation of an off-center clump with stellar mass ${\rm M}_{\star}\sim 10^{7}\,{\rm M}_{\odot}$, effective radius ${\rm R}_{\rm 1/2\,\rm Clump}\sim 20\,{\rm pc}$, and surface density $Σ_{\star} \sim 10^{4}\,{\rm M}_{\odot}\,{\rm pc}^{-2}$. The clump progenitor gas cloud is originally not star-forming, but strong ram pressure gradients driven by the quasar winds (orders of magnitude stronger than experienced in the absence of winds) lead to rapid compression and subsequent conversion of gas into stars at densities much higher than the average density of star-forming gas. The AGN-triggered star-forming clump reaches ${\rm SFR} \sim 50\,{\rm M}_{\odot}\,{\rm yr}^{-1}$ and $Σ_{\rm SFR} \sim 10^{4}\,{\rm M}_{\odot}\,{\rm yr}^{-1}\,{\rm kpc}^{-2}$, converting most of the progenitor gas cloud into stars in $\sim$2\,Myr, significantly faster than its initial free-fall time and with stellar feedback unable to stop star formation. In contrast, the same gas cloud in the absence of quasar winds forms stars over a much longer period of time ($\sim$35\,Myr), at lower densities, and losing spatial coherency. The presence of young, ultra-dense, gravitationally bound stellar clumps in recently quenched galaxies could thus indicate local positive feedback acting alongside the strong negative impact of powerful quasar winds, providing a plausible formation scenario for globular clusters. △ Less

Submitted 17 April, 2024; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: 14 pages, 12 figures. Accepted for publication in MNRAS; doi:10.1093/mnras/stae1021

Unraveling Jet Quenching Criteria Across L* Galaxies and Massive Cluster Ellipticals

Authors: Kung-Yi Su, Greg L. Bryan, Christopher C. Hayward, Rachel S. Somerville, Philip F. Hopkins, Razieh Emami, Claude-André Faucher-Giguère, Eliot Quataert, Sam B. Ponnada, Drummond Fielding, Dušan Kereš

Abstract: In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. AGN jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work (Su et al. 2021) concludes from simulations of a halo with $10^{14} M_\odot$ that a successful… ▽ More In the absence of supplementary heat, the radiative cooling of halo gas around massive galaxies (Milky Way mass and above) leads to an excess of cold gas or stars beyond observed levels. AGN jet-induced heating is likely essential, but the specific properties of the jets remain unclear. Our previous work (Su et al. 2021) concludes from simulations of a halo with $10^{14} M_\odot$ that a successful jet model should have an energy flux comparable to the free-fall energy flux at the cooling radius and should inflate a sufficiently wide cocoon with a long enough cooling time. In this paper, we investigate three jet modes with constant fluxes satisfying the criteria, including high-temperature thermal jets, cosmic ray (CR)-dominant jets, and widely precessing kinetic jets in $10^{12}-10^{15}\,{\rm M}_{\odot}$ halos using high-resolution, non-cosmological MHD simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We find that scaling the jet energy according to the free-fall energy at the cooling radius can successfully suppress the cooling flows and quench galaxies without obviously violating observational constraints. We investigate an alternative scaling method in which we adjust the energy flux based on the total cooling rate within the cooling radius. However, we observe that the strong interstellar medium (ISM) cooling dominates the total cooling rate in this scaling approach, resulting in a jet flux that exceeds the amount needed to suppress the cooling flows. With the same energy flux, the CR-dominant jet is most effective in suppressing the cooling flow across all the surveyed halo masses due to the enhanced CR pressure support. We confirm that the criteria for a successful jet model, which we proposed in Su et al. (2021), work across a much wider range, encompassing halo masses of $10^{12}-10^{15} {\rm M_\odot}$. △ Less

Submitted 2 November, 2023; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: 15 pages, 10 figures, Submitted to MNRAS; comments welcome!

Effects of multi-channel AGN feedback in FIRE cosmological simulations of massive galaxies

Authors: Lindsey Byrne, Claude-André Faucher-Giguère, Sarah Wellons, Philip F. Hopkins, Daniel Anglés-Alcázar, Imran Sultan, Nastasha Wijers, Jorge Moreno, Sam Ponnada

Abstract: Feedback from supermassive black holes is believed to be a critical driver of the observed color bimodality of galaxies above the Milky Way mass scale. AGN feedback has been modeled in many galaxy formation simulations, but most implementations have involved simplified prescriptions or a coarse-grained interstellar medium (ISM). We present the first set of FIRE-3 cosmological zoom-in simulations w… ▽ More Feedback from supermassive black holes is believed to be a critical driver of the observed color bimodality of galaxies above the Milky Way mass scale. AGN feedback has been modeled in many galaxy formation simulations, but most implementations have involved simplified prescriptions or a coarse-grained interstellar medium (ISM). We present the first set of FIRE-3 cosmological zoom-in simulations with AGN feedback evolved to $z\sim0$, examining the impact of AGN feedback on a set of galaxies with halos in the mass range $10^{12}-10^{13} M_{\odot}$. These simulations combine detailed stellar and ISM physics with multi-channel AGN feedback including radiative feedback, mechanical outflows, and in some simulations, cosmic rays (CRs). We find that massive (>L*) galaxies in these simulations can match local scaling relations including the stellar mass-halo mass relation and the $M_{\rm BH}$-$σ$ relation; in the stronger model with CRs, they also match the size-mass relation and the Faber-Jackson relation. Many of the massive galaxies in the simulations with AGN feedback have quenched star formation and elliptical morphologies, in qualitative agreement with observations. In contrast, simulations at the massive end without AGN feedback produce galaxies that are too massive and form stars too rapidly, are order-of-magnitude too compact, and have velocity dispersions well above Faber-Jackson. Despite these successes, the AGN models analyzed do not produce uniformly realistic galaxies when the feedback parameters are held constant: while the stronger model produces the most realistic massive galaxies, it tends to over-quench the lower-mass galaxies. This indicates that further refinements of the AGN modeling are needed. △ Less

Submitted 7 May, 2024; v1 submitted 24 October, 2023; originally announced October 2023.

Comments: 18 pages, 8 figures, submitted to ApJ

An Analytic Model For Magnetically-Dominated Accretion Disks

Authors: Philip F. Hopkins, Jonathan Squire, Eliot Quataert, Norman Murray, Kung-Yi Su, Ulrich P. Steinwandel, Kyle Kremer, Claude-Andre Faucher-Giguere, Sarah Wellons

Abstract: Recent numerical cosmological radiation-magnetohydrodynamic-thermochemical-star formation simulations have resolved the formation of quasar accretion disks with Eddington or super-Eddington accretion rates onto supermassive black holes (SMBHs) down to a few hundred gravitational radii. These 'flux-frozen' and hyper-magnetized disks appear to be qualitatively distinct from classical $α$ disks and m… ▽ More Recent numerical cosmological radiation-magnetohydrodynamic-thermochemical-star formation simulations have resolved the formation of quasar accretion disks with Eddington or super-Eddington accretion rates onto supermassive black holes (SMBHs) down to a few hundred gravitational radii. These 'flux-frozen' and hyper-magnetized disks appear to be qualitatively distinct from classical $α$ disks and magnetically-arrested disks: the midplane pressure is dominated by toroidal magnetic fields with plasma $β\ll 1$ powered by advection of magnetic flux from the interstellar medium (ISM), and they are super-sonically and trans-Alfvenically turbulent with cooling times short compared to dynamical times yet remain gravitationally stable owing to magnetic support. In this paper, we present a simple analytic similarity model for such disks. For reasonable assumptions, the model is entirely specified by the boundary conditions (inflow rate at the BH radius of influence [BHROI]). We show that the scalings from this model are robust to various detailed assumptions, agree remarkably well with the simulations (given their simplicity), and demonstrate the self-consistency and gravitational stability of such disks even in the outer accretion disk (approaching the BHROI) at hyper-Eddington accretion rates. △ Less

Submitted 12 March, 2024; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: 12 pages, 1 figure. Replaced with accepted version. Companion paper to 'FORGE'D IN FIRE II' (arXiv:2310.04506, part of a series with arXiv:2309.13115 -- animations of the simulations referred to here can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html)

FORGE'd in FIRE II: The Formation of Magnetically-Dominated Quasar Accretion Disks from Cosmological Initial Conditions

Authors: Philip F. Hopkins, Jonathan Squire, Kung-Yi Su, Ulrich P. Steinwandel, Kyle Kremer, Yanlong Shi, Michael Y. Grudic, Sarah Wellons, Claude-Andre Faucher-Giguere, Daniel Angles-Alcazar, Norman Murray, Eliot Quataert

Abstract: In a companion paper, we reported the self-consistent formation of quasar accretion disks with inflow rates $\sim 10\,{\rm M_{\odot}\,yr^{-1}}$ down to <300 Schwarzschild radii from cosmological radiation-magneto-thermochemical-hydrodynamical galaxy and star formation simulations. We see the formation of a well-defined, steady-state accretion disk which is stable against star formation at sub-pc s… ▽ More In a companion paper, we reported the self-consistent formation of quasar accretion disks with inflow rates $\sim 10\,{\rm M_{\odot}\,yr^{-1}}$ down to <300 Schwarzschild radii from cosmological radiation-magneto-thermochemical-hydrodynamical galaxy and star formation simulations. We see the formation of a well-defined, steady-state accretion disk which is stable against star formation at sub-pc scales. The disks are optically thick, with radiative cooling balancing accretion, but with properties that are distinct from those assumed in most previous accretion disk models. The pressure is strongly dominated by (primarily toroidal) magnetic fields, with a plasma $β\sim 10^{-4}$ even in the disk midplane. They are qualitatively distinct from magnetically elevated or arrested disks. The disks are strongly turbulent, with trans-Alfvenic and highly super-sonic turbulence, and balance this via a cooling time that is short compared to the disk dynamical time, and can sustain highly super-Eddington accretion rates. Their surface and 3D densities at $\sim 10^{3}-10^{5}$ gravitational radii are much lower than in a Shakura-Sunyaev disk, with important implications for their thermo-chemistry and stability. We show how the magnetic field strengths and geometries arise from rapid advection of flux with the inflow from much weaker galaxy-scale fields in these 'flux-frozen' disks, and how this stabilizes the disk and gives rise to efficient torques. Re-simulating without magnetic fields produces catastrophic fragmentation with a vastly smaller, lower-$\dot{M}$ Shakura-Sunyaev-like disk. △ Less

Submitted 18 January, 2024; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: 46 pages, 30 figures. Updated to match accepted version from The Open Journal of Astrophysics. Part of a series with arXiv:2309.13115 . Animations of the simulations can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html

Seeding the CGM: How Satellites Populate the Cold Phase of Milky Way Halos

Authors: Manami Roy, Kung-Yi Su, Stephanie Tonnesen, Drummond B. Fielding, Claude-André Faucher-Giguère

Abstract: The origin of the cold phase in the CGM is a highly debated question. We investigate the contribution of satellite galaxies to the cold gas budget in the circumgalactic medium (CGM)of a Milky Way-like host galaxy. We perform controlled experiments with three different satellite mass distributions and identify several mechanisms by which satellites can add cold gas to the CGM, including ram pressur… ▽ More The origin of the cold phase in the CGM is a highly debated question. We investigate the contribution of satellite galaxies to the cold gas budget in the circumgalactic medium (CGM)of a Milky Way-like host galaxy. We perform controlled experiments with three different satellite mass distributions and identify several mechanisms by which satellites can add cold gas to the CGM, including ram pressure stripping and induced cooling in the mixing layer of the stripped cold gas. These two mechanisms contribute a comparable amount of cold gas to the host CGM. We find that the less massive satellites ($\leq 10^9 M_\odot$) not only lose all of their cold gas in a short period ($\sim$ 0.5-1 Gyr), but their stripped cold clouds also mix with the hot CGM gas and get heated up quickly. However, stellar feedback from these less massive satellites can hugely alter the fate of their stripped gas. Feedback speeds up the destruction of the stripped cold clouds from these satellites by making them more diffuse with more surface area. On the other hand, the more massive satellites (LMC or SMC-like $\sim 10^{10} M_\odot$) can add cold gas to the total gas budget of the host CGM for several Gyrs. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: 17 pages, 12 figures, 2 tables Accepted for publication in MNRAS on Oct 6, 2023

Synchrotron Signatures of Cosmic Ray Transport Physics in Galaxies

Authors: Sam B. Ponnada, Iryna S. Butsky, Raphael Skalidis, Philip F. Hopkins, Georgia V. Panopoulou, Cameron Hummels, Dušan Kereš, Eliot Quataert, Claude-André Faucher-Giguère, Kung-Yi Su

Abstract: Cosmic rays (CRs) may drive outflows and alter the phase structure of the circumgalactic medium, with potentially important implications on galaxy formation. However, these effects ultimately depend on the dominant mode of transport of CRs within and around galaxies, which remains highly uncertain. To explore potential observable constraints on CR transport, we investigate a set of cosmological FI… ▽ More Cosmic rays (CRs) may drive outflows and alter the phase structure of the circumgalactic medium, with potentially important implications on galaxy formation. However, these effects ultimately depend on the dominant mode of transport of CRs within and around galaxies, which remains highly uncertain. To explore potential observable constraints on CR transport, we investigate a set of cosmological FIRE-2 CR-MHD simulations of L$_{\ast}$ galaxies which evolve CRs with transport models motivated by self-confinement (SC) and extrinsic turbulence (ET) paradigms. To first order, the synchrotron properties diverge between SC and ET models due to a CR physics driven hysteresis. SC models show a higher tendency to undergo `ejective' feedback events due to a runaway buildup of CR pressure in dense gas due to the behavior of SC transport scalings at extremal CR energy densities. The corresponding CR wind-driven hysteresis results in brighter, smoother, and more extended synchrotron emission in SC runs relative to ET and constant diffusion runs. The differences in synchrotron arise from different morphology, ISM gas and \textbf{B} properties, potentially ruling out SC as the dominant mode of CR transport in typical star-forming L$_{\ast}$ galaxies, and indicating the potential for non-thermal radio continuum observations to constrain CR transport physics. △ Less

Submitted 28 February, 2024; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: 5 pages, 3 figures. Accepted version

FORGE'd in FIRE: Resolving the End of Star Formation and Structure of AGN Accretion Disks from Cosmological Initial Conditions

Authors: Philip F. Hopkins, Michael Y. Grudic, Kung-Yi Su, Sarah Wellons, Daniel Angles-Alcazar, Ulrich P. Steinwandel, David Guszejnov, Norman Murray, Claude-Andre Faucher-Giguere, Eliot Quataert, Dusan Keres

Abstract: It has recently become possible to zoom-in from cosmological to sub-pc scales in galaxy simulations to follow accretion onto supermassive black holes (SMBHs). However, at some point the approximations used on ISM scales (e.g. optically-thin cooling and stellar-population-integrated star formation [SF] and feedback [FB]) break down. We therefore present the first cosmological radiation-magnetohydro… ▽ More It has recently become possible to zoom-in from cosmological to sub-pc scales in galaxy simulations to follow accretion onto supermassive black holes (SMBHs). However, at some point the approximations used on ISM scales (e.g. optically-thin cooling and stellar-population-integrated star formation [SF] and feedback [FB]) break down. We therefore present the first cosmological radiation-magnetohydrodynamic (RMHD) simulation which self-consistently combines the FIRE physics (relevant on galactic/ISM scales where SF/FB are ensemble-averaged) and STARFORGE physics (relevant on small scales where we track individual (proto)stellar formation and evolution), together with explicit RMHD (including non-ideal MHD and multi-band M1-RHD) which self-consistently treats both optically-thick and thin regimes. This allows us to span scales from ~100 Mpc down to <100 au (~300 Schwarzschild radii) around a SMBH at a time where it accretes as a bright quasar, in a single simulation. We show that accretion rates up to $\sim 10-100\,{\rm M_{\odot}\,yr^{-1}}$ can be sustained into the accretion disk at $\ll 10^{3}\,R_{\rm schw}$, with gravitational torques between stars and gas dominating on sub-kpc scales until star formation is shut down on sub-pc scales by a combination of optical depth to cooling and strong magnetic fields. There is an intermediate-scale, flux-frozen disk which is gravitoturbulent and stabilized by magnetic pressure sustaining strong turbulence and inflow with persistent spiral modes. In this paper we focus on how gas gets into the small-scale disk, and how star formation is efficiently suppressed. △ Less

Submitted 12 March, 2024; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: 38 pages, 18 figures. Replaced with accepted version

Testing the near-far connection with FIRE simulations: inferring the stellar mass function of the proto-Local Group at z > 6 using the fossil record of present-day galaxies

Authors: Pratik J. Gandhi, Andrew Wetzel, Michael Boylan-Kolchin, Robyn E. Sanderson, Alessandro Savino, Daniel R. Weisz, Erik J. Tollerud, Guochao Sun, Claude-Andre Faucher-Giguere

Abstract: The shape of the low-mass (faint) end of the galaxy stellar mass function (SMF) or ultraviolet luminosity function (UVLF) at z > 6 is an open question for understanding which galaxies primarily drove cosmic reionisation. Resolved photometry of Local Group low-mass galaxies allows us to reconstruct their star formation histories, stellar masses, and UV luminosities at early times, and this fossil r… ▽ More The shape of the low-mass (faint) end of the galaxy stellar mass function (SMF) or ultraviolet luminosity function (UVLF) at z > 6 is an open question for understanding which galaxies primarily drove cosmic reionisation. Resolved photometry of Local Group low-mass galaxies allows us to reconstruct their star formation histories, stellar masses, and UV luminosities at early times, and this fossil record provides a powerful `near-far' technique for studying the reionisation-era SMF/UVLF, probing orders of magnitude lower in mass than direct HST/JWST observations. Using 882 low-mass (Mstar < 10^9 Msun) galaxies across 11 Milky Way- and Local Group-analogue environments from the FIRE-2 cosmological baryonic zoom-in simulations, we characterise their progenitors at z ~ 6 - 9, the mergers/disruption of those progenitors over time, and how well their present-day fossil record traces the high-redshift SMF. A present-day galaxy with Mstar ~ 10^5 Msun (10^9 Msun) had ~1 (~30) progenitors at z ~ 7, and its main progenitor comprised ~100% (~50%) of the total stellar mass of all its progenitors at z ~ 7. We show that although only ~ 15% of the early population of low-mass galaxies survives to present day, the fossil record of surviving Local Group galaxies accurately traces the low-mass slope of the SMF at z ~ 6 - 9. We find no obvious mass dependence to the mergers and accretion, and show that applying this reconstruction technique to just the low-mass galaxies at z = 0 and not the MW/M31 hosts correctly recovers the slope of the SMF down to Mstar ~ 10^4.5 Msun at z > 6. Thus, we validate the `near-far' approach as an unbiased tool for probing low-mass reionisation-era galaxies. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 16 pages, 6 figures, 3 appendix figures, and 1 table. Submitted to MNRAS

The Cosmic Ultraviolet Baryon Survey: Empirical Characterization of Turbulence in the Cool Circumgalactic Medium

Authors: Hsiao-Wen Chen, Zhijie Qu, Michael Rauch, Mandy C. Chen, Fakhri S. Zahedy, Sean D. Johnson, Joop Schaye, Gwen C. Rudie, Erin Boettcher, Sebastiano Cantalupo, Claude-Andre Faucher-Giguere, Jenny E. Greene, Sebastian Lopez, Robert A. Simcoe

Abstract: This paper reports the first measurement of the relationship between turbulent velocity and cloud size in the diffuse circumgalactic medium (CGM) in typical galaxy halos at redshift z~0.4-1. Through spectrally-resolved absorption profiles of a suite of ionic transitions paired with careful ionization analyses of individual components, cool clumps of size as small as l_cl~1 pc and density lower tha… ▽ More This paper reports the first measurement of the relationship between turbulent velocity and cloud size in the diffuse circumgalactic medium (CGM) in typical galaxy halos at redshift z~0.4-1. Through spectrally-resolved absorption profiles of a suite of ionic transitions paired with careful ionization analyses of individual components, cool clumps of size as small as l_cl~1 pc and density lower than nH = 0.001 cm^-3 are identified in galaxy halos. In addition, comparing the line widths between different elements for kinematically matched components provides robust empirical constraints on the thermal temperature T and the non-thermal motions bNT, independent of the ionization models. On average, bNT is found to increase with l_cl following bNT \propto l_cl^0.3 over three decades in spatial scale from l_cl~1 pc to l_cl~1 kpc. Attributing the observed bNT to turbulent motions internal to the clumps, the best-fit bNT-l_cl relation shows that the turbulence is consistent with Kolmogorov at <1 kpc with a roughly constant energy transfer rate per unit mass of epsilon~0.003 cm^2 s^-3 and a dissipation time scale of <~ 100 Myr. No significant difference is found between massive quiescent and star-forming halos in the sample on scales less than 1 kpc. While the inferred epsilon is comparable to what is found in CIV absorbers at high redshift, it is considerably smaller than observed in star-forming gas or in extended line-emitting nebulae around distant quasars. A brief discussion of possible sources to drive the observed turbulence in the cool CGM is presented. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: Accepted for publication in ApJ Letters; 10 pages, 4 figures

Synchrotron Emission on FIRE: Equipartition Estimators of Magnetic Fields in Simulated Galaxies with Spectrally-Resolved Cosmic Rays

Authors: Sam B. Ponnada, Georgia V. Panopoulou, Iryna S. Butsky, Philip F. Hopkins, Raphael Skalidis, Cameron Hummels, Eliot Quataert, Dušan Kereš, Claude-André Faucher-Giguère, Kung-Yi Su

Abstract: Synchrotron emission is one of few observable tracers of galactic magnetic fields (\textbf{B}) and cosmic rays (CRs). Much of our understanding of \textbf{B} in galaxies comes from utilizing synchrotron observations in conjunction with several simplifying assumptions of equipartition models, however it remains unclear how well these assumptions hold, and what \textbf{B} these estimates physically… ▽ More Synchrotron emission is one of few observable tracers of galactic magnetic fields (\textbf{B}) and cosmic rays (CRs). Much of our understanding of \textbf{B} in galaxies comes from utilizing synchrotron observations in conjunction with several simplifying assumptions of equipartition models, however it remains unclear how well these assumptions hold, and what \textbf{B} these estimates physically represent. Using FIRE simulations which self consistently evolve CR proton, electron, and positron spectra from MeV to TeV energies, we present the first synthetic synchrotron emission predictions from simulated L$_{*}$ galaxies with "live" spectrally-resolved CR-MHD. We find that synchrotron emission can be dominated by relatively cool and dense gas, resulting in equipartition estimates of \textbf{B} with fiducial assumptions underestimating the "true" \textbf{B} in the gas that contributes the most emission by factors of 2-3 due to small volume filling factors. Motivated by our results, we present an analytic framework that expands upon equipartition models for estimating \textbf{B} in a multi-phase medium. Comparing our spectrally-resolved synchrotron predictions to simpler spectral assumptions used in galaxy simulations with CRs, we find that spectral evolution can be crucial for accurate synchrotron calculations towards galactic centers, where loss terms are large. △ Less

Submitted 9 January, 2024; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 10 pages, 5 figures (with 3 additional figures in the Appendix). Accepted version (MNRAS December 21, 2023)

The proto-galaxy of Milky Way-mass haloes in the FIRE simulations

Authors: Danny Horta, Emily C. Cunningham, Robyn Sanderson, Kathryn V. Johnston, Alis Deason, Andrew Wetzel, Fiona McCluskey, Nicolás Garavito-Camargo, Lina Necib, Claude-André Faucher-Giguère, Arpit Arora, Pratik J. Gandhi

Abstract: Observational studies are finding stars believed to be relics of the earliest stages of hierarchical mass assembly of the Milky Way (i.e., proto-Galaxy). In this work, we contextualize these findings by studying the masses, ages, spatial distributions, morphology, kinematics, and chemical compositions of proto-galaxy populations from the 13 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological… ▽ More Observational studies are finding stars believed to be relics of the earliest stages of hierarchical mass assembly of the Milky Way (i.e., proto-Galaxy). In this work, we contextualize these findings by studying the masses, ages, spatial distributions, morphology, kinematics, and chemical compositions of proto-galaxy populations from the 13 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. Our findings indicate that proto-Milky Way populations: i) can have a stellar mass range between $1\times10^{8}<\mathrm{M}_{\star}<2\times10^{10}[\mathrm{M}_{\odot}]$, a virial mass range between $3\times10^{10}<\mathrm{M}_{\star}<6\times10^{11}[\mathrm{M}_{\odot}]$, and be as young as $8 \lesssim \mathrm{Age} \lesssim 12.8$ [Gyr] ($1\lesssim z \lesssim 6$); ii) are predominantly centrally concentrated, with $\sim50\%$ of the stars contained within $5-10$ kpc; iii) on average show weak but systematic net rotation in the plane of the host's disc at $z=0$ (i.e., $0.25\lesssim\langleκ/κ_{\mathrm{disc}}\rangle\lesssim0.8$); iv) present [$α$/Fe]-[Fe/H] compositions that overlap with the metal-poor tail of the host's old disc; v) tend to assemble slightly earlier in Local Group-like environments than in systems in isolation. Interestingly, we find that ~60% of the proto-Milky Way galaxies are comprised by 1 dominant system ($1/5\lesssim$M$_{\star}$/M$_{\star,\mathrm{proto-Milky Way}}$$\lesssim4/5$) and 4-5 lower mass systems (M$_{\star}$/M$_{\star,\mathrm{proto-Milky Way}}$$\lesssim1/10$); the other ~40% are comprised by 2 dominant systems and 3-4 lower mass systems. These massive/dominant proto-Milky Way fragments can be distinguished from the lower mass ones in chemical-kinematic samples, but appear (qualitatively) indistinguishable from one another. Our results could help observational studies disentangle if the Milky Way formed from one or two dominant systems. △ Less

Submitted 11 December, 2023; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: Accepted in MNRAS. 16 pages, 12 figures and 3 tables

Bursty Star Formation Naturally Explains the Abundance of Bright Galaxies at Cosmic Dawn

Authors: Guochao Sun, Claude-André Faucher-Giguère, Christopher C. Hayward, Xuejian Shen, Andrew Wetzel, Rachel K. Cochrane

Abstract: Recent discoveries of a significant population of bright galaxies at cosmic dawn $\left(z \gtrsim 10\right)$ have enabled critical tests of cosmological galaxy formation models. In particular, the bright end of the galaxy UV luminosity function (UVLF) appears higher than predicted by many models. Using approximately 25,000 galaxy snapshots at $8 \leq z \leq 12$ in a suite of FIRE-2 cosmological "z… ▽ More Recent discoveries of a significant population of bright galaxies at cosmic dawn $\left(z \gtrsim 10\right)$ have enabled critical tests of cosmological galaxy formation models. In particular, the bright end of the galaxy UV luminosity function (UVLF) appears higher than predicted by many models. Using approximately 25,000 galaxy snapshots at $8 \leq z \leq 12$ in a suite of FIRE-2 cosmological "zoom-in'' simulations from the Feedback in Realistic Environments (FIRE) project, we show that the observed abundance of UV-bright galaxies at cosmic dawn is reproduced in these simulations with a multi-channel implementation of standard stellar feedback processes, without any fine-tuning. Notably, we find no need to invoke previously suggested modifications such as a non-standard cosmology, a top-heavy stellar initial mass function, or a strongly enhanced star formation efficiency. We contrast the UVLFs predicted by bursty star formation in these original simulations to those derived from star formation histories (SFHs) smoothed over prescribed timescales (e.g., 100 Myr). The comparison demonstrates that the strongly time-variable SFHs predicted by the FIRE simulations play a key role in correctly reproducing the observed, bright-end UVLFs at cosmic dawn: the bursty SFHs induce order-or-magnitude changes in the abundance of UV-bright ($M_\mathrm{UV} \lesssim -20$) galaxies at $z \gtrsim 10$. The predicted bright-end UVLFs are consistent with both the spectroscopically confirmed population and the photometrically selected candidates. We also find good agreement between the predicted and observationally inferred integrated UV luminosity densities, which evolve more weakly with redshift in FIRE than suggested by some other models. △ Less

Submitted 10 September, 2023; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: 12 pages, 4 figures + 1 table, ApJL accepted

Hooks & Bends in the Radial Acceleration Relation: Discriminatory Tests for Dark Matter and MOND

Authors: Francisco J. Mercado, James S. Bullock, Jorge Moreno, Michael Boylan-Kolchin, Philip F. Hopkins, Andrew Wetzel, Claude-André Faucher-Giguère, Jenna Samuel

Abstract: The Radial Acceleration Relation (RAR) connects the total gravitational acceleration of a galaxy at a given radius, $a_{\rm tot}(r)$, with that accounted for by baryons at the same radius, $a_{\rm bar}(r)$. The shape and tightness of the RAR for rotationally-supported galaxies have characteristics in line with MOdified Newtonian Dynamics (MOND) and can also arise within the Cosmological Constant +… ▽ More The Radial Acceleration Relation (RAR) connects the total gravitational acceleration of a galaxy at a given radius, $a_{\rm tot}(r)$, with that accounted for by baryons at the same radius, $a_{\rm bar}(r)$. The shape and tightness of the RAR for rotationally-supported galaxies have characteristics in line with MOdified Newtonian Dynamics (MOND) and can also arise within the Cosmological Constant + Cold Dark Matter ($Λ$CDM) paradigm. We use zoom simulations of 20 galaxies with stellar masses of $M_{\star} \, \simeq 10^{7-11} \, M_{\odot}$ to study the RAR in the \texttt{FIRE-2} simulations. We highlight the existence of simulated galaxies with non-monotonic RAR tracks that ``hook'' down from the average relation. These hooks are challenging to explain in Modified Inertia theories of MOND, but naturally arise in all of our \lcdm-simulated galaxies that are dark-matter dominated at small radii and have feedback-induced cores in their dark matter haloes. We show, analytically and numerically, that downward hooks are expected in such cored haloes because they have non-monotonic acceleration profiles. We also extend the relation to accelerations below those traced by disc galaxy rotation curves. In this regime, our simulations exhibit ``bends'' off of the MOND-inspired extrapolation of the RAR, which, at large radii, approach $a_{\rm tot} \, \approx \, a_{\rm bar} \, /f_{\rm b}$, where $f_{\rm b}$ is the cosmic baryon fraction. Future efforts to search for these hooks and bends in real galaxies will provide interesting tests for MOND and $Λ$CDM. △ Less

Submitted 19 March, 2024; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: 15 pages, 9 figures; Accepted for publication in MNRAS

Does God play dice with star clusters?

Authors: Michael Y. Grudić, Stella S. R. Offner, Dávid Guszejnov, Claude-André Faucher-Giguère, Philip F. Hopkins

Abstract: When a detailed model of a stellar population is unavailable, it is most common to assume that stellar masses are independently and identically distributed according to some distribution: the universal initial mass function (IMF). However, stellar masses resulting from causal, long-ranged physics cannot be truly random and independent, and the IMF may vary with environment. To compare stochastic s… ▽ More When a detailed model of a stellar population is unavailable, it is most common to assume that stellar masses are independently and identically distributed according to some distribution: the universal initial mass function (IMF). However, stellar masses resulting from causal, long-ranged physics cannot be truly random and independent, and the IMF may vary with environment. To compare stochastic sampling with a physical model, we run a suite of 100 STARFORGE radiation magnetohydrodynamics simulations of low-mass star cluster formation in $2000M_\odot$ clouds that form $\sim 200$ stars each on average. The stacked IMF from the simulated clouds has a sharp truncation at $\sim 28 M_\odot$, well below the typically-assumed maximum stellar mass $M_{\rm up} \sim 100-150M_\odot$ and the total cluster mass. The sequence of star formation is not totally random: massive stars tend to start accreting sooner and finish later than the average star. However, final cluster properties such as maximum stellar mass and total luminosity have a similar amount of cloud-to-cloud scatter to random sampling. Therefore stochastic sampling does not generally model the stellar demographics of a star cluster as it is forming, but may describe the end result fairly well, if the correct IMF -- and its environment-dependent upper cutoff -- are known. △ Less

Submitted 7 December, 2023; v1 submitted 30 June, 2023; originally announced July 2023.

Comments: Accepted for publication in Open Journal of Astrophysics

The Cosmic Ultraviolet Baryon Survey (CUBS) VI: Connecting Physical Properties of the Cool Circumgalactic Medium to Galaxies at $z\approx 1$

Authors: Zhijie Qu, Hsiao-Wen Chen, Gwen C. Rudie, Sean D. Johnson, Fakhri S. Zahedy, David DePalma, Erin Boettcher, Sebastiano Cantalupo, Mandy C. Chen, Kathy L. Cooksey, Claude-André Faucher-Giguère, Jennifer I-Hsiu Li, Sebastian Lopez, Joop Schaye, Robert A. Simcoe

Abstract: This paper presents a new sample of 19 unique galaxies and galaxy groups at $z\approx1$ from the CUBS program, which is designated as the CUBSz1 sample. In this CUBSz1 sample, nine galaxies or galaxy groups show absorption features, while ten systems do not have detectable absorption with 2-$σ$ upper limits of log$N$(HeI)/cm$^{-2}\lesssim 13.5$ and log$N$(OV)/cm$^{-2}\lesssim 13.3$. Environmental… ▽ More This paper presents a new sample of 19 unique galaxies and galaxy groups at $z\approx1$ from the CUBS program, which is designated as the CUBSz1 sample. In this CUBSz1 sample, nine galaxies or galaxy groups show absorption features, while ten systems do not have detectable absorption with 2-$σ$ upper limits of log$N$(HeI)/cm$^{-2}\lesssim 13.5$ and log$N$(OV)/cm$^{-2}\lesssim 13.3$. Environmental properties of the galaxies, including galaxy overdensities, the total stellar mass and gravitational potential summed over all nearby neighbors, and the presence of local ionizing sources, are found to have a significant impact on the observed CGM absorption properties. Specifically, massive galaxies and galaxies in overdense regions exhibit a higher rate of incidence of absorption. At the same time, the observed CGM absorption properties in galaxy groups appear to be driven by the galaxy closest to the QSO sightline, rather than by the most massive galaxy or by mass-weighted properties. We introduce a total projected gravitational potential $ψ$, defined as $-ψ/G =\sum M_{\rm halo}/d_{\rm proj}$ summed over all group members, to characterize the overall galaxy environment. This projected gravitational potential correlates linearly with the maximum density detected in each sightline, consistent with higher-pressure gas being confined in deeper gravitational potential wells. In addition, we find that the radial profile of cool gas density exhibits a general decline from the inner regions to the outskirts, being in pressure balance with the hot halo. Finally, we note that the ionizing flux from nearby galaxies can generate an elevated $N$(HI)/$N$(HeI) ratio, which in turn provides a unique diagnostic of possible local sources contributing to the ionizing radiation field. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: 17 pages, 11 figures and 1 table in the main text; 40 figures and 1 table in the appendix; accepted for publication in MNRAS

Accretion onto disk galaxies via hot and rotating CGM inflows

Authors: Jonathan Stern, Drummond Fielding, Zachary Hafen, Kung-Yi Su, Nadav Naor, Claude-André Faucher-Giguère, Eliot Quataert, James Bullock

Abstract: Observed accretion rates onto the Milky-Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is an inflow in the volume-filling hot phase ($\sim10^6$ K) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution va… ▽ More Observed accretion rates onto the Milky-Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is an inflow in the volume-filling hot phase ($\sim10^6$ K) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution valid in the slow-rotation limit, we show that a hot inflow spins up as it approaches the galaxy, while remaining hot, subsonic and quasi-spherical. At the radius of angular momentum support ($\approx15$ kpc for the Milky-Way) the hot flow flattens into a disk geometry and then cools from $\sim10^6$ K to $\sim10^4$ K at the disk-halo interface. Cooling affects all hot gas, rather than just a subset of individual gas clouds, implying that accretion via hot inflows does not rely on local thermal instability in contrast with 'precipitation' models for galaxy accretion. Prior to cooling and accretion the inflow completes $\sim t_{\rm cool}/t_{\rm ff}$ radians of rotation, where $t_{\rm cool}/t_{\rm ff}$ is the cooling time to free-fall time ratio in hot gas immediately outside the galaxy. The ratio $t_{\rm cool}/t_{\rm ff}$ may thus govern the development of turbulence and enhancement of magnetic fields in gas accreting onto low-redshift spirals. We argue that accretion via hot inflows can explain the observed truncation of nearby thin stellar disks at $\approx4$ disk radii. We also show that if rotating hot inflows are common in Milky-Way size disk galaxies, as predicted, then signatures should be observable with X-ray telescopes and FRB surveys. △ Less

Submitted 15 February, 2024; v1 submitted 31 May, 2023; originally announced June 2023.

Comments: 21 pages, 12 figures, resubmitted to MNRAS following referee report

Probing bursty star formation by cross-correlating extragalactic background light and galaxy surveys

Authors: Guochao Sun, Adam Lidz, Andreas L. Faisst, Claude-André Faucher-Giguère

Abstract: Understanding the star formation rate (SFR) variability and how it depends on physical properties of galaxies is important for developing and testing the theory of galaxy formation. We investigate how statistical measurements of the extragalactic background light (EBL) can shed light on this topic and complement traditional methods based on observations of individual galaxies. Using semi-empirical… ▽ More Understanding the star formation rate (SFR) variability and how it depends on physical properties of galaxies is important for developing and testing the theory of galaxy formation. We investigate how statistical measurements of the extragalactic background light (EBL) can shed light on this topic and complement traditional methods based on observations of individual galaxies. Using semi-empirical models of galaxy evolution and SFR indicators sensitive to different star formation timescales (e.g., H$α$ and UV continuum luminosities), we show that the SFR variability, quantified by the joint probability distribution of the SFR indicators (i.e., the bivariate conditional luminosity function), can be characterized as a function of galaxy mass and redshift through the cross-correlation between deep, near-infrared maps of the EBL and galaxy distributions. As an example, we consider combining upcoming SPHEREx maps of the EBL with galaxy samples from Rubin/LSST. We demonstrate that their cross-correlation over a sky fraction of $f_\mathrm{sky}\sim0.5$ can constrain the joint SFR indicator distribution at high significance up to $z\sim2.5$ for mass-complete samples of galaxies down to $M_{*}\sim10^9\,M_{\odot}$. These constraints not only allow models of different SFR variability to be distinguished, but also provide unique opportunities to investigate physical mechanisms that require large number statistics such as environmental effects. The cross-correlations investigated illustrate the power of combining cosmological surveys to extract information inaccessible from each data set alone, while the large galaxy populations probed capture ensemble-averaged properties beyond the reach of targeted observations towards individual galaxies. △ Less

Submitted 30 June, 2023; v1 submitted 15 May, 2023; originally announced May 2023.

Comments: 12 pages, 7 figures, MNRAS accepted

Seen and unseen: bursty star formation and its implications for observations of high-redshift galaxies with JWST

Authors: Guochao Sun, Claude-André Faucher-Giguère, Christopher C. Hayward, Xuejian Shen

Abstract: Both observations and simulations have shown strong evidence for highly time-variable star formation in low-mass and/or high-redshift galaxies, which has important observational implications because high-redshift galaxy samples are rest-UV selected and therefore particularly sensitive to the recent star formation. Using a suite of cosmological "zoom-in" simulations at $z>5$ from the Feedback in Re… ▽ More Both observations and simulations have shown strong evidence for highly time-variable star formation in low-mass and/or high-redshift galaxies, which has important observational implications because high-redshift galaxy samples are rest-UV selected and therefore particularly sensitive to the recent star formation. Using a suite of cosmological "zoom-in" simulations at $z>5$ from the Feedback in Realistic Environments (FIRE) project, we examine the implications of bursty star formation histories for observations of high-redshift galaxies with JWST. We characterize how the galaxy observability depends on the star formation history. We also investigate selection effects due to bursty star formation on the physical properties measured, such as the gas fraction, specific star formation rate, and metallicity. We find the observability to be highly time-dependent for galaxies near the survey's limiting flux due to the SFR variability: as the star formation rate fluctuates, the same galaxy oscillates in and out of the observable sample. The observable fraction $f_\mathrm{obs} = 50\%$ at $z \sim 7$ and $M_{\star} \sim 10^{8.5}$ to $10^{9}\,M_{\odot}$ for a JWST/NIRCam survey reaching a limiting magnitude of $m^\mathrm{lim}_\mathrm{AB} \sim 29$-$30$, representative of surveys such as JADES and CEERS. JWST-detectable galaxies near the survey limit tend to have properties characteristic of galaxies in the bursty phase: on average, they show approximately 2.5 times higher cold, dense gas fractions and 20 times higher specific star formation rates at a given stellar mass than galaxies below the rest-UV detection threshold. Our study represents a first step in quantifying selection effects and the associated biases due to bursty star formation in studying high-redshift galaxy properties. △ Less

Submitted 18 September, 2023; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: 8 pages, 4 figures, resubmitted after incorporating referee's comments; analysis expanded to include more galaxies and some quantitative results corrected

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

Authors: Fiona McCluskey, Andrew Wetzel, Sarah R. Loebman, Jorge Moreno, Claude-Andre Faucher-Giguere, Philip F. Hopkins

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 disk settling and post-formation dynamical heating. We identify three eras: a Pre-Disk Era (typically >8 Gyr ago), when stars formed on dispersion-dominated orbits; an Early-Disk Era (~ 8 - 4… ▽ More 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 disk settling and post-formation dynamical heating. We identify three eras: a Pre-Disk Era (typically >8 Gyr ago), when stars formed on dispersion-dominated orbits; an Early-Disk Era (~ 8 - 4 Gyr ago), when stars started to form on rotation-dominated orbits but with high velocity dispersion, sigma_form; and a Late-Disk Era (< 4 Gyr ago), when stars formed with low sigma_form. sigma_form increased with time during the Pre-Disk Era, peaking ~ 8 Gyr ago, then decreased throughout the Early-Disk Era as the disk settled and remained low throughout the Late-Disk Era. By contrast, the velocity dispersion measured today, sigma_now, increases monotonically with age because of stronger post-formation heating for Pre-Disk stars. Importantly, most of sigma_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, sigma_R,form > sigma_phi,form > sigma_Z,form. Post-formation heating primarily increased sigma_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 disk began to settle correlates with its dynamical state today: earlier-settling galaxies currently form colder disks. Including stellar cosmic-ray feedback does not significantly change disk rotational support at fixed stellar mass. △ Less

Submitted 21 November, 2023; v1 submitted 24 March, 2023; originally announced March 2023.

Comments: 23 pages, accepted for publication in MNRAS

The impact of AGN-driven winds on physical and observable galaxy sizes

Authors: R. K. Cochrane, D. Anglés-Alcázar, J. Mercedes-Feliz, C. C. Hayward, C. -A. Faucher-Giguère, S. Wellons, B. A. Terrazas, A. Wetzel, P. F. Hopkins, J. Moreno, K. -Y. Su, R. S. Somerville

Abstract: Without AGN feedback, simulated massive, star-forming galaxies become too compact relative to observed galaxies at z<2. In this paper, we perform high-resolution re-simulations of a massive (M_star~10^11 M_sol) galaxy at z~2.3, drawn from the Feedback in Realistic Environments (FIRE) project. In the simulation without AGN feedback, the galaxy experiences a rapid starburst and shrinking of its half… ▽ More Without AGN feedback, simulated massive, star-forming galaxies become too compact relative to observed galaxies at z<2. In this paper, we perform high-resolution re-simulations of a massive (M_star~10^11 M_sol) galaxy at z~2.3, drawn from the Feedback in Realistic Environments (FIRE) project. In the simulation without AGN feedback, the galaxy experiences a rapid starburst and shrinking of its half-mass radius. We experiment with driving mechanical AGN winds, using a state-of-the-art hyper-Lagrangian refinement technique to increase particle resolution. These winds reduce the gas surface density in the inner regions of the galaxy, suppressing the compact starburst and maintaining an approximately constant half-mass radius. Using radiative transfer, we study the impact of AGN feedback on the magnitude and extent of the multi-wavelength continuum emission. When AGN winds are included, the suppression of the compact, dusty starburst results in lowered flux at FIR wavelengths (due to decreased star formation) but increased flux at optical-to-near-IR wavelengths (due to decreased dust attenuation, in spite of the lowered star formation rate), relative to the case without AGN winds. The FIR half-light radius decreases from ~1 kpc to ~0.1 kpc in <40 Myr when AGN winds are not included, but increases to ~2 kpc when they are. Interestingly, the half-light radius at optical-NIR wavelengths remains approximately constant over 35 Myr, for simulations with and without AGN winds. In the case without winds, this occurs despite the rapid compaction, and is due to heavy dust obscuration in the inner regions of the galaxy. This work highlights the importance of forward-modelling when comparing simulated and observed galaxy populations. △ Less

Submitted 17 May, 2023; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: 14 pages, 10 figures. Accepted for publication in MNRAS

Key Physical Processes in the Circumgalactic Medium

Authors: Claude-Andre Faucher-Giguere, S. Peng Oh

Abstract: Spurred by rich, multi-wavelength observations and enabled by new simulations, ranging from cosmological to sub-pc scales, the last decade has seen major theoretical progress in our understanding of the circumgalactic medium. We review key physical processes in the CGM. Our conclusions include: (1) The properties of the CGM depend on a competition between gravity-driven infall and gas cooling. Whe… ▽ More Spurred by rich, multi-wavelength observations and enabled by new simulations, ranging from cosmological to sub-pc scales, the last decade has seen major theoretical progress in our understanding of the circumgalactic medium. We review key physical processes in the CGM. Our conclusions include: (1) The properties of the CGM depend on a competition between gravity-driven infall and gas cooling. When cooling is slow relative to free fall, the gas is hot (roughly virial temperature) whereas the gas is cold (T~10^4 K) when cooling is rapid. (2) Gas inflows and outflows play crucial roles, as does the cosmological environment. Large-scale structure collimates cold streams and provides angular momentum. Satellite galaxies contribute to the CGM through winds and gas stripping. (3) In multiphase gas, the hot and cold phases continuously exchange mass, energy and momentum. The interaction between turbulent mixing and radiative cooling is critical. A broad spectrum of cold gas structures, going down to sub-pc scales, arises from fragmentation, coagulation, and condensation onto gas clouds. (4) Magnetic fields, thermal conduction and cosmic rays can substantially modify how the cold and hot phases interact, although microphysical uncertainties are presently large. Key open questions for future work include the mutual interplay between small-scale structure and large-scale dynamics, and how the CGM affects the evolution of galaxies. △ Less

Submitted 24 January, 2023; originally announced January 2023.

Comments: 69 pages, 13 figures. Accepted for publication in Annual Review of Astronomy and Astrophysics. Authors' draft. Edited version will appear in the next volume

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

The imprint of bursty star formation on alpha-element abundance patterns in Milky Way-like galaxies

Authors: Hanna Parul, Jeremy Bailin, Andrew Wetzel, Alexander B. Gurvich, Claude-André Faucher-Giguère, Zachary Hafen, Jonathan Stern, Owain Snaith

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

Submitted 18 January, 2023; originally announced January 2023.

Comments: 15 pages; Accepted 17 January 2023 in MNRAS

$\rm [C_{II}]$ 158 $\rm μm$ emission as an indicator of galaxy star formation rate

Authors: Lichen Liang, Robert Feldmann, Norman Murray, Desika Narayanan, Christopher C. Hayward, Daniel Anglés-Alcázar, Luigi Bassini, Alexander J. Richings, Claude-André Faucher-Giguère, Dongwoo T. Chung, Jennifer Y. H. Chan, Doǧa Tolgay, Onur Çatmabacak, Dušan Kereš, Philip F. Hopkins

Abstract: Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C_{II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C_{II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C_{II}]}{}/{}\rm SFR$ than the local SFGs, including the infrared-lumin… ▽ More Observations of local star-forming galaxies (SFGs) show a tight correlation between their singly ionized carbon line luminosity ($L_{\rm [C_{II}]}$) and star formation rate (SFR), suggesting that $L_{\rm [C_{II}]}$ may be a useful SFR tracer for galaxies. Some other galaxy populations, however, are found to have lower $L_{\rm [C_{II}]}{}/{}\rm SFR$ than the local SFGs, including the infrared-luminous, starburst galaxies at low and high redshifts, as well as some moderately star-forming galaxies at the epoch of re-ionization (EoR). The origin of this `$\rm [C_{II}]$ deficit' is unclear. In this work, we study the $L_{\rm [C_{II}]}$-SFR relation of galaxies using a sample of $z=0-8$ galaxies with $M_*\approx10^7-5\times10^{11}\,M_\odot$ extracted from cosmological volume and zoom-in simulations from the Feedback in Realistic Environments (FIRE) project. We find a simple analytic expression for $L_{\rm [C_{II}]}$/SFR of galaxies in terms of the following parameters: mass fraction of $\rm [C_{II}]$-emitting gas ($f_{\rm [C_{II}]}$), gas metallicity ($Z_{\rm gas}$), gas density ($n_{\rm gas}$) and gas depletion time ($t_{\rm dep}{}={}M_{\rm gas}{}/{}\rm SFR$). We find two distinct physical regimes, where $t_{\rm dep}$ ($Z_{\rm gas}$) is the main driver of the $\rm [C_{II}]$ deficit in $\rm H_2$-rich ($\rm H_2$-poor) galaxies. The observed $\rm [C_{II}]$ deficit of IR-luminous galaxies and early EoR galaxies, corresponding to the two different regimes, is due to short gas depletion time and low gas metallicity, respectively. Our result indicates that $\rm [C_{II}]$ deficit is a common phenomenon of galaxies, and caution needs to be taken when applying a constant $L_{\rm [C_{II}]}$-to-SFR conversion factor derived from local SFGs to estimate cosmic SFR density at high redshifts and interpret data from upcoming $\rm [C_{II}]$ line intensity mapping experiments. △ Less

Submitted 6 December, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

Comments: 44 page, 24 figures, 8 tables and 8 appendices, accepted for publication in MNRAS. Key figures: fig. 16 & 17. Comments are welcome

Local positive feedback in the overall negative: the impact of quasar winds on star formation in the FIRE cosmological simulations

Authors: Jonathan Mercedes-Feliz, Daniel Anglés-Alcázar, Christopher C. Hayward, Rachel K. Cochrane, Bryan A. Terrazas, Sarah Wellons, Alexander J. Richings, Claude-André Faucher-Giguère, Jorge Moreno, Kung Yi Su, Philip F. Hopkins, Eliot Quataert, Dušan Kereš

Abstract: Negative feedback from accreting supermassive black holes is regarded as a key ingredient in suppressing star formation and quenching massive galaxies. However, several models and observations suggest that black hole feedback may have a positive effect, triggering star formation by compressing interstellar medium gas to higher densities. We investigate the dual role of black hole feedback using co… ▽ More Negative feedback from accreting supermassive black holes is regarded as a key ingredient in suppressing star formation and quenching massive galaxies. However, several models and observations suggest that black hole feedback may have a positive effect, triggering star formation by compressing interstellar medium gas to higher densities. We investigate the dual role of black hole feedback using cosmological hydrodynamic simulations from the Feedback In Realistic Environments (FIRE) project, including a novel implementation of hyper-refined accretion-disc winds. Focusing on a massive, star-forming galaxy at $z \sim 2$ ($M_{\rm halo} \sim 10^{12.5} \, {\rm M}_{\odot}$), we show that strong quasar winds with kinetic power $\sim$10$^{46}$ erg/s acting for $>$20$\,$Myr drive the formation of a central gas cavity and can dramatically reduce the star formation rate surface density across the galaxy disc. The suppression of star formation is primarily driven by reducing the amount of gas that can become star-forming, compared to directly evacuating the pre-existing star-forming gas reservoir (preventive feedback dominates over ejective feedback). Despite the global negative impact of quasar winds, we identify several plausible signatures of local positive feedback, including: (1) spatial anti-correlation of wind-dominated regions and star-forming clumps, (2) higher local star formation efficiency in compressed gas near the edge of the cavity, and (3) increased local contribution of outflowing material to star formation. Stars forming under the presence of quasar winds tend to do so at larger radial distances. Our results suggest that positive and negative AGN feedback can coexist in galaxies, but local positive triggering of star formation plays a minor role in global galaxy growth. △ Less

Submitted 1 August, 2023; v1 submitted 4 January, 2023; originally announced January 2023.

Comments: 17 pages, 12 figures. Accepted for publication in MNRAS

A jolt to the system: ram pressure on low-mass galaxies in simulations of the Local Group

Authors: Jenna Samuel, Bhavya Pardasani, Andrew Wetzel, Isaiah Santistevan, Michael Boylan-Kolchin, Jorge Moreno, Claude-Andre Faucher-Giguere

Abstract: Low-mass galaxies are highly susceptible to environmental effects that can efficiently quench star formation. We explore the role of ram pressure in quenching low-mass galaxies ($M_{*}\sim10^{5-9}\,\rm{M}_{\odot}$) within 2 Mpc of Milky Way (MW) hosts using the FIRE-2 simulations. Ram pressure is highly variable across different environments, within individual MW haloes, and for individual low-mas… ▽ More Low-mass galaxies are highly susceptible to environmental effects that can efficiently quench star formation. We explore the role of ram pressure in quenching low-mass galaxies ($M_{*}\sim10^{5-9}\,\rm{M}_{\odot}$) within 2 Mpc of Milky Way (MW) hosts using the FIRE-2 simulations. Ram pressure is highly variable across different environments, within individual MW haloes, and for individual low-mass galaxies over time. The impulsiveness of ram pressure -- the maximum ram pressure scaled to the integrated ram pressure prior to quenching -- correlates with whether a galaxy is quiescent or star-forming. The time-scale between maximum ram pressure and quenching is anticorrelated with impulsiveness, such that high impulsiveness corresponds to quenching time-scales $<1$ Gyr. Galaxies in low-mass groups ($M_\mathrm{*,host}\sim10^{7-9}\,\rm{M}_{\odot}$) outside of MW haloes experience typical ram pressure only slightly lower than ram pressure on MW satellites, helping to explain effective quenching via group pre-processing. Ram pressure on MW satellites rises sharply with decreasing distance to the host, and, at a fixed physical distance, more recent pericentre passages are typically associated with higher ram pressure because of greater gas density in the inner host halo at late times. Furthermore, the ram pressure and gas density in the inner regions of Local Group-like paired host haloes is higher at small angles off the host galaxy disc compared to isolated hosts. The quiescent fraction of satellites within these low-latitude regions is also elevated in the simulations and observations, signaling possible anisotropic quenching via ram pressure around MW-mass hosts. △ Less

Submitted 25 August, 2023; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 13 pages, 10 figures, + appendices. Accepted for publication in MNRAS

Orientations of DM Halos in FIRE-2 Milky Way-mass Galaxies

Authors: Jay Baptista, Robyn Sanderson, Dan Huber, Andrew Wetzel, Omid Sameie, Michael Boylan-Kolchin, Jeremy Bailin, Philip F. Hopkins, Claude-Andre Faucher-Giguere, Sukanya Chakrabarti, Drona Vargya, Nondh Panithanpaisal, Arpit Arora, Emily Cunningham

Abstract: The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way's DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FI… ▽ More The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way's DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FIRE-2 Milky Way-mass galaxies to explore the evolution of the DM halo's orientation with radius and time, with or without a major merger with a Large Magellanic Cloud (LMC) analog, and when varying the DM model. In three of the four CDM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20 degrees beyond about 30 galactocentric kpc, reaching a maximum of 30--90 degrees depending on the individual halo's formation history. In identical simulations using a model of self-interacting DM with $σ= 1 \, \mathrm{cm}^2 \, \mathrm{g}^{-1}$, the halo remains aligned with the stellar disk out to $\sim$200--400 kpc. Interactions with massive satellites ($M \gtrsim 4 \times 10^{10} \, \rm{M_\odot}$ at pericenter; $M \gtrsim 3.3 \times 10^{10} \, \rm{M_\odot}$ at infall) affect the orientation of the halo significantly, aligning the halo's major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of SIDM if the effects of massive satellites can be accounted for. △ Less

Submitted 29 November, 2022; originally announced November 2022.

Comments: Submitted to ApJ, 19 pages, 12 figures, 3 tables

Active galactic nucleus jet feedback in hydrostatic halos

Authors: Rainer Weinberger, Kung-Yi Su, Kristian Ehlert, Christoph Pfrommer, Lars Hernquist, Greg L. Bryan, Volker Springel, Yuan Li, Blakesley Burkhart, Ena Choi, Claude-André Faucher-Giguère

Abstract: Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution funct… ▽ More Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution function only require resolutions sufficient to keep the jet-inflated cavities stable. Comparing different model variations, as well as an independent jet model using a different hydrodynamics code, we show that the dominant uncertainties are the choices of jet properties within a given model. Independent of implementation, we find that light, thermal jets with low momentum flux tend to delay the onset of a cooling flow more efficiently on a $50$ Myr timescale than heavy, kinetic jets. The delay of the cooling flow originates from a displacement and boost in entropy of the central gas. If the jet kinetic luminosity depends on accretion rate, collimated, light, hydrodynamic jets are able to reduce cooling flows in halos, without a need for jet precession or wide opening angles. Comparing the jet feedback with a `kinetic wind' implementation shows that equal amounts of star formation rate reduction can be achieved by different interactions with the halo gas: the jet has a larger effect on the hot halo gas while leaving the denser, star forming phase in place, while the wind acts more locally on the star forming phase, which manifests itself in different time-variability properties. △ Less

Submitted 5 May, 2023; v1 submitted 21 November, 2022; originally announced November 2022.

Comments: 22 pages, 20 figures, accepted for publication in MNRAS, comments welcome

A unified model for the co-evolution of galaxies and their circumgalactic medium: the relative roles of turbulence and atomic cooling physics

Authors: Viraj Pandya, Drummond B. Fielding, Greg L. Bryan, Christopher Carr, Rachel S. Somerville, Jonathan Stern, Claude-Andre Faucher-Giguere, Zachary Hafen, Daniel Angles-Alcazar, John C. Forbes

Abstract: The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM co-evolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and… ▽ More The circumgalactic medium (CGM) plays a pivotal role in regulating gas flows around galaxies and thus shapes their evolution. However, the details of how galaxies and their CGM co-evolve remain poorly understood. We present a new time-dependent two-zone model that self-consistently tracks not just mass and metal flows between galaxies and their CGM but also the evolution of the global thermal and turbulent kinetic energy of the CGM. Our model accounts for heating and turbulence driven by both supernova winds and cosmic accretion as well as radiative cooling, turbulence dissipation, and halo outflows due to CGM overpressurization. We demonstrate that, depending on parameters, the CGM can undergo a phase transition (``thermalization'') from a cool, turbulence-supported phase to a virial-temperature, thermally-supported phase. This CGM phase transition is largely determined by the ability of radiative cooling to balance heating from supernova winds and turbulence dissipation. We perform an initial calibration of our model to the FIRE-2 cosmological hydrodynamical simulations and show that it can approximately reproduce the baryon cycles of the simulated halos. In particular, we find that, for these parameters, the phase transition occurs at high-redshift in ultrafaint progenitors and at low redshift in classical $M_{\rm vir}\sim10^{11}M_{\odot}$ dwarfs, while Milky Way-mass halos undergo the transition at $z\approx0.5$. We see a similar transition in the simulations though it is more gradual, likely reflecting radial dependence and multi-phase gas not captured by our model. We discuss these and other limitations of the model and possible future extensions. △ Less

Submitted 27 August, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: Updated to reflect version accepted by ApJ; minor text edits to clarify limitations and caveats of model

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

Authors: L. Bassini, R. Feldmann, J. Gensior, C. C. Hayward, C. -A. Faucher-Giguère, E. Cenci, L. Liang, M. Bernardini

Abstract: Recent observations indicate that galactic outflows are ubiquitous in high redshift 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 (SFH) and galactic outflow properties of galaxies in massive haloes (… ▽ More Recent observations indicate that galactic outflows are ubiquitous in high redshift 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 (SFH) and galactic outflow properties of galaxies in massive haloes ($10^{12}M_{\odot}<M_{\rm vir} <5\times 10^{12}M_{\odot}$) at $z\gtrsim5.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 (AGN). The simulated galaxies resemble $z>4$ DSFGs, with SFRs of $\sim 1000\ M_{\odot}\rm yr^{-1}$ and molecular gas masses of $M_{\rm mol}\sim 10^{10}\ M_{\odot}$. However, the simulated galaxies are characterised by higher circular velocities than those observed in high-z DSFGs. The mass loading factors from stellar feedback are of the order of $\sim 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 (ISM). 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\%$. Finally, compared to FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high redshift sample. We argue that both effects originate from the higher total and gas surface densities that characterise high$-z$ massive systems. △ Less

Submitted 15 November, 2022; originally announced November 2022.

Comments: 17 pages, 12 figures. Comments welcome

The observable properties of galaxy accretion events in Milky Way-like galaxies in the FIRE-2 cosmological simulations

Authors: Danny Horta, Emily C. Cunningham, Robyn E. Sanderson, Kathryn V. Johnston, Nondh Panithanpaisal, Arpit Arora, Lina Necib, Andrew Wetzel, Jeremy Bailin, Claude-André Faucher-Giguère

Abstract: In the $Λ$-Cold Dark Matter model of the Universe, galaxies form in part through accreting satellite systems. Previous work have built an understanding of the signatures of these processes contained within galactic stellar halos. This work revisits that picture using seven Milky Way-like galaxies in the \textit{Latte} suite of FIRE-2 cosmological simulations. The resolution of these simulations al… ▽ More In the $Λ$-Cold Dark Matter model of the Universe, galaxies form in part through accreting satellite systems. Previous work have built an understanding of the signatures of these processes contained within galactic stellar halos. This work revisits that picture using seven Milky Way-like galaxies in the \textit{Latte} suite of FIRE-2 cosmological simulations. The resolution of these simulations allows a comparison of contributions from satellites above M$_{*}$$\gtrsim$10$\times$$^{7}$M$_{\odot}$, enabling the analysis of observable properties for disrupted satellites in a fully self-consistent and cosmological context. Our results show that, the time of accretion and the stellar mass of an accreted satellite are fundamental parameters that in partnership dictate the resulting spatial distribution, orbital energy, and [$α$/Fe]-[Fe/H] compositions of the stellar debris of such mergers $at$ $present$ $day$. These parameters also govern the resulting dynamical state of an accreted galaxy at $z=0$, leading to the expectation that the inner regions of the stellar halo (R$_{\mathrm{GC}}$ $\lesssim$30 kpc) should contain fully phase-mixed debris from both lower and higher mass satellites. In addition, we find that a significant fraction of the lower mass satellites accreted at early times deposit debris in the outer halo (R$_{\mathrm{GC}}$ $>$50 kpc) that are $not$ fully phased-mixed, indicating that they could be identified in kinematic surveys. Our results suggest that, as future surveys become increasingly able to map the outer halo of our Galaxy, they may reveal the remnants of long-dead dwarf galaxies whose counterparts are too faint to be seen $in$ $situ$ in higher redshift surveys. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: Submitted for publication in ApJ

Stellar feedback-regulated black hole growth: driving factors from nuclear to halo scales

Authors: Lindsey Byrne, Claude-André Faucher-Giguère, Jonathan Stern, Daniel Anglés-Alcázar, Sarah Wellons, Alexander B. Gurvich, Philip F. Hopkins

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

Submitted 17 October, 2022; originally announced October 2022.

Comments: 21 pages, 9 figures, submitted to MNRAS

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

Authors: Sijie Yu, James S. Bullock, Alexander B. Gurvich, Zachary Hafen, Jonathan Stern, Michael Boylan-Kolchin, Claude-André Faucher-Giguère, Andrew Wetzel, Philip F. Hopkins, Jorge Moreno

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

Submitted 19 June, 2023; v1 submitted 7 October, 2022; originally announced October 2022.

Comments: 16 pages, 10 figures; this version was accepted to MNRAS on June 9, 2023

The Cosmic Ultraviolet Baryon Survey (CUBS) V: On the Thermodynamic Properties of the Cool Circumgalactic Medium at $z < 1$

Authors: Zhijie Qu, Hsiao-Wen Chen, Gwen C. Rudie, Fakhri S. Zahedy, Sean D. Johnson, Erin Boettcher, Sebastiano Cantalupo, Mandy C. Chen, Kathy L. Cooksey, David DePalma, Claude-André Faucher-Giguère, Michael Rauch, Joop Schaye, Robert A. Simcoe

Abstract: This paper presents a systematic study of the photoionization and thermodynamic properties of the cool circumgalactic medium (CGM) as traced by rest-frame ultraviolet absorption lines around 26 galaxies at redshift $z\lesssim1$. The study utilizes both high-quality far-ultraviolet and optical spectra of background QSOs and deep galaxy redshift surveys to characterize the gas density, temperature,… ▽ More This paper presents a systematic study of the photoionization and thermodynamic properties of the cool circumgalactic medium (CGM) as traced by rest-frame ultraviolet absorption lines around 26 galaxies at redshift $z\lesssim1$. The study utilizes both high-quality far-ultraviolet and optical spectra of background QSOs and deep galaxy redshift surveys to characterize the gas density, temperature, and pressure of individual absorbing components and to resolve their internal non-thermal motions. The derived gas density spans more than three decades, from $\log (n_{\rm H}/{\rm cm^{-3}}) \approx -4$ to $-1$, while the temperature of the gas is confined in a narrow range of $\log (T/{\rm K})\approx 4.3\pm 0.3$. In addition, a weak anti-correlation between gas density and temperature is observed, consistent with the expectation of the gas being in photoionization equilibrium. Furthermore, decomposing the observed line widths into thermal and non-thermal contributions reveals that more than 30% of the components at $z\lesssim 1$ exhibit line widths driven by non-thermal motions, in comparison to $<20$% found at $z\approx 2$-3. Attributing the observed non-thermal line widths to intra-clump turbulence, we find that massive quenched galaxies on average exhibit higher non-thermal broadening/turbulent energy in their CGM compared to star-forming galaxies at $z\lesssim 1$. Finally, strong absorption features from multiple ions covering a wide range of ionization energy (e.g., from Mg II to O IV) can be present simultaneously in a single absorption system with kinematically aligned component structure, but the inferred pressure in different phases may differ by a factor of $\approx 10$. △ Less

Submitted 2 September, 2022; originally announced September 2022.

Comments: 16 pages, 4 figures, 1 table, accepted for publication in MNRAS

A New Discrete Dynamical Friction Estimator Based on $N$-body Simulations

Authors: Linhao Ma, Philip F. Hopkins, Luke Zoltan Kelley, Claude-André Faucher-Giguère

Abstract: A longstanding problem in galactic simulations is to resolve the dynamical friction (DF) force acting on massive black hole particles when their masses are comparable to or less than the background simulation particles. Many sub-grid models based on the traditional Chandrasekhar DF formula have been proposed, yet they suffer from fundamental ambiguities in the definition of some terms in Chandrase… ▽ More A longstanding problem in galactic simulations is to resolve the dynamical friction (DF) force acting on massive black hole particles when their masses are comparable to or less than the background simulation particles. Many sub-grid models based on the traditional Chandrasekhar DF formula have been proposed, yet they suffer from fundamental ambiguities in the definition of some terms in Chandrasekhar's formula when applied to real galaxies, as well as difficulty in evaluating continuous quantities from (spatially) discrete simulation data. In this work we present a new sub-grid dynamical friction estimator based on the discrete nature of $N$-body simulations, which also avoids the ambiguously-defined quantities in Chandrasekhar's formula. We test our estimator in the GIZMO code and find that it agrees well with high-resolution simulations where DF is fully captured, with negligible additional computational cost. We also compare it with a Chandrasekhar estimator and discuss its applications in real galactic simulations. △ Less

Submitted 15 May, 2023; v1 submitted 25 August, 2022; originally announced August 2022.

Comments: Accepted to MNRAS

Effects of the environment on the multiplicity properties of stars in the STARFORGE simulations

Authors: Dávid Guszejnov, Aman N. Raju, Stella S. R. Offner, Michael Y. Grudić, Claude-André Faucher-Giguère, Philip F. Hopkins, Anna L. Rosen

Abstract: Most observed stars are part of a multiple star system, but the formation of such systems and the role of environment and various physical processes is still poorly understood. We present a suite of radiation-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project that include stellar feedback with varied initial surface density, magnetic fields, level of turbul… ▽ More Most observed stars are part of a multiple star system, but the formation of such systems and the role of environment and various physical processes is still poorly understood. We present a suite of radiation-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project that include stellar feedback with varied initial surface density, magnetic fields, level of turbulence, metallicity, interstellar radiation field, simulation geometry and turbulent driving. In our fiducial cloud the raw simulation data reproduces the observed multiplicity fractions for Solar-type and higher mass stars, similar to previous works. However, after correcting for observational incompleteness the simulation under-predicts these values. The discrepancy is likely due to the lack of disk fragmentation, as the simulation only resolves multiples that form either through capture or core fragmentation. The raw mass distribution of companions is consistent with randomly drawing from the initial mass function for the companions of $>1\,\mathrm{M_\odot}$ stars, however, accounting for observational incompleteness produces a flatter distribution similar to observations. We show that stellar multiplicity changes as the cloud evolves and anti-correlates with stellar density. This relationship also explains most multiplicity variations between runs, i.e., variations in the initial conditions that increase stellar density (increased surface density, reduced turbulence) decrease multiplicity. While other parameters, such as metallicity, interstellar radiation, and geometry significantly affect the star formation history or the IMF, varying them produces no clear trend in stellar multiplicity properties. △ Less

Submitted 28 September, 2022; v1 submitted 4 August, 2022; originally announced August 2022.

Comments: 20 pages, 21 figures, submitted to MNRAS

The effects of local stellar radiation and dust depletion on non-equilibrium interstellar chemistry

Authors: Alexander J. Richings, Claude-Andre Faucher-Giguere, Alexander B. Gurvich, Joop Schaye, Christopher C. Hayward

Abstract: 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 meta… ▽ More 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$μ$m, [OI] 63$μ$m, [OIII] 88$μ$m, [NII] 122$μ$m and H$α$ 6563A. Our simulations with a uniform radiation field predict fainter luminosities, by up to an order of magnitude for [OIII] 88$μ$m and H$α$ 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. △ Less

Submitted 30 August, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 29 pages, 24 figures (including appendices). Accepted for publication in MNRAS. Minor changes relative to v1