FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation. Hopkins, P. F, Wetzel, A., Keres, D., Faucher-Giguere, C., Quataert, E., Boylan-Kolchin, M., Murray, N., Hayward, C. C., Garrison-Kimmel, S., Hummels, C., Feldmann, R., Torrey, P., Ma, X., Angles-Alcazar, D., Su, K., Orr, M., Schmitz, D., Escala, I., Sanderson, R., Grudic, M. Y., Hafen, Z., Kim, J., Fitts, A., Bullock, J. S., Wheeler, C., Chan, T. K., Elbert, O. D., & Narananan, D. ArXiv e-prints, 1702:arXiv:1702.06148, February, 2017.
FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation [link]Paper  abstract   bibtex   
The Feedback In Realistic Environments (FIRE) project explores the role of feedback in cosmological simulations of galaxy formation. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Now, motivated by the development of more accurate numerics (hydrodynamic solvers, gravitational softening, supernova coupling) and the 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 show FIRE-2 improvements do not qualitatively change galaxy-scale properties relative to FIRE-1. We then pursue an extensive study of numerics versus physics in galaxy simulations. Details of the star-formation (SF) algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and SF occurs at higher-than-mean densities. We present several new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are remarkably robust to the numerics that we test, provided that: (1) Toomre masses (cold disk scale heights) are resolved; (2) feedback coupling ensures conservation and isotropy, and (3) individual supernovae are time-resolved. As resolution increases, stellar masses and profiles converge first, followed by metal abundances and visual morphologies, then properties of winds and the circumgalactic medium. The central (\textasciitildekpc) mass concentration of massive (L*) galaxies is sensitive to numerics, particularly how winds ejected into hot halos are trapped, mixed, and recycled into the galaxy. Multiple feedback mechanisms are required to reproduce observations: SNe regulate stellar masses; OB/AGB mass loss fuels late-time SF; radiative feedback suppresses instantaneous SFRs and accretion onto dwarfs. We provide tables, initial conditions, and the numerical algorithms required to reproduce our simulations.
@article{hopkins_fire-2_2017,
	title = {{FIRE}-2 {Simulations}: {Physics} versus {Numerics} in {Galaxy} {Formation}},
	volume = {1702},
	shorttitle = {{FIRE}-2 {Simulations}},
	url = {http://adsabs.harvard.edu/abs/2017arXiv170206148H},
	abstract = {The Feedback In Realistic Environments (FIRE) project explores the role of feedback in cosmological simulations of galaxy formation. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Now, motivated by the development of more accurate numerics
(hydrodynamic solvers, gravitational softening, supernova coupling) and the 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 show FIRE-2 improvements do not qualitatively change galaxy-scale properties relative to FIRE-1. We then pursue an extensive study of numerics versus physics in galaxy simulations. Details of the star-formation (SF) algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and SF occurs at higher-than-mean densities. We present several new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are remarkably robust to the numerics that we test, provided that: (1) Toomre masses (cold disk scale heights) are resolved; (2) feedback coupling ensures conservation and isotropy, and (3) individual supernovae are time-resolved. As resolution increases, stellar masses and profiles converge first, followed by metal abundances and visual morphologies, then properties of winds and the circumgalactic medium. The central ({\textasciitilde}kpc) mass concentration of massive (L*) galaxies is sensitive to numerics, particularly how winds ejected into hot halos are trapped, mixed, and recycled into the galaxy. Multiple feedback mechanisms are required to reproduce observations: SNe regulate stellar masses; OB/AGB mass loss fuels late-time SF; radiative feedback suppresses instantaneous SFRs and accretion onto dwarfs. We provide tables, initial conditions, and the numerical algorithms required to reproduce our simulations.},
	urldate = {2017-03-09},
	journal = {ArXiv e-prints},
	author = {Hopkins, Philip F and Wetzel, Andrew and Keres, Dusan and Faucher-Giguere, Claude-Andre and Quataert, Eliot and Boylan-Kolchin, Michael and Murray, Norman and Hayward, Christopher C. and Garrison-Kimmel, Shea and Hummels, Cameron and Feldmann, Robert and Torrey, Paul and Ma, Xiangcheng and Angles-Alcazar, Daniel and Su, Kung-Yi and Orr, Matthew and Schmitz, Denise and Escala, Ivanna and Sanderson, Robyn and Grudic, Michael Y. and Hafen, Zachary and Kim, Ji-Hoon and Fitts, Alex and Bullock, James S. and Wheeler, Coral and Chan, T. K. and Elbert, Oliver D. and Narananan, Desika},
	month = feb,
	year = {2017},
	keywords = {Astrophysics - Astrophysics of Galaxies, Astrophysics - Cosmology and Nongalactic Astrophysics, Astrophysics - Instrumentation and Methods for Astrophysics},
	pages = {arXiv:1702.06148},
}

Downloads: 0