The robustness of cosmological hydrodynamic simulation predictions to changes in numerics and cooling physics

Shuiyao Huang, Neal Katz, Romeel S Dave, Mark Fardal, Juna Kollmeier, Benjamin D. Oppenheimer, Molly S. Peeples, Shawn Roberts, David H. Weinberg, Philip F. Hopkins, Robert Thompson

Research output: Contribution to journalArticle

Abstract

We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined time-step criterion, and metal-line cooling that accounts for photoionization in the presence of a recently refined Haardt & Madau model of the ionizing background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations, it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Lyα forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low-mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII.

Original languageEnglish (US)
Pages (from-to)2021-2046
Number of pages26
JournalMonthly Notices of the Royal Astronomical Society
Volume484
Issue number2
DOIs
StatePublished - Jan 1 2019
Externally publishedYes

Fingerprint

physics
hydrodynamics
cooling
galaxies
prediction
predictions
simulation
metal
accretion
gas
statistics
metals
formulations
Kelvin-Helmholtz instability
cold gas
surface tension
high temperature gases
clumps
recycling
stellar mass

Keywords

  • Galaxies: evolution
  • Galaxies: general
  • Hydrodynamics
  • Methods: numerical

ASJC Scopus subject areas

  • Astronomy and Astrophysics
  • Space and Planetary Science

Cite this

The robustness of cosmological hydrodynamic simulation predictions to changes in numerics and cooling physics. / Huang, Shuiyao; Katz, Neal; Dave, Romeel S; Fardal, Mark; Kollmeier, Juna; Oppenheimer, Benjamin D.; Peeples, Molly S.; Roberts, Shawn; Weinberg, David H.; Hopkins, Philip F.; Thompson, Robert.

In: Monthly Notices of the Royal Astronomical Society, Vol. 484, No. 2, 01.01.2019, p. 2021-2046.

Research output: Contribution to journalArticle

Huang, S, Katz, N, Dave, RS, Fardal, M, Kollmeier, J, Oppenheimer, BD, Peeples, MS, Roberts, S, Weinberg, DH, Hopkins, PF & Thompson, R 2019, 'The robustness of cosmological hydrodynamic simulation predictions to changes in numerics and cooling physics', Monthly Notices of the Royal Astronomical Society, vol. 484, no. 2, pp. 2021-2046. https://doi.org/10.1093/mnras/stz057
Huang, Shuiyao ; Katz, Neal ; Dave, Romeel S ; Fardal, Mark ; Kollmeier, Juna ; Oppenheimer, Benjamin D. ; Peeples, Molly S. ; Roberts, Shawn ; Weinberg, David H. ; Hopkins, Philip F. ; Thompson, Robert. / The robustness of cosmological hydrodynamic simulation predictions to changes in numerics and cooling physics. In: Monthly Notices of the Royal Astronomical Society. 2019 ; Vol. 484, No. 2. pp. 2021-2046.
@article{fa7ecf0dccc84388aaf1f2482e2c8e08,
title = "The robustness of cosmological hydrodynamic simulation predictions to changes in numerics and cooling physics",
abstract = "We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined time-step criterion, and metal-line cooling that accounts for photoionization in the presence of a recently refined Haardt & Madau model of the ionizing background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations, it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Lyα forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low-mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII.",
keywords = "Galaxies: evolution, Galaxies: general, Hydrodynamics, Methods: numerical",
author = "Shuiyao Huang and Neal Katz and Dave, {Romeel S} and Mark Fardal and Juna Kollmeier and Oppenheimer, {Benjamin D.} and Peeples, {Molly S.} and Shawn Roberts and Weinberg, {David H.} and Hopkins, {Philip F.} and Robert Thompson",
year = "2019",
month = "1",
day = "1",
doi = "10.1093/mnras/stz057",
language = "English (US)",
volume = "484",
pages = "2021--2046",
journal = "Monthly Notices of the Royal Astronomical Society",
issn = "0035-8711",
publisher = "Oxford University Press",
number = "2",

}

TY - JOUR

T1 - The robustness of cosmological hydrodynamic simulation predictions to changes in numerics and cooling physics

AU - Huang, Shuiyao

AU - Katz, Neal

AU - Dave, Romeel S

AU - Fardal, Mark

AU - Kollmeier, Juna

AU - Oppenheimer, Benjamin D.

AU - Peeples, Molly S.

AU - Roberts, Shawn

AU - Weinberg, David H.

AU - Hopkins, Philip F.

AU - Thompson, Robert

PY - 2019/1/1

Y1 - 2019/1/1

N2 - We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined time-step criterion, and metal-line cooling that accounts for photoionization in the presence of a recently refined Haardt & Madau model of the ionizing background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations, it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Lyα forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low-mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII.

AB - We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined time-step criterion, and metal-line cooling that accounts for photoionization in the presence of a recently refined Haardt & Madau model of the ionizing background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations, it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Lyα forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low-mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII.

KW - Galaxies: evolution

KW - Galaxies: general

KW - Hydrodynamics

KW - Methods: numerical

UR - http://www.scopus.com/inward/record.url?scp=85063374935&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85063374935&partnerID=8YFLogxK

U2 - 10.1093/mnras/stz057

DO - 10.1093/mnras/stz057

M3 - Article

AN - SCOPUS:85063374935

VL - 484

SP - 2021

EP - 2046

JO - Monthly Notices of the Royal Astronomical Society

JF - Monthly Notices of the Royal Astronomical Society

SN - 0035-8711

IS - 2

ER -