-
Why does the Milky Way have a metallicity floor?
Authors:
Britton D. Smith,
Brian W. O'Shea,
Sadegh Khochfar,
Matthew J. Turk,
John H. Wise,
Michael L. Norman
Abstract:
The prevalence of light element enhancement in the most metal-poor stars is potentially an indication that the Milky Way has a metallicity floor for star formation around $\sim$10$^{-3.5}$ Z$_{\odot}$. We propose that this metallicity floor has its origins in metal-enriched star formation in the minihalos present during the Galaxy's initial formation. To arrive at this conclusion, we analyze a cos…
▽ More
The prevalence of light element enhancement in the most metal-poor stars is potentially an indication that the Milky Way has a metallicity floor for star formation around $\sim$10$^{-3.5}$ Z$_{\odot}$. We propose that this metallicity floor has its origins in metal-enriched star formation in the minihalos present during the Galaxy's initial formation. To arrive at this conclusion, we analyze a cosmological radiation hydrodynamics simulation that follows the concurrent evolution of multiple Population III star-forming minihalos. The main driver for the central gas within minihalos is the steady increase in hydrostatic pressure as the halos grow. We incorporate this insight into a hybrid one-zone model that switches between pressure-confined and modified free-fall modes to evolve the gas density with time according to the ratio of the free-fall and sound-crossing timescales. This model is able to accurately reproduce the density and chemo-thermal evolution of the gas in each of the simulated minihalos up to the point of runaway collapse. We then use this model to investigate how the gas responds to the absence of H$_{2}$. Without metals, the central gas becomes increasingly stable against collapse as it grows to the atomic cooling limit. When metals are present in the halo at a level of $\sim$10$^{-3.7}$ Z$_{\odot}$, however, the gas is able to achieve gravitational instability while still in the minihalo regime. Thus, we conclude that the Galaxy's metallicity floor is set by the balance within minihalos of gas-phase metal cooling and the radiation background associated with its early formation environment.
△ Less
Submitted 10 July, 2024; v1 submitted 12 June, 2024;
originally announced June 2024.
-
The Role of Radiation and Halo Mergers in Pop III Star Formation
Authors:
Lilia Correa Magnus,
Britton D. Smith,
Sadegh Khochfar,
Brian W. O'Shea,
John H. Wise,
Michael L. Norman,
Matthew J. Turk
Abstract:
We present a study of the co-evolution of a population of primordial star-forming minihalos at Cosmic Dawn. In this study, we highlight the influence of individual Population III stars on the ability of nearby minihalos to form sufficient molecular hydrogen to undergo star formation. In the absence of radiation, we find the minimum halo mass required to bring about collapse to be ~10^5 Msun, this…
▽ More
We present a study of the co-evolution of a population of primordial star-forming minihalos at Cosmic Dawn. In this study, we highlight the influence of individual Population III stars on the ability of nearby minihalos to form sufficient molecular hydrogen to undergo star formation. In the absence of radiation, we find the minimum halo mass required to bring about collapse to be ~10^5 Msun, this increases to ~10^6 Msun after two stars have formed. We find an inverse relationship between halo mass and the time required for it to recover its molecular gas after being disrupted by radiation from a nearby star. We also take advantage of the extremely high resolution to investigate the effects of major and minor mergers on the gas content of star-forming minihalos. Contrary to previous claims of fallback of supernova ejecta, we find minihalos evacuated after hosting Pop III stars primarily recover gas through mergers with undisturbed halos. We identify an intriguing type of major merger between recently evacuated halos and gas-rich ones, finding that these 'mixed' mergers accelerate star formation instead of suppressing it like their low redshift counterparts. We attribute this to the gas-poor nature of one of the merging halos resulting in no significant rise in temperature or turbulence and instead inducing a rapid increase in central density and hydrostatic pressure. This constitutes a novel formation pathway for Pop III stars and establishes major mergers as potentially the primary source of gas, thus redefining the role of major mergers at this epoch.
△ Less
Submitted 16 October, 2023; v1 submitted 7 July, 2023;
originally announced July 2023.
-
Powderday: Dust Radiative Transfer for Galaxy Simulations
Authors:
Desika Narayanan,
Matthew J. Turk,
Thomas Robitaille,
Ashley J. Kelly,
B. Connor McClellan,
Ray S. Sharma,
Prerak Garg,
Matthew Abruzzo,
Ena Choi,
Charlie Conroy,
Benjamin D. Johnson,
Benjamin Kimock,
Qi Li,
Christopher C. Lovell,
Sidney Lower,
George C. Privon,
Jonathan Roberts,
Snigdaa Sethuram,
Gregory F. Snyder,
Robert Thompson,
John H. Wise
Abstract:
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 fle…
▽ More
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.
△ Less
Submitted 18 June, 2020;
originally announced June 2020.
-
Cyberinfrastructure Requirements to Enhance Multi-messenger Astrophysics
Authors:
Philip Chang,
Gabrielle Allen,
Warren Anderson,
Federica B. Bianco,
Joshua S. Bloom,
Patrick R. Brady,
Adam Brazier,
S. Bradley Cenko,
Sean M. Couch,
Tyce DeYoung,
Ewa Deelman,
Zachariah B Etienne,
Ryan J. Foley,
Derek B Fox,
V. Zach Golkhou,
Darren R Grant,
Chad Hanna,
Kelly Holley-Bockelmann,
D. Andrew Howell,
E. A. Huerta,
Margaret W. G. Johnson,
Mario Juric,
David L. Kaplan,
Daniel S. Katz,
Azadeh Keivani
, et al. (17 additional authors not shown)
Abstract:
The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up,…
▽ More
The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
unyt: Handle, manipulate, and convert data with units in Python
Authors:
Nathan J. Goldbaum,
John A. ZuHone,
Matthew J. Turk,
Kacper Kowalik,
Anna L. Rosen
Abstract:
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 bas…
▽ More
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.
△ Less
Submitted 1 October, 2018; v1 submitted 6 June, 2018;
originally announced June 2018.
-
GAMER-2: a GPU-accelerated adaptive mesh refinement code -- accuracy, performance, and scalability
Authors:
Hsi-Yu Schive,
John A. ZuHone,
Nathan J. Goldbaum,
Matthew J. Turk,
Massimo Gaspari,
Chin-Yu Cheng
Abstract:
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…
▽ More
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 ($ψ$DM) simulations. GAMER-2 is open source (available at https://github.com/gamer-project/gamer) and new contributions are welcome.
△ Less
Submitted 4 October, 2018; v1 submitted 19 December, 2017;
originally announced December 2017.
-
Grackle: a Chemistry and Cooling Library for Astrophysics
Authors:
Britton D. Smith,
Greg L. Bryan,
Simon C. O. Glover,
Nathan J. Goldbaum,
Matthew J. Turk,
John Regan,
John H. Wise,
Hsi-Yu Schive,
Tom Abel,
Andrew Emerick,
Brian W. O'Shea,
Peter Anninos,
Cameron B. Hummels,
Sadegh Khochfar
Abstract:
We present the Grackle chemistry and cooling library for astrophysical simulations and models. Grackle provides a treatment of non-equilibrium primordial chemistry and cooling for H, D, and He species, including H2 formation on dust grains; tabulated primordial and metal cooling; multiple UV background models; and support for radiation transfer and arbitrary heat sources. The library has an easily…
▽ More
We present the Grackle chemistry and cooling library for astrophysical simulations and models. Grackle provides a treatment of non-equilibrium primordial chemistry and cooling for H, D, and He species, including H2 formation on dust grains; tabulated primordial and metal cooling; multiple UV background models; and support for radiation transfer and arbitrary heat sources. The library has an easily implementable interface for simulation codes written in C, C++, and Fortran as well as a Python interface with added convenience functions for semi-analytical models. As an open-source project, Grackle provides a community resource for accessing and disseminating astrochemical data and numerical methods. We present the full details of the core functionality, the simulation and Python interfaces, testing infrastructure, performance, and range of applicability. Grackle is a fully open-source project and new contributions are welcome.
△ Less
Submitted 14 December, 2016; v1 submitted 29 October, 2016;
originally announced October 2016.
-
The AGORA High-Resolution Galaxy Simulations Comparison Project. II: Isolated Disk Test
Authors:
Ji-hoon Kim,
Oscar Agertz,
Romain Teyssier,
Michael J. Butler,
Daniel Ceverino,
Jun-Hwan Choi,
Robert Feldmann,
Ben W. Keller,
Alessandro Lupi,
Thomas Quinn,
Yves Revaz,
Spencer Wallace,
Nickolay Y. Gnedin,
Samuel N. Leitner,
Sijing Shen,
Britton D. Smith,
Robert Thompson,
Matthew J. Turk,
Tom Abel,
Kenza S. Arraki,
Samantha M. Benincasa,
Sukanya Chakrabarti,
Colin DeGraf,
Avishai Dekel,
Nathan J. Goldbaum
, et al. (18 additional authors not shown)
Abstract:
Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the st…
▽ More
Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.
△ Less
Submitted 8 March, 2018; v1 submitted 10 October, 2016;
originally announced October 2016.
-
Machine Learning and Cosmological Simulations II: Hydrodynamical Simulations
Authors:
Harshil M. Kamdar,
Matthew J. Turk,
Robert J. Brunner
Abstract:
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 algorit…
▽ More
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.
△ Less
Submitted 12 January, 2016; v1 submitted 26 October, 2015;
originally announced October 2015.
-
Machine Learning and Cosmological Simulations I: Semi-Analytical Models
Authors:
Harshil M. Kamdar,
Matthew J. Turk,
Robert J. Brunner
Abstract:
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 th…
▽ More
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.
△ Less
Submitted 21 October, 2015;
originally announced October 2015.
-
Dark Sky Simulations: Early Data Release
Authors:
Samuel W. Skillman,
Michael S. Warren,
Matthew J. Turk,
Risa H. Wechsler,
Daniel E. Holz,
P. M. Sutter
Abstract:
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…
▽ More
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.
△ Less
Submitted 9 July, 2014;
originally announced July 2014.
-
The Birth of a Galaxy - III. Propelling reionisation with the faintest galaxies
Authors:
John H. Wise,
Vasiliy G. Demchenko,
Martin T. Halicek,
Michael L. Norman,
Matthew J. Turk,
Tom Abel,
Britton D. Smith
Abstract:
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 atom…
▽ More
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$α$ forest observations at z = 4-6.
△ Less
Submitted 19 May, 2014; v1 submitted 24 March, 2014;
originally announced March 2014.
-
The AGORA High-Resolution Galaxy Simulations Comparison Project
Authors:
Ji-hoon Kim,
Tom Abel,
Oscar Agertz,
Greg L. Bryan,
Daniel Ceverino,
Charlotte Christensen,
Charlie Conroy,
Avishai Dekel,
Nickolay Y. Gnedin,
Nathan J. Goldbaum,
Javiera Guedes,
Oliver Hahn,
Alexander Hobbs,
Philip F. Hopkins,
Cameron B. Hummels,
Francesca Iannuzzi,
Dusan Keres,
Anatoly Klypin,
Andrey V. Kravtsov,
Mark R. Krumholz,
Michael Kuhlen,
Samuel N. Leitner,
Piero Madau,
Lucio Mayer,
Christopher E. Moody
, et al. (21 additional authors not shown)
Abstract:
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…
▽ More
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.
△ Less
Submitted 24 December, 2013; v1 submitted 12 August, 2013;
originally announced August 2013.
-
Enzo: An Adaptive Mesh Refinement Code for Astrophysics
Authors:
The Enzo Collaboration,
Greg L. Bryan,
Michael L. Norman,
Brian W. O'Shea,
Tom Abel,
John H. Wise,
Matthew J. Turk,
Daniel R. Reynolds,
David C. Collins,
Peng Wang,
Samuel W. Skillman,
Britton Smith,
Robert P. Harkness,
James Bordner,
Ji-hoon Kim,
Michael Kuhlen,
Hao Xu,
Nathan Goldbaum,
Cameron Hummels,
Alexei G. Kritsuk,
Elizabeth Tasker,
Stephen Skory,
Christine M. Simpson,
Oliver Hahn,
Jeffrey S. Oishi
, et al. (4 additional authors not shown)
Abstract:
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-gra…
▽ More
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.
△ Less
Submitted 8 July, 2013;
originally announced July 2013.
-
Population III Star Formation In Large Cosmological Simulations I. Halo Temporal and Physical Environment
Authors:
Brian D. Crosby,
Brian W. O'Shea,
Britton D. Smith,
Matthew J. Turk,
Oliver Hahn
Abstract:
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 reflec…
▽ More
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.
△ Less
Submitted 19 June, 2013;
originally announced June 2013.
-
How to Scale a Code in the Human Dimension
Authors:
Matthew J. Turk
Abstract:
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…
▽ More
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.
△ Less
Submitted 29 January, 2013;
originally announced January 2013.
-
Dwarf Galaxies with Ionizing Radiation Feedback. II: Spatially-resolved Star Formation Relation
Authors:
Ji-hoon Kim,
Mark R. Krumholz,
John H. Wise,
Matthew J. Turk,
Nathan J. Goldbaum,
Tom Abel
Abstract:
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 stabili…
▽ More
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.
△ Less
Submitted 8 March, 2018; v1 submitted 25 October, 2012;
originally announced October 2012.
-
Dwarf Galaxies with Ionizing Radiation Feedback. I: Escape of Ionizing Photons
Authors:
Ji-hoon Kim,
Mark R. Krumholz,
John H. Wise,
Matthew J. Turk,
Nathan J. Goldbaum,
Tom Abel
Abstract:
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…
▽ More
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.
△ Less
Submitted 8 March, 2018; v1 submitted 11 October, 2012;
originally announced October 2012.
-
The Birth of a Galaxy. II. The Role of Radiation Pressure
Authors:
John H. Wise,
Tom Abel,
Matthew J. Turk,
Michael L. Norman,
Britton D. Smith
Abstract:
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 s…
▽ More
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.
△ Less
Submitted 26 July, 2012; v1 submitted 5 June, 2012;
originally announced June 2012.
-
High-Performance Astrophysical Simulations and Analysis with Python
Authors:
Matthew J. Turk,
Britton D. Smith
Abstract:
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 unde…
▽ More
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.
△ Less
Submitted 19 December, 2011;
originally announced December 2011.
-
Magnetic Fields in Population III Star Formation
Authors:
Matthew J. Turk,
Jeffrey S. Oishi,
Tom Abel,
Greg Bryan
Abstract:
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 amplifi…
▽ More
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.
△ Less
Submitted 19 December, 2011;
originally announced December 2011.
-
A Multi-Code Analysis Toolkit for Astrophysical Simulation Data
Authors:
Matthew J. Turk,
Britton D. Smith,
Jeffrey S. Oishi,
Stephen Skory,
Samuel W. Skillman,
Tom Abel,
Michael L. Norman
Abstract:
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.enzoto…
▽ More
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.
△ Less
Submitted 15 November, 2010;
originally announced November 2010.
-
The Birth of a Galaxy: Primordial Metal Enrichment and Stellar Populations
Authors:
John H. Wise,
Matthew J. Turk,
Michael L. Norman,
Tom Abel
Abstract:
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 simulati…
▽ More
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.
△ Less
Submitted 30 November, 2011; v1 submitted 11 November, 2010;
originally announced November 2010.
-
Effects of Varying the Three-Body Molecular Hydrogen Formation Rate in Primordial Star Formation
Authors:
Matthew J. Turk,
Paul Clark,
S. C. O. Glover,
T. H. Greif,
Tom Abel,
Ralf Klessen,
Volker Bromm
Abstract:
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…
▽ More
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.
△ Less
Submitted 8 November, 2010;
originally announced November 2010.
-
High-Entropy Polar Regions Around the First Protostars
Authors:
Matthew J. Turk,
Michael L. Norman,
Tom Abel
Abstract:
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…
▽ More
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.
△ Less
Submitted 28 October, 2010;
originally announced October 2010.
-
Galaxy Cluster Radio Relics in Adaptive Mesh Refinement Cosmological Simulations: Relic Properties and Scaling Relationships
Authors:
Samuel W. Skillman,
Eric J. Hallman,
Brian W. O'Shea,
Jack O. Burns,
Britton D. Smith,
Matthew J. Turk
Abstract:
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 vici…
▽ More
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.
△ Less
Submitted 9 May, 2011; v1 submitted 17 June, 2010;
originally announced June 2010.
-
Parallel HOP: A Scalable Halo Finder for Massive Cosmological Data Sets
Authors:
Stephen Skory,
Matthew J. Turk,
Michael L. Norman,
Alison L. Coil
Abstract:
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 estab…
▽ More
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.
△ Less
Submitted 17 June, 2010; v1 submitted 19 January, 2010;
originally announced January 2010.
-
The Formation of Population III Binaries from Cosmological Initial Conditions
Authors:
Matthew J. Turk,
Tom Abel,
Brian W. O'Shea
Abstract:
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 as…
▽ More
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.
△ Less
Submitted 16 July, 2009;
originally announced July 2009.
-
Three Modes of Metal-Enriched Star Formation in the Early Universe
Authors:
Britton D. Smith,
Matthew J. Turk,
Steinn Sigurdsson,
Brian W. O'Shea,
Michael L. Norman
Abstract:
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…
▽ More
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)
△ Less
Submitted 28 September, 2008; v1 submitted 10 June, 2008;
originally announced June 2008.
-
Resolving the Formation of Protogalaxies. II. Central Gravitational Collapse
Authors:
John H. Wise,
Matthew J. Turk,
Tom Abel
Abstract:
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…
▽ More
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.
△ Less
Submitted 26 March, 2008; v1 submitted 8 October, 2007;
originally announced October 2007.