We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.

Without an existing large scale coherent magnetic field in the early Universe, Population III (PopIII) stars would likely rotate at or near break-up speed. In this work, focusing on the accretion phase of PopIII stars, we investigate the possibility of generating a coherent magnetic field through large scale dynamo processes, as well as the corresponding field saturation level. Using results from hydrodynamic simulations, we demonstrate that primordial accretion disks are turbulent with a Shakura-Sunyaev disk parameter $\alpha_{ss} \gtrsim 10^{-3}$, and evidence helical turbulence with a dynamo number $\vert D_{\alpha \Omega} \vert \gg 10$. The presence of helical turbulence at these levels allows large scale dynamo modes to grow, and the saturation level is determined by the amount of net helicity remaining in the dynamo-active regions (a.k.a. the quenching problem). %We demonstrate that the magnetic field can reach approximate equipartition, with $B/B_{\rm eq} \sim 3$, indicating that the dynamo quenching problem could be alleviated through an accretion flow. We demonstrate that, if the accretion could successfully alleviate the quenching problem, the magnetic field can reach approximate equipartition with $B/B_{\rm eq} \sim 3$.

Under rapid cooling from molecular hydrogen, the accretion disks around Population III (PopIII) stars are believed to fragment, resulting in multiple accreting cores. In this paper, we build a theoretical framework for calculating the optical depth of H$_2$ ro-vibrational line cooling based on the vertical structure in these accretion disks. Applying this physically motivated prescription for the optical depth, we find that cooling in the inner disk with $r \lesssim 10 {\rm\ AU}$ is attenuated significantly as a result of high surface density; $PdV$ heating becomes more efficient than cooling, which prevents fragmentation in the inner disk. Despite this, cooling becomes dynamically important in the outer disk, favoring fragmentation. We argue that most of the resultant fragments are initially at the outer disk, and that any surviving fragment has to migrate slower than the disk-scale photo-evaporation process. Since type I migration is fast, migration slows down as a result of gap-opening in the disk structure. Two possible processes for gap-opening are studied: (1) through a massive, densely-cored ($\rho \gtrsim 10^{-8} {\rm\ g\ cm^{-3}}$) clump able to radiate away the excess gravitational potential energy, and (2) through a fast-growing central star, with $\dot{M} \gtrsim 2 \times 10^{-3} \, M_\odot {\rm\ yr^{-1}}$, whose gravity dominates the star-disk system and favors gap opening.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

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.

We present GAMER-2, a GPU-accelerated adaptive mesh refinement (AMR) code for astrophysics. It provides a rich set of features, including adaptive time-stepping, several hydrodynamic schemes, magnetohydrodynamics, self-gravity, particles, star formation, chemistry and radiative processes with GRACKLE, data analysis with yt, and memory pool for efficient object allocation. GAMER-2 is fully bitwise reproducible. For the performance optimization, it adopts hybrid OpenMP/MPI/GPU parallelization and utilizes overlapping CPU computation, GPU computation, and CPU-GPU communication. Load balancing is achieved using a Hilbert space-filling curve on a level-by-level basis without the need to duplicate the entire AMR hierarchy on each MPI process. To provide convincing demonstrations of the accuracy and performance of GAMER-2, we directly compare with Enzo on isolated disk galaxy simulations and with FLASH on galaxy cluster merger simulations. We show that the physical results obtained by different codes are in very good agreement, and GAMER-2 outperforms Enzo and FLASH by nearly one and two orders of magnitude, respectively, on the Blue Waters supercomputers using $1-256$ nodes. More importantly, GAMER-2 exhibits similar or even better parallel scalability compared to the other two codes. We also demonstrate good weak and strong scaling using up to 4096 GPUs and 65,536 CPU cores, and achieve a uniform resolution as high as $10{,}240^3$ cells. Furthermore, GAMER-2 can be adopted as an AMR+GPUs framework and has been extensively used for the wave dark matter ($\psi$DM) simulations. GAMER-2 is open source (available at https://github.com/gamer-project/gamer) and new contributions are welcome.

Software that processes real-world data or that models a physical system must have some way of managing units. While simple approaches like the understood convention that all data are in a unit system (such as the MKS SI unit system) do work in practice, they are fraught with possible sources of error both by developers and users of the software. In this paper we present unyt, a Python library based on NumPy and SymPy for handling data that has units. It is designed both to aid quick interactive calculations and to be tightly integrated into a larger Python application or library. We compare unyt with two other Python libraries for handling units, Pint and astropy.units, and find that unyt is faster, has higher test coverage, and has fewer lines of code.

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.

We describe a new method for simulating ionizing radiation and supernova feedback in the analogues of low-redshift galactic disks. In this method, which we call star-forming molecular cloud (SFMC) particles, we use a ray-tracing technique to solve the radiative transfer equation for ultraviolet photons emitted by thousands of distinct particles on the fly. Joined with high numerical resolution of 3.8 pc, the realistic description of stellar feedback helps to self-regulate star formation. This new feedback scheme also enables us to study the escape of ionizing photons from star-forming clumps and from a galaxy, and to examine the evolving environment of star-forming gas clumps. By simulating a galactic disk in a halo of 2.3e11 Msun, we find that the average escape fraction from all radiating sources on the spiral arms (excluding the central 2.5 kpc) fluctuates between 0.08% and 5.9% during a ~20 Myr period with a mean value of 1.1%. The flux of escaped photons from these sources is not strongly beamed, but manifests a large opening angle of more than 60 degree from the galactic pole. Further, we investigate the escape fraction per SFMC particle, f_esc(i), and how it evolves as the particle ages. We discover that the average escape fraction f_esc is dominated by a small number of SFMC particles with high f_esc(i). On average, the escape fraction from a SFMC particle rises from 0.27% at its birth to 2.1% at the end of a particle lifetime, 6 Myrs. This is because SFMC particles drift away from the dense gas clumps in which they were born, and because the gas around the star-forming clumps is dispersed by ionizing radiation and supernova feedback. The framework established in this study brings deeper insight into the physics of photon escape fraction from an individual star-forming clump, and from a galactic disk.

We investigate the spatially-resolved star formation relation using a galactic disk formed in a comprehensive high-resolution (3.8 pc) simulation. Our new implementation of stellar feedback includes ionizing radiation as well as supernova explosions, and we handle ionizing radiation by solving the radiative transfer equation rather than by a subgrid model. Photoheating by stellar radiation stabilizes gas against Jeans fragmentation, reducing the star formation rate. Because we have self-consistently calculated the location of ionized gas, we are able to make spatially-resolved mock observations of star formation tracers, such as H-alpha emission. We can also observe how stellar feedback manifests itself in the correlation between ionized and molecular gas. Applying our techniques to the disk in a galactic halo of 2.3e11 Msun, we find that the correlation between star formation rate density (estimated from mock H-alpha emission) and molecular hydrogen density shows large scatter, especially at high resolutions of <~ 75 pc that are comparable to the size of giant molecular clouds (GMCs). This is because an aperture of GMC size captures only particular stages of GMC evolution, and because H-alpha traces hot gas around star-forming regions and is displaced from the molecular hydrogen peaks themselves. By examining the evolving environment around star clusters, we speculate that the breakdown of the traditional star formation laws of the Kennicutt-Schmidt type at small scales is further aided by a combination of stars drifting from their birthplaces, and molecular clouds being dispersed via stellar feedback.

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 extend a machine learning (ML) framework presented previously to model galaxy formation and evolution in a hierarchical universe using N-body + hydrodynamical simulations. In this work, we show that ML is a promising technique to study galaxy formation in the backdrop of a hydrodynamical simulation. We use the Illustris Simulation to train and test various sophisticated machine learning algorithms. By using only essential dark matter halo physical properties and no merger history, our model predicts the gas mass, stellar mass, black hole mass, star formation rate, $g-r$ color, and stellar metallicity fairly robustly. Our results provide a unique and powerful phenomenological framework to explore the galaxy-halo connection that is built upon a solid hydrodynamical simulation. The promising reproduction of the listed galaxy properties demonstrably place ML as a promising and a significantly more computationally efficient tool to study small-scale structure formation. We find that ML mimics a full-blown hydrodynamical simulation surprisingly well in a computation time of mere minutes. The population of galaxies simulated by ML, while not numerically identical to Illustris, is statistically and physically robust and follows the same fundamental observational constraints. Machine learning offers an intriguing and promising technique to create quick mock galaxy catalogs in the future.

We present a new exploratory framework to model galaxy formation and evolution in a hierarchical universe by using machine learning (ML). Our motivations are two-fold: (1) presenting a new, promising technique to study galaxy formation, and (2) quantitatively analyzing the extent of the influence of dark matter halo properties on galaxies in the backdrop of semi-analytical models (SAMs). We use the influential Millennium Simulation and the corresponding Munich SAM to train and test various sophisticated machine learning algorithms (k-Nearest Neighbors, decision trees, random forests and extremely randomized trees). By using only essential dark matter halo physical properties for haloes of $M>10^{12} M_{\odot}$ and a partial merger tree, our model predicts the hot gas mass, cold gas mass, bulge mass, total stellar mass, black hole mass and cooling radius at z = 0 for each central galaxy in a dark matter halo for the Millennium run. Our results provide a unique and powerful phenomenological framework to explore the galaxy-halo connection that is built upon SAMs and demonstrably place ML as a promising and a computationally efficient tool to study small-scale structure formation.

Submillimetre-luminous galaxies at high-redshift are the most luminous, heavily star-forming galaxies in the Universe, and are characterised by prodigious emission in the far-infrared at 850 microns (S850 > 5 mJy). They reside in halos ~ 10^13Msun, have low gas fractions compared to main sequence disks at a comparable redshift, trace complex environments, and are not easily observable at optical wavelengths. Their physical origin remains unclear. Simulations have been able to form galaxies with the requisite luminosities, but have otherwise been unable to simultaneously match the stellar masses, star formation rates, gas fractions and environments. Here we report a cosmological hydrodynamic galaxy formation simulation that is able to form a submillimetre galaxy which simultaneously satisfies the broad range of observed physical constraints. We find that groups of galaxies residing in massive dark matter halos have rising star formation histories that peak at collective rates ~ 500-1000 Msun/yr at z=2-3, by which time the interstellar medium is sufficiently enriched with metals that the region may be observed as a submillimetre-selected system. The intense star formation rates are fueled in part by a reservoir gas supply enabled by stellar feedback at earlier times, not through major mergers. With a duty cycle of nearly a gigayear, our simulations show that the submillimetre-luminous phase of high-z galaxies is a drawn out one that is associated with significant mass buildup in early Universe proto-clusters, and that many submillimetre-luminous galaxies are actually composed of numerous unresolved components (for which there is some observational evidence).

The Dark Sky Simulations are an ongoing series of cosmological N-body simulations designed to provide a quantitative and accessible model of the evolution of the large-scale Universe. Such models are essential for many aspects of the study of dark matter and dark energy, since we lack a sufficiently accurate analytic model of non-linear gravitational clustering. In July 2014, we made available to the general community our early data release, consisting of over 55 Terabytes of simulation data products, including our largest simulation to date, which used $1.07 \times 10^{12}~(10240^3)$ particles in a volume $8h^{-1}\mathrm{Gpc}$ across. Our simulations were performed with 2HOT, a purely tree-based adaptive N-body method, running on 200,000 processors of the Titan supercomputer, with data analysis enabled by yt. We provide an overview of the derived halo catalogs, mass function, power spectra and light cone data. We show self-consistency in the mass function and mass power spectrum at the 1% level over a range of more than 1000 in particle mass. We also present a novel method to distribute and access very large datasets, based on an abstraction of the World Wide Web (WWW) as a file system, remote memory-mapped file access semantics, and a space-filling curve index. This method has been implemented for our data release, and provides a means to not only query stored results such as halo catalogs, but also to design and deploy new analysis techniques on large distributed datasets.

Starlight from galaxies plays a pivotal role throughout the process of cosmic reionisation. We present the statistics of dwarf galaxy properties at z > 7 in haloes with masses up to 10^9 solar masses, using a cosmological radiation hydrodynamics simulation that follows their buildup starting with their Population III progenitors. We find that metal-enriched star formation is not restricted to atomic cooling ($T_{\rm vir} \ge 10^4$ K) haloes, but can occur in haloes down to masses ~10^6 solar masses, especially in neutral regions. Even though these smallest galaxies only host up to 10^4 solar masses of stars, they provide nearly 30 per cent of the ionising photon budget. We find that the galaxy luminosity function flattens above M_UV ~ -12 with a number density that is unchanged at z < 10. The fraction of ionising radiation escaping into the intergalactic medium is inversely dependent on halo mass, decreasing from 50 to 5 per cent in the mass range $\log M/M_\odot = 7.0-8.5$. Using our galaxy statistics in a semi-analytic reionisation model, we find a Thomson scattering optical depth consistent with the latest Planck results, while still being consistent with the UV emissivity constraints provided by Ly$\alpha$ forest observations at z = 4-6.

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

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.

As scientists' needs for computational techniques and tools grow, they cease to be supportable by software developed in isolation. In many cases, these needs are being met by communities of practice, where software is developed by domain scientists to reach pragmatic goals and satisfy distinct and enumerable scientific goals. We present techniques that have been successful in growing and engaging communities of practice, specifically in the yt and Enzo communities.

Cosmological hydrodynamical simulations of galaxy evolution are increasingly able to produce realistic galaxies, but the largest hurdle remaining is in constructing subgrid models that accurately describe the behavior of stellar feedback. As an alternate way to test and calibrate such models, we propose to focus on the circumgalactic medium. To do so, we generate a suite of adaptive-mesh refinement (AMR) simulations for a Milky-Way-massed galaxy run to z=0, systematically varying the feedback implementation. We then post-process the simulation data to compute the absorbing column density for a wide range of common atomic absorbers throughout the galactic halo, including H I, Mg II, Si II, Si III, Si IV, C IV, N V, O VI, and O VII. The radial profiles of these atomic column densities are compared against several quasar absorption line studies, to determine if one feedback prescription is favored. We find that although our models match some of the observations (specifically those ions with lower ionization strengths), it is particularly difficult to match O VI observations. There is some indication that the models with increased feedback intensity are better matches. We demonstrate that sufficient metals exist in these halos to reproduce the observed column density distribution in principle, but the simulated circumgalactic medium lacks significant multiphase substructure and is generally too hot. Furthermore, we demonstrate the failings of inflow-only models (without energetic feedback) at populating the CGM with adequate metals to match observations even in the presence of multiphase structure. Additionally, we briefly investigate the evolution of the CGM from z=3 to present. Overall, we find that quasar absorption line observations of the gas around galaxies provide a new and important constraint on feedback models.

Massive stars provide feedback that shapes the interstellar medium of galaxies at all redshifts and their resulting stellar populations. Here we present three adaptive mesh refinement radiation hydrodynamics simulations that illustrate the impact of momentum transfer from ionising radiation to the absorbing gas on star formation in high-redshift dwarf galaxies. Momentum transfer is calculated by solving the radiative transfer equation with a ray tracing algorithm that is adaptive in spatial and angular coordinates. We find that momentum input partially affects star formation by increasing the turbulent support to a three-dimensional rms velocity equal to the circular velocity of early haloes. Compared to a calculation that neglects radiation pressure, the star formation rate is decreased by a factor of five to 1.8 x 10^{-2} Msun/yr in a dwarf galaxy with a dark matter and stellar mass of 2.0 x 10^8 and 4.5 x 10^5 solar masses, respectively, when radiation pressure is included. Its mean metallicity of 10^{-2.1} Z_sun is consistent with the observed dwarf galaxy luminosity-metallicity relation. However, what one may naively expect from the calculation without radiation pressure, the central region of the galaxy overcools and produces a compact, metal-rich stellar population with an average metallicity of 0.3 Z_sun, indicative of an incorrect physical recipe. In addition to photo-heating in HII regions, radiation pressure further drives dense gas from star forming regions, so supernovae feedback occurs in a warmer and more diffuse medium, launching metal-rich outflows. Capturing this aspect and a temporal separation between the start of radiative and supernova feedback are numerically important in the modeling of galaxies to avoid the "overcooling problem". We estimate that dust in early low-mass galaxies is unlikely to aid in momentum transfer from radiation to the gas.

The usage of the high-level scripting language Python has enabled new mechanisms for data interrogation, discovery and visualization of scientific data. We present yt, an open source, community-developed astrophysical analysis and visualization toolkit for data generated by high-performance computing (HPC) simulations of astrophysical phenomena. Through a separation of responsibilities in the underlying Python code, yt allows data generated by incompatible, and sometimes even directly competing, astrophysical simulation platforms to be analyzed in a consistent manner, focusing on physically relevant quantities rather than quantities native to astrophysical simulation codes. We present on its mechanisms for data access, capabilities for MPI-parallel analysis, and its implementation as an in situ analysis and visualization tool.

We study the buildup of magnetic fields during the formation of Population III star-forming regions, by conducting cosmological simulations from realistic initial conditions and varying the Jeans resolution. To investigate this in detail, we start simulations from identical initial conditions, mandating 16, 32 and 64 zones per Jeans length, and studied the variation in their magnetic field amplification. We find that, while compression results in some amplification, turbulent velocity fluctuations driven by the collapse can further amplify an initially weak seed field via dynamo action, provided there is sufficient numerical resolution to capture vortical motions (we find this requirement to be 64 zones per Jeans length, slightly larger than, but consistent with previous work run with more idealized collapse scenarios). We explore saturation of amplification of the magnetic field, which could potentially become dynamically important in subsequent, fully-resolved calculations. We have also identified a relatively surprising phenomena that is purely hydrodynamic: the higher-resolved simulations possess substantially different characteristics, including higher infall-velocity, increased temperatures inside 1000 AU, and decreased molecular hydrogen content in the innermost region. Furthermore, we find that disk formation is suppressed in higher-resolution calculations, at least at the times that we can follow the calculation. We discuss the effect this may have on the buildup of disks over the accretion history of the first clump to form as well as the potential for gravitational instabilities to develop and induce fragmentation.

By definition, Population III stars are metal-free, and their protostellar collapse is driven by molecular hydrogen cooling in the gas-phase, leading to large characteristic masses. Population II stars with lower characteristic masses form when the star-forming gas reaches a critical metallicity of 10^{-6} - 10^{-3.5} Z_\odot. We present an adaptive mesh refinement radiation hydrodynamics simulation that follows the transition from Population III to II star formation. The maximum spatial resolution of 1 comoving parsec allows for individual molecular clouds to be well-resolved and their stellar associations to be studied in detail. We model stellar radiative feedback with adaptive ray tracing. A top-heavy initial mass function for the Population III stars is considered, resulting in a plausible distribution of pair-instability supernovae and associated metal enrichment. We find that the gas fraction recovers from 5 percent to nearly the cosmic fraction in halos with merger histories rich in halos above 10^7 solar masses. A single pair-instability supernova is sufficient to enrich the host halo to a metallicity floor of 10^{-3} Z_\odot and to transition to Population II star formation. This provides a natural explanation for the observed floor on damped Lyman alpha (DLA) systems metallicities reported in the literature, which is of this order. We find that stellar metallicities do not necessarily trace stellar ages, as mergers of halos with established stellar populations can create superpositions of t-Z evolutionary tracks. A bimodal metallicity distribution is created after a starburst occurs when the halo can cool efficiently through atomic line cooling.

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.

The analysis of complex multiphysics astrophysical simulations presents a unique and rapidly growing set of challenges: reproducibility, parallelization, and vast increases in data size and complexity chief among them. In order to meet these challenges, and in order to open up new avenues for collaboration between users of multiple simulation platforms, we present yt (available at http://yt.enzotools.org/), an open source, community-developed astrophysical analysis and visualization toolkit. Analysis and visualization with yt are oriented around physically relevant quantities rather than quantities native to astrophysical simulation codes. While originally designed for handling Enzo's structure adaptive mesh refinement (AMR) data, yt has been extended to work with several different simulation methods and simulation codes including Orion, RAMSES, and FLASH. We report on its methods for reading, handling, and visualizing data, including projections, multivariate volume rendering, multi-dimensional histograms, halo finding, light cone generation and topologically-connected isocontour identification. Furthermore, we discuss the underlying algorithms yt uses for processing and visualizing data, and its mechanisms for parallelization of analysis tasks.

The transformation of atomic hydrogen to molecular hydrogen through three-body reactions is a crucial stage in the collapse of primordial, metal-free halos, where the first generation of stars (Population III stars) in the Universe are formed. However, in the published literature, the rate coefficient for this reaction is uncertain by nearly an order of magnitude. We report on the results of both adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of the collapse of metal-free halos as a function of the value of this rate coefficient. For each simulation method, we have simulated a single halo three times, using three different values of the rate coefficient. We find that while variation between halo realizations may be greater than that caused by the three-body rate coefficient being used, both the accretion physics onto Population III protostars as well as the long-term stability of the disk and any potential fragmentation may depend strongly on this rate coefficient.

We report on simulations of the formation of the first stars in the Universe, where we identify regions of hot atomic gas (fH2 < 10-6) at densities above 10-14 g/cc, heated to temperatures ranging between 3000 and 8000 K. Within this temperature range atomic hydrogen is unable to cool effectively. We describe the kinetic and thermal characteristics of these regions and investigate their origin. We find that these regions, while small in total mass fraction of the cloud, may be dynamically important over the accretion timescale for the central clump in the cloud, particularly as a chemical, rather than radiative, mechanism for clearing the polar regions of the accretion disk of material and terminating accretion along these directions. These inherently three-dimensional effects stress the need for multi-dimensional calculations of protostellar accretion for reliable predictions of the masses of the very first stars.

Modern N-body cosmological simulations contain billions ($10^9$) of dark matter particles. These simulations require hundreds to thousands of gigabytes of memory, and employ hundreds to tens of thousands of processing cores on many compute nodes. In order to study the distribution of dark matter in a cosmological simulation, the dark matter halos must be identified using a halo finder, which establishes the halo membership of every particle in the simulation. The resources required for halo finding are similar to the requirements for the simulation itself. In particular, simulations have become too extensive to use commonly-employed halo finders, such that the computational requirements to identify halos must now be spread across multiple nodes and cores. Here we present a scalable-parallel halo finding method called Parallel HOP for large-scale cosmological simulation data. Based on the halo finder HOP, it utilizes MPI and domain decomposition to distribute the halo finding workload across multiple compute nodes, enabling analysis of much larger datasets than is possible with the strictly serial or previous parallel implementations of HOP. We provide a reference implementation of this method as a part of the toolkit yt, an analysis toolkit for Adaptive Mesh Refinement (AMR) data that includes complementary analysis modules. Additionally, we discuss a suite of benchmarks that demonstrate that this method scales well up to several hundred tasks and datasets in excess of $2000^3$ particles. The Parallel HOP method and our implementation can be readily applied to any kind of N-body simulation data and is therefore widely applicable.

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.

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)

Numerous cosmological hydrodynamic studies have addressed the formation of galaxies. Here we choose to study the first stages of galaxy formation, including non-equilibrium atomic primordial gas cooling, gravity and hydrodynamics. Using initial conditions appropriate for the concordance cosmological model of structure formation, we perform two adaptive mesh refinement simulations of ~10^8 M_sun galaxies at high redshift. The calculations resolve the Jeans length at all times with more than 16 cells and capture over 14 orders of magnitude in length scales. In both cases, the dense, 10^5 solar mass, one parsec central regions are found to contract rapidly and have turbulent Mach numbers up to 4. Despite the ever decreasing Jeans length of the isothermal gas, we only find one site of fragmentation during the collapse. However, rotational secular bar instabilities transport angular momentum outwards in the central parsec as the gas continues to collapse and lead to multiple nested unstable fragments with decreasing masses down to sub-Jupiter mass scales. Although these numerical experiments neglect star formation and feedback, they clearly highlight the physics of turbulence in gravitationally collapsing gas. The angular momentum segregation seen in our calculations plays an important role in theories that form supermassive black holes from gaseous collapse.