Showing 1–50 of 52 results for author: Smith, B D

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.

Comments: 11 pages, 5 figures, accepted for publication in MNRAS

Figuring Out Gas & Galaxies In Enzo (FOGGIE) VIII: Complex and Stochastic Metallicity Gradients at z > 2

Authors: Ayan Acharyya, Molly S. Peeples, Jason Tumlinson, Brian W. O Shea, Cassandra Lochhaas, Anna C. Wright, Raymond C. Simons, Ramona Augustin, Britton D. Smith, Eugene Hyeonmin Lee

Abstract: Gas-phase metallicity gradients are a crucial element in understanding the chemical evolution of galaxies. We use the FOGGIE simulations to study the metallicity gradients ($\nabla Z$) of six Milky Way-like galaxies throughout their evolution. FOGGIE galaxies generally exhibit steep negative gradients for most of their history, with only a few short-lived instances reaching positive slopes that ap… ▽ More Gas-phase metallicity gradients are a crucial element in understanding the chemical evolution of galaxies. We use the FOGGIE simulations to study the metallicity gradients ($\nabla Z$) of six Milky Way-like galaxies throughout their evolution. FOGGIE galaxies generally exhibit steep negative gradients for most of their history, with only a few short-lived instances reaching positive slopes that appear to arise mainly from interactions with other galaxies. FOGGIE concurs with other simulation results but disagrees with the robust observational finding that flat and positive gradients are common at $z>1$. By tracking the metallicity gradient at a rapid cadence of simulation outputs ($\sim 5$--10 Myr), we find that theoretical gradients are highly stochastic: the FOGGIE galaxies spend $\sim 30-50$\% of their time far away from a smoothed trajectory inferred from analytic models or other, less high-cadence simulations. This rapid variation makes instantaneous gradients from observations more difficult to interpret in terms of physical processes. Because of these geometric and stochastic complications, we explore non-parametric methods of quantifying the evolving metallicity distribution at $z > 1$. We investigate how efficiently non-parametric measures of the 2-D metallicity distribution respond to metal production and mixing. Our results suggest that new methods of quantifying and interpreting gas-phase metallicity will be needed to relate trends in upcoming high-$z$ {\it JWST} observations with the underlying physics of gas accretion, expulsion, and recycling in early galaxies. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ

Hungry or Not: How Stellar-Mass Black Holes Grow (or Don't) in Dark Matter Mini-Haloes at High-Resolution

Authors: Simone T Gordon, Britton D Smith, Sadegh Khochfar, John Anthony Regan

Abstract: We compare the performance of the popular Bondi-Hoyle-Lyttleton (BHL) accretion scheme with a simple mass-flux scheme applied to stellar-mass black holes (BHs) across six levels of increasing spatial resolution. Simulating the formation of black holes within cosmological mini-haloes at $z \sim 20$, we investigate scenarios both with and without supernova events, which result in BHs of initial mass… ▽ More We compare the performance of the popular Bondi-Hoyle-Lyttleton (BHL) accretion scheme with a simple mass-flux scheme applied to stellar-mass black holes (BHs) across six levels of increasing spatial resolution. Simulating the formation of black holes within cosmological mini-haloes at $z \sim 20$, we investigate scenarios both with and without supernova events, which result in BHs of initial mass $10.8 \, \text{M}_\odot$ and $270 \, \text{M}_\odot$ respectively. Our explicit focus on the stellar-mass range pushes the maximum resolution down to sub-$10^{-3} \, \text{pc}$ regimes, where more complicated gas dynamics are resolved. We observe efficient growth and rotationally supported, $\sim$$10^{-1} \, \text{pc}$-scale discs around all $270 \, \text{M}_\odot$ BHs independent of resolution and accretion scheme, though clumps, bars, and spiral arm structures impact stability at high resolution. We analyse the effect of these instabilities on the accretion cycle. In contrast, all bar one of the $10.8 \, \text{M}_\odot$ BHs fail to attract a disc and experience modest growth, even when characteristic scales of accretion and dynamical friction are reasonably resolved. While the two accretion schemes somewhat converge in mass growth for the $270 \, \text{M}_\odot$ case over $1 \, \text{Myr}$, the greater degree of gas fragmentation induces more randomness in the evolution of the $10.8 \, \text{M}_\odot$ BHs. We conclude that early universe black holes of $M_{\text{BH}} \sim 10^1 \, \text{M}_\odot$ struggle to grow even in gas-rich environments without feedback in comparison to seeds of $M_{\text{BH}} \sim 10^2 \, \text{M}_\odot$, and the latter exhibit convergent growth histories across accretion schemes below a spatial resolution of $1 \times 10^{-3} \, \text{pc}$. △ Less

Submitted 18 March, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: 24 pages, 13 figures. Accepted by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 529, Issue 1, March 2024, Pages 604-627

Figuring Out Gas & Galaxies In Enzo (FOGGIE) VII: The (Dis)Assembly of Stellar Halos

Authors: Anna C. Wright, Jason Tumlinson, Molly S. Peeples, Brian W. O'Shea, Cassandra Lochhaas, Lauren Corlies, Britton D. Smith, Nguyen Binh, Ramona Augustin, Raymond C. Simons

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

Submitted 12 June, 2024; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: 33 pages, 15 figures, accepted to ApJ

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.

Comments: 15 pages, 12 figures, accepted for publication in MNRAS, revised with moderate changes

Developing Digital Twins for Earth Systems: Purpose, Requisites, and Benefits

Authors: Yuhan Rao, Rob Redmon, Kirstine Dale, Sue E. Haupt, Aaron Hopkinson, Ann Bostrom, Sid Boukabara, Thomas Geenen, David M. Hall, Benjamin D. Smith, Dev Niyogi, V. Ramaswamy, Eric A. Kihn

Abstract: The accelerated change in our planet due to human activities has led to grand societal challenges including health crises, intensified extreme weather events, food security, environmental injustice, etc. Digital twin systems combined with emerging technologies such as artificial intelligence and edge computing provide opportunities to support planning and decision-making to address these challenge… ▽ More The accelerated change in our planet due to human activities has led to grand societal challenges including health crises, intensified extreme weather events, food security, environmental injustice, etc. Digital twin systems combined with emerging technologies such as artificial intelligence and edge computing provide opportunities to support planning and decision-making to address these challenges. Digital twins for Earth systems (DT4ESs) are defined as the digital representation of the complex integrated Earth system including both natural processes and human activities. They have the potential to enable a diverse range of users to explore what-if scenarios across spatial and temporal scales to improve our understanding, prediction, mitigation, and adaptation to grand societal challenges. The 4th NOAA AI Workshop convened around 100 members who are developing or interested in participating in the development of DT4ES to discuss a shared community vision and path forward on fostering a future ecosystem of interoperable DT4ES. This paper summarizes the workshop discussions around DT4ES. We first defined the foundational features of a viable digital twins for Earth system that can be used to guide the development of various use cases of DT4ES. Finally, we made practical recommendations for the community on different aspects of collaboration in order to enable a future ecosystem of interoperable DT4ES, including equity-centered use case development, community-driven investigation of interoperability for DT4ES, trust-oriented co-development, and developing a community of practice. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: This whitepaper is an outcome of the 4th NOAA AI Workshop

Microwave-to-optical conversion in a room-temperature $^{87}$Rb vapor with frequency-division multiplexing control

Authors: Benjamin D. Smith, Bahar Babaei, Andal Narayanan, Lindsay J. LeBlanc

Abstract: Coherent microwave-to-optical conversion is crucial for transferring quantum information generated in the microwave domain to optical frequencies, where propagation losses can be minimised. Among the various physical platforms that have realized coherent microwave-to-optical transduction, those that use atoms as transducers have shown rapid progress in recent years. In this paper we report an expe… ▽ More Coherent microwave-to-optical conversion is crucial for transferring quantum information generated in the microwave domain to optical frequencies, where propagation losses can be minimised. Among the various physical platforms that have realized coherent microwave-to-optical transduction, those that use atoms as transducers have shown rapid progress in recent years. In this paper we report an experimental demonstration of coherent microwave-to-optical conversion that maps a microwave signal to a large, tunable 550(30) MHz range of optical frequencies using room-temperature $^{87}$Rb atoms. The inhomogeneous Doppler broadening of the atomic vapor advantageously supports the tunability of an input microwave channel to any optical frequency channel within the Doppler width, along with simultaneous conversion of a multi-channel input microwave field to corresponding optical channels. In addition, we demonstrate phase-correlated amplitude control of select channels, resulting in complete extinction of one of the channels, providing an analog to a frequency domain beam splitter across five orders of magnitude in frequency. With frequency-division multiplexing capability, multi-channel conversion, and amplitude control of frequency channels, neutral atomic systems may be effective quantum processors for quantum information encoded in frequency-bin qubits. △ Less

Submitted 4 December, 2023; v1 submitted 30 May, 2023; originally announced May 2023.

Journal ref: Commun Phys 6, 338 (2023)

The fate of baryons in counterfactual universes

Authors: Boon Kiat Oh, John A. Peacock, Sadegh Khochfar, Britton D. Smith

Abstract: We present results from nine simulations that compare the standard $Λ$ Cold Dark Matter cosmology ($Λ$CDM) with counterfactual universes, for approximately $100\,{\rm Gyr}$ using the Enzo simulation code. We vary the value of $Λ$ and the fluctuation amplitude to explore the effect on the evolution of the halo mass function (HMF), the intergalactic medium (IGM) and the star formation history (SFH).… ▽ More We present results from nine simulations that compare the standard $Λ$ Cold Dark Matter cosmology ($Λ$CDM) with counterfactual universes, for approximately $100\,{\rm Gyr}$ using the Enzo simulation code. We vary the value of $Λ$ and the fluctuation amplitude to explore the effect on the evolution of the halo mass function (HMF), the intergalactic medium (IGM) and the star formation history (SFH). The distinct peak in star formation rate density (SFRD) and its subsequent decline are both affected by the interplay between gravitational attraction and the accelerating effects of $Λ$. The IGM cools down more rapidly in models with a larger $Λ$ and also with a lower $σ_8$, reflecting the reduced SFRD associated with these changes -- although changing $σ_8$ is not degenerate with changing $Λ$, either regarding the thermal history of the IGM or the SFH. However, these induced changes to the IGM or ionizing background have little impact on the calculated SFRD. We provide fits for the evolution of the SFRD in these different universes, which we integrate over time to derive an asymptotic star formation efficiency. Together with Weinberg's uniform prior on $Λ$, the estimated probability of observers experiencing a value of $Λ$ no greater than the observed value is 13%, substantially larger than some alternative estimates. Within the Enzo model framework, then, observer selection within a multiverse is able to account statistically for the small value of the cosmological constant, although $Λ$ in our universe does appear to be at the low end of the predicted range. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 18 pages, 10 figures, MNRAS: MN-22-1662-MJ.R1

Figuring Out Gas & Galaxies In Enzo (FOGGIE) VI: The Circumgalactic Medium of $L^*$ Galaxies is Supported in an Emergent, Non-Hydrostatic Equilibrium

Authors: Cassandra Lochhaas, Jason Tumlinson, Molly S. Peeples, Brian W. O'Shea, Jessica K. Werk, Raymond C. Simons, James Juno, Claire E. Kopenhafer, Ramona Augustin, Anna C. Wright, Ayan Acharyya, Britton D. Smith

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

Submitted 7 March, 2023; v1 submitted 20 June, 2022; originally announced June 2022.

Comments: 36 pages, 24 figures. Accepted to ApJ. See https://foggie.science/cgm_forces.html for animated Figures 6 and 7

The AGORA High-resolution Galaxy Simulations Comparison Project. III: Cosmological zoom-in simulation of a Milky Way-mass halo

Authors: Santi Roca-Fàbrega, Ji-hoon Kim, Loic Hausammann, Kentaro Nagamine, Johnny W. Powell, Ikkoh Shimizu, Daniel Ceverino, Alessandro Lupi, Joel R. Primack, Thomas Quinn, Yves Revaz, Héctor Velázquez, Tom Abel, Michael Buehlmann, Avishai Dekel, Bili Dong, Oliver Hahn, Cameron B. Hummels, Ki-won Kim, Britton D. Smith, Clayton J. Strawn, Romain Teyssier, Matthew Turk

Abstract: We present a suite of high-resolution cosmological zoom-in simulations to $z=4$ of a $10^{12}\,{\rm M}_{\odot}$ halo at $z=0$, obtained using seven contemporary astrophysical simulation codes widely used in the numerical galaxy formation community. Physics prescriptions for gas cooling, heating, and star formation, are similar to the ones used in our previous {\it AGORA} disk comparison but now ac… ▽ More We present a suite of high-resolution cosmological zoom-in simulations to $z=4$ of a $10^{12}\,{\rm M}_{\odot}$ halo at $z=0$, obtained using seven contemporary astrophysical simulation codes widely used in the numerical galaxy formation community. Physics prescriptions for gas cooling, heating, and star formation, are similar to the ones used in our previous {\it AGORA} disk comparison but now account for the effects of cosmological processes. In this work, we introduce the most careful comparison yet of galaxy formation simulations run by different code groups, together with a series of four calibration steps each of which is designed to reduce the number of tunable simulation parameters adopted in the final run. After all the participating code groups successfully completed the calibration steps, we reach a suite of cosmological simulations with similar mass assembly histories down to $z=4$. With numerical accuracy that resolves the internal structure of a target halo, we find that the codes overall agree well with one another in e.g., gas and stellar properties, but also show differences in e.g., circumgalactic medium properties. We argue that, if adequately tested in accordance with our proposed calibration steps and common parameters, the results of high-resolution cosmological zoom-in simulations can be robust and reproducible. New code groups are invited to join this comparison by generating equivalent models by adopting the common initial conditions, the common easy-to-implement physics package, and the proposed calibration steps. Further analyses of the simulations presented here will be in forthcoming reports from our Collaboration. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Comments: Accepted for publication to the ApJ. The simulation snapshots used in this paper will be made publicly available in a user-friendly interface. The announcement of the data release will be posted on the AGORA Project webpage (www.agorasimulations.org). Early access to the data can be provided upon request to the authors, or to the project coordinator (sroca01@ucm.es, santacruzgalaxy@gmail.com)

Evolving beyond z=0: insights about the future of stars and the intergalactic medium

Authors: Boon Kiat Oh, John A. Peacock, Sadegh Khochfar, Britton D. Smith

Abstract: We present results from seven cosmological simulations that have been extended beyond the present era as far as redshift $z=-0.995$ or $t\approx96\,{\rm Gyr}$, using the Enzo simulation code. We adopt the calibrated star formation and feedback prescriptions from our previous work on reproducing the Milky Way with Enzo with modifications to the simulation code, chemistry and cooling library. We the… ▽ More We present results from seven cosmological simulations that have been extended beyond the present era as far as redshift $z=-0.995$ or $t\approx96\,{\rm Gyr}$, using the Enzo simulation code. We adopt the calibrated star formation and feedback prescriptions from our previous work on reproducing the Milky Way with Enzo with modifications to the simulation code, chemistry and cooling library. We then consider the future behaviour of the halo mass function (HMF), the equation of state (EOS) of the IGM, and the cosmic star formation history (SFH). Consistent with previous work, we find a freeze-out in the HMF at $z\approx-0.6$. The evolution of the EOS of the IGM presents an interesting case study of the cosmological coincidence problem, where there is a sharp decline in the IGM temperature immediately after $z=0$. For the SFH, the simulations produce a peak and a subsequent decline into the future. However, we do find a turnaround in the SFH after $z\approx-0.98$ in some simulations, probably due to the limitations of the criteria used for star formation. By integrating the SFH in time up to $z=-0.92$, the simulation with the best spatial resolution predicts an asymptotic total stellar mass that is very close to that obtained from extrapolating the fit of the observed SFR. Lastly, we investigate the future evolution of the partition of baryons within a Milky Way-sized galaxy, using both a zoom and a box simulation. Despite vastly different resolutions, these simulations predict individual haloes containing an equal fraction of baryons in stars and gas at the time of freeze-out ($t\approx30\,{\rm Gyr}$). △ Less

Submitted 3 March, 2021; originally announced March 2021.

Figuring Out Gas & Galaxies In Enzo (FOGGIE) V: The Virial Temperature Does Not Describe Gas in a Virialized Galaxy Halo

Authors: Cassandra Lochhaas, Jason Tumlinson, Brian W. O'Shea, Molly S. Peeples, Britton D. Smith, Jessica K. Werk, Ramona Augustin, Raymond C. Simons

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

Submitted 7 September, 2021; v1 submitted 16 February, 2021; originally announced February 2021.

Comments: 30 pages, 14 figures, accepted to ApJ

Storing short single-photon-level optical pulses in Bose-Einstein condensates for high-performance quantum memory

Authors: Erhan Saglamyurek, Taras Hrushevskyi, Anindya Rastogi, Logan W. Cooke, Benjamin D. Smith, Lindsay J. LeBlanc

Abstract: Large-scale quantum networks require quantum memories featuring long-lived storage of non-classical light together with efficient, high-speed and reliable operation. The concurrent realization of these features is challenging due to inherent limitations of matter platforms and light-matter interaction protocols. Here, we propose an approach to overcome this obstacle, based on the implementation of… ▽ More Large-scale quantum networks require quantum memories featuring long-lived storage of non-classical light together with efficient, high-speed and reliable operation. The concurrent realization of these features is challenging due to inherent limitations of matter platforms and light-matter interaction protocols. Here, we propose an approach to overcome this obstacle, based on the implementation of the Autler-Townes-splitting (ATS) quantum-memory protocol on a Bose-Einstein condensate (BEC) platform. We demonstrate a proof-of-principle of this approach by storing short pulses of single-photon-level light as a collective spin-excitation in a rubidium BEC. For 20 ns long-pulses, we achieve an ultra-low-noise memory with an efficiency of 30% and lifetime of 15 $μ$s. The non-adiabatic character of the ATS protocol (leading to high-speed and low-noise operation) in combination with the intrinsically large atomic densities and ultra-low temperatures of the BEC platform (offering highly efficient and long-lived storage) opens up a new avenue towards high-performance quantum memories. △ Less

Submitted 29 October, 2020; originally announced October 2020.

Comments: 7 pages and 7 figures, include Supplementary Info

GPU-accelerated solutions of the nonlinear Schrödinger equation for simulating 2D spinor BECs

Authors: Benjamin D. Smith, Logan W. Cooke, Lindsay J. LeBlanc

Abstract: As a first approximation beyond linearity, the nonlinear Schrödinger equation (NLSE) reliably describes a broad class of physical systems. Though numerical solutions of this model are well-established, these methods can be computationally complex. In this paper, we showcase a code development approach, demonstrating how computational time can be significantly reduced with readily available graphic… ▽ More As a first approximation beyond linearity, the nonlinear Schrödinger equation (NLSE) reliably describes a broad class of physical systems. Though numerical solutions of this model are well-established, these methods can be computationally complex. In this paper, we showcase a code development approach, demonstrating how computational time can be significantly reduced with readily available graphics processing unit (GPU) hardware and a straightforward code migration using open-source libraries. This process shows how CPU computations with power-law scaling in computation time with grid size can be made linear using GPUs. As a specific case study, we investigate the Gross-Pitaevskii equation, a specific version of the nonlinear Schrödinger model, as it describes in two dimensions a trapped, interacting, two-component Bose-Einstein condensate (BEC) subject to a spatially dependent interspin coupling, resulting in an analog to a spin-Hall system. This computational approach lets us probe high-resolution spatial features - revealing an interaction-dependent phase transition - all in a reasonable amount of time. Our computational approach is particularly relevant for research groups looking to easily accelerate straightforward numerical simulation of physical phenomena. △ Less

Submitted 23 March, 2022; v1 submitted 28 October, 2020; originally announced October 2020.

Comments: 29 pages, 4 figures

Journal ref: Comput. Phys. Commun. 275, 108314 (2022)

External Enrichment of Minihalos by the First Supernovae

Authors: William Hicks, Azton Wells, Michael L. Norman, John H. Wise, Britton D. Smith, Brian W. O'Shea

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

Submitted 10 January, 2021; v1 submitted 11 September, 2020; originally announced September 2020.

Comments: 15 pages, 12 figures, Accepted for publication in ApJ

Figuring Out Gas & Galaxies In Enzo (FOGGIE). IV. The Stochasticity of Ram Pressure Stripping in Galactic Halos

Authors: Raymond C. Simons, Molly S. Peeples, Jason Tumlinson, Brian W. O'Shea, Britton D. Smith, Lauren Corlies, Cassandra Lochhaas, Yong Zheng, Ramona Augustin, Deovrat Prasad, Gregory F. Snyder, Erik Tollerud

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

Submitted 29 April, 2020; originally announced April 2020.

Comments: 18 pages, 10 figures. Submitted to ApJ

Calibration of a star formation and feedback model for cosmological simulations with Enzo

Authors: Boon Kiat Oh, Britton D. Smith, John A. Peacock, Sadegh Khochfar

Abstract: We present results from seventy-one zoom simulations of a Milky Way-sized (MW) halo, exploring the parameter space for a widely-used star formation and feedback model in the {\tt Enzo} simulation code. We propose a novel way to match observations, using functional fits to the observed baryon makeup over a wide range of halo masses. The model MW galaxy is calibrated using three parameters: the star… ▽ More We present results from seventy-one zoom simulations of a Milky Way-sized (MW) halo, exploring the parameter space for a widely-used star formation and feedback model in the {\tt Enzo} simulation code. We propose a novel way to match observations, using functional fits to the observed baryon makeup over a wide range of halo masses. The model MW galaxy is calibrated using three parameters: the star formation efficiency $\left(f_*\right)$, the efficiency of thermal energy from stellar feedback $\left(ε\right)$ and the region into which feedback is injected $\left(r\ {\rm and}\ s\right)$. We find that changing the amount of feedback energy affects the baryon content most significantly. We then identify two sets of feedback parameter values that are both able to reproduce the baryonic properties for haloes between $10^{10}\,\mathrm{M_\odot}$ and $10^{12}\,\mathrm{M_\odot}$. We can potentially improve the agreement by incorporating more parameters or physics. If we choose to focus on one property at a time, we can obtain a more realistic halo baryon makeup. We show that the employed feedback prescription is insensitive to dark matter mass resolution between $10^5\,{\rm M_\odot}$ and $10^7\,{\rm M_\odot}$. Contrasting both star formation criteria and the corresponding combination of optimal feedback parameters, we also highlight that feedback is self-consistent: to match the same baryonic properties, with a relatively higher gas to stars conversion efficiency, the feedback strength required is lower, and vice versa. Lastly, we demonstrate that chaotic variance in the code can cause deviations of approximately 10\% and 25\% in the stellar and baryon mass in simulations evolved from identical initial conditions. △ Less

Submitted 7 February, 2020; originally announced February 2020.

Figuring Out Gas & Galaxies in Enzo (FOGGIE). III. The Mocky Way: Investigating Biases in Observing the Milky Way's Circumgalactic Medium

Authors: Yong Zheng, Molly S. Peeples, Brian W. O'Shea, Raymond C. Simons, Cassandra Lochhaas, Lauren Corlies, Jason Tumlinson, Britton D. Smith, Ramona Augustin

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

Submitted 19 June, 2020; v1 submitted 21 January, 2020; originally announced January 2020.

Comments: ApJ in press

The Impact of Enhanced Halo Resolution on the Simulated Circumgalactic Medium

Authors: Cameron B. Hummels, Britton D. Smith, Philip F. Hopkins, Brian W. O'Shea, Devin W. Silvia, Jessica K. Werk, Nicolas Lehner, John H. Wise, David C. Collins, Iryna S. Butsky

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

Submitted 29 November, 2018; originally announced November 2018.

Comments: 19 pages, 13 figures, submitted to ApJ. Simulation movies available at http://chummels.org/tempest

Imprints of the first billion years: Lyman limit systems at $z \sim 5$

Authors: Neil H. M. Crighton, J. Xavier Prochaska, Michael T. Murphy, John M. O'Meara, Gabor Worseck, Britton D. Smith

Abstract: Lyman Limit systems (LLSs) trace the low-density circumgalactic medium and the most dense regions of the intergalactic medium, so their number density and evolution at high redshift, just after reionisation, are important to constrain. We present a survey for LLSs at high redshifts, $z_{\rm LLS} =3.5$--5.4, in the homogeneous dataset of 153 optical quasar spectra at $z \sim 5$ from the Giant Gemin… ▽ More Lyman Limit systems (LLSs) trace the low-density circumgalactic medium and the most dense regions of the intergalactic medium, so their number density and evolution at high redshift, just after reionisation, are important to constrain. We present a survey for LLSs at high redshifts, $z_{\rm LLS} =3.5$--5.4, in the homogeneous dataset of 153 optical quasar spectra at $z \sim 5$ from the Giant Gemini GMOS survey. Our analysis includes detailed investigation of survey biases using mock spectra which provide important corrections to the raw measurements. We estimate the incidence of LLSs per unit redshift at $z \approx 4.4$ to be $\ell(z) = 2.6 \pm 0.4$. Combining our results with previous surveys at $z_{\rm LLS} <4$, the best-fit power-law evolution is $\ell(z) = \ell_* [(1+z)/4]^α$ with $\ell_* = 1.46 \pm 0.11$ and $α= 1.70 \pm 0.22$ (68\% confidence intervals). Despite hints in previous $z_{\rm LLS} <4$ results, there is no indication for a deviation from this single power-law soon after reionization. Finally, we integrate our new results with previous surveys of the intergalactic and circumgalactic media to constrain the hydrogen column density distribution function, $f(N_{\rm HI},X)$, over 10 orders of magnitude. The data at $z \sim 5$ are not well described by the $f(N_{\rm HI},X)$ model previously reported for $z \sim 2$--3 (after re-scaling) and a 7-pivot model fitting the full $z \sim 2$--5 dataset is statistically unacceptable. We conclude that there is significant evolution in the shape of $f(N_{\rm HI},X)$ over this $\sim$2 billion year period. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 16 pages, 19 figures, accepted to MNRAS

Figuring Out Gas & Galaxies in Enzo (FOGGIE). I. Resolving Simulated Circumgalactic Absorption at 2 < z < 2.5

Authors: Molly S. Peeples, Lauren Corlies, Jason Tumlinson, Brian W. O'Shea, Nicolas Lehner, John M. O'Meara, J. Christopher Howk, Britton D. Smith, John H. Wise, Cameron B. Hummels

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

Submitted 14 March, 2019; v1 submitted 15 October, 2018; originally announced October 2018.

Comments: ApJ, in press. New appendix with results from simulation with 78pc resolution CGM at z=2.5 now included. Appendix C with additional velocity phase diagrams has been omitted to fit within the arXiv space restraints

Journal ref: ApJ (2019), 873, 2

Formation of First Galaxies inside Density Peaks and Voids under the Influence of Dark Matter-Baryon Streaming Velocity, I: Initial Condition and Simulation Scheme

Authors: Kyungjin Ahn, Britton D. Smith

Abstract: We present a systematic study of the cosmic variance that existed in the formation of first stars and galaxies. We focus on the cosmic variance induced by the large-scale density and velocity environment engraved at the epoch of recombination. The density environment is predominantly determined by the dark-matter overdensity, and the velocity environment by the dark matter-baryon streaming velocit… ▽ More We present a systematic study of the cosmic variance that existed in the formation of first stars and galaxies. We focus on the cosmic variance induced by the large-scale density and velocity environment engraved at the epoch of recombination. The density environment is predominantly determined by the dark-matter overdensity, and the velocity environment by the dark matter-baryon streaming velocity. Toward this end, we introduce a new cosmological initial condition generator BCCOMICS, which solves the quasi-linear evolution of small-scale perturbations under the large-scale density and streaming-velocity environment and generates the initial condition for dark matter and baryons, as either particles or grid data at a specific redshift. We also describe a scheme to simulate the formation of first galaxies inside density peaks and voids, where a local environment is treated as a separate universe. The resulting cosmic variance in the minihalo number density and the amount of cooling mass are presented as an application. Density peaks become a site for enhanced formation of first galaxies, which compete with the negative effect from the dark matter-baryon streaming velocity on structure formation. △ Less

Submitted 28 November, 2018; v1 submitted 11 July, 2018; originally announced July 2018.

Comments: in press for publication in ApJ; 16 pages, 7 figures, 1 table; BCCOMICS repository specified

Probing the Dependence of the Intergalactic Medium on Large Scale Environment Using the Low Redshift Lyman Alpha Forest

Authors: Stephanie Tonnesen, Britton D. Smith, Juna Kollmeier, Renyue Cen

Abstract: We examine the statistics of the low-redshift Ly-alpha forest in an adaptive mesh refinement hydrodynamic cosmological simulation of sufficient volume to include distinct large-scale environments. We compare our HI column density distribution of absorbers both with recent work and between two highly-refined regions of our simulation: a large-scale overdensity and a large-scale underdensity (on sca… ▽ More We examine the statistics of the low-redshift Ly-alpha forest in an adaptive mesh refinement hydrodynamic cosmological simulation of sufficient volume to include distinct large-scale environments. We compare our HI column density distribution of absorbers both with recent work and between two highly-refined regions of our simulation: a large-scale overdensity and a large-scale underdensity (on scales of approximately 20 Mpc). We recover the average results presented in Kollmeier et al. (2014) using different simulation methods. We further break down these results as a function of environment to examine the detailed dependence of absorber statistics on large-scale density. We find that the slope of the HI column density distribution in the 10$^{12.5}$ $\le$ N$_{HI}$/cm$^{-2}$ $\le$ 10$^{14.5}$ range depends on environment such that the slope becomes steeper for higher environmental density, and this difference reflects distinct physical conditions of the intergalactic medium on these scales. We track this difference to the different temperature structures of filaments in varying environments. Specifically, filaments in the overdensity are hotter and, correspondingly, are composed of gas with lower HI fractions than those in underdense environments. Our results highlight that in order to understand the physics driving the HI CDD, we need not only improved accounting of the sources of ionizing UV photons, but also of the physical conditions of the IGM and how this may vary as a function of large-scale environment. △ Less

Submitted 11 July, 2017; originally announced July 2017.

Comments: 9 Figures, Accepted in ApJ

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.

Comments: 20 pages, 8 figures, accepted for publication in MNRAS. For more info, visit grackle.readthedocs.io

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.

Comments: 28 pages, 35 figures, Accepted for publication in the Astrophysical Journal, Image resolution greatly reduced, High-resolution version of this article is available at http://www.jihoonkim.org/agora/AGORA_Paper4_draft.pdf, The first paper of the AGORA Initiative is at http://adsabs.harvard.edu/abs/2014ApJS..210...14K, More information on AGORA is at http://www.AGORAsimulations.org/

Journal ref: Astrophysical Journal, 833, 202, 2016

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.

Comments: 21 pages, 15 figures, 4 tables. Accepted in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 2014 442 (2): 2560-2579

An HST/COS Survey of the Low-Redshift IGM. I. Survey, Methodology, & Overall Results

Authors: Charles W. Danforth, Brian A. Keeney, Evan M. Tilton, J. Michael Shull, Matthew Stevans, Matthew M. Pieri, John T. Stocke, Blair D. Savage, Kevin France, David Syphers, Britton D. Smith, James C. Green, Cynthia Froning, Steven V. Penton, Steven N. Osterman

Abstract: We use high-quality, medium-resolution {\it Hubble Space Telescope}/Cosmic Origins Spectrograph (\HST/COS) observations of 82 UV-bright AGN at redshifts $z_{AGN}<0.85$ to construct the largest survey of the low-redshift intergalactic medium (IGM) to date: 5343 individual extragalactic absorption lines in HI and 25 different metal-ion species grouped into 2610 distinct redshift systems at… ▽ More We use high-quality, medium-resolution {\it Hubble Space Telescope}/Cosmic Origins Spectrograph (\HST/COS) observations of 82 UV-bright AGN at redshifts $z_{AGN}<0.85$ to construct the largest survey of the low-redshift intergalactic medium (IGM) to date: 5343 individual extragalactic absorption lines in HI and 25 different metal-ion species grouped into 2610 distinct redshift systems at $z_{abs}<0.75$ covering total redshift pathlengths $Δz_{HI}=21.7$ and $Δz_{OVI}=14.5$. Our semi-automated line-finding and measurement technique renders the catalog as objectively-defined as possible. The cumulative column-density distribution of HI systems can be parametrized $dN(>N)/dz=C_{14}(N/10^{14} cm^{-2})^{-(β-1)}$, with $C_{14}=25\pm1$ and $β=1.65\pm0.02$. This distribution is seen to evolve both in amplitude, $C_{14}\sim(1+z)^{2.0\pm0.1}$, and slope $β(z)=1.73-0.26 z$ for $z<0.47$. We observe metal lines in 427 systems, and find that the fraction of IGM absorbers detected in metals is strongly dependent on N_{HI}. The distribution of OVI absorbers appear to evolve in the same sense as the Lya forest. We calculate contributions to $Ω_b$ from different components of the low-$z$ IGM and determine the Lya decrement as a function of redshift. IGM absorbers are analyzed via a two-point correlation function (TPCF) in velocity space. We find substantial clustering of \HI\ absorbers on scales of $Δv=50-300$ km/s with no significant clustering at $Δv>1000$ km/s. Splitting the sample into strong and weak absorbers, we see that most of the clustering occurs in strong, $N_{HI}>10^{13.5} cm^{-2}$, metal-bearing IGM systems. The full catalog of absorption lines and fully-reduced spectra is available via MAST as a high-level science product at http://archive.stsci.edu/prepds/igm/. △ Less

Submitted 21 December, 2015; v1 submitted 11 February, 2014; originally announced February 2014.

Comments: This is the accepted version (v3) of the paper. Previous versions (July 2015 and Feb. 2014) should be replaced by this one. In particular, please note that the associated MAST high-level-science product has been updated to reflect the of the final state of the paper. It is available at: http://archive.stsci.edu/prepds/igm/

Fragmentation in Dusty Low-Metallicity Star Forming Halos

Authors: Gregory R. Meece, Britton D. Smith, Brian W. O'Shea

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

Submitted 16 September, 2013; originally announced September 2013.

Comments: 28 pages, 25 figures, submitted to ApJ

Journal ref: Gregory R. Meece et al. 2014 ApJ 783 75

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.

Comments: 21 pages, 6 figures, Accepted for publication in the Astrophysical Journal Supplement, Image resolution greatly reduced, High-resolution version of this article and more information about the AGORA Project including the science goals of the 13 Working Groups are available at http://www.AGORAsimulations.org/ and at http://sites.google.com/site/santacruzcomparisonproject/details/

Journal ref: Astrophysical Journal Supplement Series, 210, 14, 2014

Bringing Simulation and Observation Together to Better Understand the Intergalactic Medium

Authors: Hilary Egan, Britton D. Smith, Brian W. O'Shea, J. Michael Shull

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

Submitted 18 July, 2014; v1 submitted 8 July, 2013; originally announced July 2013.

Comments: 22 pages, 21 figures. Accepted by The Astrophysical Journal

Cosmological Simulations of Isotropic Conduction in Galaxy Clusters

Authors: Britton D. Smith, Brian W. O'Shea, G. Mark Voit, David Ventimiglia, Samuel W. Skillman

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

Submitted 24 June, 2013; originally announced June 2013.

Comments: 13 pages, 13 figures. Submitted to The Astrophysical Journal

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.

Comments: 17 pages, 13 figures. Submitted to ApJ

On The Road To More Realistic Galaxy Cluster Simulations: The Effects of Radiative Cooling and Thermal Feedback Prescriptions on the Observational Properties of Simulated Galaxy Clusters

Authors: Stephen Skory, Eric Hallman, Jack O. Burns, Samuel W. Skillman, Brian W. O'Shea, Britton D. Smith

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

Submitted 13 November, 2012; originally announced November 2012.

Comments: 20 pages, 16 figures, 2 tables

The effect of feedback and reionization on star formation in low-mass dwarf galaxy haloes

Authors: Christine M. Simpson, Greg L. Bryan, Kathryn V. Johnston, Britton D. Smith, Mordecai-Mark Mac Low, Sanjib Sharma, Jason Tumlinson

Abstract: We simulate the evolution of a 10^9 Msun dark matter halo in a cosmological setting with an adaptive-mesh refinement code as an analogue to local low luminosity dwarf irregular and dwarf spheroidal galaxies. The primary goal of our study is to investigate the roles of reionization and supernova feedback in determining the star formation histories of low mass dwarf galaxies. We include a wide range… ▽ More We simulate the evolution of a 10^9 Msun dark matter halo in a cosmological setting with an adaptive-mesh refinement code as an analogue to local low luminosity dwarf irregular and dwarf spheroidal galaxies. The primary goal of our study is to investigate the roles of reionization and supernova feedback in determining the star formation histories of low mass dwarf galaxies. We include a wide range of physical effects, including metal cooling, molecular hydrogen formation and cooling, photoionization and photodissociation from a metagalactic background, a simple prescription for self-shielding, star formation, and a simple model for supernova driven energetic feedback. We carry out simulations excluding each major effect in turn. We find that reionization is primarily responsible for expelling most of the gas in our simulations, but that supernova feedback is required to disperse the dense, cold gas in the core of the halo. Moreover, we show that the timing of reionization can produce an order of magnitude difference in the final stellar mass of the system. For our full physics run with reionization at z=9, we find a stellar mass of about 10^5 Msun at z=0, and a mass-to-light ratio within the half-light radius of approximately 130 Msun/Lsun, consistent with observed low-luminosity dwarfs. However, the resulting median stellar metallicity is 0.06 Zsun, considerably larger than observed systems. In addition, we find star formation is truncated between redshifts 4 and 7, at odds with the observed late time star formation in isolated dwarf systems but in agreement with Milky Way ultrafaint dwarf spheroidals. We investigate the efficacy of energetic feedback in our simple thermal-energy driven feedback scheme, and suggest that it may still suffer from excessive radiative losses, despite reaching stellar particle masses of about 100 Msun, and a comoving spatial resolution of 11 pc. △ Less

Submitted 18 March, 2013; v1 submitted 5 November, 2012; originally announced November 2012.

Comments: 26 pages, 14 figures, accepted for publication in MNRAS; additional explanation added to sections 2,3 and 4; Figure 4 updated; results unchanged

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.

Comments: 18 pages, 11 figures, replaced with accepted version, MNRAS. Minor changes with the conclusions unaffected

Journal ref: MNRAS, 2012, vol. 427, 311

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.

Comments: Published in the Proceedings of the PyHPC2011 workshop at Supercomputing 2011. yt can be found at http://yt-project.org/

The Baryon Census in a Multiphase Intergalactic Medium: 30% of the Baryons May Still Be Missing

Authors: J. Michael Shull, Britton D. Smith, Charles W. Danforth

Abstract: For low-redshift cosmology and galaxy formation rates, it is important to account for all the baryons synthesized in the Big Bang. Although galaxies and clusters contain 10% of the baryons, many more reside in the photoionized Lyman-alpha forest and shocked-heated warm-hot intergalactic medium (WHIM) at T = 10^5 to 10^7 K. Current tracers of WHIM at 10^5 to 10^6 K include the O VI 1032, 1038 absor… ▽ More For low-redshift cosmology and galaxy formation rates, it is important to account for all the baryons synthesized in the Big Bang. Although galaxies and clusters contain 10% of the baryons, many more reside in the photoionized Lyman-alpha forest and shocked-heated warm-hot intergalactic medium (WHIM) at T = 10^5 to 10^7 K. Current tracers of WHIM at 10^5 to 10^6 K include the O VI 1032, 1038 absorption lines, together with broad Lyman-alpha absorbers (BLAs) and EUV/X-ray absorption lines from Ne VIII, O VII, and O VIII. We improve the O VI baryon surveys with corrections for oxygen metallicity (Z/Zsun) and O VI ionization fraction (f_OVI) using cosmological simulations of heating, cooling, and metal transport in a density-temperature structured medium. Statistically, their product correlates with column density, (Z/Zsun)(f_OVI) = (0.015)(N_OVI/10^{14} cm^-2)^0.70. The N_OVI-weighted mean is 0.01, which doubles previous estimates of WHIM baryon content. We also reanalyze H I data from the Hubble Space Telescope, applying redshift corrections for absorber density, photoionizing background, and proper length, dl/dz. We find substantial baryon fractions in the photoionized Lya forest (28 +/- 11%), O VI/BLA-traced WHIM (25 +/- 8%), and collapsed phase (18 +/- 4%) in galaxies, groups, clusters, and circumgalactic gas. The baryon shortfall is 29 +/- 13%, which may be detected in X-ray absorbers from hotter WHIM or in weaker Lya and O VI absorbers. Further progress will require higher-precision baryon surveys of weak absorbers at column densities N_HI > 10^{12.0} cm^-2, N_OVI > 10^{12.5} cm^-2, and N_OVII > 10^{14.5} cm^-2, with moderate-resolution UV and X-ray spectrographs. △ Less

Submitted 14 May, 2012; v1 submitted 12 December, 2011; originally announced December 2011.

Comments: 14 pages, 9 figures, revised submission to ApJ (new Fig 2, Appendices A,B)

Numerical Simulations of Supernova Dust Destruction. II. Metal-Enriched Ejecta Knots

Authors: D. W. Silvia, B. D. Smith, J. M. Shull

Abstract: Following our previous work, we investigate through hydrodynamic simulations the destruction of newly-formed dust grains by sputtering in the reverse shocks of supernova remnants. Using an idealized setup of a planar shock impacting a dense, spherical clump, we implant a population of Lagrangian particles into the clump to represent a distribution of dust grains in size and composition. We vary th… ▽ More Following our previous work, we investigate through hydrodynamic simulations the destruction of newly-formed dust grains by sputtering in the reverse shocks of supernova remnants. Using an idealized setup of a planar shock impacting a dense, spherical clump, we implant a population of Lagrangian particles into the clump to represent a distribution of dust grains in size and composition. We vary the relative velocity between the reverse shock and ejecta clump to explore the effects of shock-heating and cloud compression. Because supernova ejecta will be metal-enriched, we consider gas metallicities from Z/Zsun = 1 to 100 and their influence on cooling properties of the cloud and the thermal sputtering rates of embedded dust grains. We post-process the simulation output to calculate grain sputtering for a variety of species and size distributions. In the metallicity regime considered in this paper, the balance between increased radiative cooling and increased grain erosion depends on the impact velocity of the reverse shock. For slow shocks (velocity less than or equal to 3000 km/s), the amount of dust destruction is comparable across metallicities, or in some cases is decreased with increased metallicity. For higher shock velocities (velocity greater than or equal to 5000 km/s), an increase in metallicity from Z/Zsun = 10 to 100 can lead to an additional 24% destruction of the initial dust mass. While the total dust destruction varies widely across grain species and simulation parameters, our most extreme cases result in complete destruction for some grain species and only 44% dust mass survival for the most robust species. These survival rates are important in understanding how early supernovae contribute to the observed dust masses in high-redshift galaxies. △ Less

Submitted 19 January, 2012; v1 submitted 1 November, 2011; originally announced November 2011.

Comments: 10 pages, 6 figures, 2 tables, changes made to the text and figures as suggested by the anonymous referee, accepted by the Astrophysical Journal

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.

Comments: 18 pages, 6 figures, emulateapj format. Resubmitted to Astrophysical Journal Supplement Series with revisions from referee. yt can be found at http://yt.enzotools.org/

The Properties of X-ray Cold Fronts in a Statistical Sample of Simulated Galaxy Clusters

Authors: Eric J. Hallman, Samuel W. Skillman, Tesla E. Jeltema, Britton D. Smith, Brian W. O'Shea, Jack O. Burns, Michael L. Norman

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

Submitted 15 October, 2010; originally announced October 2010.

Comments: 16 pages, 21 figures, Accepted to ApJ

The Nature of the Warm/Hot Intergalactic Medium I. Numerical Methods, Convergence, and OVI Absorption

Authors: Britton D. Smith, Eric J. Hallman, J. Michael Shull, Brian W. O'Shea

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

Submitted 4 February, 2011; v1 submitted 1 September, 2010; originally announced September 2010.

Comments: 22 pages, 21 figures, emulateapj, accepted for publication in ApJ

Protostellar Feedback Processes and the Mass of the First Stars

Authors: Jonathan C. Tan, Britton D. Smith, Brian W. O'Shea

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

Submitted 18 August, 2010; originally announced August 2010.

Comments: 6 pages, Proceedings of 'The First Stars and Galaxies: Challenges for the Next Decade", Austin, TX, March 8-11, 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.

Comments: 20 pages, 15 figures, further discussion and appendix added, accepted to ApJ

How well do cosmological simulations reproduce individual-halo properties?

Authors: M. Trenti, B. D. Smith, E. J. Hallman, S. W. Skillman, J. M. Shull

Abstract: Cosmological simulations of galaxy formation often rely on prescriptions for star formation and feedback that depend on halo properties such as halo mass, central over-density, and virial temperature. In this paper we address the convergence of individual halo properties, based on their number of particles N, focusing in particular on the mass of halos near the resolution limit of a simulation.… ▽ More Cosmological simulations of galaxy formation often rely on prescriptions for star formation and feedback that depend on halo properties such as halo mass, central over-density, and virial temperature. In this paper we address the convergence of individual halo properties, based on their number of particles N, focusing in particular on the mass of halos near the resolution limit of a simulation. While it has been established that the halo mass function is sampled on average down to N~30 particles, we show that individual halo properties exhibit significant scatter, and some systematic biases, as one approaches the resolution limit. We carry out a series of cosmological simulations using the Gadget2 and Enzo codes with N_p=64^3 to N_p=1024^3 total particles, keeping the same large-scale structure in the simulation box. We consider boxes from l_{box} = 8 Mpc/h to l_{box} = 512 Mpc/h to probe different halo masses and formation redshifts. We cross-identify dark matter halos in boxes at different resolutions and measure the scatter in their properties. The uncertainty in the mass of single halos depends on the number of particles (scaling approximately as N^{-1/3}), but the rarer the density peak, the more robust its identification. The virial radius of halos is very stable and can be measured without bias for halos with N>30. In contrast, the average density within a sphere containing 25% of the total halo mass is severely underestimated (by more than a factor 2) and the halo spin is moderately overestimated for N<100. If sub-grid physics is implemented upon a cosmological simulation, we recommend that rare halos (~3sigma peaks) be resolved with N>100 particles and common halos (~1sigma peaks) with N>400 particles to avoid excessive numerical noise and possible systematic biases in the results. △ Less

Submitted 27 January, 2010; originally announced January 2010.

Comments: 25 pages, 9 figures, ApJ accepted

Numerical Simulations of Supernova Dust Destruction. I. Cloud-crushing and Post-processed Grain Sputtering

Authors: D. W. Silvia, B. D. Smith, J. M. Shull

Abstract: We investigate through hydrodynamic simulations the destruction of newly-formed dust grains by sputtering in the reverse shocks of supernova remnants. Using an idealized setup of a planar shock impacting a dense, spherical clump, we implant a population of Lagrangian particles into the clump to represent a distribution of dust grains in size and composition. We then post-process the simulation ou… ▽ More We investigate through hydrodynamic simulations the destruction of newly-formed dust grains by sputtering in the reverse shocks of supernova remnants. Using an idealized setup of a planar shock impacting a dense, spherical clump, we implant a population of Lagrangian particles into the clump to represent a distribution of dust grains in size and composition. We then post-process the simulation output to calculate the grain sputtering for a variety of species and size distributions. We explore the parameter space appropriate for this problem by altering the over-density of the ejecta clumps and the speed of the reverse shocks. Since radiative cooling could lower the temperature of the medium in which the dust is embedded and potentially protect the dust by slowing or halting grain sputtering, we study the effects of different cooling methods over the time scale of the simulations. In general, our results indicate that grains with radii less than 0.1 microns are sputtered to much smaller radii and often destroyed completely, while larger grains survive their interaction with the reverse shock. We also find that, for high ejecta densities, the percentage of dust that survives is strongly dependent on the relative velocity between the clump and the reverse shock, causing up to 50% more destruction for the highest velocity shocks. The fraction of dust destroyed varies widely across grain species, ranging from total destruction of Al2O3 grains to minimal destruction of Fe grains (only 20% destruction in the most extreme cases). C and SiO2 grains show moderate to strong sputtering as well, with 38% and 80% mass loss. The survival rate of grains formed by early supernovae is crucial in determining whether or not they can act as the "dust factories" needed to explain high-redshift dust. △ Less

Submitted 20 April, 2010; v1 submitted 26 January, 2010; originally announced January 2010.

Comments: 18 pages, 16 figures, 2 tables, changes make to the text and figures as suggested by the anonymous referee, accepted by the Astrophysical Journal

The Santa Fe Light Cone Simulation Project: II. The Prospects for Direct Detection of the WHIM with SZE Surveys

Authors: Eric J. Hallman, Brian W. O'Shea, Britton D. Smith, Jack O. Burns, Michael L. Norman

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

Submitted 18 March, 2009; originally announced March 2009.

Comments: 9 pages, submitted to Astrophysical Journal

Journal ref: Astrophys.J.698:1795-1802,2009

The effect of photo-ionization on the cooling rates of enriched, astrophysical plasmas

Authors: Robert P. C. Wiersma, Joop Schaye, Britton D. Smith

Abstract: Radiative cooling is central to a wide range of astrophysical problems. Despite its importance, cooling rates are generally computed using very restrictive assumptions, such as collisional ionization equilibrium and solar relative abundances. We simultaneously relax both assumptions and investigate the effects of photo-ionization of heavy elements by the meta-galactic UV/X-ray background and of… ▽ More Radiative cooling is central to a wide range of astrophysical problems. Despite its importance, cooling rates are generally computed using very restrictive assumptions, such as collisional ionization equilibrium and solar relative abundances. We simultaneously relax both assumptions and investigate the effects of photo-ionization of heavy elements by the meta-galactic UV/X-ray background and of variations in relative abundances on the cooling rates of optically thin gas in ionization equilibrium. We find that photo-ionization by the meta-galactic background radiation reduces the net cooling rates by up to an order of magnitude for gas densities and temperatures typical of the shock-heated intergalactic medium and proto-galaxies. In addition, photo-ionization changes the relative contributions of different elements to the cooling rates. We conclude that photo-ionization by the ionizing background and heavy elements both need to be taken into account in order for the cooling rates to be correct to order of magnitude. Moreover, if the rates need to be known to better than a factor of a few, then departures of the relative abundances from solar need to be taken into account. We propose a method to compute cooling rates on an element-by-element basis by interpolating pre-computed tables that take photo-ionization into account. We provide such tables for a popular model of the evolving UV/X-ray background radiation, computed using the photo-ionization package CLOUDY. △ Less

Submitted 4 November, 2008; v1 submitted 23 July, 2008; originally announced July 2008.

Comments: 10 pages, 6 figures, Accepted for publication in MNRAS. One figure added and minor textual changes made to first version. Downloadable tables and videos available at http://www.strw.leidenuniv.nl/WSS08/

Journal ref: MNRAS 393 (2009) 99-107

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.

Comments: Significant revisions made, emulateapj, 15 pages, 3 tables, 8 figures, accepted for publication in ApJ. High-res version available at: http://solo.colorado.edu/~brittons/papers/threemodes.pdf

Journal ref: Astrophys.J.691:441-451,2009

Modeling of Emission Signatures of Massive Black Hole Binaries: I Methods

Authors: Tamara Bogdanovic, Britton D. Smith, Steinn Sigurdsson, Michael Eracleous

Abstract: We model the electromagnetic signatures of massive black hole binaries (MBHBs) with an associated gas component. The method comprises numerical simulations of relativistic binaries and gas coupled with calculations of the physical properties of the emitting gas. We calculate the UV/X-ray and the Halpha light curves and the Halpha emission profiles. The simulations are carried out with a modified… ▽ More We model the electromagnetic signatures of massive black hole binaries (MBHBs) with an associated gas component. The method comprises numerical simulations of relativistic binaries and gas coupled with calculations of the physical properties of the emitting gas. We calculate the UV/X-ray and the Halpha light curves and the Halpha emission profiles. The simulations are carried out with a modified version of the parallel tree SPH code Gadget. The heating, cooling, and radiative processes are calculated for two different physical scenarios, where the gas is approximated as a black-body or a solar metallicity gas. The calculation for the solar metallicity scenario is carried out with the photoionization code Cloudy. We focus on sub-parsec binaries which have not yet entered the gravitational radiation phase. The results from the first set of calculations, carried out for a coplanar binary and gas disk, suggest that there are pronounced outbursts in the X-ray light curve during pericentric passages. If such outbursts persist for a large fraction of the lifetime of the system, they can serve as an indicator of this type of binary. The predicted Halpha emission line profiles may be used as a criterion for selection of MBHB candidates from existing archival data. The orbital period and mass ratio of a binary may be inferred after carefully monitoring the evolution of the Halpha profiles of the candidates. The discovery of sub-parsec binaries is an important step in understanding of the merger rates of MBHBs and their evolution towards the detectable gravitational wave window. △ Less

Submitted 2 August, 2007; originally announced August 2007.

Comments: 28 pages, 29 figures, to appear in ApJS

Metal Cooling in Simulations of Cosmic Structure Formation

Authors: Britton D. Smith, Steinn Sigurdsson, Tom Abel

Abstract: The addition of metals to any gas can significantly alter its evolution by increasing the rate of radiative cooling. In star-forming environments, enhanced cooling can potentially lead to fragmentation and the formation of low-mass stars, where metal-free gas-clouds have been shown not to fragment. Adding metal cooling to numerical simulations has traditionally required a choice between speed an… ▽ More The addition of metals to any gas can significantly alter its evolution by increasing the rate of radiative cooling. In star-forming environments, enhanced cooling can potentially lead to fragmentation and the formation of low-mass stars, where metal-free gas-clouds have been shown not to fragment. Adding metal cooling to numerical simulations has traditionally required a choice between speed and accuracy. We introduce a method that uses the sophisticated chemical network of the photoionization software, Cloudy, to include radiative cooling from a complete set of metals up to atomic number 30 (Zn) that can be used with large-scale three-dimensional hydrodynamic simulations. Our method is valid over an extremely large temperature range (10 K < T < 10^8 K), up to hydrogen number densities of 10^12 cm^-3. At this density, a sphere of 1 Msun has a radius of roughly 40 AU. We implement our method in the adaptive mesh refinement (AMR) hydrodynamic/N-body code, Enzo. Using cooling rates generated with this method, we study the physical conditions that led to the transition from Population III to Population II star formation. While C, O, Fe, and Si have been previously shown to make the strongest contribution to the cooling in low-metallicity gas, we find that up to 40% of the metal cooling comes from fine-structure emission by S, when solar abundance patterns are present. At metallicities, Z > 10^-4 Zsun, regions of density and temperature exist where gas is both thermally unstable and has a cooling time less than its dynamical time. We identify these doubly unstable regions as the most inducive to fragmentation. At high redshifts, the CMB inhibits efficient cooling at low temperatures and, thus, reduces the size of the doubly unstable regions, making fragmentation more difficult. △ Less

Submitted 3 January, 2008; v1 submitted 5 June, 2007; originally announced June 2007.

Comments: 19 pages, 12 figures, significant revision, including new figures