-
Why does the Milky Way have a metallicity floor?
Authors:
Britton D. Smith,
Brian W. O'Shea,
Sadegh Khochfar,
Matthew J. Turk,
John H. Wise,
Michael L. Norman
Abstract:
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 cos…
▽ More
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.
△ Less
Submitted 10 July, 2024; v1 submitted 12 June, 2024;
originally announced June 2024.
-
The formation and survival of cold gas in a magnetized cool-core galaxy cluster
Authors:
Martin Fournier,
Philipp Grete,
Marcus Brüggen,
Forrest W. Glines,
Brian W. O'Shea
Abstract:
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…
▽ More
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.
△ Less
Submitted 10 June, 2024; v1 submitted 7 June, 2024;
originally announced June 2024.
-
The Case for Hot-Mode Accretion in Abell 2029
Authors:
Deovrat Prasad,
G. Mark Voit,
Brian W. O'Shea
Abstract:
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 clust…
▽ More
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.
△ Less
Submitted 9 November, 2023;
originally announced November 2023.
-
Figuring Out Gas & Galaxies In Enzo (FOGGIE) VII: The (Dis)Assembly of Stellar Halos
Authors:
Anna C. Wright,
Jason Tumlinson,
Molly S. Peeples,
Brian W. O'Shea,
Cassandra Lochhaas,
Lauren Corlies,
Britton D. Smith,
Nguyen Binh,
Ramona Augustin,
Raymond C. Simons
Abstract:
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…
▽ More
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.
△ Less
Submitted 12 June, 2024; v1 submitted 18 September, 2023;
originally announced September 2023.
-
The Role of Radiation and Halo Mergers in Pop III Star Formation
Authors:
Lilia Correa Magnus,
Britton D. Smith,
Sadegh Khochfar,
Brian W. O'Shea,
John H. Wise,
Michael L. Norman,
Matthew J. Turk
Abstract:
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…
▽ More
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.
△ Less
Submitted 16 October, 2023; v1 submitted 7 July, 2023;
originally announced July 2023.
-
As a matter of dynamical range -- scale dependent energy dynamics in MHD turbulence
Authors:
Philipp Grete,
Brian W. O'Shea,
Kris Beckwith
Abstract:
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énic, high plasma-beta MHD turbulence simulations that only vary in their dynamical range, i.e., in their separation between the large-scal…
▽ More
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énic, 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.
△ Less
Submitted 17 November, 2022;
originally announced November 2022.
-
Figuring Out Gas & Galaxies In Enzo (FOGGIE) VI: The Circumgalactic Medium of $L^*$ Galaxies is Supported in an Emergent, Non-Hydrostatic Equilibrium
Authors:
Cassandra Lochhaas,
Jason Tumlinson,
Molly S. Peeples,
Brian W. O'Shea,
Jessica K. Werk,
Raymond C. Simons,
James Juno,
Claire E. Kopenhafer,
Ramona Augustin,
Anna C. Wright,
Ayan Acharyya,
Britton D. Smith
Abstract:
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. T…
▽ More
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 ($Δ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.
△ Less
Submitted 7 March, 2023; v1 submitted 20 June, 2022;
originally announced June 2022.
-
Possibilities and Limitations of Kinematically Identifying Stars from Accreted Ultra-Faint Dwarf Galaxies
Authors:
Kaley Brauer,
Hillary Diane Andales,
Alexander P. Ji,
Anna Frebel,
Mohammad K. Mardini,
Facundo A. Gomez,
Brian W. O'Shea
Abstract:
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 curre…
▽ More
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.
△ Less
Submitted 1 August, 2022; v1 submitted 14 June, 2022;
originally announced June 2022.
-
Seeking Self-Regulating Simulations of Idealized Milky Way-Like Galaxies
Authors:
Claire Kopenhafer,
Brian W. O'Shea,
G. Mark Voit
Abstract:
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 feedba…
▽ More
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.
△ Less
Submitted 9 June, 2022;
originally announced June 2022.
-
Horizons: Nuclear Astrophysics in the 2020s and Beyond
Authors:
H. Schatz,
A. D. Becerril Reyes,
A. Best,
E. F. Brown,
K. Chatziioannou,
K. A. Chipps,
C. M. Deibel,
R. Ezzeddine,
D. K. Galloway,
C. J. Hansen,
F. Herwig,
A. P. Ji,
M. Lugaro,
Z. Meisel,
D. Norman,
J. S. Read,
L. F. Roberts,
A. Spyrou,
I. Tews,
F. X. Timmes,
C. Travaglio,
N. Vassh,
C. Abia,
P. Adsley,
S. Agarwal
, et al. (140 additional authors not shown)
Abstract:
Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit…
▽ More
Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities.
△ Less
Submitted 16 May, 2022;
originally announced May 2022.
-
Atmospheric Circulation in Simulations of the AGN-CGM Connection at Halo Masses $\sim 10^{13.5}, M_\odot$
Authors:
Deovrat Prasad,
G. Mark Voit,
Brian W. O'Shea
Abstract:
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 sc…
▽ More
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$ ).
△ Less
Submitted 2 May, 2022; v1 submitted 28 October, 2021;
originally announced October 2021.
-
Some First Stars Were Red: Detecting Signatures of Massive Population III Formation Through Long-Term Stochastic Color Variations
Authors:
Tyrone E. Woods,
Chris J. Willott,
John A. Regan,
John H. Wise,
Turlough P. Downes,
Michael L. Norman,
Brian W. O'Shea
Abstract:
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…
▽ More
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.
△ Less
Submitted 27 September, 2021;
originally announced September 2021.
-
Figuring Out Gas & Galaxies In Enzo (FOGGIE) V: The Virial Temperature Does Not Describe Gas in a Virialized Galaxy Halo
Authors:
Cassandra Lochhaas,
Jason Tumlinson,
Brian W. O'Shea,
Molly S. Peeples,
Britton D. Smith,
Jessica K. Werk,
Ramona Augustin,
Raymond C. Simons
Abstract:
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 m…
▽ More
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.
△ Less
Submitted 7 September, 2021; v1 submitted 16 February, 2021;
originally announced February 2021.
-
External Enrichment of Minihalos by the First Supernovae
Authors:
William Hicks,
Azton Wells,
Michael L. Norman,
John H. Wise,
Britton D. Smith,
Brian W. O'Shea
Abstract:
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…
▽ More
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.
△ Less
Submitted 10 January, 2021; v1 submitted 11 September, 2020;
originally announced September 2020.
-
As a matter of tension -- kinetic energy spectra in MHD turbulence
Authors:
Philipp Grete,
Brian W. O'Shea,
Kris Beckwith
Abstract:
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 dyn…
▽ More
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énic 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.
△ Less
Submitted 7 September, 2020;
originally announced September 2020.
-
Magnetized Decaying Turbulence in the Weakly Compressible Taylor-Green Vortex
Authors:
Forrest W. Glines,
Philipp Grete,
Brian W. O'Shea
Abstract:
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…
▽ More
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.
△ Less
Submitted 26 April, 2021; v1 submitted 2 September, 2020;
originally announced September 2020.
-
The Formation of Very Massive Stars in Early Galaxies and Implications for Intermediate Mass Black Holes
Authors:
John A. Regan,
John H. Wise,
Tyrone E. Woods,
Turlough P. Downes,
Brian W. O'Shea,
Michael L. Norman
Abstract:
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 massi…
▽ More
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.
△ Less
Submitted 17 December, 2020; v1 submitted 18 August, 2020;
originally announced August 2020.
-
Massive Star Formation in Metal-Enriched Haloes at High Redshift
Authors:
John A. Regan,
Zoltán Haiman,
John H. Wise,
Brian W. O'Shea,
Michael L. Norman
Abstract:
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…
▽ More
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.
△ Less
Submitted 20 August, 2020; v1 submitted 25 June, 2020;
originally announced June 2020.
-
Environmental Dependence of Self-Regulating Black-hole Feedback in Massive Galaxies
Authors:
Deovrat Prasad,
G. Mark Voit,
Brian W. O'shea,
Forrest Glines
Abstract:
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 t…
▽ More
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.
△ Less
Submitted 23 November, 2020; v1 submitted 18 June, 2020;
originally announced June 2020.
-
A Black-Hole Feedback Valve in Massive Galaxies
Authors:
G. M. Voit,
G. L. Bryan,
D. Prasad,
R. Frisbie,
Y. Li,
M. Donahue,
B. W. O'Shea,
M. Sun,
N. Werner
Abstract:
Star formation in the universe's most massive galaxies proceeds furiously early in time but then nearly ceases. Plenty of hot gas remains available but does not cool and condense into star-forming clouds. Active galactic nuclei (AGN) release enough energy to inhibit cooling of the hot gas, but energetic arguments alone do not explain why quenching of star formation is most effective in high-mass g…
▽ More
Star formation in the universe's most massive galaxies proceeds furiously early in time but then nearly ceases. Plenty of hot gas remains available but does not cool and condense into star-forming clouds. Active galactic nuclei (AGN) release enough energy to inhibit cooling of the hot gas, but energetic arguments alone do not explain why quenching of star formation is most effective in high-mass galaxies. In fact, optical observations show that quenching is more closely related to a galaxy's central stellar velocity dispersion ($σ_v$) than to any other characteristic. Here, we show that high $σ_v$ is critical to quenching because a deep central potential well maximizes the efficacy of AGN feedback. In order to remain quenched, a galaxy must continually sweep out the gas ejected from its aging stars. Supernova heating can accomplish this task as long as the AGN sufficiently reduces the gas pressure of the surrounding circumgalactic medium (CGM). We find that CGM pressure acts as the control knob on a valve that regulates AGN feedback and suggest that feedback power self-adjusts so that it suffices to lift the CGM out of the galaxy's potential well. Supernova heating then drives a galactic outflow that remains homogeneous if $σ_v \gtrsim 240 \, {\rm km \, s^{-1}}$. AGN feedback can effectively quench galaxies with a comparable velocity dispersion, but feedback in galaxies with a much lower velocity dispersion tends to result in convective circulation and accumulation of multiphase gas within the galaxy.
△ Less
Submitted 23 October, 2020; v1 submitted 16 June, 2020;
originally announced June 2020.
-
Figuring Out Gas & Galaxies In Enzo (FOGGIE). IV. The Stochasticity of Ram Pressure Stripping in Galactic Halos
Authors:
Raymond C. Simons,
Molly S. Peeples,
Jason Tumlinson,
Brian W. O'Shea,
Britton D. Smith,
Lauren Corlies,
Cassandra Lochhaas,
Yong Zheng,
Ramona Augustin,
Deovrat Prasad,
Gregory F. Snyder,
Erik Tollerud
Abstract:
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---…
▽ More
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.
△ Less
Submitted 29 April, 2020;
originally announced April 2020.
-
Tests of AGN Feedback Kernels in Simulated Galaxy Clusters
Authors:
Forrest W. Glines,
Brian W. O'Shea,
G. Mark Voit
Abstract:
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 times…
▽ More
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.
△ Less
Submitted 9 November, 2021; v1 submitted 31 March, 2020;
originally announced April 2020.
-
Figuring Out Gas & Galaxies in Enzo (FOGGIE). III. The Mocky Way: Investigating Biases in Observing the Milky Way's Circumgalactic Medium
Authors:
Yong Zheng,
Molly S. Peeples,
Brian W. O'Shea,
Raymond C. Simons,
Cassandra Lochhaas,
Lauren Corlies,
Jason Tumlinson,
Britton D. Smith,
Ramona Augustin
Abstract:
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…
▽ More
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.
△ Less
Submitted 19 June, 2020; v1 submitted 21 January, 2020;
originally announced January 2020.
-
Cool-Core Cycles and Phoenix
Authors:
Deovrat Prasad,
Prateek Sharma,
Arif Babul,
G. Mark Voit,
Brian W. O'Shea
Abstract:
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…
▽ More
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}.
△ Less
Submitted 28 April, 2020; v1 submitted 27 September, 2019;
originally announced September 2019.
-
As A Matter of State: The role of thermodynamics in magnetohydrodynamic turbulence
Authors:
Philipp Grete,
Brian W. O'Shea,
Kris Beckwith
Abstract:
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 t…
▽ More
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, $γ$, 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énic, 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 $γ$. Similarly, with the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing $γ$. 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.
△ Less
Submitted 5 December, 2019; v1 submitted 11 August, 2019;
originally announced August 2019.
-
The Emergence of the First Star-free Atomic Cooling Haloes in the Universe
Authors:
John A. Regan,
John H. Wise,
Brian W. O'Shea,
Michael L. Norman
Abstract:
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 phys…
▽ More
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.
△ Less
Submitted 8 January, 2020; v1 submitted 7 August, 2019;
originally announced August 2019.
-
K-Athena: a performance portable structured grid finite volume magnetohydrodynamics code
Authors:
Philipp Grete,
Forrest W. Glines,
Brian W. O'Shea
Abstract:
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…
▽ More
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 .
△ Less
Submitted 15 July, 2020; v1 submitted 10 May, 2019;
originally announced May 2019.
-
Circumgalactic Gas and the Precipitation Limit
Authors:
G. M. Voit,
A. Babul,
Iu. Babyk,
G. L. Bryan,
H. -W. Chen,
M. Donahue,
D. Fielding,
M. Gaspari,
Y. Li,
M. McDonald,
B. W. O'Shea,
D. Prasad,
P. Sharma,
M. Sun,
G. Tremblay,
J. Werk,
N. Werner,
F. Zahedy
Abstract:
During the last decade, numerous and varied observations, along with increasingly sophisticated numerical simulations, have awakened astronomers to the central role the circumgalactic medium (CGM) plays in regulating galaxy evolution. It contains the majority of the baryonic matter associated with a galaxy, along with most of the metals, and must continually replenish the star forming gas in galax…
▽ More
During the last decade, numerous and varied observations, along with increasingly sophisticated numerical simulations, have awakened astronomers to the central role the circumgalactic medium (CGM) plays in regulating galaxy evolution. It contains the majority of the baryonic matter associated with a galaxy, along with most of the metals, and must continually replenish the star forming gas in galaxies that continue to sustain star formation. And while the CGM is complex, containing gas ranging over orders of magnitude in temperature and density, a simple emergent property may be governing its structure and role. Observations increasingly suggest that the ambient CGM pressure cannot exceed the limit at which cold clouds start to condense out and precipitate toward the center of the potential well. If feedback fueled by those clouds then heats the CGM and causes it to expand, the pressure will drop and the "rain" will diminish. Such a feedback loop tends to suspend the CGM at the threshold pressure for precipitation. The coming decade will offer many opportunities to test this potentially fundamental principle of galaxy evolution.
△ Less
Submitted 29 March, 2019; v1 submitted 26 March, 2019;
originally announced March 2019.
-
Understanding the circumgalactic medium is critical for understanding galaxy evolution
Authors:
Molly S. Peeples,
Peter Behroozi,
Rongmon Bordoloi,
Alyson Brooks,
James S. Bullock,
Joseph N. Burchett,
Hsiao-Wen Chen,
John Chisholm,
Charlotte Christensen,
Alison Coil,
Lauren Corlies,
Aleksandar Diamond-Stanic,
Megan Donahue,
Claude-André Faucher-Giguère,
Henry Ferguson,
Drummond Fielding,
Andrew J. Fox,
David M. French,
Steven R. Furlanetto,
Mario Gennaro,
Karoline M. Gilbert,
Erika Hamden,
Nimish Hathi,
Matthew Hayes,
Alaina Henry
, et al. (47 additional authors not shown)
Abstract:
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…
▽ More
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.
△ Less
Submitted 13 March, 2019;
originally announced March 2019.
-
Cyberinfrastructure Requirements to Enhance Multi-messenger Astrophysics
Authors:
Philip Chang,
Gabrielle Allen,
Warren Anderson,
Federica B. Bianco,
Joshua S. Bloom,
Patrick R. Brady,
Adam Brazier,
S. Bradley Cenko,
Sean M. Couch,
Tyce DeYoung,
Ewa Deelman,
Zachariah B Etienne,
Ryan J. Foley,
Derek B Fox,
V. Zach Golkhou,
Darren R Grant,
Chad Hanna,
Kelly Holley-Bockelmann,
D. Andrew Howell,
E. A. Huerta,
Margaret W. G. Johnson,
Mario Juric,
David L. Kaplan,
Daniel S. Katz,
Azadeh Keivani
, et al. (17 additional authors not shown)
Abstract:
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,…
▽ More
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.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
Circumgalactic Pressure Profiles Indicate Precipitation-Limited Atmospheres for $M_* \sim 10^9$$-$$10^{11.5}\,M_\odot$
Authors:
G. M. Voit,
M. Donahue,
F. Zahedy,
H. -W. Chen,
J. Werk,
G. L. Bryan,
B. W. O'Shea
Abstract:
Cosmic gas cycles in and out of galaxies, but outside of galaxies it is difficult to observe except for the absorption lines that circumgalactic clouds leave in the spectra of background quasars. Using photoionization modeling of those lines to determine cloud pressures, we find that galaxies are surrounded by extended atmospheres that confine the clouds and have a radial pressure profile that dep…
▽ More
Cosmic gas cycles in and out of galaxies, but outside of galaxies it is difficult to observe except for the absorption lines that circumgalactic clouds leave in the spectra of background quasars. Using photoionization modeling of those lines to determine cloud pressures, we find that galaxies are surrounded by extended atmospheres that confine the clouds and have a radial pressure profile that depends on galaxy mass. Motivated by observations of the universe's most massive galaxies, we compare those pressure measurements with models predicting the critical pressure at which cooler clouds start to precipitate out of the hot atmosphere and rain toward the center. We find excellent agreement, implying that the precipitation limit applies to galaxies over a wide mass range.
△ Less
Submitted 14 June, 2019; v1 submitted 28 February, 2019;
originally announced March 2019.
-
Formation of massive black holes in rapidly growing pre-galactic gas clouds
Authors:
John H. Wise,
John A. Regan,
Brian W. O'Shea,
Michael L. Norman,
Turlough P. Downes,
Hao Xu
Abstract:
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 f…
▽ More
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.
△ Less
Submitted 22 January, 2019;
originally announced January 2019.
-
The Impact of Enhanced Halo Resolution on the Simulated Circumgalactic Medium
Authors:
Cameron B. Hummels,
Britton D. Smith,
Philip F. Hopkins,
Brian W. O'Shea,
Devin W. Silvia,
Jessica K. Werk,
Nicolas Lehner,
John H. Wise,
David C. Collins,
Iryna S. Butsky
Abstract:
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 h…
▽ More
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.
△ Less
Submitted 29 November, 2018;
originally announced November 2018.
-
Figuring Out Gas & Galaxies in Enzo (FOGGIE). II. Emission from the z=3 Circumgalactic Medium
Authors:
Lauren Corlies,
Molly S. Peeples,
Jason Tumlinson,
Brian W. O'Shea,
Nicolas Lehner,
J. Christopher Howk,
John M. O'Meara
Abstract:
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…
▽ More
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.
△ Less
Submitted 12 November, 2018;
originally announced November 2018.
-
Figuring Out Gas & Galaxies in Enzo (FOGGIE). I. Resolving Simulated Circumgalactic Absorption at 2 < z < 2.5
Authors:
Molly S. Peeples,
Lauren Corlies,
Jason Tumlinson,
Brian W. O'Shea,
Nicolas Lehner,
John M. O'Meara,
J. Christopher Howk,
Britton D. Smith,
John H. Wise,
Cameron B. Hummels
Abstract:
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. U…
▽ More
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.
△ Less
Submitted 14 March, 2019; v1 submitted 15 October, 2018;
originally announced October 2018.
-
The Origin of r-process Enhanced Metal-Poor Halo Stars In Now-Destroyed Ultra-Faint Dwarf Galaxies
Authors:
Kaley Brauer,
Alexander P. Ji,
Anna Frebel,
Gregory A. Dooley,
Facundo A. Gomez,
Brian W. O'Shea
Abstract:
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-destro…
▽ More
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.
△ Less
Submitted 6 February, 2019; v1 submitted 14 September, 2018;
originally announced September 2018.
-
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.
-
FRIB and the GW170817 Kilonova
Authors:
A. Aprahamian,
R. Surman,
A. Frebel,
G. C. McLaughlin,
A. Arcones,
A. B. Balantekin,
J. Barnes,
Timothy C. Beers,
Erika M. Holmbeck,
Jinmi Yoon,
Maxime Brodeur,
T. M. Sprouse,
Nicole Vassh,
Jolie A. Cizewski,
Jason A. Clark,
Benoit Cote,
Sean M. Couch,
M. Eichler,
Jonathan Engel,
Rana Ezzeddine,
George M. Fuller,
Samuel A. Giuliani,
Robert Grzywacz,
Sophia Han,
C. J. Horowitz
, et al. (23 additional authors not shown)
Abstract:
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 experim…
▽ More
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.
△ Less
Submitted 3 September, 2018;
originally announced September 2018.
-
r-Process Nucleosynthesis: Connecting Rare-Isotope Beam Facilities with the Cosmos
Authors:
C. J. Horowitz,
A. Arcones,
B. Côté,
I. Dillmann,
W. Nazarewicz,
I. U. Roederer,
H. Schatz,
A. Aprahamian,
D. Atanasov,
A. Bauswein,
J. Bliss,
M. Brodeur,
J. A. Clark,
A. Frebel,
F. Foucart,
C. J. Hansen,
O. Just,
A. Kankainen,
G. C. McLaughlin,
J. M. Kelly,
S. N. Liddick,
D. M. Lee,
J. Lippuner,
D. Martin,
J. Mendoza-Temis
, et al. (13 additional authors not shown)
Abstract:
This is an exciting time for the study of r-process nucleosynthesis. Recently, a neutron star merger GW170817 was observed in extraordinary detail with gravitational waves and electromagnetic radiation from radio to gamma rays. The very red color of the associated kilonova suggests that neutron star mergers are an important r-process site. Astrophysical simulations of neutron star mergers and core…
▽ More
This is an exciting time for the study of r-process nucleosynthesis. Recently, a neutron star merger GW170817 was observed in extraordinary detail with gravitational waves and electromagnetic radiation from radio to gamma rays. The very red color of the associated kilonova suggests that neutron star mergers are an important r-process site. Astrophysical simulations of neutron star mergers and core collapse supernovae are making rapid progress. Detection of both, electron neutrinos and antineutrinos from the next galactic supernova will constrain the composition of neutrino-driven winds and provide unique nucleosynthesis information. Finally FRIB and other rare-isotope beam facilities will soon have dramatic new capabilities to synthesize many neutron-rich nuclei that are involved in the r-process. The new capabilities can significantly improve our understanding of the r-process and likely resolve one of the main outstanding problems in classical nuclear astrophysics. However, to make best use of the new experimental capabilities and to fully interpret the results, a great deal of infrastructure is needed in many related areas of astrophysics, astronomy, and nuclear theory. We will place these experiments in context by discussing astrophysical simulations and observations of r-process sites, observations of stellar abundances, galactic chemical evolution, and nuclear theory for the structure and reactions of very neutron-rich nuclei. This review paper was initiated at a three-week International Collaborations in Nuclear Theory program in June 2016 where we explored promising r-process experiments and discussed their likely impact, and their astrophysical, astronomical, and nuclear theory context.
△ Less
Submitted 11 May, 2018;
originally announced May 2018.
-
As a matter of force - Systematic biases in idealized turbulence simulations
Authors:
Philipp Grete,
Brian W. O'Shea,
Kris Beckwith
Abstract:
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 quanti…
▽ More
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 $δ$-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 $δ$-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.
△ Less
Submitted 6 May, 2018; v1 submitted 14 March, 2018;
originally announced March 2018.
-
Descendants of the first stars: the distinct chemical signature of second generation stars
Authors:
Tilman Hartwig,
Naoki Yoshida,
Mattis Magg,
Anna Frebel,
Simon C. O. Glover,
Facundo A. Gómez,
Brendan Griffen,
Miho N. Ishigaki,
Alexander P. Ji,
Ralf S. Klessen,
Brian W. O'Shea,
Nozomu Tominaga
Abstract:
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…
▽ More
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.
△ Less
Submitted 16 July, 2018; v1 submitted 15 January, 2018;
originally announced January 2018.
-
Selecting ultra-faint dwarf candidate progenitors in cosmological N-body simulations at high redshifts
Authors:
Mohammadtaher Safarzadeh,
Alexander P. Ji,
Gregory A. Dooley,
Anna Frebel,
Evan Scannapieco,
Facundo A. Gómez,
Brian W. O'Shea
Abstract:
The smallest satellites of the Milky Way ceased forming stars during the epoch of reionization and thus provide archaeological access to galaxy formation at $z>6$. Numerical studies of these ultra-faint dwarf galaxies (UFDs) require expensive cosmological simulations with high mass resolution that are carried out down to $z=0$. However, if we are able to statistically identify UFD host progenitors…
▽ More
The smallest satellites of the Milky Way ceased forming stars during the epoch of reionization and thus provide archaeological access to galaxy formation at $z>6$. Numerical studies of these ultra-faint dwarf galaxies (UFDs) require expensive cosmological simulations with high mass resolution that are carried out down to $z=0$. However, if we are able to statistically identify UFD host progenitors at high redshifts \emph{with relatively high probabilities}, we can avoid this high computational cost. To find such candidates, we analyze the merger trees of Milky Way type halos from the high-resolution ${\it Caterpillar}$ suite of dark matter only simulations. Satellite UFD hosts at $z=0 $ are identified based on four different abundance matching (AM) techniques. All the halos at high redshifts are traced forward in time in order to compute the probability of surviving as satellite UFDs today. Our results show that selecting potential UFD progenitors based solely on their mass at z=12 (8) results in a 10\% (20\%) chance of obtaining a surviving UFD at $z=0$ in three of the AM techniques we adopted. We find that the progenitors of surviving satellite UFDs have lower virial ratios ($η$), and are preferentially located at large distances from the main MW progenitor, while they show no correlation with concentration parameter. Halos with favorable locations and virial ratios are $\approx 3$ times more likely to survive as satellite UFD candidates at $z=0.$
△ Less
Submitted 1 March, 2018; v1 submitted 11 December, 2017;
originally announced December 2017.
-
Validating Semi-Analytic Models of High-Redshift Galaxy Formation using Radiation Hydrodynamical Simulations
Authors:
Benoit Côté,
Devin Silvia,
Brian W. O'Shea,
Britton Smith,
John H. Wise
Abstract:
We use a cosmological hydrodynamic simulation calculated with Enzo and the semi-analytic galaxy formation model (SAM) GAMMA to address the chemical evolution of dwarf galaxies in the early universe. The long-term goal of the project is to better understand the origin of metal-poor stars and the formation of dwarf galaxies and the Milky Way halo by cross-validating these theoretical approaches. We…
▽ More
We use a cosmological hydrodynamic simulation calculated with Enzo and the semi-analytic galaxy formation model (SAM) GAMMA to address the chemical evolution of dwarf galaxies in the early universe. The long-term goal of the project is to better understand the origin of metal-poor stars and the formation of dwarf galaxies and the Milky Way halo by cross-validating these theoretical approaches. We combine GAMMA with the merger tree of the most massive galaxy found in the hydrodynamic simulation and compare the star formation rate, the metallicity distribution function (MDF), and the age-metallicity relationship predicted by the two approaches. We found that the SAM can reproduce the global trends of the hydrodynamic simulation. However, there are degeneracies between the model parameters and more constraints (e.g., star formation efficiency, gas flows) need to be extracted from the simulation to isolate the correct semi-analytic solution. Stochastic processes such as bursty star formation histories and star formation triggered by supernova explosions cannot be reproduced by the current version of GAMMA. Non-uniform mixing in the galaxy's interstellar medium, coming primarily from self-enrichment by local supernovae, causes a broadening in the MDF that can be emulated in the SAM by convolving its predicted MDF with a Gaussian function having a standard deviation of ~0.2 dex. We found that the most massive galaxy in the simulation retains nearby 100% of its baryonic mass within its virial radius, which is in agreement with what is needed in GAMMA to reproduce the global trends of the simulation.
△ Less
Submitted 28 March, 2018; v1 submitted 17 October, 2017;
originally announced October 2017.
-
First Light II: Emission Line Extinction, Population III Stars, and X-ray Binaries
Authors:
Kirk S. S. Barrow,
John H. Wise,
Aycin Aykutalp,
Brian W. O'Shea,
Michael L. Norman,
Hao Xu
Abstract:
We produce synthetic spectra and observations for metal-free stellar populations and high mass X-ray binaries in the Renaissance Simulations at a redshift of 15. We extend our methodology from the first paper in the series by modelling the production and extinction of emission lines throughout a dusty and metal-enriched interstellar and circum-galactic media extracted from the simulation, using a…
▽ More
We produce synthetic spectra and observations for metal-free stellar populations and high mass X-ray binaries in the Renaissance Simulations at a redshift of 15. We extend our methodology from the first paper in the series by modelling the production and extinction of emission lines throughout a dusty and metal-enriched interstellar and circum-galactic media extracted from the simulation, using a Monte Carlo calculation. To capture the impact of high-energy photons, we include all frequencies from hard X-ray to far infrared with enough frequency resolution to discern line emission and absorption profiles. The most common lines in our sample in order of their rate of occurrence are Ly$α$, the C IV $λ\lambda1548,1551$ doublet, H-$α$, and the Ca II $λλ\lambda8498,8542,8662$ triplet. The best scenario for a direct observation of a metal-free stellar population is a merger between two Population III galaxies. In mergers between metal-enriched and metal-free stellar populations, some characteristics may be inferred indirectly. Single Population III galaxies are too dim to be observed photometrically at $z = 15$. Ly$α$ emission is discernible by $JWST$ as an increase in $\rm{J_{200w} - J_{277w}}$ colour off the intrinsic stellar tracks. Observations of metal-free stars will be difficult, though not impossible, with the next generation of space telescopes.
△ Less
Submitted 8 December, 2017; v1 submitted 13 September, 2017;
originally announced September 2017.
-
The Second Data Release of the KODIAQ Survey
Authors:
John M. O'Meara,
Nicolas Lehner,
J. Christopher Howk,
J. Xavier Prochaska,
Andrew J. Fox,
Molly S. Peeples,
Jason Tumlinson,
Brian W. O'Shea
Abstract:
We present and make publicly available the second data release (DR2) of the Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ) survey. KODIAQ DR2 consists of a fully-reduced sample of 300 quasars at 0.07 < z_em < 5.29 observed with HIRES at high resolution (36,000 <= R <= 103,000). DR2 contains 831 spectra available in continuum normalized form, representing a sum total exposu…
▽ More
We present and make publicly available the second data release (DR2) of the Keck Observatory Database of Ionized Absorption toward Quasars (KODIAQ) survey. KODIAQ DR2 consists of a fully-reduced sample of 300 quasars at 0.07 < z_em < 5.29 observed with HIRES at high resolution (36,000 <= R <= 103,000). DR2 contains 831 spectra available in continuum normalized form, representing a sum total exposure time of ~4.9 megaseconds on source. These co-added spectra arise from a total of 1577 individual exposures of quasars taken from the Keck Observatory Archive (KOA) in raw form and uniformly processed. DR2 extends DR1 by adding 130 new quasars to the sample, including additional observations of QSOs in DR1. All new data in DR2 were obtained with the single-chip Tektronix TK2048 CCD configuration of HIRES in operation between 1995 and 2004. DR2 is publicly available to the community, housed as a higher level science product at the KOA and in the igmspec database (v03).
△ Less
Submitted 25 July, 2017;
originally announced July 2017.
-
Energy transfer in compressible magnetohydrodynamic turbulence
Authors:
Philipp Grete,
Brian W. O'Shea,
Kris Beckwith,
Wolfram Schmidt,
Andrew Christlieb
Abstract:
Magnetic fields, compressibility and turbulence are important factors in many terrestrial and astrophysical processes. While energy dynamics, i.e. how energy is transferred within and between kinetic and magnetic reservoirs, has been previously studied in the context of incompressible magnetohydrodynamic (MHD) turbulence, we extend shell-to-shell energy transfer analysis to the compressible regime…
▽ More
Magnetic fields, compressibility and turbulence are important factors in many terrestrial and astrophysical processes. While energy dynamics, i.e. how energy is transferred within and between kinetic and magnetic reservoirs, has been previously studied in the context of incompressible magnetohydrodynamic (MHD) turbulence, we extend shell-to-shell energy transfer analysis to the compressible regime. We derive four new transfer functions specifically capturing compressibility effects in the kinetic and magnetic cascade, and capturing energy exchange via magnetic pressure. To illustrate their viability, we perform and analyze four simulations of driven isothermal MHD turbulence in the sub- and supersonic regime with two different codes. On the one hand, our analysis reveals robust characteristics across regime and numerical method. For example, energy transfer between individual scales is local and forward for both cascades with the magnetic cascade being stronger than the kinetic one. Magnetic tension and magnetic pressure related transfers are less local and weaker than the cascades. We find no evidence for significant nonlocal transfer. On the other hand, we show that certain functions, e.g., the compressive component of the magnetic energy cascade, exhibit a more complex behavior that varies both with regime and numerical method. Having established a basis for the analysis in the compressible regime, the method can now be applied to study a broader parameter space.
△ Less
Submitted 3 October, 2017; v1 submitted 20 June, 2017;
originally announced June 2017.
-
First Light: exploring the Spectra of High-Redshift Galaxies in the Renaissance Simulations
Authors:
Kirk S. S. Barrow,
John H. Wise,
Michael L. Norman,
Brian W. O'Shea,
Hao Xu
Abstract:
We present synthetic observations for the first generations of galaxies in the Universe and make predictions for future deep field observations for redshifts greater than 6. Due to the strong impact of nebular emission lines and the relatively compact scale of HII regions, high resolution cosmological simulations and a robust suite of analysis tools are required to properly simulate spectra. We cr…
▽ More
We present synthetic observations for the first generations of galaxies in the Universe and make predictions for future deep field observations for redshifts greater than 6. Due to the strong impact of nebular emission lines and the relatively compact scale of HII regions, high resolution cosmological simulations and a robust suite of analysis tools are required to properly simulate spectra. We created a software pipeline consisting of FSPS, Hyperion, Cloudy and our own tools to generate synthetic IR observations from a fully three-dimensional arrangement of gas, dust, and stars. Our prescription allows us to include emission lines for a complete chemical network and tackle the effect of dust extinction and scattering in the various lines of sight. We provide spectra, 2-D binned photon imagery for both HST and JWST IR filters, luminosity relationships, and emission line strengths for a large sample of high redshift galaxies in the Renaissance Simulations. Our resulting synthetic spectra show high variability between galactic halos with a strong dependence on stellar mass, metallicity, gas mass fraction, and formation history. Halos with the lowest stellar mass have the greatest variability in [OIII]/H$β$, [OIII] and CIII] while halos with higher masses are seen to show consistency in their spectra and [OIII] equivalent widths (EW) between 1Å and 10Å. Viewing angle accounted for three-fold difference in flux due to the presence of ionized gas channels in a halo. Furthermore, JWST color plots show a discernible relationship between redshift, color, and mean stellar age.
△ Less
Submitted 19 June, 2017; v1 submitted 10 January, 2017;
originally announced January 2017.
-
Tracing the first stars and galaxies of the Milky Way
Authors:
Brendan F. Griffen,
Gregory A. Dooley,
Alexander P. Ji,
Brian W. O'Shea,
Facundo A. Gómez,
Anna Frebel
Abstract:
We use 30 high-resolution dark matter halos of the $Caterpillar$ simulation suite to probe the first stars and galaxies of Milky Way-mass systems. We quantify the environment of the high-$z$ progenitors of the Milky Way and connect them to the properties of the host and satellites today. We identify the formation sites of the first generation of Population III (Pop III) stars ($z$ ~ 25) and first…
▽ More
We use 30 high-resolution dark matter halos of the $Caterpillar$ simulation suite to probe the first stars and galaxies of Milky Way-mass systems. We quantify the environment of the high-$z$ progenitors of the Milky Way and connect them to the properties of the host and satellites today. We identify the formation sites of the first generation of Population III (Pop III) stars ($z$ ~ 25) and first galaxies ($z$ ~ 22) with several different models based on a minimum halo mass including a simple model for Lyman-Werner feedback. Through this method we find approximately 23,000 $\pm$ 5,000 Pop III potentially star-forming sites per Milky Way-mass host, though this number is drastically reduced to ~550 star-forming sites when Lyman-Werner feedback is included, as it has critical effects at these length scales. The majority of these halos identified form in isolation (96% at $z$ = 15) and are not subject to external enrichment by neighboring halos (median separation ~1 pkpc at $z$ = 15), though half merge with a system larger than themselves within 1.5 Gyrs. Approximately 55% of the entire population has merged into the host halo by $z$ = 0. Using particle tagging, we additionally trace the Pop III remnant population to $z$ = 0 and find an order of magnitude scatter in their number density at small (i.e. r $<$ 5 kpc) and large (i.e. r $>$ 50 kpc) galactocentric radii at $z$ = 0. Using our large number of realizations, we provide fitting functions for determining the number of progenitor minihalo and atomic cooling halo systems that present-day dwarf galaxies and the Magellanic cloud system might have accreted since their formation. We demonstrate that observed dwarf galaxies with stellar masses below 10$^{4.6}$ M$_{\odot}$ are unlikely to have merged with any other star-forming systems.
△ Less
Submitted 21 November, 2016; v1 submitted 2 November, 2016;
originally announced November 2016.
-
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.
-
The AGORA High-Resolution Galaxy Simulations Comparison Project. II: Isolated Disk Test
Authors:
Ji-hoon Kim,
Oscar Agertz,
Romain Teyssier,
Michael J. Butler,
Daniel Ceverino,
Jun-Hwan Choi,
Robert Feldmann,
Ben W. Keller,
Alessandro Lupi,
Thomas Quinn,
Yves Revaz,
Spencer Wallace,
Nickolay Y. Gnedin,
Samuel N. Leitner,
Sijing Shen,
Britton D. Smith,
Robert Thompson,
Matthew J. Turk,
Tom Abel,
Kenza S. Arraki,
Samantha M. Benincasa,
Sukanya Chakrabarti,
Colin DeGraf,
Avishai Dekel,
Nathan J. Goldbaum
, et al. (18 additional authors not shown)
Abstract:
Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the st…
▽ More
Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.
△ Less
Submitted 8 March, 2018; v1 submitted 10 October, 2016;
originally announced October 2016.