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.

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.

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.

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 ($\sigma_v$) than to any other characteristic. Here, we show that high $\sigma_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 $\sigma_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.

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.

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

Observations since the 1970's have revealed the existence of cool-core (CC) clusters, which are galaxy clusters with a central cooling time much shorter than the age of the universe. Both observations and theory suggest that the ambient gas at the centers of galaxy clusters is thermally regulated by a central heating mechanism that suppresses condensation (most likely an active galactic nucleus, or AGN). Previous analytical work has suggested specific configurations of heating kernels that may result in thermal balance and a steady state. To test this hypothesis, we simulated idealized galaxy clusters using the ENZO cosmology code with a spatial heat-input kernel meant to mimic feedback from a central AGN. Thermal heating as a function of radius was injected according to a range of kernels, with global thermal balance enforced at all times. We compare our simulation results with observed entropy profiles from the ACCEPT cluster dataset. Although some heating kernels produced thermally steady galaxy clusters, no kernel was able to produce a steady cluster with a central entropy as low as the central entropies typically observed among CC clusters. The general behavior of the simulations depended on the amount 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 heating leading to a flat but elevated entropy core. The simulated clusters enter an unsteady multiphase state on a timescale proportional to the square of the cooling time of the lowest entropy gas in the simulation, with centrally concentrated heating resulting in a steady state lasting for a longer period of time.

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.

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

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.

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 83.1%. Finally, we present the strategies used for implementation and the challenges encountered 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 .

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.

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.

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.

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.

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

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.

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.

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.

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.

Using a high resolution simulation of an isolated dwarf galaxy, accounting for multi-channel stellar feedback and chemical evolution on a star-by-star basis, we investigate how each of 15 metal species are distributed within our multi-phase interstellar medium (ISM) and ejected from our galaxy by galactic winds. For the first time, we demonstrate that the mass fraction probability distribution functions (PDFs) of individual metal species in the ISM are well described by a piecewise log-normal and power-law distribution. The PDF properties vary within each ISM phase. Hot gas is dominated by recent enrichment, with a significant power-law tail to high metal fractions, while cold gas is predominately log-normal. In addition, elements dominated by asymptotic giant branch (AGB) wind enrichment (e.g. N and Ba) mix less efficiently than elements dominated by supernova enrichment (e.g. $\alpha$ elements and Fe). This result is driven by the differences in source energetics and source locations, particularly the higher chance compared to massive stars for AGB stars to eject material into cold gas. Nearly all of the produced metals are ejected from the galaxy (only 4% are retained), but over 20% of metals dominated by AGB enrichment are retained. In dwarf galaxies, therefore, elements synthesized predominately through AGB winds should be both overabundant and have a larger spread compared to elements synthesized in either core collapse or Type Ia supernovae. We discuss the observational implications of these results, their potential use in developing improved models of galactic chemical evolution, and their generalization to more massive galaxies.

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.

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

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.

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.

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

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.

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.

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 ($\eta$), 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.$

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$\alpha$, the C IV $\lambda\lambda1548,1551$ doublet, H-$\alpha$, and the Ca II $\lambda\lambda\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$\alpha$ 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.

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.

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

We provide an analytic framework for interpreting observations of multiphase circumgalactic gas that is heavily informed by recent numerical simulations of thermal instability and precipitation in cool-core galaxy clusters. We start by considering the local conditions required for the formation of multiphase gas via two different modes: (1) uplift of ambient gas by galactic outflows, and (2) condensation in a stratified stationary medium in which thermal balance is explicitly maintained. Analytic exploration of these two modes provides insights into the relationships between the local ratio of the cooling and freefall time scales (i.e., t_cool / t_ff), the large-scale gradient of specific entropy, and development of precipitation and multiphase media in circumgalactic gas. We then use these analytic findings to interpret recent simulations of circumgalactic gas in which global thermal balance is maintained. We show that long-lasting configurations of gas with 5 < t_cool / t_ff < 20 and radial entropy profiles similar to observations of local cool-core galaxy cluster cores are a natural outcome of precipitation-regulated feedback. We conclude with some observational predictions that follow from these models. This work focuses primarily on precipitation and AGN feedback in galaxy cluster cores, because that is where the observations of multiphase gas around galaxies are most complete. However, many of the physical principles that govern condensation in those environments apply to circumgalactic gas around galaxies of all masses.

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$\beta$, [OIII] and CIII] while halos with higher masses are seen to show consistency in their spectra and [OIII] equivalent widths (EW) between 1\AA\ and 10\AA. 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.

The role of compact binary mergers as the main production site of r-process elements is investigated by combining stellar abundances of Eu observed in the Milky Way, galactic chemical evolution (GCE) simulations, binary population synthesis models, and Advanced LIGO gravitational wave measurements. We compiled and reviewed seven recent GCE studies to extract the frequency of neutron star - neutron star (NS-NS) mergers that is needed in order to reproduce the observed [Eu/Fe] vs [Fe/H] relationship. We used our simple chemical evolution code to explore the impact of different analytical delay-time distribution (DTD) functions for NS-NS mergers. We then combined our metallicity-dependent population synthesis models with our chemical evolution code to bring their predictions, for both NS-NS mergers and black hole - neutron star mergers, into a GCE context. Finally, we convolved our results with the cosmic star formation history to provide a direct comparison with current and upcoming Advanced LIGO measurements. When assuming that NS-NS mergers are the exclusive r-process sites, and that the ejected r-process mass per merger event is 0.01 Msun, the number of NS-NS mergers needed in GCE studies is about 10 times larger than what is predicted by standard population synthesis models. These two distinct fields can only be consistent with each other when assuming optimistic rates, massive NS-NS merger ejecta, and low Fe yields for massive stars. Population synthesis models and GCE simulations are in agreement with the current upper limit (O1) established by Advanced LIGO during their first run of observations. Upcoming measurements will provide important constraints on the actual local NS-NS merger rate, will provide insights on the plausibility of the GCE requirement, and will help to define whether or not compact binary mergers can be the dominant source of r-process elements in the Universe.

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.

We use the OMEGA galactic chemical evolution code to investigate how the assumptions used for the treatment of galactic inflows and outflows impact numerical predictions. The goal is to determine how our capacity to reproduce the chemical evolution trends of a galaxy is affected by the choice of implementation used to include those physical processes. In pursuit of this goal, we experiment with three different prescriptions for galactic inflows and outflows and use OMEGA within a Markov Chain Monte Carlo code to recover the set of input parameters that best reproduces the chemical evolution of nine elements in the dwarf spheroidal galaxy Sculptor. This provides a consistent framework for comparing the best-fit solutions generated by our different models. Despite their different degrees of intended physical realism, we found that all three prescriptions can reproduce in an almost identical way the stellar abundance trends observed in Sculptor. While the three models have the same capacity to fit the data, the best values recovered for the parameters controlling the number of Type Ia supernovae and the strength of galactic outflows, are substantially different and in fact mutually exclusive from one model to another. For the purpose of understanding how a galaxy evolves, we conclude that only reproducing the evolution of a limited number of elements is insufficient and can lead to misleading conclusions. More elements or additional constraints such as the galaxy's star formation efficiency and the gas fraction are needed in order to break the degeneracy between the different modeling assumptions. Our results show that the successes and failures of chemical evolution models are predominantly driven by the input stellar yields, rather than by the complexity of the galaxy model itself. Simple models such as OMEGA are therefore sufficient to test and validate stellar yields.

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.

Over the last decade block-structured adaptive mesh refinement (SAMR) has found increasing use in large, publicly available codes and frameworks. SAMR frameworks have evolved along different paths. Some have stayed focused on specific domain areas, others have pursued a more general functionality, providing the building blocks for a larger variety of applications. In this survey paper we examine a representative set of SAMR packages and SAMR-based codes that have been in existence for half a decade or more, have a reasonably sized and active user base outside of their home institutions, and are publicly available. The set consists of a mix of SAMR packages and application codes that cover a broad range of scientific domains. We look at their high-level frameworks, and their approach to dealing with the advent of radical changes in hardware architecture. The codes included in this survey are BoxLib, Cactus, Chombo, Enzo, FLASH, and Uintah.

Due to their long mean free paths, X-rays are expected to have many significant impacts globally on the properties of the intergalactic medium (IGM) by their heating and ionizing processes on large scales. At high redshifts, X-rays from Population (Pop) III binaries might have important effects on cosmic reionization and the Lyman alpha forest. As a continuation of our previous work on Pop III binary X-rays (Xu et al. 2014), we use the Pop III distribution and evolution from the Renaissance Simulations, a suite of self-consistent cosmological radiation hydrodynamics simulations of the formation of the first galaxies, to calculate the X-ray luminosity density and background over the redshift range 20 > z > 7.6. As we find that Pop III star formation continues at a low, nearly constant rate to the end of reionization, X-rays are being continuously produced at significant rates compared to other possible X-ray sources, such as AGNs and normal X-ray binaries during the same period of time. We estimate that Pop III binaries produce approximately 6 eV of energy in the X-rays per hydrogen atom. We calculate the X-ray background for different monochromatic photon energies. KeV X-rays redshift and accumulate to produce a strong X-ray background spectrum extending to roughly 500 eV. The X-ray background is strong enough to heat the IGM to ~ 1000 K and to ionize a few percent of the neutral hydrogen. These effects are important for an understanding of the neutral hydrogen hyperfine transition 21-cm line signatures, the Ly alpha forest, and optical depth of the CMB to Thomson scattering.

Cosmic reionization is thought to be primarily fueled by the first generations of galaxies. We examine their stellar and gaseous properties, focusing on the star formation rates and the escape of ionizing photons, as a function of halo mass, redshift, and environment using the full suite of the {\it Renaissance Simulations} with an eye to provide better inputs to global reionization simulations. This suite, carried out with the adaptive mesh refinement code Enzo, is unprecedented in terms of their size and physical ingredients. The simulations probe overdense, average, and underdense regions of the universe of several hundred comoving Mpc$^3$, each yielding a sample of over 3,000 halos in the mass range $10^7 - 10^{9.5}~\Ms$ at their final redshifts of 15, 12.5, and 8, respectively. In the process, we simulate the effects of radiative and supernova feedback from 5,000 to 10,000 metal-free (Population III) stars in each simulation. We find that halos as small as $10^7~\Ms$ are able to form stars due to metal-line cooling from earlier enrichment by massive Population III stars. However, we find such halos do not form stars continuously. Using our large sample, we find that the galaxy-halo occupation fraction drops from unity at virial masses above $10^{8.5}~\Ms$ to $\sim$50\% at $10^8 ~\Ms$ and $\sim$10\% at $10^7~\Ms$, quite independent of redshift and region. Their average ionizing escape fraction is $\sim$5\% in the mass range $10^8 - 10^9~\Ms$ and increases with decreasing halo mass below this range, reaching 40--60\% at $10^7~\Ms$. Interestingly, we find that the escape fraction varies between 10--20\% in halos with virial masses $\sim 3 \times 10^9~\Ms$. Taken together, our results confirm the importance of the smallest galaxies as sources of ionizing radiation contributing to the reionization of the universe.

Galactic chemical evolution is a topic that involves nuclear physics, stellar evolution, galaxy evolution, observation, and cosmology. Continuous communication and feedback between these fields is a key component in improving our understanding of how galaxies form and how elements are created and recycled in galaxies and intergalactic space. In this proceedings, we present the current state of the JINA-NuGrid chemical evolution pipeline. It is designed to probe the impact of nuclear astrophysics uncertainties on galactic chemical evolution, to improve our knowledges regarding the origin of the elements in a cosmological context, and to create the required interdisciplinary connections.

We used a one-zone chemical evolution model to address the question of how many masses and metallicities are required in grids of massive stellar models in order to ensure reliable galactic chemical evolution predictions. We used a set of yields that includes seven masses between 13 and 30 Msun, 15 metallicities between 0 and 0.03 in mass fraction, and two different remnant mass prescriptions. We ran several simulations where we sampled subsets of stellar models to explore the impact of different grid resolutions. Stellar yields from low- and intermediate-mass stars and from Type Ia supernovae have been included in our simulations, but with a fixed grid resolution. We compared our results with the stellar abundances observed in the Milky Way for O, Na, Mg, Si, Ca, Ti, and Mn. Our results suggest that the range of metallicity considered is more important than the number of metallicities within that range, which only affects our numerical predictions by about 0.1 dex. We found that our predictions at [Fe/H] < -2 are very sensitive to the metallicity range and the mass sampling used for the lowest metallicity included in the set of yields. Variations between results can be as high as 0.8 dex, for any remnant mass prescription. At higher [Fe/H], we found that the required number of masses depends on the element of interest and on the remnant mass prescription. With a monotonic remnant mass prescription where every model explodes as a core-collapse supernova, the mass resolution induces variations of 0.2 dex on average. But with a remnant mass prescription that includes islands of non-explodability, the mass resolution can cause variations of about 0.2 to 0.7 dex depending on the choice of metallicity range. With such a prescription, explosive or non-explosive models can be missed if not enough masses are selected, resulting in over- or under-estimations of the mass ejected by massive stars.

Additional physics beyond standard hydrodynamics is needed to fully model the intracluster medium (ICM); however, as we move to more sophisticated models, it is important to consider the role of magnetic fields and the way the fluid approximation breaks down. This paper represents a first step towards developing a self-consistent model of the ICM by characterizing the statistical properties of magnetic fields in cosmological simulations of galaxy clusters. We find that plasma conditions are largely homogeneous across a range of cluster masses and relaxation states. We also find that the magnetic field length scales are resolution dependent and not based on any particular physical process. Energy transfer mechanisms and scales are also identified, and imply the existence of small scale dynamo action. The scales of the small scale dynamo are resolution limited and driven by numerical resistivity and viscosity.

Computational skills are required across all astronomy disciplines. Many students enter degree programs without sufficient skills to solve computational problems in their core classes or contribute immediately to research. We recommend advocacy for computational literacy, familiarity with fundamental software carpentry skills, and mastery of basic numerical methods by the completion of an undergraduate degree in Astronomy. We recommend the AAS Education Task Force advocate for a significant increase in computational literacy. We encourage the AAS to modestly fund efforts aimed at providing Open Education Resources (OER) that will significantly impact computational literacy in astronomy education.

We use a simple one-zone galactic chemical evolution model to quantify the uncertainties generated by the input parameters in numerical predictions, for a galaxy with properties similar to those of the Milky Way. We compiled several studies from the literature to gather the current constraints for our simulations regarding the typical value and uncertainty of seven basic parameters, which are: the lower and upper mass limits of the stellar initial mass function (IMF), the slope of the high-mass end of the stellar IMF, the slope of the delay-time distribution function of Type Ia supernovae (SNe Ia), the number of SNeIa per Msun formed, the total stellar mass formed, and the final mass of gas. We derived a probability distribution function to express the range of likely values for every parameter, which were then included in a Monte Carlo code to run several hundred simulations with randomly selected input parameters. This approach enables us to analyze the predicted chemical evolution of 16 elements in a statistical way by identifying the most probable solutions along with their 68% and 95% confidence levels. Our results show that the overall uncertainties are shaped by several input parameters that individually contribute at different metallicities, and thus at different galactic ages. The level of uncertainty then depends on the metallicity and is different from one element to another. Among the seven input parameters considered in this work, the slope of the IMF and the number of SNe Ia are currently the two main sources of uncertainty.

We present results on the formation of Pop III stars at redshift 7.6 from the Renaissance Simulations, a suite of extremely high-resolution and physics-rich radiation transport hydrodynamics cosmological adaptive-mesh refinement simulations of high redshift galaxy formation performed on the Blue Waters supercomputer. In a survey volume of about 220 comoving Mpc$^3$, we found 14 Pop III galaxies with recent star formation. The surprisingly late formation of Pop III stars is possible due to two factors: (i) the metal enrichment process is local and slow, leaving plenty of pristine gas to exist in the vast volume; and (ii) strong Lyman-Werner radiation from vigorous metal-enriched star formation in early galaxies suppresses Pop III formation in ("not so") small primordial halos with mass less than $\sim$ 3 $\times$ 10$^7$ M$_\odot$. We quantify the properties of these Pop III galaxies and their Pop III star formation environments. We look for analogues to the recently discovered luminous Ly $\alpha$ emitter CR7 (Sobral et al. 2015), which has been interpreted as a Pop III star cluster within or near a metal-enriched star forming galaxy. We find and discuss a system similar to this in some respects, however the Pop III star cluster is far less massive and luminous than CR7 is inferred to be.

This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 timescale is at least two orders of magnitude -- and in some cases greater -- than that available currently. 2) The growth rate of data produced by simulations is overwhelming the current ability, of both facilities and researchers, to store and analyze it. Additional resources and new techniques for data analysis are urgently needed. 3) Data rates and volumes from HEP experimental facilities are also straining the ability to store and analyze large and complex data volumes. Appropriately configured leadership-class facilities can play a transformational role in enabling scientific discovery from these datasets. 4) A close integration of HPC simulation and data analysis will aid greatly in interpreting results from HEP experiments. Such an integration will minimize data movement and facilitate interdependent workflows. 5) Long-range planning between HEP and ASCR will be required to meet HEP's research needs. To best use ASCR HPC resources the experimental HEP program needs a) an established long-term plan for access to ASCR computational and data resources, b) an ability to map workflows onto HPC resources, c) the ability for ASCR facilities to accommodate workflows run by collaborations that can have thousands of individual members, d) to transition codes to the next-generation HPC platforms that will be available at ASCR facilities, e) to build up and train a workforce capable of developing and using simulations and analysis to support HEP scientific research on next-generation systems.

This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade.

We compare several common sub-grid implementations of AGN feedback, focusing on the effects of different triggering mechanisms and the differences between thermal and kinetic feedback. Our main result is that pure thermal feedback that is centrally injected behaves differently from feedback with even a small kinetic component. Specifically, pure thermal feedback results in excessive condensation and smothering of the AGN by cold gas because the feedback energy does not propagate to large enough radii. We do not see large differences between implementations of different triggering mechanisms, as long as the spatial resolution is sufficiently high, probably because all of the implementations tested here trigger strong AGN feedback under similar conditions. In order to assess the role of resolution, we vary the size of the "accretion zone" in which properties are measured to determine the AGN accretion rate and resulting feedback power. We find that a larger accretion zone results in steadier jets but can also allow too much cold-gas condensation in simulations with a Bondi-like triggering algorithm. We also vary the opening angle of jet precession and find that a larger precession angle causes more of the jet energy to thermalize closer to the AGN, thereby producing results similar to pure thermal feedback. Our simulations confirm that AGN can regulate the thermal state of cool-core galaxy clusters and maintain the core in a state that is marginally susceptable to thermal instability and precipitation.

We present the largest number of Milky Way sized dark matter halos simulated at very high mass ($\sim$$10^4$ M$_\odot$/particle) and temporal resolution ($\sim$5 Myrs/snapshot) done to date, quadrupling what is currently available in the literature. This initial suite consists of the first 24 halos of the $Caterpillar$ $Project$ (www.caterpillarproject.org) whose project goal of 60 - 70 halos will be made public when complete. We resolve $\sim$20,000 gravitationally bound subhalos within the virial radius of each host halo. Over the ranges set by our spatial resolution our convergence is excellent and improvements were made upon current state-of-the-art halo finders to better identify substructure at such high resolutions (e.g., on average we recover $\sim$4 subhalos in each host halo above 10$^8$ M$_\odot$ which would have otherwise not been found using conventional methods). For our relaxed halos, the inner profiles are reasonably fit by Einasto profiles ($\alpha$ = 0.169 $\pm$ 0.023) though this depends on the relaxed nature and assembly history of a given halo. Averaging over all halos, the substructure mass fraction is $f_{m,subs} = 0.121 \pm 0.041$, and mass function slope is d$N$/d$M\propto M^{-1.88 \pm 0.10}$ though we find scatter in the normalizations for fixed halo mass due to more concentrated hosts having less subhalos at fixed subhalo mass. There are no biases stemming from Lagrangian volume selection as all Lagrangian volume types are included in our sample. Our detailed contamination study of 264 low resolution halos has resulted in obtaining very large and unprecedented, high-resolution regions around our host halos for our target resolution (sphere of radius $\sim$$1.4 \pm 0.4$ Mpc) allowing for accurate studies of low mass dwarf galaxies at large galactocentric radii and the very first stellar systems at high redshift ($z \geq$ 10).

This paper presents the first results from a model for chemical evolution that can be applied to N-body cosmological simulations and quantitatively compared to measured stellar abundances from large astronomical surveys. This model convolves the chemical yield sets from a range of stellar nucleosynthesis calculations (including AGB stars, Type Ia and II supernovae, and stellar wind models) with a user-specified stellar initial mass function (IMF) and metallicity to calculate the time-dependent chemical evolution model for a "simple stellar population" of uniform metallicity and formation time. These simple stellar population models are combined with a semi-analytic model for galaxy formation and evolution that uses merger trees from N-body cosmological simulations to track several $\alpha$- and iron-peak elements for the stellar and multiphase interstellar medium components of several thousand galaxies in the early ($z \geq 6$) universe. The simulated galaxy population is then quantitatively compared to two complementary datasets of abundances in the Milky Way stellar halo, and is capable of reproducing many of the observed abundance trends. The observed abundance ratio distributions are qualitatively well matched by our model, and the observational data is best reproduced with a Chabrier IMF, a chemically-enriched star formation efficiency of $0.2$, and a redshift of reionization of $7$.

We show that diffusion due to chaotic mixing in the Neighbourhood of the Sun may not be as relevant as previously suggested in erasing phase space signatures of past Galactic accretion events. For this purpose, we analyse Solar Neighbourhood-like volumes extracted from cosmological simulations that naturally account for chaotic orbital behaviour induced by the strongly triaxial and cuspy shape of the resulting dark matter haloes, among other factors. In the approximation of an analytical static triaxial model, our results show that a large fraction of stellar halo particles in such local volumes have chaos onset times (i.e., the timescale at which stars commonly associated with chaotic orbits will exhibit their chaotic behaviour) significantly larger than a Hubble time. Furthermore, particles that do present a chaotic behaviour within a Hubble time do not exhibit significant diffusion in phase space.

We present a simulation of the formation of the earliest Population II stars, starting from cosmological initial conditions and ending when metals created in the first supernovae are incorporated into a collapsing gas-cloud. This occurs after a supernova blast-wave collides with a nearby mini-halo, inducing further turbulence that efficiently mixes metals into the dense gas in the center of the halo. The gas that first collapses has been enriched to a metallicity of Z ~ 2e-5 Zsun. Due to the extremely low metallicity, collapse proceeds similarly to metal-free gas until dust cooling becomes efficient at high densities, causing the cloud to fragment into a large number of low mass objects. This external enrichment mechanism provides a plausible origin for the most metal-poor stars observed, such as SMSS J031300.36-670839.3, that appear to have formed out of gas enriched by a single supernova. This mechanism operates on shorter timescales than the time for low-mass mini-halos (M < 5e5 Msun) to recover their gas after experiencing a supernova. As such, metal-enriched stars will likely form first via this channel if the conditions are right for it to occur. We identify a number of other externally enriched halos that may form stars in this manner. These halos have metallicities as high as 0.01 Zsun, suggesting that some members of the first generation of metal-enriched stars may be hiding in plain sight in current stellar surveys.

In this paper, we present the first results from the Renaissance Simulations, a suite of extremely high-resolution and physics-rich AMR calculations of high redshift galaxy formation performed on the Blue Waters supercomputer. These simulations contain hundreds of well-resolved galaxies at $z \sim 25-8$, and make several novel, testable predictions. Most critically, we show that the ultraviolet luminosity function of our simulated galaxies is consistent with observations of high-z galaxy populations at the bright end of the luminosity function (M$_{1600} \leq -17$), but at lower luminosities is essentially flat rather than rising steeply, as has been inferred by Schechter function fits to high-z observations, and has a clearly-defined lower limit in UV luminosity. This behavior of the luminosity function is due to two factors: (i) the strong dependence of the star formation rate on halo virial mass in our simulated galaxy population, with lower-mass halos having systematically lower star formation rates and thus lower UV luminosities; and (ii) the fact that halos with virial masses below $\simeq 2 \times 10^8$ M$_\odot$ do not universally contain stars, with the fraction of halos containing stars dropping to zero at $\simeq 7 \times 10^6$ M$_\odot$. Finally, we show that the brightest of our simulated galaxies may be visible to current and future ultra-deep space-based surveys, particularly if lensed regions are chosen for observation.

Galaxy growth depends critically on the interplay between radiative cooling of cosmic gas and the resulting energetic feedback that cooling triggers. This interplay has proven exceedingly difficult to model, even with large supercomputer simulations, because of its complexity. Nevertheless, real galaxies are observed to obey simple scaling relations among their primary observable characteristics. Here we show that a generic emergent property of the interplay between cooling and feedback can explain the observed scaling relationships between a galaxy's stellar mass, its total mass, and its chemical enrichment level, as well as the relationship between the average orbital velocity of its stars and the mass of its central black hole. These relationships naturally result from any feedback mechanism that strongly heats a galaxy's circumgalactic gas in response to precipitation of colder clouds out of that gas, because feedback then suspends the gas in a marginally precipitating state.

Most of the massive elliptical galaxies in the universe stopped forming stars billions of years ago, even though plenty of hot gas remains available for star formation. Here we present compelling evidence indicating that quenching of star formation depends on both black-hole feedback and Type Ia supernova heating. We analyze Chandra X-ray observations of ten massive ellipticals, five with extended, potentially star-forming multiphase gas and five single-phase ellipticals with no star formation. The ratio of cooling time to freefall time at 1--10 kpc in the multiphase galaxies is tc/tff ~10, indicating that precipitation-driven feedback limits cooling but does not eliminate condensation. In the same region of the single-phase galaxies, the radial profiles of gas entropy are consistent with a thermally stable (tc/tff > 20) supernova-driven outflow that sweeps stellar ejecta out of the galaxy. However, in one of those single-phase ellipticals (NGC 4261) we find tc/tff < 10 at < 300 pc. Notably, its jets are ~50 times more powerful than in the other nine ellipticals, in agreement with models indicating that precipitation near the black hole should switch its fueling mode from Bondi-like accretion to cold chaotic accretion. We conclude by hypothesizing that particularly strong black-hole outbursts can shut off star formation in massive elliptical galaxies by boosting the entropy of the hot gas and flipping the system into the supernova-sweeping state.

Numerical simulations of active galactic nuclei (AGN) feedback in cool-core galaxy clusters have successfully avoided classical cooling flows, but often produce too much cold gas. We perform adaptive mesh simulations that include momentum-driven AGN feedback, self-gravity, star formation and stellar feedback, focusing on the interplay between cooling, AGN heating and star formation in an isolated cool-core cluster. Cold clumps triggered by AGN jets and turbulence form filamentary structures tens of kpc long. This cold gas feeds both star formation and the supermassive black hole (SMBH), triggering an AGN outburst that increases the entropy of the ICM and reduces its cooling rate. Within 1-2 Gyr, star formation completely consumes the cold gas, leading to a brief shutoff of the AGN. The ICM quickly cools and redevelops multiphase gas, followed by another cycle of star formation/AGN outburst. Within 6.5 Gyr, we observe three such cycles. There is good agreement between our simulated cluster and the observations of cool-core clusters. ICM cooling is dynamically balanced by AGN heating, and a cool-core appearance is preserved. The minimum cooling time to free-fall time ratio typically varies between a few and $\gtrsim 20$. The star formation rate (SFR) covers a wide range, from 0 to a few hundred $\rm M_{\odot}\, yr^{-1}$, with an average of $\sim 40 \,\rm M_{\odot}\, yr^{-1}$. The instantaneous SMBH accretion rate shows large variations on short timescales, but the average value correlates well with the SFR. Simulations without stellar feedback or self-gravity produce qualitatively similar results, but a lower SMBH feedback efficiency (0.1% compared to 1%) results in too many stars.

Heat input roughly balances radiative cooling in the gaseous cores of galaxy clusters even when the central cooling time is short, implying that cooling triggers a feedback loop that maintains thermal balance. Furthermore, cores with short cooling times tend to have multiphase structure, suggesting that the intracluster medium (ICM) becomes locally thermally unstable for cooling times < 1 Gyr. In this work, we use 2D and 3D hydrodynamic simulations to study the onset of condensation in idealized galaxy-cluster cores. In particular, we look at how the condensation process depends on the ratio of cooling time to freefall time and on the geometry of the gravitational potential. We conclude that the ICM can always evolve to a state in which condensation occurs if given enough time, but that an initial timescale ratio tcool /tff < 10 is needed for thermal instability to grow quickly enough to affect realistic cluster cores within a timescale that is relevant for cosmological structure formation. We find that instability leads to convection and that perturbations continue to grow while the gas convects. Condensation occurs when the timescale ratio in the low-entropy tail of the perturbation distribution drops below tcool /tff < 3, even if the volume-averaged timescale ratio is substantially greater. In our simulations, the geometry of the gravitational potential does not have a strong effect on thermal stability. Finally, we find that if condensation is powering feedback, a conversion efficiency of around 10^-3 for converting the condensed mass into thermal energy is sufficient to maintain thermal balance in the ICM.

Motivated by recent studies suggesting that the Large Magellanic Cloud (LMC) could be significantly more massive than previously thought, we explore whether the approximation of an inertial Galactocentric reference frame is still valid in the presence of such a massive LMC. We find that previous estimates of the LMC's orbital period and apocentric distance derived assuming a fixed Milky Way are significantly shortened for models where the Milky Way is allowed to move freely in response to the gravitational pull of the LMC. Holding other parameters fixed, the fraction of models favoring first infall is reduced. Due to this interaction, the Milky Way center of mass within the inner 50 kpc can be significantly displaced in phase-space in a very short period of time that ranges from 0.3 to 0.5 Gyr by as much as 30 kpc and 75 km/s. Furthermore, we show that the gravitational pull of the LMC and response of the Milky Way are likely to significantly affect the orbit and phase space distribution of tidal debris from the Sagittarius dwarf galaxy (Sgr). Such effects are larger than previous estimates based on the torque of the LMC alone. As a result, Sgr deposits debris in regions of the sky that are not aligned with the present-day Sgr orbital plane. In addition, we find that properly accounting for the movement of the Milky Way around its common center of mass with the LMC significantly modifies the angular distance between apocenters and tilts its orbital pole, alleviating tensions between previous models and observations. While these models are preliminary in nature, they highlight the central importance of accounting for the mutual gravitational interaction between the MW and LMC when modeling the kinematics of objects in the Milky Way and Local Group.

The first galaxies in the Universe are the building blocks of all observed galaxies. We present scaling relations for galaxies forming at redshifts $z \ge 15$ when reionization is just beginning. We utilize the ``Rarepeak' cosmological radiation hydrodynamics simulation that captures the complete star formation history in over 3,300 galaxies, starting with massive Population III stars that form in dark matter halos as small as ~$10^6 M_\odot$. We make various correlations between the bulk halo quantities, such as virial, gas, and stellar masses and metallicities and their respective accretion rates, quantifying a variety of properties of the first galaxies up to halo masses of $10^9 M_\odot$. Galaxy formation is not solely relegated to atomic cooling halos with virial temperatures greater than $10^4$ K, where we find a dichotomy in galaxy properties between halos above and below this critical mass scale. Halos below the atomic cooling limit have a stellar mass -- halo mass relationship $\log M_\star \simeq 3.5 + 1.3\log(M_{\rm vir} / 10^7 M_\odot)$. We find a non-monotonic relationship between metallicity and halo mass for the smallest galaxies. Their initial star formation events enrich the interstellar medium and subsequent star formation to a median of $10^{-2} Z_\odot$ and $10^{-1.5} Z_\odot$, respectively, in halos of total mass $10^7 M_\odot$ that is then diluted by metal-poor inflows, well beyond Population III pre-enrichment levels of $10^{-3.5} Z_\odot$. The scaling relations presented here can be employed in models of reionization, galaxy formation and chemical evolution in order to consider these galaxies forming prior to reionization.

Due to their long mean free path, X-rays are expected to have an important impact on cosmic reionization by heating and ionizing the IGM on large scales, especially after simulations have suggested that Population III stars may form in pairs at redshifts as high as 20 - 30. We use the Pop III distribution and evolution from a self-consistent cosmological radiation hydrodynamics simulation of the formation of the first galaxies and a simple Pop III X-ray binaries model to estimate their X-ray output in a high density region larger than 100 comoving (Mpc)$^3$. We then combine three different methods --- ray tracing, a one-zone model, and X-ray background modeling --- to investigate the X-ray propagation, intensity distribution, and long term effects on the IGM thermal and ionization state. The efficiency and morphology of photo-heating and photo-ionization are dependent on the photon energies. The sub-keV X-rays only impact the IGM near the sources while the keV photons contribute significantly to the X-ray background and heat and ionize the IGM smoothly. The X-rays just below 1 keV are most efficient in heating and ionizing the IGM. We find that the IGM might be heated to over 100 K by $z=10$ and the high density source region might reach 10$^4$ K, limited by atomic hydrogen cooling. This may be important for predicting the 21-cm neutral hydrogen signals. But, on the other hand, the free electrons from X-ray ionizations are not enough to contribute significantly to the optical depth of CMB to the Thomson scattering.

The methods by which one characterizes the distribution of matter in cosmological simulations is intrinsically different from how one performs the same task observationally. In this paper, we make substantial steps towards comparing simulations and observations of the intergalactic medium (IGM) in a more sensible way. We present a pipeline that generates and fits synthetic QSO absorption spectra using sight lines cast through a cosmological simulation, and simultaneously identifies structure by directly analyzing the variations in HI and OVI number density. We compare synthetic absorption spectra with a less observationally motivated, but more straightforward density threshold-based method for finding absorbers. Our efforts focus on HI and OVI to better characterize the warm/hot intergalactic medium, a subset of the IGM that is challenging to conclusively identify observationally. We find that the two methods trace roughly the same quantities of HI and OVI above observable column density limits, but the synthetic spectra typically identify more substructure in absorbers. We use both methods to characterize HI and OVI absorber properties. We find that both integrated and differential column density distributions from both methods generally agree with observations. The distribution of Doppler parameters between the two methods are similar for Lya and compare reasonably with observational results, but while the two methods agree with each other with OVI systems, they both are systematically different from observations. We find a strong correlation between OVI baryon fraction and OVI column density. We also discuss a possible bimodality in the temperature distribution of the gas traced by OVI.

[Abridged] We present an application of a statistical tool known as Sensitivity Analysis to characterize the relationship between input parameters and observational predictions of semi-analytic models of galaxy formation coupled to cosmological $N$-body simulations. We show how a sensitivity analysis can be performed on our chemo-dynamical model, ChemTreeN, to characterize and quantify its relationship between model input parameters and predicted observable properties. The result of this analysis provides the user with information about which parameters are most important and most likely to affect the prediction of a given observable. It can also be used to simplify models by identifying input parameters that have no effect on the outputs of interest. Conversely, it allows us to identify what model parameters can be most efficiently constrained by the given observational data set. We have applied this technique to real observational data sets associated with the Milky Way, such as the luminosity function of satellite galaxies. A statistical comparison of model outputs and real observables is used to obtain a "best-fitting" parameter set. We consider different Milky Way-like dark matter halos to account for the dependence of the best-fitting parameters selection process on the underlying merger history of the models. For all formation histories considered, running ChemTreeN with best-fitting parameters produced luminosity functions that tightly fit their observed counterpart. However, only one of the resulting stellar halo models was able to reproduce the observed stellar halo mass within 40 kpc of the Galactic center. On the basis of this analysis it is possible to disregard certain models, and their corresponding merger histories, as good representations of the underlying merger history of the Milky Way.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

This document presents (off-line) computing requrements and challenges for Cosmic Frontier science, covering the areas of data management, analysis, and simulations. We invite contributions to extend the range of covered topics and to enhance the current descriptions.

The first stars in the universe, termed Population III, are thought to have been very massive compared to the stars that form in the present epoch. As feedback from the first generation of stars altered the contents of the interstellar medium, the universe switched to a low-mass mode of star formation, which continues in the high metallicity stars formed in the present era. Several studies have investigated the transition between metal-free and metal-enriched star formation, with tentative evidence being found for a metallicity threshold near 10^-3.5 Z_sun due to atomic and molecular transitions and another threshold near 10^-5.5 Z_sun due to dust. In this work, we simulate the formation of stars in idealized low-metallicity halos using the AMR code Enzo. We conduct several simulations of 10^6 M_sun and 10^7 M_sun halos in which the metal content, initial rotation, and degree of turbulence are varied in order to study the effect of these properties on gas fragmentation over a range of densities. We find tentative support for the idea of a critical metallicity, but the effect of varying metallicity on the gas we observe is not as dramatic as what has been reported in earlier studies. We find no clear relation between the initial spin or the initial level of turbulence in the halo and the final properties of the gas contained therein. Additionally, we find that the degree to which the Jeans length is refined, the initial density profile of the gas, and the inclusion of deuterium chemistry each have a significant effect on the evolution and fragmentation of the gas in the halo - in particular, we find that at least 64 grid cells are needed to cover the Jeans length in order to properly resolve the fragmentation.

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

Simulations of galaxy clusters have a difficult time reproducing the radial gas-property gradients and red central galaxies observed to exist in the cores of galaxy clusters. Thermal conduction has been suggested as a mechanism that can help bring simulations of cluster cores into better alignment with observations by stabilizing the feedback processes that regulate gas cooling, but this idea has not yet been well tested with cosmological numerical simulations. Here we present cosmological simulations of ten galaxy clusters performed with five different levels of isotropic Spitzer conduction, which alters both the cores and outskirts of clusters, but not dramatically. In the cores, conduction flattens central temperature gradients, making them nearly isothermal and slightly lowering the central density but failing to prevent a cooling catastrophe there. Conduction has little effect on temperature gradients outside of cluster cores because outward conductive heat flow tends to inflate the outer parts of the intracluster medium (ICM) instead of raising its temperature. In general, conduction tends reduce temperature inhomogeneity in the ICM, but our simulations indicate that those homogenizing effects would be extremely difficult to observe in ~5 keV clusters. Outside the virial radius, our conduction implementation lowers the gas densities and temperatures because it reduces the Mach numbers of accretion shocks. We conclude that despite the numerous small ways in which conduction alters the structure of galaxy clusters, none of these effects are significant enough to make the efficiency of conduction easily measurable unless its effects are more pronounced in clusters hotter than those we have simulated.

We present a semi-analytic, computationally inexpensive model to identify halos capable of forming a Population III star in cosmological simulations across a wide range of times and environments. This allows for a much more complete and representative set of Population III star forming halos to be constructed, which will lead to Population III star formation simulations that more accurately reflect the diversity of Population III stars, both in time and halo mass. This model shows that Population III and chemically enriched stars coexist beyond the formation of the first generation of stars in a cosmological simulation until at least z~10, and likely beyond, though Population III stars form at rates that are 4-6 orders of magnitude lower than chemically enriched stars by z=10. A catalog of more than 40,000 candidate Population III forming halos were identified, with formation times temporally ranging from z=30 to z=10, and ranging in mass from 2.3x10^5 M_sun to 1.2x10^10 M_sun. At early times, the environment that Population III stars form in is very similar to that of halos hosting chemically enriched star formation. At later times Population III stars are found to form in low-density regions that are not yet chemically polluted due to a lack of previous star formation in the area. Population III star forming halos become increasingly spatially isolated from one another at later times, and are generally closer to halos hosting chemically enriched star formation than to another halo hosting Population III star formation by z~10.

Recently, Widrow and collaborators announced the discovery of vertical density waves in the Milky Way disk. Here we investigate a scenario where these waves were induced by the Sagittarius dwarf galaxy as it plunged through the Galaxy. Using numerical simulations, we find that the Sagittarius impact produces North-South asymmetries and vertical wave-like behavior that qualitatively agrees with what is observed. The extent to which vertical modes can radially penetrate into the disc, as well as their amplitudes, depend on the mass of the perturbing satellite. We show that the mean height of the disc is expected to vary more rapidly in the radial than in the azimuthal direction. If the observed vertical density asymmetry is indeed caused by vertical oscillations, we predict radial and azimuthal variations of the mean vertical velocity, correlating with the spatial structure. These variations can have amplitudes as large as 8 km/s.

Flux limited X-ray surveys of galaxy clusters show that clusters come in two roughly equally proportioned varieties: "cool core" clusters (CCs) and non-"cool core" clusters (NCCs). In previous work, we have demonstrated using cosmological $N$-body + Eulerian hydrodynamic simulations that NCCs are often consistent with early major mergers events that destroy embryonic CCs. In this paper we extend those results and conduct a series of simulationsusing different methods of gas cooling, and of energy and metal feedback from supernovae, where we attempt to produce a population of clusters with realistic central cooling times, entropies, and temperatures. We find that the use of metallicity-dependent gas cooling is essential to prevent early overcooling,and that adjusting the amount of energy and metal feedback can have a significant impact on observable X-ray quantities of the gas. We are able to produce clusters with more realistic central observable quantities than have previously been attained. However, there are still significant discrepancies between the simulated clusters and observations, which indicates that a different approach to simulating galaxies in clusters is needed. We conclude by looking towards a promising subgrid method of modeling galaxy feedback in clusters which may help to ameliorate the discrepancies between simulations and observations.

Non-thermal radio emission from cosmic ray electrons in the vicinity of merging galaxy clusters is an important tracer of cluster merger activity, and is the result of complex physical processes that involve magnetic fields, particle acceleration, gas dynamics, and radiation. In particular, objects known as radio relics are thought to be the result of shock-accelerated electrons that, when embedded in a magnetic field, emit synchrotron radiation in the radio wavelengths. In order to properly model this emission, we utilize the adaptive mesh refinement simulation of the magnetohydrodynamic evolution of a galaxy cluster from cosmological initial conditions. We locate shock fronts and apply models of cosmic ray electron acceleration that are then input into radio emission models. We have determined the thermodynamic properties of this radio-emitting plasma and constructed synthetic radio observations to compare to observed galaxy clusters. We find a significant dependence of the observed morphology and radio relic properties on the viewing angle of the cluster, raising concerns regarding the interpretation of observed radio features in clusters. We also find that a given shock should not be characterized by a single Mach number. We find that the bulk of the radio emission comes from gas with T>5x10^7, \rho~10^(-28)-10^(-27) g/cm^3, with magnetic field strengths of 0.1-1.0 \mu G and shock Mach numbers of M~3-6. We present an analysis of the radio spectral index which suggests that the spatial variation of the spectral index can mimic synchrotron aging. Finally, we examine the polarization fraction and position angle of the simulated radio features, and compare to observations.

We use the semi-analytic model ChemTreeN, coupled to cosmological N-body simulations, to explore how different galaxy formation histories can affect observational properties of Milky Way-like galaxies' stellar haloes and their satellite populations. Gaussian processes are used to generate model emulators that allow one to statistically estimate a desired set of model outputs at any location of a p-dimensional input parameter space. This enables one to explore the full input parameter space orders of magnitude faster than could be done otherwise. Using mock observational data sets generated by ChemTreeN itself, we show that it is possible to successfully recover the input parameter vectors used to generate the mock observables if the merger history of the host halo is known. However, our results indicate that for a given observational data set the determination of "best fit" parameters is highly susceptible to the particular merger history of the host. Very different halo merger histories can reproduce the same observational dataset, if the "best fit" parameters are allowed to vary from history to history. Thus, attempts to characterize the formation history of the Milky Way using these kind of techniques must be performed statistically, analyzing large samples of high resolution N-body simulations.

It is now known that minor mergers are capable of creating structure in the phase-space distribution of their host galaxy's disc. In order to search for such imprints in the Milky Way, we analyse the SEGUE F/G-dwarf and the Schuster et al. (2006) stellar samples. We find similar features in these two completely independent stellar samples, consistent with the predictions of a Milky Way minor-merger event. We next apply the same analyses to high-resolution, idealised N-body simulations of the interaction between the Sagittarius dwarf galaxy and the Milky Way. The energy distributions of stellar particle samples in small spatial regions in the host disc reveal strong variations of structure with position. We find good matches to the observations for models with a mass of Sagittarius' dark matter halo progenitor $\lessapprox 10^{11}$ M$_{\odot}$. Thus, we show that this kind of analysis could be used to provide unprecedentedly tight constraints on Sagittarius' orbital parameters, as well as place a lower limit on its mass.

By means of N-body simulations we study the response of a galactic disc to a minor merger event. We find that non-self-gravitating, spiral-like features are induced in the thick disc. As we have shown in a previous work, this "ringing" also leaves an imprint in velocity space (the u-v plane) in small spatial regions, such as the solar neighbourhood. As the disc relaxes after the event, clumps in the u-v plane get closer with time, allowing us to estimate the time of impact. In addition to confirming the possibility of this diagnostic, here we show that in a more realistic scenario, the in-fall trajectory of the perturber gives rise to an azimuthal dependence of the structure in phase-space. We also find that the space defined by the energy and angular momentum of stars is a better choice than velocity space, as clumps remain visible even in large local volumes. This makes their observational detection much easier since one need not be restricted to a small spatial volume. We show that information about the time of impact, the mass of the perturber, and its trajectory is stored in the kinematics of disc stars.

Cosmological shocks are a critical part of large-scale structure formation, and are responsible for heating the intracluster medium in galaxy clusters. In addition, they are also capable of accelerating non-thermal electrons and protons. In this work, we focus on the acceleration of electrons at shock fronts, which is thought to be responsible for radio relics - extended radio features in the vicinity of merging galaxy clusters. By combining high resolution AMR/N-body cosmological simulations with an accurate shock finding algorithm and a model for electron acceleration, we calculate the expected synchrotron emission resulting from cosmological structure formation. We produce synthetic radio maps of a large sample of galaxy clusters and present luminosity functions and scaling relationships. With upcoming long wavelength radio telescopes, we expect to see an abundance of radio emission associated with merger shocks in the intracluster medium. By producing observationally motivated statistics, we provide predictions that can be compared with observations to further improve our understanding of magnetic fields and electron shock acceleration.

We perform a series of cosmological simulations using Enzo, an Eulerian adaptive-mesh refinement, N-body + hydrodynamical code, applied to study the warm/hot intergalactic medium. The WHIM may be an important component of the baryons missing observationally at low redshift. We investigate the dependence of the global star formation rate and mass fraction in various baryonic phases on spatial resolution and methods of incorporating stellar feedback. Although both resolution and feedback significantly affect the total mass in the WHIM, all of our simulations find that the WHIM fraction peaks at z ~ 0.5, declining to 35-40% at z = 0. We construct samples of synthetic OVI absorption lines from our highest-resolution simulations, using several models of oxygen ionization balance. Models that include both collisional ionization and photoionization provide excellent fits to the observed number density of absorbers per unit redshift over the full range of column densities (10^13 cm-2 <= N_OVI <= 10^15 cm^-2). Models that include only collisional ionization provide better fits for high column density absorbers (N_OVI > 10^14 cm^-2). The distribution of OVI in density and temperature exhibits two populations: one at T ~ 10^5.5 K (collisionally ionized, 55% of total OVI) and one at T ~ 10^4.5 K (photoionized, 37%) with the remainder located in dense gas near galaxies. While not a perfect tracer of hot gas, OVI provides an important tool for a WHIM baryon census.

We examine the incidence of cold fronts in a large sample of galaxy clusters extracted from a (512h^-1 Mpc) hydrodynamic/N-body cosmological simulation with adiabatic gas physics computed with the Enzo adaptive mesh refinement code. This simulation contains a sample of roughly 4000 galaxy clusters with M > 10^14 M_sun at z=0. For each simulated galaxy cluster, we have created mock 0.3-8.0 keV X-ray observations and spectroscopic-like temperature maps. We have searched these maps with a new automated algorithm to identify the presence of cold fronts in projection. Using a threshold of a minimum of 10 cold front pixels in our images, corresponding to a total comoving length L_cf > 156h^-1 kpc, we find that roughly 10-12% of all projections in a mass-limited sample would be classified as cold front clusters. Interestingly, the fraction of clusters with extended cold front features in our synthetic maps of a mass-limited sample trends only weakly with redshift out to z=1.0. However, when using different selection functions, including a simulated flux limit, the trending with redshift changes significantly. The likelihood of finding cold fronts in the simulated clusters in our sample is a strong function of cluster mass. In clusters with M>7.5x10^14 M_sun the cold front fraction is 40-50%. We also show that the presence of cold fronts is strongly correlated with disturbed morphology as measured by quantitative structure measures. Finally, we find that the incidence of cold fronts in the simulated cluster images is strongly dependent on baryonic physics.

We present an overview of recent work that focuses on understanding the radiative feedback processes that are potentially important during Population III star formation. Specifically, we examine the effect of the Lyman-Werner (photodissociating) background on the early stages of primordial star formation, which serves to delay the onset of star formation in a given halo but never suppresses it entirely. We also examine the effect that both photodissociating and ionizing radiation in I-fronts from nearby stellar systems have on the formation of primordial protostellar clouds. Depending on the strength of the incoming radiation field and the central density of the halos, Pop III star formation can be suppressed, unaffected, or even enhanced. Understanding these and other effects is crucial to modeling Population III star formation and to building the earliest generations of galaxies in the Universe.

We review theoretical models of Population III.1 star formation, focusing on the protostellar feedback processes that are expected to terminate accretion and thus set the mass of these stars. We discuss how dark matter annihilation may modify this standard feedback scenario. Then, under the assumption that dark matter annihilation is unimportant, we predict the mass of stars forming in 12 cosmological minihalos produced in independent numerical simulations. This allows us to make a simple estimate of the Pop III.1 initial mass function and how it may evolve with redshift.

Recently, Suzaku has produced temperature and entropy profiles, along with profiles of gas density, gas fraction, and mass, for multiple galaxy clusters out to ~r_200 (~= virial radius). In this paper, we compare these novel X-ray observations with results from N-body + hydrodynamic adaptive mesh refinement cosmological simulations using the Enzo code. There is excellent agreement in the temperature, density, and entropy profiles between a sample of 27 mostly substructure-free massive clusters in the simulated volume and the observed clusters. This supports our previous contention that clusters have "universal" outer temperature profiles. Furthermore, it appears that the simplest adiabatic gas physics used in these Enzo simulations is adequate to model the outer regions of these clusters without other mechanisms (e.g., non-gravitational heating, cooling, magnetic fields, or cosmic rays). However, the outskirts of these clusters are not in hydrostatic equilibrium. There is significant bulk flow and turbulence in the outer intracluster medium created by accretion from filaments. Thus, the gas is not fully supported by thermal pressure. The implications for mass estimation from X-ray data are discussed.

GRB 080913 and GRB 090423 are the most distant gamma-ray bursts (GRBs) known to-date, with spectroscopically determined redshifts of z=6.7 and z=8.1, respectively. The detection of bursts at this early epoch of the Universe significantly constrains the nature of GRBs and their progenitors. We perform population synthesis studies of the formation and evolution of early stars, and calculate the resulting formation rates of short and long-duration GRBs at high redshift. The peak of the GRB rate from Population II stars occurs at z=7 for a model with efficient/fast mixing of metals, while it is found at z=3 for an inefficient/slow metallicity evolution model. We show that in the redshift range 6<z<10 essentially all GRBs originate from Population II stars, regardless of metallicity evolution model. These stars (having small, but non-zero metallicity) are the most likely progenitors for both long GRBs (collapsars) and short GRBs (NS-NS or BH-NS mergers) at this epoch. Although the predicted intrinsic rates of long and short GRBs are similar at these high redshifts, observational selection effects lead to higher (factor of 10) observed rates for long GRBs. We conclude that the two recently observed high-z GRB events are most likely long GRBs originating from Population II collapsars.

Previous high resolution cosmological simulations predict the first stars to appear in the early universe to be very massive and to form in isolation. Here we discuss a cosmological simulation in which the central 50 solar mass clump breaks up into two cores, having a mass ratio of two to one, with one fragment collapsing to densities of 10^{-8} g/cc. The second fragment, at a distance of 800 astronomical units, is also optically thick to its own cooling radiation from molecular hydrogen lines, but is still able to cool via collision-induced emission. The two dense peaks will continue to accrete from the surrounding cold gas reservoir over a period of 10^5 years and will likely form a binary star system.

Detection of the Warm-Hot Intergalactic Medium (WHIM) using Sunyaev-Zeldovich effect (SZE) surveys is an intriguing possibility, and one that may allow observers to quantify the amount of "missing baryons" in the WHIM phase. We estimate the necessary sensitivity for detecting low density WHIM gas with the South Pole Telescope (SPT) and Planck Surveyor for a synthetic 100 square degree sky survey. This survey is generated from a very large, high dynamic range adaptive mesh refinement cosmological simulation performed with the Enzo code. We find that for a modest increase in the SPT survey sensitivity (a factor of 2-4), the WHIM gas makes a detectable contribution to the integrated sky signal. For a Planck-like satellite, similar detections are possible with a more significant increase in sensitivity (a factor of 8-10). We point out that for the WHIM gas, the kinematic SZE signal can sometimes dominate the thermal SZE where the thermal SZE decrement is maximal (150 GHz), and that using the combination of the two increases the chance of WHIM detection using SZE surveys. However, we find no evidence of unique features in the thermal SZE angular power spectrum that may aid in its detection. Interestingly, there are differences in the power spectrum of the kinematic SZE, which may not allow us to detect the WHIM directly, but could be an important contaminant in cosmological analyses of the kSZE-derived velocity field. Corrections derived from numerical simulations may be necessary to account for this contamination.

The cosmic dark ages are the mysterious epoch during which the pristine gas began to condense and ultimately form the first stars. Although these beginnings have long been a topic of theoretical interest, technology has only recently allowed the beginnings of observational insight into this epoch. Many questions surround the formation of stars in metal-free gas and the history of the build-up of metals in the intergalactic medium: (1) What were the properties of the first stellar and galactic sources to form in pristine (metal-free) gas? (2) When did the epoch of Population III (metal-free) star formation take place and how long did it last? (3) Was the stellar initial mass function dramatically different for the first stars and galaxies? These questions are all active areas of theoretical research. However, new observational constraints via the direct detection of Population III star formation are vital to making progress in answering the broader questions surrounding how galaxies formed and how the cosmological properties of the universe have affected the objects it contains.

This white paper, directed to the Stars and Stellar Evolution panel, has three objectives: 1) to provide the Astro2010 Decadal Survey with a vista into the goals of the nuclear physics and nuclear astrophysics community; 2) to alert the astronomical community of joint opportunities for discoveries at the interface between nuclear physics and astronomy; and 3) to delineate efforts in nuclear physics and describe the observational and theoretical advances in astrophysics necessary to make progress towards answering the following questions in the Nuclear Science 2007 Long Range Plan: 1) What is the origin and distribution of the elements? 2) What are the nuclear reactions that power stars and stellar explosions? 3) What is the nature of dense matter? The scope of this white paper concerns the specific area of "low energy" nuclear astrophysics. We define this as the area of overlap between astrophysics and the study of nuclear structure and reactions. Of the questions listed above, two -- What is the origin of the elements? and What is the nature of dense matter? -- were specifically listed in the National Academies Study, "Connecting Quarks with the Cosmos".

Simulations of the formation of Population III (Pop III) stars suggest that they were much more massive than the Pop II and Pop I stars observed today. This is due to the collapse dynamics of metal-free gas, which is regulated by the radiative cooling of molecular hydrogen. We study how the collapse of gas clouds is altered by the addition of metals to the star-forming environment by performing a series of simulations of pre-enriched star formation at various metallicities. For metallicities below the critical metallicity, Z_cr, collapse proceeds similarly to the metal-free case, and only massive objects form. For metallicities well above Z_cr, efficient cooling rapidly lowers the gas temperature to the temperature of the CMB. The gas is unable to radiatively cool below the CMB temperature, and becomes thermally stable. For high metallicities, Z >= 10^-2.5 Zsun, this occurs early in the evolution of the gas cloud, when the density is still relatively low. The resulting cloud-cores show little or no fragmentation, and would most likely form massive stars. If the metallicity is not vastly above Z_cr, the cloud cools efficiently but does not reach the CMB temperature, and fragmentation into multiple objects occurs. We conclude that there were three distinct modes of star formation at high redshift (z >= 4): a `primordial' mode, producing massive stars (10s to 100s Msun) at very low metallicities (Z <= 10^-3.75 Zsun); a CMB-regulated mode, producing moderate mass (10s of Msun) stars at high metallicites (Z >= 10^-2.5 Zsun at redshift z ~ 15-20); and a low-mass (a few Msun) mode existing between those two metallicities. As the universe ages and the CMB temperature decreases, the range of the low mass mode extends to higher metallicities, eventually becoming the only mode of star formation. (Abridged)

We present new results characterizing cosmological shocks within adaptive mesh refinement N-Body/hydrodynamic simulations that are used to predict non-thermal components of large-scale structure. This represents the first study of shocks using adaptive mesh refinement. We propose a modified algorithm for finding shocks from those used on unigrid simulations that reduces the shock frequency of low Mach number shocks by a factor of ~3. We then apply our new technique to a large, (512 Mpc/h)^3, cosmological volume and study the shock Mach number (M) distribution as a function of pre-shock temperature, density, and redshift. Because of the large volume of the simulation, we have superb statistics that results from having thousands of galaxy clusters. We find that the Mach number evolution can be interpreted as a method to visualize large-scale structure formation. Shocks with Mach<5 typically trace mergers and complex flows, while 5<Mach<20 and Mach>20 generally follow accretion onto filaments and galaxy clusters, respectively. By applying results from nonlinear diffusive shock acceleration models using the first-order Fermi process, we calculate the amount of kinetic energy that is converted into cosmic ray protons. The acceleration of cosmic ray protons is large enough that in order to use galaxy clusters as cosmological probes, the dynamic response of the gas to the cosmic rays must be included in future numerical simulations.

Current numerical studies suggest that the first galaxies formed a few stars at a time and were enriched only gradually by the first heavy elements. However, the large box sizes in these models cannot resolve primordial supernova explosions or the mixing of their metals with ambient gas, which could result in intervening, prompt generations of low-mass stars. We present multiscale 1D models of Population III supernovae in cosmological halos that evolve the blast from its earliest stages as a free expansion. We find that if the star ionizes the halo, the ejecta strongly interacts with the dense shell swept up by the H II region, potentially cooling and fragmenting it into clumps that are unstable to gravitational collapse. If the star fails to ionize the halo, the explosion propagates metals out to 20 - 40 pc and then collapses, enriching tens of thousands of solar masses of primordial gas, in contrast to previous models that suggest such explosions 'fizzle'. Rapid formation of low-mass stars trapped in the gravitational potential well of the halo appears to be inevitable in such circumstances. Consequently, it is possible that far more stars were swept up into the first galaxies, at earlier times and with distinct chemical signatures, than in present models. Upcoming measurements by JWST and ALMA may discriminate between these two paradigms.

We investigate the effects of weakly-interacting massive particle (WIMP) dark matter annihilation on the formation of Population III.1 stars, which are theorized to form from the collapse of gas cores at the centers of dark matter minihalos. We consider the relative importance of cooling due to baryonic radiative processes and heating due to WIMP annihilation. We analyze the dark matter and gas profiles of several halos formed in cosmological-scale numerical simulations. The heating rate depends sensitively on the dark matter density profile, which we approximate with a power law rho_chi ~ r^{-alpha_chi}, in the numerically unresolved inner regions of the halo. If we assume a self-similar structure so that alpha_chi ~= 1.5 as measured on the resolved scales ~1pc, then for a fiducial WIMP mass of 100GeV, the heating rate is typically much smaller (<10^{-3}) than the cooling rate for densities up to n_H=10^{17}cm^{-3}. In one case, where alpha_chi=1.65, the heating rate becomes similar to the cooling rate by a density of n_H=10^{15}cm^{-3}. The dark matter density profile is expected to steepen in the central baryon-dominated region <~1pc due to adiabatic contraction, and we observe this effect (though with relatively low resolution) in our numerical models. From these we estimate alpha_chi~=2.0. The heating now dominates cooling above n_H~=10^{14}cm^{-3}, in agreement with the previous study of Spolyar, Freese & Gondolo. We expect this leads to the formation of an equilibrium structure with a baryonic and dark matter density distribution exhibiting a flattened central core. Examining such equilibria, we find total luminosities due to WIMP annihilation are relatively constant and ~10^3 L_sun, set by the radiative luminosity of the baryonic core. We discuss the implications for Pop III.1 star formation... (abridged)

We report the results of the first self-consistent three-dimensional adaptive mesh refinement magnetohydrodynamical simulations of Population III star formation including the Biermann Battery effect. We find that the Population III stars formed including this effect are both qualitatively and quantitatively similar to those from hydrodynamics-only (non-MHD) cosmological simulations. We observe peak magnetic fields of ~10^-9 G in the center of our star-forming halo at z ~ 17.55. The magnetic fields created by the Biermann Battery effect are predominantly formed early in the evolution of the primordial halo at low density and large spatial scales, and then grow through compression and by shear flows. The fields seen in this calculation are never large enough to be dynamically important (with beta >= 10^{15} at all times), and should be considered the minimum possible fields in existence during Population III star formation, and may be seed fields for the stellar dynamo or the magnetorotational instability at higher densities and smaller spatial scales.

We present numerical simulations of primordial supernovae in cosmological minihalos at $z \sim$ 20. We consider Type II supernovae, hypernovae, and pair instability supernovae (PISN) in halos from 6.9 $\times$ 10$^5$ - 1.2 $\times$ 10$^7$ $\Ms$, those in which Population III stars are expected to form via H$_2$ cooling. The supernovae evolve along two evolutionary paths according to whether they explode in \ion{H}{2} regions or neutral halos. Those in \ion{H}{2} regions first expand adiabatically and then radiate strongly upon collision with baryons ejected from the halo during its photoevaporation by the progenitor. Explosions in neutral halos promptly emit most of their kinetic energy as x-rays, but retain enough momentum to seriously disrupt the halo. We find that the least energetic of the supernovae are capable of destroying halos $\lesssim$ 10$^7$ $\Ms$, while a single PISN can destroy even more massive halos. Blasts in \ion{H}{2} regions disperse heavy elements into the IGM, but neutral halos confine the explosion and its metals. In \ion{H}{2} regions, a prompt second generation of stars may form in the remnant at radii of 100 - 200 pc in the halo. Explosions confined by large halos instead recollapse, with infall rates in excess of 10$^{-2}$ $\Ms$ yr$^{-1}$ that heavily contaminate their interior. This fallback may either fuel massive black hole growth at very high redshifts or create the first globular cluster with a radius of 10 - 20 pc at the center of the halo. Our findings allow the possibility that the first primitive galaxies formed sooner, with greater numbers of stars and distinct chemical abundance patterns, than in current models.

We present numerical simulations of how a 120 M$_\odot$ primordial star regulates star formation in nearby cosmological halos at $z \sim$ 20 by photoevaporation. Our models include nine-species primordial chemistry and self-consistent multifrequency conservative transfer of UV photons with all relevant radiative processes. Whether or not new stars form in halos clustered around a Population III star ultimately depends on their core densities and proximity to the star. Diffuse halos with central densities below 2 - 3 cm$^{-3}$ are completely ionized and evaporated anywhere in the cluster. Evolved halos with core densities above 2000 cm$^{-3}$ are impervious to both ionizing and Lyman-Werner flux at most distances from the star and collapse as quickly as they would in its absence. Star formation in halos of intermediate density can be either promoted or suppressed depending on how the I-front remnant shock compresses, deforms and enriches the core with H$_2$. We find that the 120 M$_\odot$ star photodissociates H$_2$ in most halos in the cluster but that catalysis by H- restores it a few hundred kyr after the death of the star, with little effect on star formation. Our models exhibit significant departures from previous one-dimensional spherically-symmetric simulations, which are prone to serious errors due to unphysical geometric focusing effects.

The understanding of the formation, life, and death of Population III stars, as well as the impact that these objects had on later generations of structure formation, is one of the foremost issues in modern cosmological research and has been an active area of research during the past several years. We summarize the results presented at "First Stars III," a conference sponsored by Los Alamos National Laboratory, the Kavli Institute for Particle Astrophysics and Cosmology, and the Joint Institute for Nuclear Astrophysics. This conference, the third in a series, took place in July 2007 at the La Fonda Hotel in Santa Fe, New Mexico, U.S.A.

We examine aspects of primordial star formation in the presence of a molecular hydrogen-dissociating ultraviolet background. We compare a set of AMR hydrodynamic cosmological simulations using a single cosmological realization but with a range of ultraviolet background strengths in the Lyman-Werner band. This allows us to study the effects of Lyman-Werner radiation on suppressing H2 cooling at low densities as well as the high-density evolution of the collapsing core in a self-consistent cosmological framework. We find that the addition of a photodissociating background results in a delay of the collapse of high density gas at the center of the most massive halo in the simulation and, as a result, an increase in the virial mass of this halo at the onset of baryon collapse. We find that, contrary to previous results, Population III star formation is not suppressed for J$_{21} \geq 0.1$, but occurs even with backgrounds as high as J$_{21} = 1$. We find that H2 cooling leads to collapse despite the depressed core molecular hydrogen fractions due to the elevated H2 cooling rates at $T=2-5 \times 10^3$ K. We observe a relationship between the strength of the photodissociating background and the rate of accretion onto the evolving protostellar cloud core, with higher LW background fluxes resulting in higher accretion rates. Finally, we find that the collapsing halo cores in our simulations do not fragment at densities below $n \sim 10^{10}$ cm$^{-3}$ regardless of the strength of the LW background, suggesting that Population III stars forming in halos with T$_{vir} \sim 10^4$ K may still form in isolation.

We present the first results from a new generation of simulated large sky coverage (~100 square degrees) Sunyaev-Zeldovich effect (SZE) cluster surveys using the cosmological adaptive mesh refinement N-body/hydro code Enzo. We have simulated a very large (512^3h^{-3}Mpc^3) volume with unprecedented dynamic range. We have generated simulated light cones to match the resolution and sensitivity of current and future SZE instruments. Unlike many previous studies of this type, our simulation includes unbound gas, where an appreciable fraction of the baryons in the universe reside. We have found that cluster line-of-sight overlap may be a significant issue in upcoming single-dish SZE surveys. Smaller beam surveys (~1 arcmin) have more than one massive cluster within a beam diameter 5-10% of the time, and a larger beam experiment like Planck has multiple clusters per beam 60% of the time. We explore the contribution of unresolved halos and unbound gas to the SZE signature at the maximum decrement. We find that there is a contribution from gas outside clusters of ~16% per object on average for upcoming surveys. This adds both bias and scatter to the deduced value of the integrated SZE, adding difficulty in accurately calibrating a cluster Y-M relationship. Finally, we find that in images where objects with M > 5x10^{13} M_{\odot} have had their SZE signatures removed, roughly a third of the total SZE flux still remains. This gas exists at least partially in the Warm Hot Intergalactic Medium (WHIM), and will possibly be detectable with the upcoming generation of SZE surveys.

We present numerical simulations of the photoevaporation of cosmological halos clustered around a 120 M$_\odot$ primordial star, confining our study to structures capable of hosting Population III star formation. The calculations include self-consistent multifrequency conservative transfer of UV photons together with nine-species primordial chemistry and all relevant radiative processes. The ultimate fates of these halos varies with central density and proximity to the central source but generally fall into one of four categories. Diffuse halos with central densities below 2 - 3 cm$^{-3}$ are completely ionized and evaporated by the central star anywhere in the cluster. More evolved halo cores at densities above 2000 cm$^{-3}$ are impervious to both ionizing and Lyman-Werner flux at most distances from the star and collapse of their cores proceeds without delay. Radiative feedback in halos of intermediate density can be either positive or negative, depending on how the I-front remnant shock both compresses and deforms the core and enriches it with H$_2$. We find that the 120 M$_\odot$ star photodissociates H$_2$ in most halos within the cluster but that catalysis by H- rapidly restores molecular hydrogen within a few hundred Kyr after the death of the star, with little delay in star formation. Our models exhibit significant departures from previous one-dimensional spherically-symmetric simulations, which are prone to serious errors due to unphysical geometric focusing effects.

Current and upcoming cosmological observations allow us to probe structures on smaller and smaller scales, entering highly nonlinear regimes. In order to obtain theoretical predictions in these regimes, large cosmological simulations have to be carried out. The promised high accuracy from observations make the simulation task very demanding: the simulations have to be at least as accurate as the observations. This requirement can only be fulfilled by carrying out an extensive code validation program. The first step of such a program is the comparison of different cosmology codes including gravitation interactions only. In this paper we extend a recently carried out code comparison project to include five more simulation codes. We restrict our analysis to a small cosmological volume which allows us to investigate properties of halos. For the matter power spectrum and the mass function, the previous results hold, with the codes agreeing at the 10% level over wide dynamic ranges. We extend our analysis to the comparison of halo profiles and investigate the halo count as a function of local density. We introduce and discuss ParaView as a flexible analysis tool for cosmological simulations, the use of which immensely simplifies the code comparison task.

In this paper we describe our massively parallel version of Enzo, a multiphysics, parallel, AMR application for simulating cosmological structure formation developed at UCSD and Columbia. We describe its physics, numerical algorithms, implementation, and performance on current terascale platforms. We also discuss our future plans and some of the challenges we face as we move to the petascale.

We present preliminary results from simulated large sky coverage (~100 square degrees) Sunyaev-Zeldovich effect (SZE) cluster surveys using the cosmological adaptive mesh refinement N-body/hydro code Enzo. We have generated simulated light cones to match the resolution and sensitivity of current and future SZE instruments. These simulations are the most advanced calculations of their kind. The simulated sky surveys allow a direct comparison of large N-body/hydro cosmological simulations to current and pending sky surveys. Our synthetic surveys provide an indispensable guide for observers in the interpretation of large area sky surveys, and will develop the tools necessary to discriminate between models for cluster baryonic physics, and to accurately determine cosmological parameters.

(abridged) We perform 12 extremely high resolution adaptive mesh refinement cosmological hydrodynamic simulations of Population III star formation in a Lambda CDM universe, varying the box size and large-scale structure, to understand systematic effects in the formation of primordial protostellar cores. We find results that are qualitatively similar to those observed previously. We observe that the threshold halo mass for formation of a Population III protostar does not evolve significantly with time in the redshift range studied (33 > z > 19) but exhibits substantial scatter due to different halo assembly histories: Halos which assembled more slowly develop cooling cores at lower mass than those that assemble more rapidly, in agreement with Yoshida et al. (2003). We do, however, observe significant evolution in the accretion rates of Population III protostars with redshift, with objects that form later having higher maximum accretion rates, with a variation of two orders of magnitude (10^-4 - 10^-2 Msolar/year). This can be explained by considering the evolving virial properties of the halos with redshift and the physics of molecular hydrogen formation at low densities. Our result implies that the mass distribution of Population III stars inferred from their accretion rates may be broader than previously thought, and may evolve with redshift. Finally, we observe that our collapsing protostellar cloud cores do not fragment, consistent with previous results, which suggests that Population III stars which form in halos of mass 10^5 - 10^6 Msun always form in isolation.

In this paper we examine aspects of primordial star formation in a gravitino warm dark matter universe with a cosmological constant. We compare a set of simulations using a single cosmological realization but with a wide range of warm dark matter particle masses which have not yet been conclusively ruled out by observations. The addition of a warm dark matter component to the initial power spectrum results in a delay in the collapse of high density gas at the center of the most massive halo in the simulation and, as a result, an increase in the virial mass of this halo at the onset of baryon collapse. Both of these effects become more pronounced as the warm dark matter particle mass becomes smaller. A cosmology using a gravitino warm dark matter power spectrum assuming a particle mass of m_{WDM} ~ 40keV is effectively indistinguishable from the cold dark matter case, whereas the m_{WDM} ~ 15 keV case delays star formation by approx. 10^8 years. There is remarkably little scatter between simulations in the final properties of the primordial protostar which forms at the center of the halo, possibly due to the overall low rate of halo mergers which is a result of the WDM power spectrum. The detailed evolution of the collapsing halo core in two representative WDM cosmologies is described. At low densities (n_{b} <= 10^5 cm^{-3}), the evolution of the two calculations is qualitatively similar, but occurs on significantly different timescales, with the halo in the lower particle mass calculation taking much longer to evolve over the same density range and reach runaway collapse. Once the gas in the center of the halo reaches relatively high densities (n_{b} >= 10^5 cm^{-3}) the overall evolution is essentially identical in the two calculations.

Scattering of the cosmic microwave background (CMB) in clusters of galaxies polarizes the radiation. We explore several polarization components which have their origin in the kinematic quadrupole moments induced by the motion of the scattering electrons, either directed or random. Polarization levels and patterns are determined in a cluster simulated by the hydrodynamical Enzo code. We find that polarization signals can be as high as $\sim 1 \mu$K, a level that may be detectable by upcoming CMB experiments.

There has been considerable theoretical debate over whether photoionization and supernova feedback from the first Population III stars facilitate or suppress the formation of the next generation of stars. We present results from an Eulerian adaptive mesh refinement simulation demonstrating the formation of a primordial star within a region ionized by an earlier nearby star. Despite the higher temperatures of the ionized gas and its flow out of the dark matter potential wells, this second star formed within 23 million years of its neighbor's death. The enhanced electron fraction within the HII region catalyzes rapid molecular hydrogen formation that leads to faster cooling in the subsequent star forming halos than in the first halos. This "second generation" primordial protostar has a much lower accretion rate because, unlike the first protostar, it forms in a rotationally supported disk of approx. 10-100 solar masses. This is primarily due to the much higher angular momentum of the halo in which the second star forms. In contrast to previously published scenarios, such configurations may allow binaries or multiple systems of lower mass stars to form. These first high resolution calculations offer insight into the impact of feedback upon subsequent populations of stars and clearly demonstrate how primordial chemistry promotes the formation of subsequent generations of stars even in the presence of the entropy injected by the first stars into the IGM.

We compare two cosmological hydrodynamic simulation codes in the context of hierarchical galaxy formation: The SPH code GADGET, and the Eulerian AMR code ENZO. Both codes represent dark matter with the N-body method, but use different gravity solvers and fundamentally different approaches to hydrodynamics. We compare the GADGET `entropy conserving' SPH formulation with two ENZO methods: The piecewise parabolic method (PPM), and the artificial viscosity-based scheme used in the ZEUS code. In this paper we focus on a comparison of cosmological simulations that follow either only dark matter, or also adiabatic baryonic gas. The dark matter-only runs agree generally quite well, provided ENZO is run with a comparatively fine root grid and a low overdensity threshold for mesh refinement, otherwise the abundance of low-mass halos is suppressed. This is due to the hierarchical particle-mesh method used to compute gravitational forces in ENZO, which tends to deliver lower force resolution than the tree algorithm of GADGET. At comparable force resolution, we find that the latter offers substantially better performance and lower memory consumption than the present gravity solver in ENZO. In simulations that include adiabatic gas dynamics, we find general agreement in the distribution functions of temperature, entropy, and density for gas of moderate to high overdensity, as found inside dark matter halos. However, there are some significant differences at lower overdensities. We argue that these discrepancies are presumably owing to differences in the shock-capturing abilities of the different methods. In particular, ZEUS hydro leads to some unphysical heating at early times in preshock regions. Overall, the GADGET hydro results are bracketed by those for ENZO/ZEUS and ENZO/PPM. (abridged)

We present 40 fully hydrodynamical numerical simulations of the intergalactic gas that gives rise to the Ly-alpha forest. We make artificial spectra from each and measure three output parameters: the mean flux, a measure of the most common Ly-alpha line width (b-value), and the 1D power spectrum of the flux on scales from 0.01 - 0.1 s/km. We investigate how these output parameters are correlated with the input parameters that describe the simulations. These input parameters include the main cosmological parameters; and two astrophysical parameters gammah and gammahe. The parameter gammah controls the rate of ionization of H I, He I and He II while gammahe controls the rate of heating. We compare the output from the simulations to data from HIRES and UVES spectra and we give a new measurement of the flux power from these spectra at z=1.95. We find that simulations with a wide variety of sigma_8 values, from at least 0.8 - 1.1, can fit the small scale flux power and b-values when we adjust gammahe to compensate for the sigma_8 change. We can also use gammah to adjust the H I ionization rate to simultaneously match the mean flux. We derive scaling relations that give the output parameter values expected for a variety of input parameters with errors of 2% for optical depth and 1.4% for b-balues. We can break the degeneracy between gammahe and sigma_8 using large scale power or other data to fix sigma_8.

In this paper we introduce Enzo, a 3D MPI-parallel Eulerian block-structured adaptive mesh refinement cosmology code. Enzo is designed to simulate cosmological structure formation, but can also be used to simulate a wide range of astrophysical situations. Enzo solves dark matter N-body dynamics using the particle-mesh technique. The Poisson equation is solved using a combination of fast fourier transform (on a periodic root grid) and multigrid techniques (on non-periodic subgrids). Euler's equations of hydrodynamics are solved using a modified version of the piecewise parabolic method. Several additional physics packages are implemented in the code, including several varieties of radiative cooling, a metagalactic ultraviolet background, and prescriptions for star formation and feedback. We also show results illustrating properties of the adaptive mesh portion of the code. Enzo is publicly available and can be downloaded at http://cosmos.ucsd.edu/enzo/ .

We examine the dynamical evolution and statistical properties of the supernova ejecta of massive primordial stars in a cosmological framework to determine whether this first population of stars could have enriched the universe to the levels and dispersions seen by the most recent observations of the Lyman-Alpha forest. We evolve a lambda CDM model in a 1 Mpc^3 volume to a redshift of z = 15 and add ``bubbles'' of metal corresponding to the supernova ejecta of the first generation of massive stars in all dark matter halos with masses greater than 5 times 10^5 solar masses. These initial conditions are then evolved to z = 3 and the distribution and levels of metals are compared to observations. In the absence of further star formation the primordial metal is initially contained in halos and filaments. Photoevaporation of metal-enriched gas due to the metagalactic ultraviolet background radiation at the epoch of reionization (z ~ 6) causes a sharp increase of the metal volume filling factor. At z = 3, ~ 2.5% of the simulation volume (approx. 20% of the total gas mass) is filled with gas enriched above a metallicity of 10^-4 Z_solar, and less than 0.6% of the volume is enriched above a metallicity of 10^-3 Z_solar. This suggests that, even with the most optimistic prescription for placement of primordial supernova and the amount of metals produced by each supernova, this population of stars cannot entirely be responsible for the enrichment of the Lyman-$\alpha$ forest to the levels and dispersions seen by current observations unless we have severely underestimated the duration of the Pop III epoch. However, comparison to observations show that Pop III supernovae can be significant contributors to the very low overdensity Lyman-Alpha forest.

Approximately 30-40% of all baryons in the present day universe reside in a warm-hot intergalactic medium (WHIM), with temperatures between 10^5<T<10^7 K. This is a generic prediction from six hydrodynamic simulations of currently favored structure formation models having a wide variety of numerical methods, input physics, volumes, and spatial resolutions. Most of these warm-hot baryons reside in diffuse large-scale structures with a median overdensity around 10-30, not in virialized objects such as galaxy groups or galactic halos. The evolution of the WHIM is primarily driven by shock heating from gravitational perturbations breaking on mildly nonlinear, non-equilibrium structures such as filaments. Supernova feedback energy and radiative cooling play lesser roles in its evolution. WHIM gas is consistent with observations of the 0.25 keV X-ray background without being significantly heated by non-gravitational processes because the emitting gas is very diffuse. Our results confirm and extend previous work by Cen & Ostriker and Dave' et al.