Search

Precipitation of cold gas due to thermal instability in both galaxy clusters and the circumgalactic medium may regulate AGN feedback. We investigate thermal instability in idealized simulations of the circumgalactic medium with a parameter study of over 600 three-dimensional hydrodynamic simulations of stratified turbulence with cooling, each evolved for 10 Gyr. The entropy profiles are maintained in a steady state via an idealized `thermostat' process, consistent with galaxy cluster entropy profiles. In the presence of external turbulent driving, we find cold gas precipitates, with a strong dependence whether the turbulent driving mechanism is solenoidal, compressive, or purely vertical. In the purely-vertical turbulent driving regime, we find that significant cold gas may form when the cooling time to free-fall time $t_{\rm cool} / t_{\text{ff}} \lesssim 5$. Our simulations with a ratio of $t_{\rm cool} / t_{\text{ff}} \sim 10$ do not precipitate under any circumstances, perhaps because the thermostat mechanism we use maintains a significant non-zero entropy gradient., Comment: 13 pages, 9 figures. Submitted to MNRAS. Simulation movies are available at this URL: https://benwibking.github.io/precipitation.html

We investigate the co-evolution of the angular momentum of Milky Way-like galaxies, their circumgalactic gas, and their dark matter halos using zoom-in simulations from the Figuring Out Gas & Galaxies in Enzo (FOGGIE) suite. We examine how the magnitude and orientation of the angular momentum varies over time within the halo and between the components of mass. From z~2 to today, and in general across the simulated halos, the specific angular momenta of the central galaxies and the cool gas in their circumgalactic media (T < 10^5 K) increase together. Over that same period, the specific angular momenta of the hot (>10^6 K) and dark components of the halo change minimally. By z~1, the central galaxies have generally lost association with the angular momentum of their full dark matter halo -- both in magnitude and orientation. We find a wide distribution of angular momentum orientations in the halo, varying by up to 180 degrees over small (~tens of kpc) scales and between the different components of mass. The net angular momenta of the galaxies, their circumgalactic gas, and their dark matter halos are generally misaligned with one another at all cosmic times. The present-day orientation of the central galaxies are established at late times (after z=1), after the rates of cosmic accretion and mergers decline and the disks are able to settle and stabilize their orientation., Comment: 35 pages, 14 figures, submitted to ApJ

The advent of new and near-future observatories probing the earliest epochs of the Universe has opened the opportunity to investigate the formation and growth of the first massive black holes (MBHs). Additionally, the use of high resolution cosmological simulations to investigate these high-redshift environments is needed to predict the dark matter halos in which these MBH seeds will form. We use the $\textit{Renaissance}$ simulations to analyse the formation and growth of so-called heavy seed black holes. Other past work has investigated the formation and growth of light (black hole) seeds with $\textit{Renaissance}$ and found that these black holes do not grow in the environments in which they reside. In this work we seed MBHs, in post-processing, and track accretion onto the MBHs as well as mergers with other MBHs at high-redshift. We show that the heavy seeds struggle to achieve high accretion rates with only the most massive black holes ($\gtrsim 10^5 \text{M}_\odot$) growing at close to the Eddington limit under optimistic conditions. Despite the lack of significant growth for these early MBHs, the signals from their merger events will be sufficiently strong (SNR $\sim 10^2$) to be probed by the next generation of gravitational wave observatories, such as $\textit{LISA}$. We predict that $\textit{LISA}$ will observe of the order of $10$ MBH merger events per year where the mergers occur at z $\gtrsim$ 10 or at least begin their early inspiral phase at z $\gtrsim$ 10., Comment: 15 pages, 10 figures, Submitted to the Open Journal of Astrophysics for review

The prevalence of light element enhancement in the most metal-poor stars is potentially an indication that the Milky Way has a metallicity floor for star formation around $\sim$10$^{-3.5}$ Z$_{\odot}$. We propose that this metallicity floor has its origins in metal-enriched star formation in the minihalos present during the Galaxy's initial formation. To arrive at this conclusion, we analyze a cosmological radiation hydrodynamics simulation that follows the concurrent evolution of multiple Population III star-forming minihalos. The main driver for the central gas within minihalos is the steady increase in hydrostatic pressure as the halos grow. We incorporate this insight into a hybrid one-zone model that switches between pressure-confined and modified free-fall modes to evolve the gas density with time according to the ratio of the free-fall and sound-crossing timescales. This model is able to accurately reproduce the density and chemo-thermal evolution of the gas in each of the simulated minihalos up to the point of runaway collapse. We then use this model to investigate how the gas responds to the absence of H$_{2}$. Without metals, the central gas becomes increasingly stable against collapse as it grows to the atomic cooling limit. When metals are present in the halo at a level of $\sim$10$^{-3.7}$ Z$_{\odot}$, however, the gas is able to achieve gravitational instability while still in the minihalo regime. Thus, we conclude that the Galaxy's metallicity floor is set by the balance within minihalos of gas-phase metal cooling and the radiation background associated with its early formation environment., Comment: 11 pages, 5 figures, accepted for publication in MNRAS

Filaments of cold gas ($T\leq 10^{4}$ K) are found in the inner regions of many cool-core clusters. These structures are thought to play a major role in the regulation of feedback from active galactic nuclei (AGN). We study the morphology of the filaments, their formation, and their impact on the propagation of the outflowing AGN jets. We present a set of GPU-accelerated 3D (magneto)hydrodynamical simulations of an idealized Perseus-like cluster using the performance portable code AthenaPK. We include radiative cooling, and a self-regulated AGN feedback model that redistributes accreted material through kinetic, thermal and magnetic feedback. We confirm that magnetic fields play an important role in both the formation and evolution of the cold material. These suppress the formation of massive cold disks and favor magnetically supported filaments over clumpy structures. Achieving resolutions of $25-50$ pc, we find that filaments are not monolithic as they contain numerous and complex magnetically--supported substructures. We find that the mass distribution of these clumps follows a power law for all investigated filaments, consistent with previous cloud-crushing simulations of individual clumps. Studying the evolution of individual filaments, we find that their formation pathways can be diverse. We find examples of filaments forming through a combination of cold gas uplifting and condensation, as well as systems of purely infalling clumps condensing out of the intracluster medium. The density contrast between the cold gas and the outflowing hot material leads to recurring deflections of the jets, favoring inflation of bubbles., Comment: 19 pages, 19 figures, submitted to Astronomy & Astrophysics

Data-driven methods of model identification are able to discern governing dynamics of a system from data. Such methods are well suited to help us learn about systems with unpredictable evolution or systems with ambiguous governing dynamics given our current understanding. Many plasma problems of interest fall into these categories as there are a wide range of models that exist, however each model is only useful in a certain regime and often limited by computational complexity. To ensure data-driven methods align with theory, they must be consistent and predictable when acting on data whose governing dynamics are known. Weak Sparse Identification of Nonlinear Dynamics (WSINDy) is a recently developed data-driven method that has shown promise in learning governing dynamics from data with high noise levels [1]. This work examines how WSINDy acts on ideal MHD test problems as the initial conditions are varied and specifies limiting requirements for successful equation identification. It is hard to recover the governing dynamics from data that emphasize a single dominant behavior. In these low information cases, Shannon information entropy is able to pick up on the redundancies in the data that affect recoverability., Comment: 25 pages, 15 figures. Submitted to Journal of Computational Physics. Feedback welcome. Source code and data available upon request

Radiative cooling and AGN heating are thought to form a feedback loop that regulates the evolution of low redshift cool-core galaxy clusters. Numerical simulations suggest that formation of multiphase gas in the cluster core imposes a floor on the ratio of cooling time ($t_{\rm cool}$) to free-fall time ($t_{\rm ff}$) at $\min ( t_{\rm cool} / t_{\rm ff} ) \approx 10$. Observations of galaxy clusters show evidence for such a floor, and usually the cluster cores with $\min ( t_{\rm cool} / t_{\rm ff} ) \lesssim 30$ contain abundant multiphase gas. However, there are important outliers. One of them is Abell 2029, a massive galaxy cluster ($M_{200} \gtrsim 10^{15}$ M$_\odot$) with $\min( t_{\rm cool}/t_{\rm ff}) \sim 20$, but little apparent multiphase gas. In this paper, we present high resolution 3D hydrodynamic AMR simulations of a cluster similar to A2029 and study how it evolves over a period of 1-2 Gyr. Those simulations suggest that Abell 2029 self-regulates without producing multiphase gas because the mass of its central black hole ($\sim 5\times 10^{10} \, M_\odot$) is great enough for Bondi accretion of hot ambient gas to produce enough feedback energy to compensate for radiative cooling., Comment: 8 pages, 5 figures, submitted to MNRAS

Over the next decade, the astronomical community will be commissioning multiple wide-field observatories well-suited for studying stellar halos in both integrated light and resolved stars. In preparation for this, we use five high-resolution cosmological simulations of Milky Way-like galaxies from the FOGGIE suite to explore the properties and components of stellar halos. These simulations are run with high time (5 Myr) and stellar mass (1000 M$_\odot$) resolution to better model the properties and origins of low density regions like stellar halos. We find that the FOGGIE stellar halos have masses, metallicity gradients, and surface brightness profiles that are consistent with observations. In agreement with other simulations, the FOGGIE stellar halos receive 30-40% of their mass from in situ stars. However, this population is more centrally concentrated in the FOGGIE simulations and therefore does not contribute excess light to the halo outskirts. The remaining stars are accreted from 10-50 other galaxies, with the majority of the accreted mass originating in 2-4 galaxies. While the inner halo ($r<50$ kpc) of each FOGGIE galaxy has a large number of contributors, the halo outskirts of three of the five galaxies are primarily made up of stars from only a few contributors. We predict that upcoming wide-field observatories, like the Nancy Grace Roman Space Telescope, will probe stellar halos around Milky Way-like galaxies out to ~100 kpc in integrated light and will be able to distinguish the debris of dwarf galaxies with extended star formation histories from the underlying halo with resolved color-magnitude diagrams., Comment: 33 pages, 15 figures, accepted to ApJ

We present a study of the co-evolution of a population of primordial star-forming minihalos at Cosmic Dawn. In this study, we highlight the influence of individual Population III stars on the ability of nearby minihalos to form sufficient molecular hydrogen to undergo star formation. In the absence of radiation, we find the minimum halo mass required to bring about collapse to be ~10^5 Msun, this increases to ~10^6 Msun after two stars have formed. We find an inverse relationship between halo mass and the time required for it to recover its molecular gas after being disrupted by radiation from a nearby star. We also take advantage of the extremely high resolution to investigate the effects of major and minor mergers on the gas content of star-forming minihalos. Contrary to previous claims of fallback of supernova ejecta, we find minihalos evacuated after hosting Pop III stars primarily recover gas through mergers with undisturbed halos. We identify an intriguing type of major merger between recently evacuated halos and gas-rich ones, finding that these 'mixed' mergers accelerate star formation instead of suppressing it like their low redshift counterparts. We attribute this to the gas-poor nature of one of the merging halos resulting in no significant rise in temperature or turbulence and instead inducing a rapid increase in central density and hydrostatic pressure. This constitutes a novel formation pathway for Pop III stars and establishes major mergers as potentially the primary source of gas, thus redefining the role of major mergers at this epoch., Comment: 15 pages, 12 figures, accepted for publication in MNRAS, revised with moderate changes

Magnetized turbulence is ubiquitous in many astrophysical and terrestrial plasmas but no universal theory exists. Even the detailed energy dynamics in magnetohydrodynamic (MHD) turbulence are still not well understood. We present a suite of subsonic, super-Alfv\'enic, high plasma-beta MHD turbulence simulations that only vary in their dynamical range, i.e., in their separation between the large-scale forcing and dissipation scales, and their dissipation mechanism (implicit large eddy simulation, ILES, versus and direct numerical simulation, DNS). Using an energy transfer analysis framework we calculate the effective, numerical viscosities and resistivities and demonstrate and that all ILES calculations of MHD turbulence are resolved and correspond to an equivalent visco-resistive MHD turbulence calculation. Increasing the number of grid points used in an ILES corresponds to lowering the dissipation coefficients, i.e., larger (kinetic and magnetic) Reynolds numbers for a constant forcing scale. Independently, we use this same framework to demonstrate that -- contrary to hydrodynamic turbulence -- the cross-scale energy fluxes are not constant in MHD turbulence. This applies both to different mediators (such as cascade processes or magnetic tension) for a given dynamical range as well as to a dependence on the dynamical range itself, which determines the physical properties of the flow. We do not observe any indication of convergence even at the highest resolution (largest Reynolds numbers) simulation at $2{,}048^3$ cells, calling into question whether an asymptotic regime in MHD turbulence exists, and, if so, what it looks like., Comment: under review, comments welcome

The circumgalactic medium (CGM) is often assumed to exist in or near hydrostatic equilibrium with the regulation of accretion and the effects of feedback treated as perturbations to a stable balance between gravity and thermal pressure. We investigate global hydrostatic equilibrium in the CGM using four highly-resolved $L^*$ galaxies from the Figuring Out Gas & Galaxies In Enzo (FOGGIE) project. The FOGGIE simulations were specifically targeted at fine spatial and mass resolution in the CGM ($\Delta x \lesssim 1$ kpc $h^{-1}$ and $M \simeq 200M_\odot$). We develop a new analysis framework that calculates the forces provided by thermal pressure gradients, turbulent pressure gradients, ram pressure gradients of large-scale radial bulk flows, centrifugal rotation, and gravity acting on the gas in the CGM. Thermal and turbulent pressure gradients vary strongly on scales of $\lesssim5$ kpc throughout the CGM. Thermal pressure gradients provide the main supporting force only beyond $\sim 0.25R_{200}$, or $\sim50$ kpc at $z=0$. Within $\sim0.25R_{200}$, turbulent pressure gradients and rotational support provide stronger forces than thermal pressure. More generally, we find that global equilibrium models are neither appropriate nor predictive for the small scales probed by absorption line observations of the CGM. Local conditions generally cannot be derived by assuming a global equilibrium, but an emergent global equilibrium balancing radially inward and outward forces is obtained when averaging over the non-equilibrium local conditions on large scales in space and time. Approximate hydrostatic equilibrium holds only at large distances from galaxies even when averaging out small-scale variations., Comment: 36 pages, 24 figures. Accepted to ApJ. See https://foggie.science/cgm_forces.html for animated Figures 6 and 7

The Milky Way has accreted many ultra-faint dwarf galaxies (UFDs), and stars from these galaxies can be found throughout our Galaxy today. Studying these stars provides insight into galaxy formation and early chemical enrichment, but identifying them is difficult. Clustering stellar dynamics in 4D phase space ($E$, $L_z$, $J_r$, $J_z$) is one method of identifying accreted structure which is currently being utilized in the search for accreted UFDs. We produce 32 simulated stellar halos using particle tagging with the \textit{Caterpillar} simulation suite and thoroughly test the abilities of different clustering algorithms to recover tidally disrupted UFD remnants. We perform over 10,000 clustering runs, testing seven clustering algorithms, roughly twenty hyperparameter choices per algorithm, and six different types of data sets each with up to 32 simulated samples. Of the seven algorithms, HDBSCAN most consistently balances UFD recovery rates and cluster realness rates. We find that even in highly idealized cases, the vast majority of clusters found by clustering algorithms do not correspond to real accreted UFD remnants and we can generally only recover $6\%$ of UFDs remnants at best. These results focus exclusively on groups of stars from UFDs, which have weak dynamic signatures compared to the background of other stars. The recoverable UFD remnants are those that accreted recently, $z_{\text{accretion}}\lesssim 0.5$. Based on these results, we make recommendations to help guide the search for dynamically-linked clusters of UFD stars in observational data. We find that real clusters generally have higher median energy and $J_r$, providing a way to help identify real vs. fake clusters. We also recommend incorporating chemical tagging as a way to improve clustering results., Comment: 21 pages, 14 figures, accepted to ApJ

Precipitation is potentially a mechanism through which the circumgalactic medium (CGM) can regulate a galaxy's star formation. Here we present idealized simulations of isolated Milky Way-like galaxies intended to examine the ability of galaxies to self-regulate their star formation, particularly via precipitation. Our simulations are the first CGM-focused idealized models to include stellar feedback due to the explicit formation of stars. We also examine the impact of rotation in the CGM. Using six simulations, we explore variations in the initial CGM $t_{\rm cool}/t_{\rm ff}$ ratio and rotation profile. Those variations affect the amount of star formation and gas accretion within the galactic disk. Our simulations are sensitive to their initial conditions, requiring us to gradually increase the efficiency of stellar feedback to avoid destroying the CGM before its gas can be accreted. Despite this gradual increase, the resulting outflows still evacuate large, hot cavities within the CGM and even beyond $r_{200}$. Some of the CGM gas avoids interacting with the cavities and is able to feed the disk along its midplane, but the cooling of feedback-heated gas far from the midplane is too slow to supply the disk with additional gas. Our simulations illustrate the importance of physical mechanisms in the outer CGM and IGM for star formation regulation in Milky Way-scale halos., Comment: 28 pages, 13 figures. Submitted to the Astrophysical Journal

We present the KODIAQ-Z survey aimed to characterize the cool, photoionized gas at 2.2100) between absorbers separated by less than 500 km/s imply that the metals are poorly mixed in 14.6

Coupling between active galactic nuclei (AGN) and the circumgalactic medium (CGM) is critical to the interplay between radiative cooling and feedback heating in the atmospheres of the universe's most massive galaxies. This paper presents a detailed analysis of numerical simulations showing how kinetic AGN feedback with a strong momentum flux interacts with the CGM. Our analysis shows that large scale CGM circulation plays an important role in reconfiguring the galactic atmosphere and regulating the atmosphere's central entropy level. We find that most of the AGN energy output goes into lifting of circumgalactic gas rather than heating of atmospheric gas within the galaxy, consequently reconfiguring the circumgalactic medium (CGM) in our simulations. Large scale (10s of kpc) circulation of the CGM on ~ 10-100 kpc scales therefore plays a critical role in preventing over-cooling of gas in these simulated galaxies. The simulations also show that our choices of accretion efficiency and jet opening angle significantly affect the AGN-CGM coupling. Reducing the jet opening angle to quarter of the fiducial opening angle increases the jet momentum flux, enabling it to drill through to larger radii without effectively coupling with the CGM at the center ( $r < 5$ kpc). Outflows with a lower momentum flux decelerate and thermalize the bulk of their energy at smaller radii ($r \lesssim 10$ )., Comment: 26 pages, 17 figures, submitted to ApJ

Identifying stars formed in pristine environments (Pop III) within the first billion years is vital to uncovering the earliest growth and chemical evolution of galaxies. Pop III galaxies, however, are typically expected to be too faint and too few in number to be detectable by forthcoming instruments without extremely long integration times and/or extreme lensing. In an environment, however, where star formation is suppressed until a halo crosses the atomic cooling limit (e.g., by a modest Lyman-Werner flux, high baryonic streaming velocities, and/or dynamical heating effects),primordial halos can form substantially more numerous and more massive stars. Some of these stars will in-turn be accreting more rapidly than they can thermally relax at any given time. Using high resolution cosmological zoom-in simulations of massive star formation in high-z halos, we find that such rapidly accreting stars produce prominent spectral features which would be detectable by {\it JWST}. The rapid accretion episodes within the halo lead to stochastic reprocessing of 0--20\% of the total stellar emission into the rest-frame optical over long timescales, a unique signature which may allow deep observations to identify such objects out to $z \sim 10-13$ using mid- and wide-band NIRCam colors alone., Comment: 8 pages, 4 figures, Astrophysical Journal Letters accepted

The classical definition of the virial temperature of a galaxy halo excludes a fundamental contribution to the energy partition of the halo: the kinetic energy of non-thermal gas motions. Using simulations of low-redshift, $\sim L^*$ galaxies from the FOGGIE project (Figuring Out Gas & Galaxies In Enzo) that are optimized to resolve low-density gas, we show that the kinetic energy of non-thermal motions is roughly equal to the energy of thermal motions. The simulated FOGGIE halos have $\sim 2\times$ lower bulk temperatures than expected from a classical virial equilibrium, owing to significant non-thermal kinetic energy that is formally excluded from the definition of $T_\mathrm{vir}$. We derive a modified virial temperature explicitly including non-thermal gas motions that provides a more accurate description of gas temperatures for simulated halos in virial equilibrium. Strong bursts of stellar feedback drive the simulated FOGGIE halos out of virial equilibrium, but the halo gas cannot be accurately described by the standard virial temperature even when in virial equilibrium. Compared to the standard virial temperature, the cooler modified virial temperature implies other effects on halo gas: (i) the thermal gas pressure is lower, (ii) radiative cooling is more efficient, (iii) O VI absorbing gas that traces the virial temperature may be prevalent in halos of a higher mass than expected, (iv) gas mass estimates from X-ray surface brightness profiles may be incorrect, and (v) turbulent motions make an important contribution to the energy balance of a galaxy halo., Comment: 30 pages, 14 figures, accepted to ApJ

Recent high-resolution simulations of early structure formation have shown that externally enriched halos may form some of the first metal enriched stars. This study utilizes a 1 comoving Mpc$^3$ high-resolution simulation to study the enrichment process of metal-enriched halos down to $z=9.3$. Our simulation uniquely tracks the metals ejected from Population III stars, and we use this information to identify the origin of metals within metal-enriched halos. These halos show a wide range of metallicities, but we find that the source of metals for $\gtrsim$ 50\% of metal-enriched halos is supernova explosions of Population III stars occuring outside their virial radii. The results presented here indicate that external enrichment by metal-free stars dominates the enrichment process of halos with virial mass below $10^{6}\,M_\odot$ down to $z=9.3$. Despite the prevalence of external enrichment in low mass halos, Pop II stars forming due to external enrichment are rare because of the small contribution of low-mass halos to the global star formation rate combined with low metallicities towards the center of these halos resulting from metal ejecta from external sources mixing from the outside-in. The enriched stars that do form through this process have absolute metallicities below $10^{-3}\,Z_\odot$. We also find that the fraction of externally enriched halos increases with time, $\sim 90\%$ of halos that are externally enriched have $M_\mathrm{vir} < 10^6\,M_\odot$, and that pair-instability supernovae contribute the most to the enrichment of the IGM as a whole and are thus are the predominant supernova type contributing to the external enrichment of halos., Comment: 15 pages, 12 figures, Accepted for publication in ApJ

Magnetized turbulence is ubiquitous in many astrophysical and terrestrial systems but no complete, uncontested theory even in the simplest form, magnetohydrodynamics (MHD), exists. Many theories and phenomenologies focus on the joint (kinetic and magnetic) energy fluxes and spectra. We highlight the importance of treating kinetic and magnetic energies separately to shed light on MHD turbulence dynamics. We conduct an implicit large eddy simulation of subsonic, super-Alfv\'enic MHD turbulence and analyze the scale-wise energy transfer over time. Our key finding is that the kinetic energy spectrum develops a scaling of approximately $k^{-4/3}$ in the stationary regime as the kinetic energy cascade is suppressed by magnetic tension. This motivates a reevaluation of existing MHD turbulence theories with respect to a more differentiated modeling of the energy fluxes., Comment: 10 pages, 6 figures, submitted to The Astrophysical Journal

Magnetohydrodynamic turbulence affects both terrestrial and astrophysical plasmas. The properties of magnetized turbulence must be better understood to more accurately characterize these systems. This work presents ideal MHD simulations of the compressible Taylor-Green vortex under a range of initial sub-sonic Mach numbers and magnetic field strengths. We find that regardless of the initial field strength, the magnetic energy becomes dominant over the kinetic energy on all scales after at most several dynamical times. The spectral indices of the kinetic and magnetic energy spectra become shallower than $k^{-5/3}$ over time and generally fluctuate. Using a shell-to-shell energy transfer analysis framework, we find that the magnetic fields facilitate a significant amount of the energy flux and that the kinetic energy cascade is suppressed. Moreover, we observe nonlocal energy transfer from the large scale kinetic energy to intermediate and small scale magnetic energy via magnetic tension. We conclude that even in intermittently or singularly driven weakly magnetized systems, the dynamical effects of magnetic fields cannot be neglected.

We investigate the ab-initio formation of super-massive stars in a pristine atomic cooling halo. The halo is extracted from a larger self-consistent parent simulation. The halo remains metal-free and star formation is suppressed due to a combination of dynamical heating from mergers and a mild ($J_{\rm LW} \sim 2 - 10 \ J_{21}$(z)) Lyman-Werner (LW) background. We find that more than 20 very massive stars form with stellar masses greater than 1000 M$_{\odot}$. The most massive star has a stellar mass of over 6000 M$_{\odot}$. However, accretion onto all stars declines significantly after the first $\sim$ 100 kyr of evolution as the surrounding material is accreted and the turbulent nature of the gas causes the stars to move to lower density regions. We post-process the impact of ionising radiation from the stars and find that ionising radiation is not a limiting factor when considering SMS formation and growth. Rather the birth environments are highly turbulent and a steady accretion flow is not maintained within the timescale (2 Myr) of our simulations. As the massive stars end their lives as direct collapse black holes this will seed these embryonic haloes with a population of black holes with masses between approximately 300 M$_{\odot}$ and 10,000 M$_{\odot}$. Afterwards they may sink to the centre of the haloes, eventually coalescing to form larger intermediate mass black holes whose in-situ mergers will be detectable by LISA., Comment: 13 pages. Accepted for publication in the Open Journal of Astrophysics

The formation of supermassive stars has generally been studied under the assumption of rapid accretion of pristine metal-free gas. Recently it was found, however, that gas enriched to metallicities up to $Z \sim 10^{-3}$ Z$_{\odot}$ can also facilitate supermassive star formation, as long as the total mass infall rate onto the protostar remains sufficiently high. We extend the analysis further by examining how the abundance of supermassive star candidate haloes would be affected if all haloes with super-critical infall rates, regardless of metallicity were included. We investigate this scenario by identifying all atomic cooling haloes in the Renaissance simulations with central mass infall rates exceeding a fixed threshold. We find that among these haloes with central mass infall rates above 0.1 M$_{\odot}$ yr$^{-1}$ approximately two-thirds of these haloes have metallicities of $Z > 10^{-3}$ Z$_{\odot}$. If metal mixing within these haloes is inefficient early in their assembly and pockets of metal-poor gas can remain then the number of haloes hosting supermassive stars can be increased by at least a factor of four. Additionally the centres of these high infall-rate haloes provide ideal environments in which to grow pre-existing black holes. Further research into the (supermassive) star formation dynamics of rapidly collapsing haloes, with inhomogeneous metal distributions, is required to gain more insight into both supermassive star formation in early galaxies as well as early black hole growth., Comment: 11 pages. Accepted for publication in the Open Journal of Astrophysics

In the universe's most massive galaxies, kinetic feedback from a central supermassive black hole appears to limit star formation. Abundant circumstantial evidence suggests that accumulation of cold gas near the central black hole strongly boosts the feedback output, keeping the ambient medium in a state marginally unstable to condensation and formation of cold gas clouds. However, the ability of that mechanism to self-regulate may depend on numerous environmental factors, including the depth of the potential well and the pressure of the surrounding circumgalactic medium (CGM). Here we present a suite of numerical simulations that explores the dependence of cold-fuelled bipolar kinetic feedback on those environmental factors. Halo mass in this simulation suite ranges from $2 \times 10^{12} \, M_\odot$ to $8 \times 10^{14} \, M_\odot$. We include the spatially extended mass and energy input from the massive galaxy's old stellar population, which is capable of sweeping gas out of the galaxy and away from the central black hole if the confining CGM pressure is sufficiently low. Our simulations show that this feedback mechanism is tightly self-regulating in a massive galaxy with a deep central potential and low CGM pressure, permitting only small amounts of multiphase gas to accumulate and allowing almost no star formation. In a massive galaxy of similar mass but a shallower central potential and greater CGM pressure the same feedback mechanism is more episodic, producing extended multiphase gas and occasionally allowing small rates of star formation ($\sim 0.1 \, M_\odot \, {\rm yr}^{-1}$). At the low-mass end of the explored range the mechanism becomes implausibly explosive, perhaps because the ambient gas initially has no angular momentum, which would have reduced the amount of condensed gas capable of fueling feedback., Comment: 17 pages, 9 figures, Accepted for publication in ApJ

We study ram pressure stripping in simulated Milky Way-like halos at z>=2 from the Figuring Out Gas & Galaxies In Enzo (FOGGIE) project. These simulations reach exquisite resolution in their circumgalactic medium (CGM) gas owing to FOGGIE's novel refinement scheme. The CGM of each halo spans a wide dynamic range in density and velocity over its volume---roughly 6 dex and 1000 km/s, respectively---translating into a 5 dex range in ram pressure imparted to interacting satellites. The ram pressure profiles of the simulated CGM are highly stochastic, owing to kpc-scale variations of the density and velocity fields of the CGM gas. As a result, the efficacy of ram pressure stripping depends strongly on the specific path a satellite takes through the CGM. The ram-pressure history of a single satellite is generally unpredictable and not well correlated with its approach vector with respect to the host galaxy. The cumulative impact of ram pressure on the simulated satellites is dominated by only a few short strong impulses---on average, 90% of the total surface momentum gained through ram pressure is imparted in 20% or less of the total orbital time. These results reveal an erratic mode of ram pressure stripping in Milky-Way like halos at high redshift---one that is not captured by a smooth spherically-averaged model of the circumgalactic medium., Comment: 18 pages, 10 figures. Submitted to ApJ

In cool-core galaxy clusters with central cooling times much shorter than a Hubble time, condensation of the ambient central gas is regulated by a heating mechanism, probably an active galactic nucleus (AGN). Previous analytical work has suggested that certain radial distributions of heat input may result in convergence to a quasi-steady global state that does not substantively change on the timescale for radiative cooling, even if the heating and cooling are not locally in balance. To test this hypothesis, we simulate idealized galaxy cluster halos using the \ENZO code with an idealized, spherically symmetric heat-input kernel intended to emulate. Thermal energy is distributed with radius according to a range of kernels, in which total heating is updated to match total cooling every $10 ~\text{Myr}$. Some heating kernels can maintain quasi-steady global configurations, but no kernel we tested produces a quasi-steady state with central entropy as low as those observed in cool-core clusters. The general behavior of the simulations depends on the proportion of heating in the inner $10 ~\text{kpc}$, with low central heating leading to central cooling catastrophes, high central heating creating a central convective zone with an inverted entropy gradient, and intermediate central heating resulting in a flat central entropy profile that exceeds observations. The timescale on which our simulated halos fall into an unsteady multiphase state is proportional to the square of the cooling time of the lowest entropy gas, allowing more centrally concentrated heating to maintain a longer lasting steady state., Comment: 17 pages, 8 figures

The circumgalactic medium (CGM) of the Milky Way is mostly obscured by nearby gas in position-velocity space because we reside inside the Galaxy. Substantial biases exist in most studies on the Milky Way's CGM that focus on easier-to-detect high-velocity gas. With mock observations on a Milky-Way analog from the FOGGIE simulation, we investigate four observational biases related to the Milky Way's CGM. First, QSO absorption-line studies probe a limited amount of the CGM mass: only 35% of the mass is at high Galactic latitudes $|b|>20$ degrees, of which only half is moving at $|v_{\rm LSR}|\gtrsim100$ km s$^{-1}$. Second, the inflow rate ($\dot{M}$) of the cold gas observable in HI 21cm is reduced by a factor of $\sim10$ as we switch from the local standard of rest to the galaxy's rest frame; meanwhile $\dot{M}$ of the cool and warm gas does not change significantly. Third, OVI and NV are promising ions to probe the Milky Way's outer CGM ($r\gtrsim$15 kpc), but CIV may be less sensitive. Lastly, the scatter in ion column density is a factor of 2 higher if the CGM is observed from inside-out than from external views because of the gas radial density profile. Our work highlights that observations of the Milky Way's CGM, especially those using HI 21cm and QSO absorption lines, are highly biased. We demonstrate that these biases can be quantified and calibrated through synthetic observations with simulated Milky-Way analogs., Comment: ApJ in press

Recent observations show that the star formation rate (SFR) in the {\it Phoenix} cluster's central galaxy is $\sim 500$ M$_\odot$ yr$^{-1}$. Even though {\it Phoenix} is a massive cluster ($M_{200} \approx 2.0\times 10^{15}$ M$_\odot$; $z\approx 0.6$) such a high central SFR is not expected in a scenario in which feedback from an active galactic nucleus (AGN) maintains the intracluster medium (ICM) in a state of rough thermal balance. It has been argued that either AGN feedback saturates in very massive clusters or the central supermassive black hole (SMBH) is too small to produce enough kinetic feedback and hence is unable to quench the catastrophic cooling. In this work, we present an alternate scenario wherein intense short-lived cooling and star formation phases followed by strong AGN outbursts are part of the AGN feedback loop. Using results from a 3D hydrodynamic simulation of a standard cool-core cluster ($M_{200}\sim 7\times10^{14}$ M$_\odot$; $z=0$), scaled to account for differences in mass and redshift, we argue that {\it Phoenix} is at the end of a cooling phase in which an AGN outburst has begun but has not yet arrested core cooling. This state of high cooling rate and star formation is expected to last for $\lesssim$ 100 Myr in {\it Phoenix}., Comment: 6 pages, 5 figures, Accepted for publication in MNRAS

Turbulence simulations play a key role in advancing the general understanding of the physical properties turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence simulations are often conducted in the isothermal limit. Given that the majority of astrophysical systems are not governed by isothermal dynamics, we aim to quantify the impact of thermodynamics on the physics of turbulence, through varying adiabatic index, $\gamma$, combined with a range of optically thin cooling functions. In this paper, we present a suite of ideal magnetohydrodynamics simulations of thermally balanced stationary turbulence in the subsonic, super-Alfv\'enic, high beta (ratio of thermal to magnetic pressure) regime, where turbulent dissipation is balanced by two idealized cooling functions (approximating linear cooling and free-free emission) and examine the impact of the equation of state by considering cases that correspond to isothermal, monatomic and diatomic gases. We find a strong anticorrelation between thermal and magnetic pressure independent of thermodynamics, whereas the strong anticorrelation between density and magnetic field found in the isothermal case weakens with increasing $\gamma$. Similarly, with the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing $\gamma$. This suggests that there exists a degeneracy in these relations with respect to thermodynamics and Mach number in this regime, which is dominated by slow magnetosonic modes. These results have implications for attempts to infer (e.g.) Mach numbers from (e.g.) Faraday rotation measurements, without additional information regarding the thermodynamics of the plasma. However, our results suggest that this degeneracy can be broken by utilizing higher-order moments of observable distribution functions., Comment: 16 pages, 7 figures, accepted for publication in ApJ

Using the Renaissance suite of simulations we examine the emergence of pristine atomic cooling haloes that are both metal-free and star-free in the early Universe. The absence of metals prevents catastrophic cooling, suppresses fragmentation, and may allow for the formation of massive black hole seeds. Here we report on the abundance of pristine atomic cooling haloes found and on the specific physical conditions that allow for the formation of these direct-collapse-black-hole (DCBH) haloes. In total in our simulations we find that 79 DCBH haloes form before a redshift of 11.6. We find that the formation of pristine atomic haloes is driven by the rapid assembly of the atomic cooling haloes with mergers, both minor and/or major, prior to reaching the atomic cooling limit a requirement. However, the ability of assembling haloes to remain free of (external) metal enrichment is equally important and underlines the necessity of following the transport of metals in such simulations. The candidate DCBH hosting haloes we find, have been exposed to mean Lyman-Werner radiation fields of J$_{LW}$ $\sim$ 1 J$_{21}$ and typically lie at least 10 kpc (physical) from the nearest massive galaxy. Growth rates of the haloes reach values of greater than 10$^7$ M$_{\odot}$ per unit redshift, leading to significant dynamical heating and the suppression of efficient cooling until the halo crosses the atomic cooling threshold. Finally, we also find five synchronised halo candidates where pairs of pristine atomic cooling haloes emerge that are both spatially and temporally synchronised., Comment: Minor updates after acceptance by MNRAS

Large scale simulations are a key pillar of modern research and require ever-increasing computational resources. Different novel manycore architectures have emerged in recent years on the way towards the exascale era. Performance portability is required to prevent repeated non-trivial refactoring of a code for different architectures. We combine Athena++, an existing magnetohydrodynamics (MHD) CPU code, with Kokkos, a performance portable on-node parallel programming paradigm, into K-Athena to allow efficient simulations on multiple architectures using a single codebase. We present profiling and scaling results for different platforms including Intel Skylake CPUs, Intel Xeon Phis, and NVIDIA GPUs. K-Athena achieves $>10^8$ cell-updates/s on a single V100 GPU for second-order double precision MHD calculations, and a speedup of 30 on up to 24,576 GPUs on Summit (compared to 172,032 CPU cores), reaching $1.94\times10^{12}$ total cell-updates/s at 76% parallel efficiency. Using a roofline analysis we demonstrate that the overall performance is currently limited by DRAM bandwidth and calculate a performance portability metric of 62.8%. Finally, we present the implementation strategies used and the challenges encountered in maximizing performance. This will provide other research groups with a straightforward approach to prepare their own codes for the exascale era. K-Athena is available at https://gitlab.com/pgrete/kathena ., Comment: 13 pages, 6 figures, 2 tables; accepted for publication in IEEE Transactions on Parallel and Distributed Systems (TPDS)

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution., Comment: Astro2020 Decadal Science White Paper

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges., Comment: 7 pages, astro2020 white paper

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions., Comment: 14 pages including 3 figures

The origin of supermassive black holes (SMBHs) that inhabit the centers of massive galaxies is largely unconstrained. Remnants from supermassive stars (SMSs) with masses around 10,000 solar masses provide the ideal seed candidates, known as direct collapse black holes. However, their very existence and formation environment in the early Universe are still under debate, with their supposed rarity further exacerbating the problem of modeling their ab-initio formation. SMS models have shown that rapid collapse, with an infall rate above a critical value, in metal-free haloes is a requirement for the formation of a proto-stellar core which will then form an SMS. Using a radiation hydrodynamics simulation of early galaxy formation, we show the natural emergence of metal-free haloes both massive enough, and with sufficiently high infall rates, to form an SMS. We find that haloes that are exposed to both a Lyman-Werner intensity of J_LW ~ 3 J_21 and that undergo at least one period of rapid growth early in their evolution are ideal cradles for SMS formation. This rapid growth induces substantial dynamical heating, amplifying the existing Lyman-Werner suppression originating from a group of young galaxies 20 kiloparsecs away. Our results strongly indicate that structure formation dynamics, rather than a critical Lyman-Werner (LW) flux, may be the main driver of massive black hole formation in the early Universe. We find that massive black hole seeds may be much more common in overdense regions of the early Universe than previously considered with a comoving number density up to 10^-3 Mpc^-3., Comment: To appear in Nature on 23 Jan 2019 for Advanced Online Publication. Authors' latest version. 13 pages in Nature print format, 13 figures, 1 table

Traditional cosmological hydrodynamics simulations fail to spatially resolve the circumgalatic medium (CGM), the reservoir of tenuous gas surrounding a galaxy and extending to its virial radius. We introduce the technique of Enhanced Halo Resolution (EHR), enabling more realistic physical modeling of the simulated CGM by consistently forcing gas refinement to smaller scales throughout the virial halo of a simulated galaxy. We investigate the effects of EHR in the Tempest simulations, a suite of Enzo-based cosmological zoom simulations following the evolution of an L* galaxy, resolving spatial scales of 500 comoving pc out to 100 comoving kpc in galactocentric radius. Among its many effects, EHR (1) changes the thermal balance of the CGM, increasing its cool gas content and decreasing its warm/hot gas content; (2) preserves cool gas structures for longer periods; and (3) enables these cool clouds to exist at progressively smaller size scales. Observationally, this results in a boost in "low ions" like H I and a drop in "high ions" like O VI throughout the CGM. These effects of EHR do not converge in the Tempest simulations, but extrapolating these trends suggests that the CGM in reality is a mist consisting of ubiquitous, small, long-lived, cool clouds suspended in a hot medium at the virial temperature of the halo. Additionally, we explore the physical mechanisms to explain why EHR produces the above effects, proposing that it works both by (1) better sampling the distribution of CGM phases enabling runaway cooling in the denser, cooler tail of the phase distribution; and (2) preventing cool gas clouds from artificially mixing with the ambient hot halo and evaporating. Evidence is found for both EHR mechanisms occurring in the Tempest simulations., Comment: 19 pages, 13 figures, submitted to ApJ. Simulation movies available at http://chummels.org/tempest

Observing the circumgalactic medium (CGM) in emission provides 3D maps of the spatial and kinematic extent of the gas that fuels galaxies and receives their feedback. We present mock emission-line maps of highly resolved CGM gas from the FOGGIE project (Figuring Out Gas & Galaxies in Enzo) and link these maps back to physical and spatial properties of the gas. By increasing the spatial resolution alone, the total luminosity of the line emission increases by an order of magnitude. This increase arises in the abundance of dense small-scale structure resolved when the CGM gas is simulated to < 100 pc scales. Current integral field unit instruments like KCWI and MUSE should be able to detect the brightest knots and filaments of such emission, and from this to infer the bulk kinematics of the CGM gas with respect to the galaxy. We conclude that accounting for small-scale structure well below the level of instrument spatial resolution is necessary to properly interpret such observations in terms of the underlying gas structure driving observable emission., Comment: 18 pages, 10 figures. Submitted to ApJ. Comments welcome

We present simulations from the new "Figuring Out Gas & Galaxies in Enzo" (FOGGIE) project. In contrast to most extant simulations of galaxy formation, which concentrate computational resources on galactic disks and spheroids with fluid and particle elements of fixed mass, the FOGGIE simulations focus on extreme spatial and mass resolution in the circumgalactic medium (CGM) surrounding galaxies. Using the Enzo code and a new refinement scheme, FOGGIE reaches spatial resolutions of 381 comoving $h^{-1}$ pc and resolves extremely low masses ($\lesssim 1$--$100$ Msun out to 100 comoving $h^{-1}$ kpc from the central halo. At these resolutions, cloud and filament-like structures giving rise to simulated absorption are smaller, and better resolved, than the same structures simulated with standard density-dependent refinement. Most of the simulated absorption arises in identifiable and well-resolved structures with masses $\lesssim 10^4$ Msun, well below the mass resolution of typical zoom simulations. However, integrated quantities such as mass surface density and ionic covering fractions change at only the $\lesssim 30$% level as resolution is varied. This relatively small changes in projected quantities---even when the sizes and distribution of absorbing clouds change dramatically---indicate that commonly used observables provide only weak constraints on the physical structure of the underlying gas. Comparing the simulated absorption features to the KODIAQ (Keck Observatory Database of Ionized Absorption toward Quasars) survey of $z \sim2$--$3.5$ Lyman limit systems, we show that high-resolution FOGGIE runs better resolve the internal kinematic structure of detected absorption, and better match the observed distribution of absorber properties. These results indicate that CGM resolution is key in properly testing simulations of galaxy evolution with circumgalactic observations., Comment: ApJ, in press. New appendix with results from simulation with 78pc resolution CGM at z=2.5 now included. Appendix C with additional velocity phase diagrams has been omitted to fit within the arXiv space restraints

The highly r-process enhanced (r-II) metal-poor halo stars we observe today could play a key role in understanding early ultra-faint dwarf galaxies, the smallest building blocks of the Milky Way. If a significant fraction of metal-poor r-II halo stars originated in the ultra-faint dwarf galaxies that merged to help form the Milky Way, observations of r-II stars could help us study these now-destroyed systems and probe the formation history of our Galaxy. To conduct our initial investigation into this possible connection, we use high-resolution cosmological simulations of Milky-Way-mass galaxies from the Caterpillar suite in combination with a simple, empirically motivated treatment of r-process enrichment. We determine the fraction of metal-poor halo stars that could have formed from highly r-process enhanced gas in now-destroyed low-mass ultra-faint dwarf galaxies, the simulated r-II fraction, and compare it to the "as observed" r-II fraction. We find that the simulated fraction, f_{r-II,sim} ~ 1-2%, can account for around half of the "as observed" fraction, f_{r-II,obs} ~ 2-4%. The "as observed" fraction likely overrepresents the fraction of r-II stars due to incomplete sampling, though, meaning f_{r-II,sim} likely accounts for more than half of the true f_{r-II,obs}. Further considering some parameter variations and scatter between individual simulations, the simulated fraction can account for around 20-80% of the "as observed" fraction., Comment: 16 pages, 6 figures, 2 tables, published in ApJ

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. $\alpha$ 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., Comment: 18 pages, 7 figures (plus 2 page, 2 figure appendix). Accepted to ApJ

In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event., Comment: Proceedings for the FRIB Theory Alliance workshop "FRIB and the GW170817 kilonova", held 16-27 July 2018 at the Facility for Rare Isotope Beams, Michigan State University, East Lansing, MI USA; 44 pages, 14 figures

The formation of stellar mass black holes from the remnants of Population III stars provides a source of initial black hole seeds with the potential to grow into intermediate or, in rare cases, possibly supermassive black holes. We use the Renaissance simulation suite to follow the growth of over 15,000 black holes born into mini-haloes in the early Universe. We compute the evolution of the black holes by post-processing individual remnant Population III star particles in the Renaissance simulation snapshots. The black holes populate haloes from 10$^{6}$ M$_{\odot}$ up to 10$^{9}$ M$_{\odot}$. We find that all of the black holes display very inefficient growth. On average the black holes increase their initial mass by a factor 10$^{-5}$, with the most active black holes increasing their mass by approximately 10%. Only a single black hole experiences any period of super-Eddington accretion, but the duration is very short and not repeated. Furthermore, we find no correlation of black hole accretion with halo mass in the mass range sampled. Within most haloes, we identify clumps of cool, dense gas for which accretion rates would be high, but instances of black holes encountering these clumps are rare and short-lived. Star formation competes with black hole growth by consuming available gas and driving down accretion rates through feedback. We conclude that the black holes born from Population III remnants do not form a significant population of intermediate mass black holes in the early Universe and will need to wait until later times to undergo significant accretion, if at all., Comment: 12 pages, 9 figures, accepted for publication in MNRAS

Many astrophysical systems encompass very large dynamical ranges in space and time, which are not accessible by direct numerical simulations. Thus, idealized subvolumes are often used to study small-scale effects including the dynamics of turbulence. These turbulent boxes require an artificial driving in order to mimic energy injection from large-scale processes. In this Letter, we show and quantify how the autocorrelation time of the driving and its normalization systematically change properties of an isothermal compressible magnetohydrodynamic flow in the sub- and supersonic regime and affect astrophysical observations such as Faraday rotation. For example, we find that $\delta$-in-time forcing with a constant energy injection leads to a steeper slope in kinetic energy spectrum and less efficient small-scale dynamo action. In general, we show that shorter autocorrelation times require more power in the acceleration field, which results in more power in compressive modes that weaken the anticorrelation between density and magnetic field strength. Thus, derived observables, such as the line-of-sight magnetic field from rotation measures, are systematically biased by the driving mechanism. We argue that $\delta$-in-time forcing is unrealistic and numerically unresolved, and conclude that special care needs to be taken in interpreting observational results based on the use of idealized simulations., Comment: 10 pages, 4 figures, version accepted for publication in ApJL

Extremely metal-poor (EMP) stars in the Milky Way (MW) allow us to infer the properties of their progenitors by comparing their chemical composition to the metal yields of the first supernovae. This method is most powerful when applied to mono-enriched stars, i.e. stars that formed from gas that was enriched by only one previous supernova. We present a novel diagnostic to identify this subclass of EMP stars. We model the first generations of star formation semi-analytically, based on dark matter halo merger trees that yield MW-like halos at the present day. Radiative and chemical feedback are included self-consistently and we trace all elements up to zinc. Mono-enriched stars account for only $\sim 1\%$ of second generation stars in our fiducial model and we provide an analytical formula for this probability. We also present a novel analytical diagnostic to identify mono-enriched stars, based on the metal yields of the first supernovae. This new diagnostic allows us to derive our main results independently from the specific assumptions made regarding Pop III star formation, and we apply it to a set of observed EMP stars to demonstrate its strengths and limitations. Our results may provide selection criteria for current and future surveys and therefore contribute to a deeper understanding of EMP stars and their progenitors., Comment: 18 pages, 20 figures, published in MNRAS