Chemical enrichment in cosmological, smoothed particle hydrodynamics simulations. Wiersma, R. P. C., Schaye, J., Theuns, T., Dalla Vecchia, C., & Tornatore, L. Monthly Notices of the Royal Astronomical Society, 399:574–600, October, 2009.
Paper doi abstract bibtex We present an implementation of stellar evolution and chemical feedback for smoothed particle hydrodynamics simulations. We consider the timed release of individual elements by both massive (Type II supernovae and stellar winds) and intermediate-mass stars (Type Ia supernovae and asymptotic giant branch stars). We illustrate the results of our method using a suite of cosmological simulations that include new prescriptions for radiative cooling, star formation and galactic winds. Radiative cooling is implemented element-by-element, in the presence of an ionizing radiation background, and we track all 11 elements that contribute significantly to the radiative cooling. While all simulations presented here use a single set of physical parameters, we take specific care to investigate the robustness of the predictions of chemodynamical simulations with respect to the ingredients, the methods and the numerical convergence. A comparison of nucleosynthetic yields taken from the literature indicates that relative abundance ratios may only be reliable at the factor of 2 level, even for a fixed initial mass function. Abundances relative to iron are even more uncertain because the rate of Type Ia supernovae is not well known. We contrast two reasonable definitions of the metallicity of a resolution element and find that while they agree for high metallicities, there are large differences at low metallicities. We argue that the discrepancy is indicative of the lack of metal mixing caused by the fact that metals are stuck to particles. We argue that since this is a (numerical) sampling problem, solving it by using a poorly constrained physical process such as diffusion could have undesired consequences. We demonstrate that the two metallicity definitions result in redshift z = 0 stellar masses that can differ by up to a factor of 2, because of the sensitivity of the cooling rates to the elemental abundances. Finally, we use several 5123 particle simulations to investigate the evolution of the distribution of heavy elements, which we find to be in reasonably good agreement with available observational constraints. We find that by z = 0 most of the metals are locked up in stars. The gaseous metals are distributed over a very wide range of gas densities and temperatures. The shock-heated warm-hot intergalactic medium has a relatively high metallicity of \textasciitilde10-1Zsolar that evolves only weakly, and is therefore an important reservoir of metals. Any census aiming to account for most of the metal mass will have to take a wide variety of objects and structures into account.
@article{wiersma_chemical_2009,
title = {Chemical enrichment in cosmological, smoothed particle hydrodynamics simulations},
volume = {399},
issn = {0035-8711},
url = {http://adsabs.harvard.edu/abs/2009MNRAS.399..574W},
doi = {10.1111/j.1365-2966.2009.15331.x},
abstract = {We present an implementation of stellar evolution and chemical feedback
for smoothed particle hydrodynamics simulations. We consider the timed
release of individual elements by both massive (Type II supernovae and
stellar winds) and intermediate-mass stars (Type Ia supernovae and
asymptotic giant branch stars). We illustrate the results of our method
using a suite of cosmological simulations that include new prescriptions
for radiative cooling, star formation and galactic winds. Radiative
cooling is implemented element-by-element, in the presence of an
ionizing radiation background, and we track all 11 elements that
contribute significantly to the radiative cooling.
While all simulations presented here use a single set of physical
parameters, we take specific care to investigate the robustness of the
predictions of chemodynamical simulations with respect to the
ingredients, the methods and the numerical convergence. A comparison of
nucleosynthetic yields taken from the literature indicates that relative
abundance ratios may only be reliable at the factor of 2 level, even for
a fixed initial mass function. Abundances relative to iron are even more
uncertain because the rate of Type Ia supernovae is not well known. We
contrast two reasonable definitions of the metallicity of a resolution
element and find that while they agree for high metallicities, there are
large differences at low metallicities. We argue that the discrepancy is
indicative of the lack of metal mixing caused by the fact that metals
are stuck to particles. We argue that since this is a (numerical)
sampling problem, solving it by using a poorly constrained physical
process such as diffusion could have undesired consequences. We
demonstrate that the two metallicity definitions result in redshift z =
0 stellar masses that can differ by up to a factor of 2, because of the
sensitivity of the cooling rates to the elemental abundances.
Finally, we use several 5123 particle simulations to
investigate the evolution of the distribution of heavy elements, which
we find to be in reasonably good agreement with available observational
constraints. We find that by z = 0 most of the metals are locked up in
stars. The gaseous metals are distributed over a very wide range of gas
densities and temperatures. The shock-heated warm-hot intergalactic
medium has a relatively high metallicity of
{\textasciitilde}10-1Zsolar that evolves only weakly, and is
therefore an important reservoir of metals. Any census aiming to account
for most of the metal mass will have to take a wide variety of objects
and structures into account.},
urldate = {2019-05-13},
journal = {Monthly Notices of the Royal Astronomical Society},
author = {Wiersma, Robert P. C. and Schaye, Joop and Theuns, Tom and Dalla Vecchia, Claudio and Tornatore, Luca},
month = oct,
year = {2009},
keywords = {cosmology: theory, galaxies: abundances, galaxies: formation, intergalactic medium, methods: numerical},
pages = {574--600},
}
Downloads: 0
{"_id":"sNLHnS5ppfGzT8oen","bibbaseid":"wiersma-schaye-theuns-dallavecchia-tornatore-chemicalenrichmentincosmologicalsmoothedparticlehydrodynamicssimulations-2009","author_short":["Wiersma, R. P. C.","Schaye, J.","Theuns, T.","Dalla Vecchia, C.","Tornatore, L."],"bibdata":{"bibtype":"article","type":"article","title":"Chemical enrichment in cosmological, smoothed particle hydrodynamics simulations","volume":"399","issn":"0035-8711","url":"http://adsabs.harvard.edu/abs/2009MNRAS.399..574W","doi":"10.1111/j.1365-2966.2009.15331.x","abstract":"We present an implementation of stellar evolution and chemical feedback for smoothed particle hydrodynamics simulations. We consider the timed release of individual elements by both massive (Type II supernovae and stellar winds) and intermediate-mass stars (Type Ia supernovae and asymptotic giant branch stars). We illustrate the results of our method using a suite of cosmological simulations that include new prescriptions for radiative cooling, star formation and galactic winds. Radiative cooling is implemented element-by-element, in the presence of an ionizing radiation background, and we track all 11 elements that contribute significantly to the radiative cooling. While all simulations presented here use a single set of physical parameters, we take specific care to investigate the robustness of the predictions of chemodynamical simulations with respect to the ingredients, the methods and the numerical convergence. A comparison of nucleosynthetic yields taken from the literature indicates that relative abundance ratios may only be reliable at the factor of 2 level, even for a fixed initial mass function. Abundances relative to iron are even more uncertain because the rate of Type Ia supernovae is not well known. We contrast two reasonable definitions of the metallicity of a resolution element and find that while they agree for high metallicities, there are large differences at low metallicities. We argue that the discrepancy is indicative of the lack of metal mixing caused by the fact that metals are stuck to particles. We argue that since this is a (numerical) sampling problem, solving it by using a poorly constrained physical process such as diffusion could have undesired consequences. We demonstrate that the two metallicity definitions result in redshift z = 0 stellar masses that can differ by up to a factor of 2, because of the sensitivity of the cooling rates to the elemental abundances. Finally, we use several 5123 particle simulations to investigate the evolution of the distribution of heavy elements, which we find to be in reasonably good agreement with available observational constraints. We find that by z = 0 most of the metals are locked up in stars. The gaseous metals are distributed over a very wide range of gas densities and temperatures. The shock-heated warm-hot intergalactic medium has a relatively high metallicity of \\textasciitilde10-1Zsolar that evolves only weakly, and is therefore an important reservoir of metals. Any census aiming to account for most of the metal mass will have to take a wide variety of objects and structures into account.","urldate":"2019-05-13","journal":"Monthly Notices of the Royal Astronomical Society","author":[{"propositions":[],"lastnames":["Wiersma"],"firstnames":["Robert","P.","C."],"suffixes":[]},{"propositions":[],"lastnames":["Schaye"],"firstnames":["Joop"],"suffixes":[]},{"propositions":[],"lastnames":["Theuns"],"firstnames":["Tom"],"suffixes":[]},{"propositions":[],"lastnames":["Dalla","Vecchia"],"firstnames":["Claudio"],"suffixes":[]},{"propositions":[],"lastnames":["Tornatore"],"firstnames":["Luca"],"suffixes":[]}],"month":"October","year":"2009","keywords":"cosmology: theory, galaxies: abundances, galaxies: formation, intergalactic medium, methods: numerical","pages":"574–600","bibtex":"@article{wiersma_chemical_2009,\n\ttitle = {Chemical enrichment in cosmological, smoothed particle hydrodynamics simulations},\n\tvolume = {399},\n\tissn = {0035-8711},\n\turl = {http://adsabs.harvard.edu/abs/2009MNRAS.399..574W},\n\tdoi = {10.1111/j.1365-2966.2009.15331.x},\n\tabstract = {We present an implementation of stellar evolution and chemical feedback \nfor smoothed particle hydrodynamics simulations. We consider the timed\nrelease of individual elements by both massive (Type II supernovae and\nstellar winds) and intermediate-mass stars (Type Ia supernovae and\nasymptotic giant branch stars). We illustrate the results of our method\nusing a suite of cosmological simulations that include new prescriptions\nfor radiative cooling, star formation and galactic winds. Radiative\ncooling is implemented element-by-element, in the presence of an\nionizing radiation background, and we track all 11 elements that\ncontribute significantly to the radiative cooling.\n\nWhile all simulations presented here use a single set of physical\nparameters, we take specific care to investigate the robustness of the\npredictions of chemodynamical simulations with respect to the\ningredients, the methods and the numerical convergence. A comparison of\nnucleosynthetic yields taken from the literature indicates that relative\nabundance ratios may only be reliable at the factor of 2 level, even for\na fixed initial mass function. Abundances relative to iron are even more\nuncertain because the rate of Type Ia supernovae is not well known. We\ncontrast two reasonable definitions of the metallicity of a resolution\nelement and find that while they agree for high metallicities, there are\nlarge differences at low metallicities. We argue that the discrepancy is\nindicative of the lack of metal mixing caused by the fact that metals\nare stuck to particles. We argue that since this is a (numerical)\nsampling problem, solving it by using a poorly constrained physical\nprocess such as diffusion could have undesired consequences. We\ndemonstrate that the two metallicity definitions result in redshift z =\n0 stellar masses that can differ by up to a factor of 2, because of the\nsensitivity of the cooling rates to the elemental abundances.\n\nFinally, we use several 5123 particle simulations to\ninvestigate the evolution of the distribution of heavy elements, which\nwe find to be in reasonably good agreement with available observational\nconstraints. We find that by z = 0 most of the metals are locked up in\nstars. The gaseous metals are distributed over a very wide range of gas\ndensities and temperatures. The shock-heated warm-hot intergalactic\nmedium has a relatively high metallicity of\n{\\textasciitilde}10-1Zsolar that evolves only weakly, and is\ntherefore an important reservoir of metals. Any census aiming to account\nfor most of the metal mass will have to take a wide variety of objects\nand structures into account.},\n\turldate = {2019-05-13},\n\tjournal = {Monthly Notices of the Royal Astronomical Society},\n\tauthor = {Wiersma, Robert P. C. and Schaye, Joop and Theuns, Tom and Dalla Vecchia, Claudio and Tornatore, Luca},\n\tmonth = oct,\n\tyear = {2009},\n\tkeywords = {cosmology: theory, galaxies: abundances, galaxies: formation, intergalactic medium, methods: numerical},\n\tpages = {574--600},\n}\n\n","author_short":["Wiersma, R. P. C.","Schaye, J.","Theuns, T.","Dalla Vecchia, C.","Tornatore, L."],"key":"wiersma_chemical_2009","id":"wiersma_chemical_2009","bibbaseid":"wiersma-schaye-theuns-dallavecchia-tornatore-chemicalenrichmentincosmologicalsmoothedparticlehydrodynamicssimulations-2009","role":"author","urls":{"Paper":"http://adsabs.harvard.edu/abs/2009MNRAS.399..574W"},"keyword":["cosmology: theory","galaxies: abundances","galaxies: formation","intergalactic medium","methods: numerical"],"metadata":{"authorlinks":{}}},"bibtype":"article","biburl":"https://bibbase.org/zotero/polyphant","dataSources":["7gvjSdWrEu7z5vjjj"],"keywords":["cosmology: theory","galaxies: abundances","galaxies: formation","intergalactic medium","methods: numerical"],"search_terms":["chemical","enrichment","cosmological","smoothed","particle","hydrodynamics","simulations","wiersma","schaye","theuns","dalla vecchia","tornatore"],"title":"Chemical enrichment in cosmological, smoothed particle hydrodynamics simulations","year":2009}