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 H 2 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.

We study the production of bottomonium states in heavy-ion reactions at collider energies available at RHIC and LHC. We employ an earlier constructed rate equation approach which accounts for both suppression and regeneration mechanisms in the quark-gluon plasma (QGP) and hadronization phases of the evolving thermal medium. Our previous predictions utilizing two limiting cases of strong and weak bottomonium binding in the QGP are updated by (i) checking the compatibility of the pertinent spectral functions with lattice-QCD results for euclidean correlators, (ii) adapting the initial conditions of the rate equation by updating bottom-related input cross sections and the charged-particle multiplicity of the fireball, and (iii) converting our calculations into observables as recently measured by the STAR and CMS experiments. Our main findings are a preference for strong Υ binding as well as a significant regeneration component at the LHC.

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 E . 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.

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.

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 L X > 10 44 erg s −1 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 × 10 44 erg s −1 show a 6σ detection of Mpc-scale diffuse emission with a 1.4 GHz luminosity of 2.4±0.4 ×10 23 W Hz −1 . This is 1.5-2 times lower than the upper limits for radio quiescent clusters from the GMRT Radio Halo Survey , 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±0.6 ×10 23 W Hz −1 ) 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 the 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. Subject headings: word -word -word

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.

Enzo (Enzo Developers, 2019a) is a block-structured adaptive mesh refinement code that is widely used to simulate astrophysical fluid flows (primarily, but not exclusively, cosmological structure formation, star formation, and turbulence). The code is a community project with dozens of users, and has contributed to hundreds of peer-reviewed publications in astrophysics, physics, and computer science. The code utilizes a Cartesian mesh can be run in one, two, or

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 powerlaw 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.