BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Britton D. Smith

Britton D. Smith appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.IMastro-phastro-ph.SR

Catalog footprint

What is connected

22works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2026arXiv

Figuring Out Gas & Galaxies in Enzo (FOGGIE). XV. Examining the Spatial and Kinematic Relationship between Circumgalactic Mg II and O VI

Understanding the thermodynamic properties of the circumgalactic medium (CGM) is key to uncovering the baryon cycle in galaxies. Here we present spatial and kinematic relationships between Mg II and O VI as representatives for low and high ion-bearing gas, in the cosmological zoom-in galaxy simulation suite FOGGIE, a set of Milky-way-like galaxy simulations with high CGM resolution. We find the O VI-bearing gas exists as a diffuse halo around the galactic disk, while the Mg II-bearing gas is more centrally located. We investigate the covering fraction, probability of co-observation, co-kinematic correspondence of these ions using two different analysis methods. We make both mock sightlines using two-dimensional projections of our simulations treating these cells as integrated lines of sight and we create one-dimensional ray objects and use the SALSA (Boyd et al. 2020) code to investigate individual gas structures that contribute most to the line of sight column densities, which we call mock absorbers. We explore the relative kinematics of these mock absorbers and find Mg II and O VI appear to have a co-kinematic relationship when looking at absorber pairs with the closest relative velocity like in Werk et al. (2016). However, this does not necessarily correspond with a close spatial separation meaning many O VI and Mg II absorber pairs only appear to be co-kinematic but are physically unrelated. Taking a more holistic look at Mg II and O VI absorber pairs reveals a much weaker correlation between these two ions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Bringing Simulation and Observation Together to Better Understand the Intergalactic Medium

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Birth of a Galaxy - III. Propelling reionisation with the faintest galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cosmological Simulations of Isotropic Conduction in Galaxy Clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fragmentation in Dusty Low-Metallicity Star Forming Halos

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Population III Star Formation In Large Cosmological Simulations I. Halo Temporal and Physical Environment

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Birth of a Galaxy. II. The Role of Radiation Pressure

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Galaxy Cluster Radio Relics in Adaptive Mesh Refinement Cosmological Simulations: Relic Properties and Scaling Relationships

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

High-Performance Astrophysical Simulations and Analysis with Python

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Multi-Code Analysis Toolkit for Astrophysical Simulation Data

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Protostellar Feedback Processes and the Mass of the First Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Three Modes of Metal-Enriched Star Formation in the Early Universe

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)

0 citations0 reviewsNo code linkedOpen access