FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
read the original abstract
The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Motivated by the development of more accurate numerics - including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms - and exploration of new physics (e.g. magnetic fields), we introduce FIRE-2, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (~kpc) mass concentrations in massive (L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.
This paper has not been read by Pith yet.
Forward citations
Cited by 51 Pith papers
-
Lyman-alpha Pressure Strongly Enhances Pre-Supernova Feedback at Cosmic Dawn: The First Multi-Dimensional Lyman-alpha Radiation Hydrodynamics Simulations
First 2D Lyα RHD simulations show Lyman-alpha radiation pressure yields radiative forces of 2-16 times L_bol/c and force multipliers of 10-60, dominating other pre-supernova feedback in metal-poor environments.
-
An Acceleration is Worth a Hundred Thousand Phase Space Measurements
One acceleration measurement equals ~10^5 phase-space measurements for local dark matter density estimation, with acceleration outperforming Jeans modeling in both equilibrium and perturbed Milky Way simulations.
-
Strongest constraints on dark acoustic oscillations from the Lyman-alpha forest
First constraints from Ly-alpha forest limit DAO participation to at most 30% of dark matter for peaks below 50 h/Mpc at 95% CL using a deep kernel emulator of simulations.
-
The Lifecycle and Emission Properties of PAHs in Cosmological Hydrodynamic Galaxy Formation Simulations
Cosmological zoom-in simulations find that grain-grain shattering in diffuse ISM gas drives rising PAH mass fraction with time, naturally producing the observed PAH-metallicity relation and inverse qPAH-molecular gas trends.
-
Individual Star Sampling in Star Formation Simulations: A Semi-Deterministic Model
Introduces a semi-deterministic star sampling method in simulations that derives the instantaneous IMF from current cluster mass using reservoir particles and on-the-fly cluster finding.
-
TDEs on FIRE: Illuminating the Cosmic Evolution of Tidal Disruption Rates
FIRE-2 simulations show per-galaxy tidal disruption rates peak near z=2.5 at 4e-4 per year, correlate with SFR and central density, and remain high in satellite galaxies at early times.
-
Mergers Matter: Gravothermal Collapse in Dwarf Halos with Self-Interacting Dark Matter
SIDM simulations of dwarf halos show that quiescent merger histories produce gravothermal core collapse while sustained mergers prevent collapse and can yield central densities below gravothermal fluid model predictions.
-
Tunneling and tidal stripping in multifield ultralight dark matter halos
A semiclassical tunneling model shows that two-field ultralight DM halos have stability bounds that can be relaxed for some density-mass ratios but become more stringent across much of the parameter space compared to ...
-
Resolving the Unresolved Galactic Winds in Multi-phase Models. I. Methodology and Application
A new fitting methodology applied to UV absorption data recovers radial trends in galactic wind velocities and mass-loading factors by constraining initial hot and cool phase parameters in a multiphase model.
-
The LMC Corona Favors a First Passage
Idealized simulations with live gas particles show the LMC corona's present-day velocity and column density profiles match a first-passage orbit but are too low in a second-passage orbit, yielding truncation radii of ...
-
Reconstructing Galactic Gravitational Potentials from Stellar Kinematics with Physics-Informed Neural Networks
A PINN approach learns galactic gravitational potentials from acceleration data, achieving sub-percent errors on simulations while outperforming analytic models and retaining interpretability via structured priors.
-
Feedback-Free Star Formation in Clusters within a Galaxy Simulated at High Resolution in Cosmic Dawn
High-resolution zoom-in simulation of a z~10 galaxy produces over 10^5 star clusters with scale-free mass function, 90% of star formation in clusters, global SFE 0.2-0.3, and reproduces JWST super-bright galaxies via ...
-
Learning the Universe with PRFM-vol: Introducing a new subgrid model for star formation in cosmological simulations
PRFM-vol is a new subgrid star formation model for cosmological simulations that computes SFR from ambient densities via PRFM theory and a modified effective EOS, producing taller stellar scale heights, slightly highe...
-
An Inverse-Compton-Boosted Cool Core Unifies Perseus's Radio and X-ray Halos
Ancient cosmic ray halos from the central galaxy boost Perseus's cool core via inverse-Compton scattering, simultaneously explaining radio minihalo, giant halo, X-ray properties, and gamma-ray data without re-acceleration.
-
$V/\sigma$ Trends with Mass for Dwarf Galaxies from the Marvelous Massive Dwarfs Suite
Simulations of dwarf galaxies show V/σ rising with stellar mass, with HI gas and young stars more rotation-supported than old stars.
-
Efficiently emulating distribution functions in gigaparsec volumes for varying cosmological parameters
A new overdensity-conditioned emulator trained on small subvolumes from Quijote recovers the global halo mass function via integration over the overdensity distribution at 0.026% of the simulation cost.
-
Little Red Dots on FIRE: The Ability of Bursty Galaxies to Host an Abundant Population of High-Redshift AGN
FIRE-2 simulations with gravitational torque-driven and free-fall accretion models predict enough high-redshift AGN to explain little red dots, with a super-Eddington Eddington-limited scenario for M_BH >= 2e5 Msun in...
-
The DREAMS Project: Disentangling the Impact of Halo-to-Halo Variance and Baryonic Feedback on Milky Way Dark Matter Density Profiles
Milky Way-mass dark matter density profiles in IllustrisTNG are largely insensitive to astrophysics and cosmology variations, dominated by halo-to-halo variance instead.
-
Metallicity Gradients in Modern Cosmological Simulations II: The Role of Bursty Versus Smooth Feedback at High-Redshift
Bursty stellar feedback produces systematically flatter metallicity gradients than smooth feedback in high-redshift galaxies across multiple simulation suites.
-
Galaxy Metallicity Gradients in the Reionization Epoch from the FIRE-2 Simulations
FIRE-2 simulations find metallicity gradients in EoR galaxies flatten from median -0.15 dex/kpc at z~10 to -0.1 at z~6, with positive correlations to stellar mass and gas flow proxy Δv/2σ and links to central SFR density.
-
OCCAM X. Neutron Capture Abundances with Keck/HIRES & Magellan/MIKE
New high-resolution spectra yield abundances for 7 neutron-capture elements in open cluster stars, revealing flat Milky Way gradients for second-peak s- and r-process species and shallower slopes for first-peak s-process.
-
Columba: isolated dwarf galaxy populations in diverse cosmological environments simulated with a cold interstellar medium
New hydrodynamical simulations show that dwarf galaxy stellar mass-halo mass relations and star formation histories are more influenced by host halo concentration than by the 5 cMpc scale environment.
-
Investigating black hole accretion and feedback self-regulation in Seyfert galaxies using the FIRE-3 cosmological hydrodynamic simulations
FIRE-3 cosmological simulations of Seyfert galaxies produce episodic AGN feedback and gas clearing but no clear anti-correlation between nuclear gas concentration and AGN luminosity, highlighting timing mismatches wit...
-
The dependence of Circumgalactic Medium properties on halo assembly histories in the IllustrisTNG simulations
TNG50 simulation analysis finds early-forming halos show higher stellar mass, metallicity, lower CGM gas mass and sSFR at z=0, with mass-dependent CGM metallicity reversal and cold-gas kinematics differences.
-
Characterising magnetic fields at the onset of star cluster formation: From giant molecular clouds to infrared dark clumps
Polarization observations reveal scale-dependent differences in magnetic field morphology between molecular clouds and clumps, a velocity-dispersion correlation, and unreliable field-strength estimates that contradict...
-
TNG SAM: Bridging Hydrodynamical Complexity and Semi-Analytic Efficiency to Model Galaxy Formation
The TNG SAM reproduces TNG hydro simulation gas and metal flows plus galaxy and halo properties within 30% accuracy out to z=6 via five targeted updates to the Santa Cruz SAM calibrated on stellar feedback-dominated galaxies.
-
Resolving galaxy formation in the early Universe with BonFIRE and CampFIRE
BonFIRE and CampFIRE simulations show bursty clustered star formation in early galaxies and predict UV luminosity functions matching observations at faint magnitudes with a turnover at M_UV approx -14 but overpredicti...
-
Introducing the Lumina project: large-volume radiation-hydrodynamic simulations of the epochs of hydrogen and helium reionization
Lumina runs a 500 cMpc radiation-hydrodynamic simulation combining IllustrisTNG galaxy formation with six-bin M1 radiation transport to predict late stellar-driven HI reionization ending around z=4.75 and AGN-driven H...
-
The Simulated Oxygen Shortage (SOS): Mapping the Missing OVI in Simulated Dwarf Galaxies to Subgrid Physics
The paper identifies underproduction of oxygen in low-mass simulated dwarf galaxies as the likely cause of missing OVI in the CGM, based on comparisons across two simulation suites.
-
Radial redistribution of stellar orbits in FIRE simulations of Milky-Way-mass galaxies
FIRE-2 simulations show that stellar radial redistribution scatter saturates at ~2 kpc for stars older than ~3 Gyr, with net orbital changes depending on age and current radius, broadly matching Milky Way observations.
-
Environmental Quenching of High-Redshift Galaxies: Interpreting JWST Observations with Simulations
Multiple galaxy formation simulations show that low-mass quenched galaxies at z>3 are predominantly environmentally quenched satellites, often only temporarily so, and match JWST observations.
-
Dynamically cold discs in high-redshift galaxies: comparison between ALMA observations and TNG50
TNG50 shows most massive high-z star-forming galaxies are dynamically hotter than ALMA data indicate, with rare cold discs forming from aligned accretion and evolving into one-third discs and two-thirds early-type gal...
-
GOALS-JWST: Resolved multi-phase molecular gas in IRAS 20551-4250 using JWST and ALMA
Multi-phase molecular gas in IRAS20551-4250 is dominated by cold CO, shows UV-heated warm H2, tidal features from a merger, and no molecular outflows, consistent with ongoing star formation.
-
Too Big to Quench? I. Constraining ISM Stripping of Dwarf Satellites in Milky Way-like Halos
Simulations constrain the mass scale for efficient ISM stripping of dwarf satellites to M_star ≲ 10^7 M_sun in MW-like halos, 0.5-1 dex below observed values, indicating additional quenching mechanisms are needed.
-
cuRAMSES: Scalable AMR Optimizations for Large-Scale Cosmological Simulations
cuRAMSES replaces Hilbert-curve domain decomposition with recursive k-section partitioning and adds Morton-key hashing plus spatial binning to cut communication volume and accelerate feedback routines by up to 260x wh...
-
Stellar feedback drives the baryon deficiency in low-mass galaxies
Stacking DESI spectra reveals star-formation-driven Mg II outflows from low-mass galaxies that escape dark matter halos, providing indirect evidence that stellar feedback causes baryon deficiency.
-
FOGGIE: Figuring Out Gas & Galaxies In Enzo XII. The Formation and Evolution of Extended HI Galactic Disks and Warps with a Dynamic Circumgalactic medium
FOGGIE cosmological simulations find that less HI-populated CGMs produce thin, coherently rotating extended disks while more populated CGMs yield irregular morphologies, with all systems settling kinematically by z=0 ...
-
COLIBRE: calibrating subgrid feedback in cosmological simulations that include a cold gas phase
COLIBRE calibrates supernova and AGN feedback parameters in multi-phase ISM cosmological simulations using emulator-based fitting to reproduce the z=0 galaxy stellar mass function and size-stellar mass relation at thr...
-
Resolved Ages and Stellar Metallicities in Progenitors of Milky Way Analogs: A Closer Look at their Star Formation Histories since $z=5$
Resolved stellar property gradients in Milky Way analog progenitors show inside-out assembly with minor, temporary disruption from major mergers.
-
Blowing star formation away in AGN hosts (BAH) -- V: The Feeding-Feedback Cycle in local AGNs as revealed by their stellar populations
Stellar population synthesis of three AGN hosts shows recent rejuvenation and central metallicity drops consistent with gas inflows fueling both star formation and AGN activity.
-
Inside-Out vs. Outside-In Quenching of MaNGA Galaxies: Dependence on Stellar Mass and Environment
Different quenching tracers applied to MaNGA galaxies produce varying inside-out versus outside-in fractions, with inside-out rising at higher stellar mass while halo-mass trends are weaker.
-
It's Not Just Star Formation: A trend of low dark matter densities in the Andromeda dwarf galaxy system
Five of seven modeled M31 dwarf spheroidals show anomalously low central DM densities at 150 pc, with star formation heating disfavored as the sole cause.
-
Determination of the best dark matter profile for the Milky Way with Gaia DR3 using Bayesian Model Comparison
Bayesian comparison of Gaia DR3 rotation curves favors the Einasto dark matter profile over NFW, prefers cored over cuspy profiles, and finds MOND variants provide poorer fits than the best dark matter models.
-
Unmasking Stellar Feedback-Driven Bubbles: Identification and Properties Analysis
Simulations of an NGC 300-like galaxy find exponential distributions of feedback-driven bubble lifetimes and sizes that increase with galactocentric radius, plus matching Hα predictions.
-
The diverse outcomes of binary white dwarf mergers and connections to Galactic LISA sources
Mock catalogs of Milky Way white dwarf mergers are created with COSMIC and released publicly, summarizing outcomes by mass and composition while varying binary evolution uncertainties to connect with LISA sources.
-
Contribution of interstellar objects to local dark matter density
Interstellar objects may contribute enough baryonic mass to reduce the local dark matter halo density to 0.24 GeV/cm³.
-
Beyond Cloud-9: The case for discovering more HI-rich failed halos
Comparisons of three cosmological simulations show HI-rich failed halos occupy different mass regimes and predict that more can be discovered locally in HI-poor environments rather than at high redshift.
-
Effects of Varied Cosmic Ray Feedback from AGN on Massive Galaxy Properties
Varied cosmic ray feedback models from AGN in FIRE-3 simulations all quench massive galaxies consistently with observations but produce vastly different circumgalactic medium properties.
-
Blowing star formation away in AGN hosts (BAH) -- V: The Feeding-Feedback Cycle in local AGNs as revealed by their stellar populations
Spatially resolved stellar population synthesis of three local AGN hosts shows recent rejuvenation and central metallicity drops consistent with gas inflows fueling both star formation and AGN activity.
-
Cosmological Galaxy Formation Modelling in the Era of the Square Kilometre Array
Review of state-of-the-art cosmological galaxy formation models for HI, molecular gas and radio continuum in preparation for SKA, advocating coordinated multi-scale simulations, forward modelling and AI emulators.
-
Revealing Cosmic Ecosystems with the Hubble Space Telescope in 2030s and Beyond
HST UV spectroscopy is presented as the unique tool for probing multiphase gas at the disk-CGM interface to understand how galaxies acquire fuel, recycle metals, and drive feedback.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.