-
Code Comparison in Galaxy Scale Simulations with Resolved Supernova Feedback: Lagrangian vs. Eulerian Methods
Authors:
Chia-Yu Hu,
Matthew C. Smith,
Romain Teyssier,
Greg L. Bryan,
Robbert Verbeke,
Andrew Emerick,
Rachel S. Somerville,
Blakesley Burkhart,
Yuan Li,
John C. Forbes,
Tjitske Starkenburg
Abstract:
We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid mo…
▽ More
We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid models, eliminating one of the major uncertainties in cosmological simulations. We find reasonable agreement on the time-averaged star formation rates as well as the joint density-temperature distributions between all codes. However, the Lagrangian codes show significantly burstier star formation, larger supernova-driven bubbles, and stronger galactic outflows compared to the Eulerian code. This is caused by the behavior in the dense, collapsing gas clouds when the Jeans length becomes unresolved: gas in Lagrangian codes collapses to much higher densities than in Eulerian codes, as the latter is stabilized by the minimal cell size. Therefore, more of the gas cloud is converted to stars and SNe are much more clustered in the Lagrangian models, amplifying their dynamical impact. The differences between Lagrangian and Eulerian codes can be reduced by adopting a higher star formation efficiency in Eulerian codes, which significantly enhances SN clustering in the latter. Adopting a zero SN delay time reduces burstiness in all codes, resulting in vanishing outflows as SN clustering is suppressed.
△ Less
Submitted 3 May, 2023; v1 submitted 22 August, 2022;
originally announced August 2022.
-
Spatially Resolved Gas-phase Metallicity in FIRE-2 Dwarfs: Late-Time Evolution of Metallicity Relations in Simulations with Feedback and Mergers
Authors:
Lori E. Porter,
Matthew E. Orr,
Blakesley Burkhart,
Andrew Wetzel,
Xiangcheng Ma,
Philip F. Hopkins,
Andrew Emerick
Abstract:
We present an analysis of spatially resolved gas-phase metallicity relations in five dwarf galaxies ($M_{halo} \approx 10^{11} M_\odot$, $M_\star \approx 10^{8.8}-10^{9.6} M_\odot$) from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, which include an explicit model for sub-grid turbulent mixing of metals in gas, near $z\approx 0$, over a period of 1.4 Gyrs,…
▽ More
We present an analysis of spatially resolved gas-phase metallicity relations in five dwarf galaxies ($M_{halo} \approx 10^{11} M_\odot$, $M_\star \approx 10^{8.8}-10^{9.6} M_\odot$) from the FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulation suite, which include an explicit model for sub-grid turbulent mixing of metals in gas, near $z\approx 0$, over a period of 1.4 Gyrs, and compare our findings with observations. While these dwarf galaxies represent a diverse sample, we find that all simulated galaxies match the observed mass-metallicity (MZR) and mass-metallicity gradient (MZGR) relations. We note that in all five galaxies, the metallicities are effectively identical between phases of the interstellar medium (ISM), with 95$\%$ being within $\pm$0.1 dex between various ISM phases, including the cold and dense gas ($T < 500$ K and $n_{\rm H} > 1$ cm$^{-3}$), ionized gas (near the H$α$ $T \approx 10^4$ K ridge-line), and nebular regions (ionized gas where the 10 Myr-averaged star formation rate is non-zero). We find that most of the scatter in relative metallicity between cold and dense gas and ionized gas/nebular regions can be attributed to either local starburst events or metal-poor inflows. We also note the presence of a major merger in one of our galaxies, m11e, with a substantial impact on the metallicity distribution in the spatially resolved map, showing two strong metallicity peaks and triggering a starburst in the main galaxy.
△ Less
Submitted 13 April, 2022;
originally announced April 2022.
-
Characterizing the Circumgalactic Medium of the Lowest-Mass Galaxies: A Case Study of IC 1613
Authors:
Yong Zheng,
Andrew Emerick,
Mary E. Putman,
Jessica K. Werk,
Evan N. Kirby,
Joshua E. G. Peek
Abstract:
Using 10 sightlines observed with the Hubble Space Telescope/Cosmic Origins Spectrograph, we study the circumgalactic medium (CGM) and outflows of IC1613, which is a low-mass ($M_*\sim10^8~M_\odot$), dwarf irregular galaxy on the outskirts of the Local Group. Among the sightlines, 4 are pointed towards UV-bright stars in IC1613, and the other 6 sightlines are background QSOs at impact parameters f…
▽ More
Using 10 sightlines observed with the Hubble Space Telescope/Cosmic Origins Spectrograph, we study the circumgalactic medium (CGM) and outflows of IC1613, which is a low-mass ($M_*\sim10^8~M_\odot$), dwarf irregular galaxy on the outskirts of the Local Group. Among the sightlines, 4 are pointed towards UV-bright stars in IC1613, and the other 6 sightlines are background QSOs at impact parameters from 6 kpc ($<0.1R_{200}$) to 61 kpc ($0.6R_{200}$). We detect a number of Si II, Si III, Si IV, C II, and C IV absorbers, most of which have velocities less than the escape velocity of IC1613 and thus are gravitationally bound. The line strengths of these ion absorbers are consistent with the CGM absorbers detected in dwarf galaxies at low redshifts. Assuming that Si II, Si III, and Si IV comprise nearly 100% of the total silicon, we find 3% ($\sim$8$\times$10$^3~{\rm M_\odot}$), 2% ($\sim$7$\times$10$^3~{\rm M_\odot}$), and 32--42% [$\sim$(1.0--1.3)$\times$10$^5~{\rm M_\odot}$] of the silicon mass in the stars, interstellar medium, and within $0.6R_{200}$ of the CGM of IC1613. We also estimate the metal outflow rate to be ${\rm \dot{M}_{out, Z}\geq1.1\times10^{-5}~M_\odot~yr^{-1}}$ and the instantaneous metal mass loading factor to be $η_{\rm Z}\geq0.004$, which are in broad agreement with available observation and simulation values. This work is the first time a dwarf galaxy of such low mass is probed by a number of both QSO and stellar sightlines, and it shows that the CGM of low-mass gas-rich galaxies can be a large reservoir enriched with metals from past and ongoing outflows.
△ Less
Submitted 29 October, 2020;
originally announced October 2020.
-
Progenitor-mass-dependent yields amplify intrinsic scatter in dwarf-galaxy elemental abundance ratios
Authors:
Dhruv A. Muley,
Coral R. Wheeler,
Philip F. Hopkins,
Andrew Wetzel,
Andrew Emerick,
Dusan Keres
Abstract:
In hydrodynamic simulations, prevailing subgrid chemical-evolution models often use a single, "IMF-averaged" supernova yield, ignoring variations in elemental abundance ratios (particularly [$α$/Fe]) in the ejecta of higher- and lower-mass supernova progenitors within a stellar population. To understand the impact of this simplification and understand the impact of more explicit models, we run FIR…
▽ More
In hydrodynamic simulations, prevailing subgrid chemical-evolution models often use a single, "IMF-averaged" supernova yield, ignoring variations in elemental abundance ratios (particularly [$α$/Fe]) in the ejecta of higher- and lower-mass supernova progenitors within a stellar population. To understand the impact of this simplification and understand the impact of more explicit models, we run FIRE simulations of a dwarf galaxy $(M_\star($z = 0$) \sim 10^6 M_\odot)$ using nucleosynthetic yields from the NuGrid database that depend on the stellar progenitor mass and metallicity. While NuGrid exhibits lower aggregate $α$-element production than default-FIRE yields, we find that its explicit mass dependence substantially widens the intrinsic scatter in the simulated [Fe/H]-[$α$/Fe] -- a phenomenon potentially visible in recent observations of dwarf galaxies.
△ Less
Submitted 8 September, 2021; v1 submitted 11 August, 2020;
originally announced August 2020.
-
The Role of Stellar Feedback in the Chemical Evolution of a Low Mass Dwarf Galaxy
Authors:
Andrew Emerick,
Greg L. Bryan,
Mordecai-Mark Mac Low
Abstract:
We investigate how each aspect of a multi-channel stellar feedback model drives the chemodynamical evolution of a low-mass, isolated dwarf galaxy using a suite of high-resolution simulations. Our model follows individual star particles sampled randomly from an adopted initial mass function, considering independently feedback from: supernovae; stellar radiation causing photoelectric heating of dust…
▽ More
We investigate how each aspect of a multi-channel stellar feedback model drives the chemodynamical evolution of a low-mass, isolated dwarf galaxy using a suite of high-resolution simulations. Our model follows individual star particles sampled randomly from an adopted initial mass function, considering independently feedback from: supernovae; stellar radiation causing photoelectric heating of dust grains, ionization and associated heating, Lyman-Werner (LW) dissociation of H$_2$, and radiation pressure; and winds from massive main sequence (neglecting their energy input) and asymptotic giant branch (AGB) stars. Radiative transfer is done by ray tracing. We consider the effects each of these processes have on regulating the star formation rate, global properties, multi-phase interstellar medium (ISM), and driving of galactic winds. We follow individual metal species from distinct nucleosynthetic enrichment channels (AGB winds, massive star stellar winds, core collapse and Type Ia supernovae) and pay particular attention to how these feedback processes regulate metal mixing in the ISM, the metal content of outflows, and the stellar abundance patterns in our galaxy. We find that---for a low-metallicity, low-mass dwarf galaxy ---stellar radiation, particularly ionizing radiation and LW radiation, are important sources of stellar feedback whose effects dominate over photoelectric heating and HI radiation pressure. However, feedback is coupled non-linearly, and the inclusion or exclusion of each process produces non-negligible effects. We find strong variations with: the star formation history; the ejection fractions of metals, mass, and energy; and the distribution of elements from different nucleosynthetic sources in both the gas and stars.
△ Less
Submitted 7 July, 2020;
originally announced July 2020.
-
Exploring Chemical Homogeneity in Dwarf Galaxies: A VLT-MUSE study of JKB 18
Authors:
Bethan L. James,
Nimisha Kumari,
Andrew Emerick,
Sergey E. Koposov,
Kristen B. W. McQuinn,
Daniel P. Stark,
Vasily Belokurov,
Roberto Maiolino
Abstract:
Deciphering the distribution of metals throughout galaxies is fundamental in our understanding of galaxy evolution. Nearby, low-metallicity, star-forming dwarf galaxies in particular can offer detailed insight into the metal-dependent processes that may have occurred within galaxies in the early Universe. Here we present VLT/MUSE observations of one such system, JKB 18, a blue diffuse dwarf galaxy…
▽ More
Deciphering the distribution of metals throughout galaxies is fundamental in our understanding of galaxy evolution. Nearby, low-metallicity, star-forming dwarf galaxies in particular can offer detailed insight into the metal-dependent processes that may have occurred within galaxies in the early Universe. Here we present VLT/MUSE observations of one such system, JKB 18, a blue diffuse dwarf galaxy with a metallicity of only 12+log(O/H)=7.6$\pm$0.2 (~0.08 Zsol). Using high spatial-resolution integral field spectroscopy of the entire system, we calculate chemical abundances for individual HII regions using the direct method and derive oxygen abundance maps using strong-line metallicity diagnostics. With large-scale dispersions in O/H, N/H and N/O of ~0.5-0.6 dex and regions harbouring chemical abundances outside this 1$σ$ distribution, we deem JKB 18 to be chemically inhomogeneous. We explore this finding in the context of other chemically inhomogeneous dwarf galaxies and conclude that neither the accretion of metal-poor gas, short mixing timescales, or self-enrichment from Wolf-Rayet stars are accountable. Using a galaxy-scale, multi-phase, hydrodynamical simulation of a low-mass dwarf galaxy, we find that chemical inhomogeneities of this level may be attributable to the removal of gas via supernovae and the specific timing of the observations with respect to star-formation activity. This study not only draws attention to the fact that dwarf galaxies can be chemically inhomogeneous, but also that the methods used in the assessment of this characteristic can be subject to bias.
△ Less
Submitted 18 June, 2020;
originally announced June 2020.
-
Tentative Detection of the Circumgalactic Medium of the Isolated Low-Mass Dwarf Galaxy WLM
Authors:
Y. Zheng,
M. E. Putman,
A. Emerick,
K. B. W. McQuinn,
J. K. Werk,
F. J. Lockman,
B. D. Oppenheimer,
A. J. Fox,
E. N. Kirby,
J. N. Burchett
Abstract:
We report a tentative detection of the circumgalactic medium (CGM) of WLM, an isolated, low-mass (log$M_*/M_\odot\approx7.6$), dwarf irregular galaxy in the Local Group (LG). We analyze an HST/COS archival spectrum of a quasar sightline (PHL2525), which is 45 kpc (0.5 virial radius) from WLM and close to the Magellanic Stream (MS). Along this sightline, two ion absorbers are detected in Si II, Si…
▽ More
We report a tentative detection of the circumgalactic medium (CGM) of WLM, an isolated, low-mass (log$M_*/M_\odot\approx7.6$), dwarf irregular galaxy in the Local Group (LG). We analyze an HST/COS archival spectrum of a quasar sightline (PHL2525), which is 45 kpc (0.5 virial radius) from WLM and close to the Magellanic Stream (MS). Along this sightline, two ion absorbers are detected in Si II, Si III, Si IV, C II, and C IV at velocities of $\sim$-220 km s$^{-1}$ (Component v-220) and $\sim$-150 km s$^{-1}$ (Component v-150). To identify their origins, we study the position-velocity alignment of the components with WLM and the nearby MS. Near the Magellanic longitude of PHL2525, the MS-related neutral and ionized gas moves at $\lesssim-190$ km s$^{-1}$, suggesting an MS origin for Component v-220, but not for Component v-150. Because PHL2525 passes near WLM and Component v-150 is close to WLM's systemic velocity ($\sim$-132 km s$^{-1}$), it is likely that Component v-150 arises from the galaxy's CGM. This results in a total Si mass in WLM's CGM of $M_{\rm Si}^{\rm CGM}\sim(0.2-1.0)\times10^5~M_\odot$ using assumption from other COS dwarf studies. Comparing $M_{\rm Si}^{\rm CGM}$ to the total Si mass synthesized in WLM over its lifetime ($\sim$1.3$\times10^5~M_\odot$), we find $\sim$3% is locked in stars, $\sim$6% in the ISM, $\sim$15%-77% in the CGM, and the rest ($\sim$14%-76%) is likely lost beyond the virial radius. Our finding resonates with other COS dwarf galaxy studies and theoretical predictions that low-mass galaxies can easily lose metals into their CGM due to stellar feedback and shallow gravitational potential.
△ Less
Submitted 11 September, 2019;
originally announced September 2019.
-
Simulating Metal Mixing of Both Common and Rare Enrichment Sources in a Low Mass Dwarf Galaxy
Authors:
Andrew Emerick,
Greg L. Bryan,
Mordecai-Mark Mac Low
Abstract:
One-zone models constructed to match observed stellar abundance patterns have been used extensively to constrain the sites of nucleosynthesis with sophisticated libraries of stellar evolution and stellar yields. The metal mixing included in these models is usually highly simplified, although it is likely to be a significant driver of abundance evolution. In this work we use high-resolution hydrody…
▽ More
One-zone models constructed to match observed stellar abundance patterns have been used extensively to constrain the sites of nucleosynthesis with sophisticated libraries of stellar evolution and stellar yields. The metal mixing included in these models is usually highly simplified, although it is likely to be a significant driver of abundance evolution. In this work we use high-resolution hydrodynamics simulations to investigate how metals from individual enrichment events with varying source energies $E_{\rm ej}$ mix throughout the multi-phase interstellar medium (ISM) of a low-mass ($M_{\rm gas}=2\times 10^{6}$~M$_{\odot}$), low-metallicity, isolated dwarf galaxy. These events correspond to the characteristic energies of both common and exotic astrophysical sites of nucleosynthesis, including: asymptotic giant branch winds ($E_{\rm ej}\sim$10$^{46}$~erg), neutron star-neutron star mergers ($E_{\rm ej}\sim$10$^{49}$~erg), supernovae ($E_{\rm ej}\sim$10$^{51}$~erg), and hypernovae ($E_{\rm ej}\sim$10$^{52}$~erg). We find the mixing timescales for individual enrichment sources in our dwarf galaxy to be long (100~Myr--1~Gyr), with a clear trend of increasing homogeneity for the more energetic events. Given these timescales, we conclude that the spatial distribution and frequency of events are important drivers of abundance homogeneity on large scales; rare, low $E_{\rm ej}$ events should be characterized by particularly broad abundance distributions. The source energy $E_{\rm ej}$ also correlates with the fraction of metals ejected in galactic winds, ranging anywhere from 60\% at the lowest energy to 95\% for hypernovae. We conclude by examining how the radial position, local ISM density, and global star formation rate influence these results.
△ Less
Submitted 3 January, 2020; v1 submitted 10 September, 2019;
originally announced September 2019.
-
Metal Mixing and Ejection in Dwarf Galaxies is Dependent on Nucleosynthetic Source
Authors:
Andrew Emerick,
Greg L. Bryan,
Mordecai-Mark Mac Low,
Benoit Côté,
Kathryn V. Johnston,
Brian W. O'Shea
Abstract:
Using a high resolution simulation of an isolated dwarf galaxy, accounting for multi-channel stellar feedback and chemical evolution on a star-by-star basis, we investigate how each of 15 metal species are distributed within our multi-phase interstellar medium (ISM) and ejected from our galaxy by galactic winds. For the first time, we demonstrate that the mass fraction probability distribution fun…
▽ More
Using a high resolution simulation of an isolated dwarf galaxy, accounting for multi-channel stellar feedback and chemical evolution on a star-by-star basis, we investigate how each of 15 metal species are distributed within our multi-phase interstellar medium (ISM) and ejected from our galaxy by galactic winds. For the first time, we demonstrate that the mass fraction probability distribution functions (PDFs) of individual metal species in the ISM are well described by a piecewise log-normal and power-law distribution. The PDF properties vary within each ISM phase. Hot gas is dominated by recent enrichment, with a significant power-law tail to high metal fractions, while cold gas is predominately log-normal. In addition, elements dominated by asymptotic giant branch (AGB) wind enrichment (e.g. N and Ba) mix less efficiently than elements dominated by supernova enrichment (e.g. $α$ elements and Fe). This result is driven by the differences in source energetics and source locations, particularly the higher chance compared to massive stars for AGB stars to eject material into cold gas. Nearly all of the produced metals are ejected from the galaxy (only 4% are retained), but over 20% of metals dominated by AGB enrichment are retained. In dwarf galaxies, therefore, elements synthesized predominately through AGB winds should be both overabundant and have a larger spread compared to elements synthesized in either core collapse or Type Ia supernovae. We discuss the observational implications of these results, their potential use in developing improved models of galactic chemical evolution, and their generalization to more massive galaxies.
△ Less
Submitted 7 November, 2018; v1 submitted 4 September, 2018;
originally announced September 2018.
-
Stellar Radiation is Critical for Regulating Star Formation and Driving Outflows in Low Mass Dwarf Galaxies
Authors:
Andrew Emerick,
Greg L. Bryan,
Mordecai-Mark Mac Low
Abstract:
Effective stellar feedback is used in models of galaxy formation to drive realistic galaxy evolution. Models typically include energy injection from supernovae as the dominant form of stellar feedback, often in some form of sub-grid recipe. However, it has been recently suggested that pre-SN feedback (stellar winds or radiation) is necessary in high-resolution simulations of galaxy evolution to pr…
▽ More
Effective stellar feedback is used in models of galaxy formation to drive realistic galaxy evolution. Models typically include energy injection from supernovae as the dominant form of stellar feedback, often in some form of sub-grid recipe. However, it has been recently suggested that pre-SN feedback (stellar winds or radiation) is necessary in high-resolution simulations of galaxy evolution to properly regulate star formation and properties of the interstellar medium (ISM). Following these processes is computationally challenging, so many prescriptions model this feedback approximately, accounting for the local destruction of dense gas clouds around newly formed stars in lieu of a full radiative transfer calculation. In this work we examine high resolution simulations (1.8~pc) of an isolated dwarf galaxy with detailed stellar feedback tracked on a star-by-star basis. By following stellar ionizing radiation with an adaptive ray-tracing radiative transfer method, we test its importance in regulating star formation and driving outflows in this galaxy. We find that including ionizing radiation reduces the star formation rate (SFR) by over a factor of 5, and is necessary to produce the ISM conditions needed for supernovae to drive significant outflows. We find that a localized approximation for radiation feedback is sufficient to regulate the SFR on short timescales, but does not allow significant outflows. Short and long range radiation effects are both important in driving the evolution of our low-metallicity, low-mass dwarf galaxy. Generalizing these results to more massive galaxies would be a valuable avenue of future research.
△ Less
Submitted 19 September, 2018; v1 submitted 1 August, 2018;
originally announced August 2018.
-
Simulating an Isolated Dwarf Galaxy with Multi-Channel Feedback and Chemical Yields from Individual Stars
Authors:
Andrew Emerick,
Greg L. Bryan,
Mordecai-Mark Mac Low
Abstract:
In order to better understand the relationship between feedback and galactic chemical evolution, we have developed a new model for stellar feedback at grid resolutions of only a few parsecs in global disk simulations, using the adaptive mesh refinement hydrodynamics code Enzo. For the first time in galaxy scale simulations, we simulate detailed stellar feedback from individual stars including asym…
▽ More
In order to better understand the relationship between feedback and galactic chemical evolution, we have developed a new model for stellar feedback at grid resolutions of only a few parsecs in global disk simulations, using the adaptive mesh refinement hydrodynamics code Enzo. For the first time in galaxy scale simulations, we simulate detailed stellar feedback from individual stars including asymptotic giant branch winds, photoelectric heating, Lyman-Werner radiation, ionizing radiation tracked through an adaptive ray-tracing radiative transfer method, and core collapse and Type Ia supernovae. We furthermore follow the star-by-star chemical yields using tracer fields for 15 metal species: C, N, O, Na, Mg, Si, S, Ca, Mn, Fe, Ni, As, Sr, Y, and Ba. We include the yields ejected in massive stellar winds, but greatly reduce the winds' velocities due to computational constraints. We describe these methods in detail in this work and present the first results from 500~Myr of evolution of an isolated dwarf galaxy with properties similar to a Local Group, low-mass dwarf galaxy. We demonstrate that our physics and feedback model is capable of producing a dwarf galaxy whose evolution is consistent with observations in both the Kennicutt-Schmidt relationship and extended Schmidt relationship. Effective feedback drives outflows with a greater metallicity than the ISM, leading to low metal retention fractions consistent with observations. Finally, we demonstrate that these simulations yield valuable information on the variation in mixing behavior of individual metal species within the multi-phase interstellar medium.
△ Less
Submitted 1 October, 2018; v1 submitted 18 July, 2018;
originally announced July 2018.
-
Grackle: a Chemistry and Cooling Library for Astrophysics
Authors:
Britton D. Smith,
Greg L. Bryan,
Simon C. O. Glover,
Nathan J. Goldbaum,
Matthew J. Turk,
John Regan,
John H. Wise,
Hsi-Yu Schive,
Tom Abel,
Andrew Emerick,
Brian W. O'Shea,
Peter Anninos,
Cameron B. Hummels,
Sadegh Khochfar
Abstract:
We present the Grackle chemistry and cooling library for astrophysical simulations and models. Grackle provides a treatment of non-equilibrium primordial chemistry and cooling for H, D, and He species, including H2 formation on dust grains; tabulated primordial and metal cooling; multiple UV background models; and support for radiation transfer and arbitrary heat sources. The library has an easily…
▽ More
We present the Grackle chemistry and cooling library for astrophysical simulations and models. Grackle provides a treatment of non-equilibrium primordial chemistry and cooling for H, D, and He species, including H2 formation on dust grains; tabulated primordial and metal cooling; multiple UV background models; and support for radiation transfer and arbitrary heat sources. The library has an easily implementable interface for simulation codes written in C, C++, and Fortran as well as a Python interface with added convenience functions for semi-analytical models. As an open-source project, Grackle provides a community resource for accessing and disseminating astrochemical data and numerical methods. We present the full details of the core functionality, the simulation and Python interfaces, testing infrastructure, performance, and range of applicability. Grackle is a fully open-source project and new contributions are welcome.
△ Less
Submitted 14 December, 2016; v1 submitted 29 October, 2016;
originally announced October 2016.
-
Gas Loss by Ram Pressure Stripping and Internal Feedback From Low Mass Milky Way Satellites
Authors:
Andrew Emerick,
Mordecai-Mark Mac Low,
Jana Grcevich,
Andrea Gatto
Abstract:
The evolution of dwarf satellites of the Milky Way is affected by the combination of ram pressure and tidal stripping, and internal feedback from massive stars. We investigate gas loss processes in the smallest satellites of the Milky Way using three-dimensional, high resolution, idealized wind tunnel simulations, accounting for gas loss through both ram pressure stripping and expulsion by superno…
▽ More
The evolution of dwarf satellites of the Milky Way is affected by the combination of ram pressure and tidal stripping, and internal feedback from massive stars. We investigate gas loss processes in the smallest satellites of the Milky Way using three-dimensional, high resolution, idealized wind tunnel simulations, accounting for gas loss through both ram pressure stripping and expulsion by supernova feedback. Using initial conditions appropriate for a dwarf galaxy like Leo T, we investigate whether or not environmental gas stripping and internal feedback can quench these low mass galaxies on the expected timescales, shorter than 2 Gyr. We find that supernova feedback contributes negligibly to the stripping rate for these low star formation rate galaxies. However, we also find that ram pressure stripping is less efficient than expected in the stripping scenarios we consider. Our work suggests that, although ram pressure stripping can eventually completely strip these galaxies, other physics is likely at play to reconcile our computed stripping times with the rapid quenching timescales deduced from observations of low mass Milky Way dwarf galaxies. We discuss the roles additional physics may play in this scenario, including host-satellite tidal interactions, cored vs. cuspy dark matter profiles, reionization, and satellite pre-processing. We conclude that a proper accounting of these physics together is necessary to understand the quenching of low mass Milky Way satellites.
△ Less
Submitted 9 May, 2016;
originally announced May 2016.
-
Warm Gas in and Around Simulated Galaxy Clusters as Probed by Absorption Lines
Authors:
Andrew Emerick,
Greg Bryan,
Mary E. Putman
Abstract:
Understanding gas flows into and out of the most massive dark matter structures in our Universe, galaxy clusters, is fundamental to understanding their evolution. Gas in clusters is well studied in the hot ($>$ 10$^{6}$ K) and cold ($<$ 10$^{4}$ K) regimes, but the warm gas component (10$^{4}$ - 10$^{6}$ K) is poorly constrained. It is challenging to observe directly, but can be probed through Ly…
▽ More
Understanding gas flows into and out of the most massive dark matter structures in our Universe, galaxy clusters, is fundamental to understanding their evolution. Gas in clusters is well studied in the hot ($>$ 10$^{6}$ K) and cold ($<$ 10$^{4}$ K) regimes, but the warm gas component (10$^{4}$ - 10$^{6}$ K) is poorly constrained. It is challenging to observe directly, but can be probed through Ly$α$ absorption studies. We produce the first systematic study of the warm gas content of galaxy clusters through synthetic Ly$α$ absorption studies using cosmological simulations of two galaxy clusters produced with Enzo. We explore the spatial and kinematic properties of our cluster absorbers, and show that the majority of the identified absorbers are due to fast moving gas associated with cluster infall from IGM filaments. Towards the center of the clusters, however, the warm IGM filaments are no longer dominant and the absorbers tend to have higher column densities and metallicities, representing stripped galaxy material. We predict that the absorber velocity distribution should generally be bi-modal and discuss the effects of cluster size, mass, and morphology on the properties of the identified absorbers, and the overall cluster warm gas content. We find tentative evidence for a change in the well known increasing N$_{HI}$ with decreasing impact parameter for the most massive dark matter halos. Our results are compared directly to observations of Ly$α$ absorbers in the Virgo cluster, and provide predictions for future Ly$α$ absorption studies.
△ Less
Submitted 19 August, 2015; v1 submitted 28 July, 2015;
originally announced July 2015.
-
Probing the Off-State of Cluster Giant Radio Halos
Authors:
Shea Brown,
Andrew Emerick,
Lawrence Rudnick,
Gianfranco Brunetti
Abstract:
We derive the best characterization to date of the properties of radio quiescent massive galaxy clusters through a statistical analysis of their average synchrotron emissivity. We stacked 105 radio images of clusters from the 843 MHz SUMSS survey, all with X-ray luminosities greater than 1.0E+44 erg/s and redshifts z < 0.2, after removing point-source contamination and rescaling to a common physic…
▽ More
We derive the best characterization to date of the properties of radio quiescent massive galaxy clusters through a statistical analysis of their average synchrotron emissivity. We stacked 105 radio images of clusters from the 843 MHz SUMSS survey, all with X-ray luminosities greater than 1.0E+44 erg/s and redshifts z < 0.2, after removing point-source contamination and rescaling to a common physical size. Each stacked cluster individually shows no significant large-scale diffuse radio emission at current sensitivity levels. Stacking of sub-samples leads to the following results: (i) clusters with L_{X} > 3.0E+44 erg/s show a 6-sigma detection of Mpc-scale diffuse emission with a 1.4 GHz luminosity of 2.4\pm0.4 x 1.0E+23 W/Hz. This is 1.5-2 times lower than the upper limits for radio quiescent clusters from the GMRT Radio Halo Survey (Venturi et al. 2008), and is the first independent confirmation of radio halo bi-modality; (ii) clusters with low X-ray concentrations have a mean radio luminosity (2.6\pm0.6 x 1.0E+23 W/Hz) that is at least twice that of high X-ray concentration clusters, and (iii) both of these detections are likely close to the low-level "off-state" of GRHs in most or all luminous X-ray clusters, and not due to contributions from a much smaller subset of "on-state" GRHs following the radio/X-ray luminosity correlation. Upcoming deep radio surveys will conclusively distinguish between these two options. We briefly discuss possible origins for the "off-state" emission and its implications for magnetic fields in most or all luminous X-ray clusters.
△ Less
Submitted 1 September, 2011;
originally announced September 2011.
