Source author record

Michael L. Norman

Michael L. Norman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph astro-ph.IM astro-ph.SR physics.flu-dyn astro-ph.EP

Catalog footprint

What is connected

64works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2021arXiv

Connecting Primordial Star Forming Regions and Second Generation Star Formation in the Phoenix Simulations

We introduce the {\em Phoenix Simulations}, a suite of highly resolved cosmological simulations featuring hydrodynamics, primordial gas chemistry, Population III and II star formation and feedback, UV radiative transfer, and saved outputs with $Δt$=200 kyr. The suite samples 73,523 distinct primordial star formation events within \npiii distinct regions, forming \ngii second-generation enriched star clusters by $z \geq 12$ within a cumulative 156.25 Mpc$^3$ volume. The regions that lead to enriched star formation contain up to $167$ primordial stars, with 78.7 \% of regions having experienced multiple types of primordial supernovae. The extent of a primordial region, measured by its metal-rich surrounding cloud, is highly variable: the average region has radius $\sim 3$ kpc, with 95 \% confidence limit on the distribution of measured radii is $\sim 5-7$ kpc. For continuing star formation, we find that the metallicity distribution of second generation stars is similar to that of subsequent Population II star formation, with both distributions spanning hyper metal-deficient ([Z/H]$\sim-7$) to super-solar ([Z/H]$\sim0.8$). We find that the metallicity of second generation stars has no strong dependence on the configuration of progenitor supernovae, with the mean metallicity of second-generation stars having $-1.73 < $[Z/H]$<-2.15$. Finally, we create an interpretable regression model to predict the radius of metal-rich influence of \piii star systems within the first 7-18 Myr after the first light. The model predicts the radius with $R_2 \gtrsim 0.4$ and mean squared error $\leq 0.06$. The probability distribution function of predicted radii compares well to that of observed radii with Jensen-Shannon distance $\lesssim 0.2$ for all modelled times.

preprint2020arXiv

Massive Star Formation in Metal-Enriched Haloes at High Redshift

The formation of supermassive stars has generally been studied under the assumption of rapid accretion of pristine metal-free gas. Recently it was found, however, that gas enriched to metallicities up to $Z \sim 10^{-3}$ Z$_{\odot}$ can also facilitate supermassive star formation, as long as the total mass infall rate onto the protostar remains sufficiently high. We extend the analysis further by examining how the abundance of supermassive star candidate haloes would be affected if all haloes with super-critical infall rates, regardless of metallicity were included. We investigate this scenario by identifying all atomic cooling haloes in the Renaissance simulations with central mass infall rates exceeding a fixed threshold. We find that among these haloes with central mass infall rates above 0.1 M$_{\odot}$ yr$^{-1}$ approximately two-thirds of these haloes have metallicities of $Z > 10^{-3}$ Z$_{\odot}$. If metal mixing within these haloes is inefficient early in their assembly and pockets of metal-poor gas can remain then the number of haloes hosting supermassive stars can be increased by at least a factor of four. Additionally the centres of these high infall-rate haloes provide ideal environments in which to grow pre-existing black holes. Further research into the (supermassive) star formation dynamics of rapidly collapsing haloes, with inhomogeneous metal distributions, is required to gain more insight into both supermassive star formation in early galaxies as well as early black hole growth.

preprint2020arXiv

The Emergence of the First Star-free Atomic Cooling Haloes in the Universe

Using the Renaissance suite of simulations we examine the emergence of pristine atomic cooling haloes that are both metal-free and star-free in the early Universe. The absence of metals prevents catastrophic cooling, suppresses fragmentation, and may allow for the formation of massive black hole seeds. Here we report on the abundance of pristine atomic cooling haloes found and on the specific physical conditions that allow for the formation of these direct-collapse-black-hole (DCBH) haloes. In total in our simulations we find that 79 DCBH haloes form before a redshift of 11.6. We find that the formation of pristine atomic haloes is driven by the rapid assembly of the atomic cooling haloes with mergers, both minor and/or major, prior to reaching the atomic cooling limit a requirement. However, the ability of assembling haloes to remain free of (external) metal enrichment is equally important and underlines the necessity of following the transport of metals in such simulations. The candidate DCBH hosting haloes we find, have been exposed to mean Lyman-Werner radiation fields of J$_{LW}$ $\sim$ 1 J$_{21}$ and typically lie at least 10 kpc (physical) from the nearest massive galaxy. Growth rates of the haloes reach values of greater than 10$^7$ M$_{\odot}$ per unit redshift, leading to significant dynamical heating and the suppression of efficient cooling until the halo crosses the atomic cooling threshold. Finally, we also find five synchronised halo candidates where pairs of pristine atomic cooling haloes emerge that are both spatially and temporally synchronised.

preprint2020arXiv

The Formation of Very Massive Stars in Early Galaxies and Implications for Intermediate Mass Black Holes

We investigate the ab-initio formation of super-massive stars in a pristine atomic cooling halo. The halo is extracted from a larger self-consistent parent simulation. The halo remains metal-free and star formation is suppressed due to a combination of dynamical heating from mergers and a mild ($J_{\rm LW} \sim 2 - 10 \ J_{21}$(z)) Lyman-Werner (LW) background. We find that more than 20 very massive stars form with stellar masses greater than 1000 M$_{\odot}$. The most massive star has a stellar mass of over 6000 M$_{\odot}$. However, accretion onto all stars declines significantly after the first $\sim$ 100 kyr of evolution as the surrounding material is accreted and the turbulent nature of the gas causes the stars to move to lower density regions. We post-process the impact of ionising radiation from the stars and find that ionising radiation is not a limiting factor when considering SMS formation and growth. Rather the birth environments are highly turbulent and a steady accretion flow is not maintained within the timescale (2 Myr) of our simulations. As the massive stars end their lives as direct collapse black holes this will seed these embryonic haloes with a population of black holes with masses between approximately 300 M$_{\odot}$ and 10,000 M$_{\odot}$. Afterwards they may sink to the centre of the haloes, eventually coalescing to form larger intermediate mass black holes whose in-situ mergers will be detectable by LISA.

preprint2016arXiv

Galaxy Properties and UV Escape Fractions During Epoch of Reionization: Results from the Renaissance Simulations

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

preprint2016arXiv

Late Pop III Star Formation During the Epoch of Reionization: Results from the Renaissance Simulations

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

preprint2016arXiv

Length Scales and Turbulent Properties of Magnetic Fields in Simulated Galaxy Clusters

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

preprint2015arXiv

Probing The Ultraviolet Luminosity Function of the Earliest Galaxies with the Renaissance Simulations

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

preprint2015arXiv

Spatially Extended 21 cm Signal from Strongly Clustered UV and X-Ray Sources in the Early Universe

We present our prediction for the local 21 cm differential brightness temperature ($δT_{b}$) from a set of strongly clustered sources of Population III (Pop III) and II (Pop II) objects in the early Universe, by a numerical simulation of their formation and radiative feedback. These objects are located inside a highly biased environment, which is a rare, high-density peak ("Rarepeak") extending to $\sim7$ comoving Mpc. We study the impact of ultraviolet (UV) and X-ray photons on the intergalactic medium (IGM) and the resulting $δT_{b}$, when Pop III stars are assumed to emit X-ray photons by forming X-ray binaries very efficiently. We parameterize the rest-frame spectral energy distribution (SED) of X-ray photons, which regulates X-ray photon-trapping, IGM-heating, secondary Lyman-alpha pumping and the resulting morphology of $δT_{b}$. A combination of emission ($δT_{b}>0$) and absorption ($δT_{b}<0$) regions appears in varying amplitudes and angular scales. The boost of the signal by the high-density environment ($δ\sim0.64$) and on a relatively large scale combine to make Rarepeak a discernible, spatially-extended ($θ\sim10'$) object for 21 cm observation at $13\lesssim z\lesssim17$, which is found to be detectable as a single object by SKA with integration time of $\sim1000$ hours. Power spectrum analysis by some of the SKA precursors (LOFAR, MWA, PAPER) of such rare peaks is found difficult due to the rarity of these peaks, and the contribution only by these rare peaks to the total power spectrum remains subdominant compared to that by all astrophysical sources.

preprint2014arXiv

Fully-Coupled Simulation of Cosmic Reionization. I: Numerical Methods and Tests

We describe an extension of the Enzo code to enable fully-coupled radiation hydrodynamical simulation of inhomogeneous reionization in large $\sim (100 Mpc)^3$ cosmological volumes with thousands to millions of point sources. We solve all dynamical, radiative transfer, thermal, and ionization processes self-consistently on the same mesh, as opposed to a postprocessing approach which coarse-grains the radiative transfer. We do, however, employ a simple subgrid model for star formation which we calibrate to observations. Radiation transport is done in the grey flux-limited diffusion (FLD) approximation, which is solved by implicit time integration split off from the gas energy and ionization equations, which are solved separately. This results in a faster and more robust scheme for cosmological applications compared to the earlier method. The FLD equation is solved using the hypre optimally scalable geometric multigrid solver from LLNL. By treating the ionizing radiation as a grid field as opposed to rays, our method is scalable with respect to the number of ionizing sources, limited only by the parallel scaling properties of the radiation solver. We test the speed and accuracy of our approach on a number of standard verification and validation tests. We show by direct comparison with Enzo's adaptive ray tracing method Moray that the well-known inability of FLD to cast a shadow behind opaque clouds has a minor effect on the evolution of ionized volume and mass fractions in a reionization simulation validation test. We illustrate an application of our method to the problem of inhomogeneous reionization in a 80 Mpc comoving box resolved with $3200^3$ Eulerian grid cells and dark matter particles.

preprint2014arXiv

Heating the IGM by X-rays from Population III Binaries in High Redshift Galaxies

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

preprint2014arXiv

Scaling Relations for Galaxies Prior to Reionization

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

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.

preprint2013arXiv

A supersonic turbulence origin of Larson's laws

We revisit the origin of Larson's scaling laws describing the structure and kinematics of molecular clouds. Our analysis is based on recent observational measurements and data from a suite of six simulations of the interstellar medium, including effects of self-gravity, turbulence, magnetic field, and multiphase thermodynamics. Simulations of isothermal supersonic turbulence reproduce observed slopes in linewidth-size and mass-size relations. Whether or not self-gravity is included, the linewidth-size relation remains the same. The mass-size relation, instead, substantially flattens below the sonic scale, as prestellar cores start to form. Our multiphase models with magnetic field and domain size 200 pc reproduce both scaling and normalization of the first Larson law. The simulations support a turbulent interpretation of Larson's relations. This interpretation implies that: (i) the slopes of linewidth-size and mass-size correlations are determined by the inertial cascade; (ii) none of the three Larson laws is fundamental; (iii) instead, if one is known, the other two follow from scale invariance of the kinetic energy transfer rate. It does not imply that gravity is dynamically unimportant. The self-similarity of structure established by the turbulence breaks in star-forming clouds due to the development of gravitational instability in the vicinity of the sonic scale. The instability leads to the formation of prestellar cores with the characteristic mass set by the sonic scale. The high-end slope of the core mass function predicted by the scaling relations is consistent with the Salpeter power-law index.

preprint2013arXiv

Direct Numerical Simulation of Reionization II: Recombinations, Clumping Factors, and the Photon Budget for Reionization

In this first of several application papers, we investigate the mechanics of reionization from stellar sources in high-z galaxies, the utility of various clumping factors on estimating the recombination time in the IGM, and the photon budget required to achieve reionization. We test the accuracy of the static and time-dependent models of Madau et al. as predictors of reionization completion/maintenance. We simulate a WMAP7 LCDM cosmological model in a 20 Mpc comoving cube with 800^3 uniform fluid cells and dark matter particles. By tuning our star formation to approximately match the observed star formation rate density and luminosity function, we created a fully coupled radiation-hydro realization of H reionization which begins to ionize at z~10 and completes at z~5.8. We find that roughly 2 ionizing photons per H atom are required to convert the neutral IGM to a highly ionized state, which supports the "photon starved" scenario discussed by Bolton & Haehnelt. The events during reionization that lead to this number can generally be described as inside-out, but in reality the narrative depends on the level of ionization of the gas one attributes to as ionized. We find that the formula for the UV photon production rate dN/dt_ion(z) needed to maintain the IGM in an ionized state derived by Madau et al. should not be used to predict the epoch of reionization completion because it ignores history-dependent terms in the global ionization balance which are not ignorable. We find that the time-dependent model for the ionized volume fraction Q_HII is more predictive, but overestimates the redshift of reionization completion by delta_z~1. We propose a revised formulation of the time-dependent model which agrees with our simulation to O(1%). Finally, we use our simulation to estimate a global UV escape fraction due to circumgalactic gas resolved on our mesh to be <f_esc>~0.7.

preprint2013arXiv

Energy cascade and scaling in supersonic isothermal turbulence

Supersonic turbulence plays an important role in a number of extreme astrophysical and terrestrial environments, yet its understanding remains rudimentary. We use data from a three-dimensional simulation of supersonic isothermal turbulence to reconstruct an exact fourth-order relation derived analytically from the Navier-Stokes equations (Galtier and Banerjee, Phys. Rev. Lett., vol. 107, 2011, p. 134501). Our analysis supports a Kolmogorov-like inertial energy cascade in supersonic turbulence previously discussed on a phenomenological level. We show that two compressible analogues of the four-fifths law exist describing fifth- and fourth-order correlations, but only the fourth-order relation remains `universal' in a wide range of Mach numbers from incompressible to highly compressible regimes. A new approximate relation valid in the strongly supersonic regime is derived and verified. We also briefly discuss the origin of bottleneck bumps in simulations of compressible turbulence.

preprint2013arXiv

Enzo: An Adaptive Mesh Refinement Code for Astrophysics

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

preprint2013arXiv

Population III Stars and Remnants in High Redshift Galaxies

Recent simulations of Population III star formation have suggested that some fraction form in binary systems, in addition to having a characteristic mass of tens of solar masses. The deaths of metal-free stars result in the initial chemical enrichment of the universe and the production of the first stellar-mass black holes. Here we present a cosmological adaptive mesh refinement simulation of an overdense region that forms a few 10^9 Msun dark matter halos and over 13,000 Population III stars by redshift 15. We find that most halos do not form Population III stars until they reach Mvir ~ 10^7 Msun because this biased region is quickly enriched from both Population III and galaxies, which also produce high levels of ultraviolet radiation that suppress H2 formation. Nevertheless, Population III stars continue to form, albeit in more massive halos, at a rate of ~ 10^{-4} Msun yr^{-1} Mpc^{-3} at redshift 15. The most massive starless halo has a mass of 7 x 10^7 Msun, which could host massive black hole formation through the direct gaseous collapse scenario. We show that the multiplicity of the Population III remnants grows with halo mass above 10^8 Msun, culminating in 50 remnants located in 10^9 Msun halos on average. This has implications that high mass X-ray binaries and intermediate mass black holes that originate from metal-free stars may be abundant in high-redshift galaxies.

preprint2012arXiv

Comparisons of Cosmological MHD Galaxy Cluster Simulations to Radio Observations

Radio observations of galaxy clusters show that there are $μ$G magnetic fields permeating the intra-cluster medium (ICM), but it is hard to accurately constrain the strength and structure of the magnetic fields without the help of advanced computer simulations. We present qualitative comparisons of synthetic VLA observations of simulated galaxy clusters to radio observations of Faraday Rotation Measure (RM) and radio halos. The cluster formation is modeled using adaptive mesh refinement (AMR) magneto-hydrodynamic (MHD) simulations with the assumption that the initial magnetic fields are injected into the ICM by active galactic nuclei (AGNs) at high redshift. In addition to simulated clusters in Xu et al. (2010, 2011), we present a new simulation with magnetic field injections from multiple AGNs. We find that the cluster with multiple injection sources is magnetized to a similar level as in previous simulations with a single AGN. The RM profiles from simulated clusters, both $|RM|$ and the dispersion of RM ($σ_{RM}$), are consistent at a first-order with the radial distribution from observations. The correlations between the $σ_{RM}$ and X-ray surface brightness from simulations are in a broad agreement with the observations, although there is an indication that the simulated clusters could be slightly over-dense and less magnetized with respect to those in the observed sample. In addition, the simulated radio halos agree with the observed correlations between the radio power versus the cluster X-ray luminosity and between the radio power versus the radio halo size. These studies show that the cluster wide magnetic fields that originate from AGNs and are then amplified by the ICM turbulence (Xu et al. 2010) match observations of magnetic fields in galaxy clusters.

preprint2012arXiv

Cosmological MHD Simulations of Galaxy Cluster Radio Relics: Insights and Warnings for Observations

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

preprint2012arXiv

Flux Correlations in Supersonic Isothermal Turbulence

Using data from a large-scale three-dimensional simulation of supersonic isothermal turbulence, we have tested the validity of an exact flux relation derived analytically from the Navier--Stokes equation by Falkovich, Fouxon and Oz [2010 New relations for correlation functions in Navier--Stokes turbulence. J. Fluid Mech. 644, 465]. That relation, for compressible barotropic fluids, was derived assuming turbulence generated by a large-scale force. However, compressible turbulence in simulations is usually initialized and maintained by a large-scale acceleration, as in gravity-driven astrophysical flows. We present a new approximate flux relation for isothermal turbulence driven by a large-scale acceleration, and find it in reasonable agreement with the simulation results.

preprint2012arXiv

Profiles of Dark Matter Velocity Anisotropy in Simulated Clusters

We report statistical results for dark matter (DM) velocity anisotropy, β, from a sample of some 6000 cluster-size halos (at redshift zero) identified in a ΛCDM hydrodynamical adaptive mesh refinement simulation performed with the Enzo code. These include profiles of β in clusters with different masses, relaxation states, and at several redshifts, modeled both as spherical and triaxial DM configurations. Specifically, although we find a large scatter in the DM velocity anisotropy profiles of different halos (across elliptical shells extending to at least ~$1.5 r_{vir}$), universal patterns are found when these are averaged over halo mass, redshift, and relaxation stage. These are characterized by a very small velocity anisotropy at the halo center, increasing outward to about 0.27 and leveling off at about $0.2 r_{vir}$. Indirect measurements of the DM velocity anisotropy fall on the upper end of the theoretically expected range. Though measured indirectly, the estimations are derived by using two different surrogate measurements - X-ray and galaxy dynamics. Current estimates of the DM velocity anisotropy are based on very small cluster sample. Increasing this sample will allow testing theoretical predictions, including the speculation that the decay of DM particles results in a large velocity boost. We also find, in accord with previous works, that halos are triaxial and likely to be more prolate when unrelaxed, whereas relaxed halos are more likely to be oblate. Our analysis does not indicate that there is significant correlation (found in some previous studies) between the radial density slope, γ, and β at large radii, $0.3 r_{vir} < r < r_{vir}$.

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.

preprint2012arXiv

The Two States of Star Forming Clouds

We examine the effects of self-gravity and magnetic fields on supersonic turbulence in isothermal molecular clouds with high resolution simulations and adaptive mesh refinement. These simulations use large root grids (512^3) to capture turbulence and four levels of refinement to capture high density, for an effective resolution of 8,196^3. Three Mach 9 simulations are performed, two super-Alfvénic and one trans-Alfvénic. We find that gravity splits the clouds into two populations, one low density turbulent state and one high density collapsing state. The low density state exhibits properties similar to non-self-gravitating in this regime, and we examine the effects of varied magnetic field strength on statistical properties: the density probability distribution function is approximately lognormal; velocity power spectral slopes decrease with field strength; alignment between velocity and magnetic field increases with field; the magnetic field probability distribution can be fit to a stretched exponential. The high density state is characterized by self-similar spheres; the density PDF is a power-law; collapse rate decreases with increasing mean field; density power spectra have positive slopes, P(ρ,k) \propto k; thermal-to-magnetic pressure ratios are unity for all simulations; dynamic-to-magnetic pressure ratios are larger than unity for all simulations; magnetic field distribution is a power-law. The high Alfvén Mach numbers in collapsing regions explain recent observations of magnetic influence decreasing with density. We also find that the high density state is found in filaments formed by converging flows, consistent with recent Herschel observations. Possible modifications to existing star formation theories are explored.

preprint2012arXiv

What Shapes the Structure of Molecular Clouds: Turbulence or Gravity?

We revisit the origin of Larson's scaling relations, which describe the structure and kinematics of molecular clouds, based on recent observations and large-scale simulations of supersonic turbulence. Using dimensional analysis, we first show that both linewidth-size and mass-size correlations observed on scales 0.1-50 pc can be explained by a simple conceptual theory of compressible turbulence without resorting to the often assumed virial equilibrium or detailed energy balance condition. The scaling laws can be consistently interpreted as a signature of supersonic turbulence with no need to invoke gravity. We then show how self-similarity of structure established by the turbulence breaks in star-forming clouds through development of gravitational instabilities in the vicinity of the sonic scale, l_s~0.1 pc, leading to the formation of prestellar cores.

preprint2011arXiv

Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence

We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss convergence of various characteristics for the turbulence decay test and impacts of various components of numerical schemes on the accuracy of solutions. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the magnetic field using the constrained transport method, and using little to no explicit artificial viscosity. Codes which fall short in one or more of these areas are still useful, but they must compensate higher numerical dissipation with higher numerical resolution. This paper is the largest, most comprehensive MHD code comparison on an application-like test problem to date. We hope this work will help developers improve their numerical algorithms while helping users to make informed choices in picking optimal applications for their specific astrophysical problems.

preprint2011arXiv

Dispersal of Galactic Magnetic Fields into Intracluster Space

Little is known about the origin and basic properties of magnetic fields in clusters of galaxies. High conductivity in magnetized interstellar plasma suggests that galactic magnetic fields are (at least partly) ejected into intracluster (IC) space by the same processes that enrich IC gas with metals. We explore the dispersal of galactic fields by hydrodynamical simulations with our new {\em Enzo-Galcon} code, which is capable of tracking a large number galaxies during cluster assembly, and modeling the processes that disperse their interstellar media. Doing so we are able to describe the evolution of the mean strength of the field and its profile across the cluster. With the known density profile of dispersed gas and an estimated range of coherence scales, we predict the spatial distribution of Faraday rotation measure and find it to be consistent with observational data.

preprint2011arXiv

Evolution and Distribution of Magnetic Fields from AGNs in Galaxy Clusters II. The Effects of Cluster Size and Dynamical State

Theory and simulations suggest that magnetic fields from radio jets and lobes powered by their central super massive black holes can be an important source of magnetic fields in the galaxy clusters. This is paper II in a series of studies where we present self-consistent high-resolution adaptive mesh refinement cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and evolution of magnetic fields ejected by an active galactic nucleus (AGN). We studied 12 different galaxy clusters with virial masses ranging from 1 $\times$ 10$^{14}$ to 2 $\times$ 10$^{15}$ M$_{\odot}$. In this work we examine the effects of the mass and merger history on the final magnetic properties. We find that the evolution of magnetic fields is qualitatively similar to those of previous studies. In most clusters, the injected magnetic fields can be transported throughout the cluster and be further amplified by the intra-cluster medium (ICM) turbulence during the cluster formation process with hierarchical mergers, while the amplification history and the magnetic field distribution depend on the cluster formation and magnetism history. This can be very different for different clusters. The total magnetic energies in these clusters are between 4 $\times$ 10$^{57}$ and $10^{61}$ erg, which is mainly decided by the cluster mass, scaling approximately with the square of the total mass. Dynamically older relaxed clusters usually have more magnetic fields in their ICM. The dynamically very young clusters may be magnetized weakly since there is not enough time for magnetic fields to be amplified.

preprint2011arXiv

On the Density Distribution in Star-forming Interstellar Clouds

We use deep adaptive mesh refinement simulations of isothermal self-gravitating supersonic turbulence to study the imprints of gravity on the mass density distribution in molecular clouds. The simulations show that the density distribution in self-gravitating clouds develops an extended power-law tail at high densities on top of the usual lognormal. We associate the origin of the tail with self-similar collapse solutions and predict the power index values in the range from -7/4 to -3/2 that agree with both simulations and observations of star-forming molecular clouds.

preprint2011arXiv

The Birth of a Galaxy: Primordial Metal Enrichment and Stellar Populations

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

preprint2011arXiv

Turbulent Clustering of Protoplanetary Dust and Planetesimal Formation

We study clustering of inertial particles in turbulent flows and discuss its applications to dust particles in protoplanetary disks. Using numerical simulations, we compute the radial distribution function (RDF), which measures the probability of finding particle pairs at given distances, and the probability density function of the particle concentration. The clustering statistics depend on the Stokes number, $St$, defined as the ratio of the particle friction timescale, $τ_{\rm p} $, to the Kolmogorov timescale in the flow. In the dissipation range, the clustering intensity strongly peaks at $St \simeq 1$, and the RDF for $St \sim 1$ shows a fast power-law increase toward small scales, suggesting that turbulent clustering may considerably enhance the particle collision rate. Clustering at inertial-range scales is of particular interest to the problem of planetesimal formation. At these scales, the strongest clustering is from particles with $τ_{\rm p}$ in the inertial range. Clustering of these particles occurs primarily around a scale where the eddy turnover time is $\simτ_{\rm p}$. Particles of different sizes tend to cluster at different locations, leading to flat RDFs between different particles at small scales. In the presence of multiple particle sizes, the overall clustering strength decreases as the particle size distribution broadens. We discuss particle clustering in recent models for planetesimal formation. We point out that, in the model based on turbulent clustering of chondrule-size particles, the probability of finding strong clusters that can seed planetesimals may have been significantly overestimated. We discuss various clustering mechanisms in simulations of planetesimal formation by gravitational collapse of dense clumps of meter-size particles, in particular the contribution from turbulent clustering due to the limited numerical resolution.

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.

preprint2010arXiv

Evolution and Distribution of Magnetic Fields from AGNs in Galaxy Clusters. I. The Effect of Injection Energy and Redshift

We present a series of cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and evolution of magnetic fields ejected by an Active Galactic Nucleus (AGN). Specifically, we investigate the influence of both the epoch of AGN (z $\sim$ 3-0.5) and the AGN energy ($\sim$ 3 $\times$ 10$^{57}$ - 2 $\times$ 10$^{60}$ ergs)on the final magnetic field distribution in a relatively massive cluster (M$_{vir}$ $\sim$10$^{15}$ M$_\odot$). We find that as long as the AGN magnetic fields are ejected before the major mergers in the cluster formation history, magnetic fields can be transported throughout the cluster and can be further amplified by the intra-cluster medium (ICM) turbulence cause by hierarchical mergers during the cluster formation process. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs, with micro Gauss fields distributed over $\sim$ Mpc scale. The amplification of the total magnetic energy by the ICM turbulence can be significant, up to $\sim$1000 times in some cases. Therefore even weak magnetic fields from AGNs can be used to magnetize the cluster to the observed level. The final magnetic energy in the ICM is determined by the ICM turbulent energy, with a weak dependence on the AGN injection energy. We discuss the properties of magnetic fields throughout the cluster and the synthetic Faraday rotation measure maps they produce. We also show that high spatial resolution over most of the magnetic regions of the cluster is very important to capture the small scale dynamo process and maintain the magnetic field structure in our simulations.

preprint2010arXiv

High-Entropy Polar Regions Around the First Protostars

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

preprint2010arXiv

Historical perspective on astrophysical MHD simulations

This contribution contains the introductory remarks that I presented at IAU Symposium 270 on ``Computational Star Formation" held in Barcelona, Spain, May 31 -- June 4, 2010. I discuss the historical development of numerical MHD methods in astrophysics from a personal perspective. The recent advent of robust, higher order-accurate MHD algorithms and adaptive mesh refinement numerical simulations promises to greatly improve our understanding of the role of magnetic fields in star formation.

preprint2010arXiv

Hydrodynamical Simulations of Galaxy Clusters with Galcons

We present our recently developed {\em galcon} approach to hydrodynamical cosmological simulations of galaxy clusters - a subgrid model added to the {\em Enzo} adaptive mesh refinement code - which is capable of tracking galaxies within the cluster potential and following the feedback of their main baryonic processes. Galcons are physically extended galactic constructs within which baryonic processes are modeled analytically. By identifying galaxy halos and initializing galcons at high redshift ($z \sim 3$, well before most clusters virialize), we are able to follow the evolution of star formation, galactic winds, and ram-pressure stripping of interstellar media, along with their associated mass, metals and energy feedback into intracluster (IC) gas, which are deposited through a well-resolved spherical interface layer. Our approach is fully described and all results from initial simulations with the enhanced {\em Enzo-Galcon} code are presented. With a galactic star formation rate derived from the observed cosmic star formation density, our galcon simulation better reproduces the observed properties of IC gas, including the density, temperature, metallicity, and entropy profiles. By following the impact of a large number of galaxies on IC gas we explicitly demonstrate the advantages of this approach in producing a lower stellar fraction, a larger gas core radius, an isothermal temperature profile in the central cluster region, and a flatter metallicity gradient than in a standard simulation.

preprint2010arXiv

Interstellar Turbulence and Star Formation

We provide a brief overview of recent advances and outstanding issues in simulations of interstellar turbulence, including isothermal models for interior structure of molecular clouds and larger-scale multiphase models designed to simulate the formation of molecular clouds. We show how self-organization in highly compressible magnetized turbulence in the multiphase ISM can be exploited in simple numerical models to generate realistic initial conditions for star formation.

preprint2010arXiv

Mass and Magnetic distributions in Self Gravitating Super Alfvenic Turbulence with AMR

In this work, we present the mass and magnetic distributions found in a recent Adaptive Mesh Refinement (AMR) MHD simulation of supersonic, \sa, self gravitating turbulence. Powerlaw tails are found in both volume density and magnetic field probability density functions, with $P(ρ) \propto ρ^{-1.67}$ and $P(B)\propto B^{-2.74}$. A power law is also found between magnetic field strength and density, with $B\propto ρ^{0.48}$, throughout the collapsing gas. The mass distribution of gravitationally bound cores is shown to be in excellent agreement with recent observation of prestellar cores. The mass to flux distribution of cores is also found to be in excellent agreement with recent Zeeman splitting measurements.

preprint2010arXiv

Parallel HOP: A Scalable Halo Finder for Massive Cosmological Data Sets

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

preprint2010arXiv

Pop III Stellar Masses and IMF

We provide a status report on our current understanding of the mass scales for Pop III.1 and Pop III.2 stars. Since the last review (Norman 2008), substantial progress has been made both numerically and analytically on the late stages of protostellar cloud core collapse, protostar formation and accretion, and stellar evolution taking into account cloud core properties and radiative feedback effects. Based on this, there are growing indications that primordial stars forming from purely cosmological initial conditions (Pop III.1) were substantially more massive than stars forming in preionized gas (Pop III.2) where HD cooling is important. Different stellar endpoints are predicted for these two types of Pop III stars with different chemical enrichment signatures: the former die as pair instability supernovae or intermediate mass black holes, whereas the latter die as iron core-collapse supernovae, leaving behind neutron star and stellar black hole remnants. We review recent simulations which show evidence for binary fragmentation at high densities, and comment on the significance of these results. We then summarize an attempt to directly calculate the Pop III.1 IMF taking into account the latest numerical and analytical models. We conclude with suggestions for the kind of simulations needed next to continue improving our understanding of Pop III star formation, which is a necessary input to understanding high redshift galaxy formation.

preprint2010arXiv

Radiation Hydrodynamical Instabilities in Cosmological and Galactic Ionization Fronts

Ionization fronts, the sharp radiation fronts behind which H/He ionizing photons from massive stars and galaxies propagate through space, were ubiquitous in the universe from its earliest times. The cosmic dark ages ended with the formation of the first primeval stars and galaxies a few hundred Myr after the Big Bang. Numerical simulations suggest that stars in this era were very massive, 25 - 500 solar masses, with H II regions of up to 30,000 light-years in diameter. We present three-dimensional radiation hydrodynamical calculations that reveal that the I-fronts of the first stars and galaxies were prone to violent instabilities, enhancing the escape of UV photons into the early intergalactic medium (IGM) and forming clumpy media in which supernovae later exploded. The enrichment of such clumps with metals by the first supernovae may have led to the prompt formation of a second generation of low-mass stars, profoundly transforming the nature of the first protogalaxies. Cosmological radiation hydrodynamics is unique because ionizing photons coupled strongly to both gas flows and primordial chemistry at early epochs, introducing a hierarchy of disparate characteristic timescales whose relative magnitudes can vary greatly throughout a given calculation. We describe the adaptive multistep integration scheme we have developed for the self-consistent transport of both cosmological and galactic ionization fronts.

preprint2010arXiv

Simulating Galaxy Clusters

This pedagogical review of galaxy cluster simulations is based on three lectures given at the 2008 Enrico Fermi Summer School entitled "Astrophysics of Galaxy Clusters". It covers the standard cosmological framework, growth of perturbations in the linear regime, analytic models for nonlinear perturbation growth, statistics of galaxy cluster populations, virial scaling relations, overview of numerical methods, simulating gas in galaxy clusters, basic results on adiabatic clusters (Santa Barbara cluster comparison project), effect of additional physics, recent progress in galaxy clustering modeling (Galcons, turbulence, AGN jets, cluster-wide B-fields), simulating statistical samples and lightcones, and simulated SZE surveys.

preprint2010arXiv

The Effect of Projection on Derived Mass-Size and Linewidth-Size Relationships

Power law mass-size and linewidth-size correlations, two of "Larson's laws," are often studied to assess the dynamical state of clumps within molecular clouds. Using the result of a hydrodynamic simulation of a molecular cloud, we investigate how geometric projection may affect the derived Larson relationships. We find that large scale structures in the column density map have similar masses and sizes to those in the 3D simulation (PPP). Smaller scale clumps in the column density map are measured to be more massive than the PPP clumps, due to the projection of all emitting gas along lines of sight. Further, due to projection effects, structures in a synthetic spectral observation (PPV) may not necessarily correlate with physical structures in the simulation. In considering the turbulent velocities only, the linewidth-size relationship in the PPV cube is appreciably different from that measured from the simulation. Including thermal pressure in the simulated linewidths imposes a minimum linewidth, which results in a better agreement in the slopes of the linewidth-size relationships, though there are still discrepancies in the offsets, as well as considerable scatter. Employing commonly used assumptions in a virial analysis, we find similarities in the computed virial parameters of the structures in the PPV and PPP cubes. However, due to the discrepancies in the linewidth- and mass- size relationships in the PPP and PPV cubes, we caution that applying a virial analysis to observed clouds may be misleading due to geometric projection effects. We speculate that consideration of physical processes beyond kinetic and gravitational pressure would be required for accurately assessing whether complex clouds, such as those with highly filamentary structure, are bound.

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.

preprint2009arXiv

Baryon Acoustic Oscillations in the Lyman Alpha Forest

We use hydrodynamic cosmological simulations in a (600 Mpc)^3 volume to study the observability of baryon acoustic oscillations (BAO) in the intergalactic medium as probed by Lyman alpha forest (LAF) absorption. The large scale separation between the wavelength of the BAO mode (~150 Mpc) and the size of LAF absorbers (~100 kpc) makes this a numerically challenging problem. We report on several 2048^3 simulations of the LAF using the ENZO code. We adopt WMAP5 concordance cosmological parameters and power spectrum including BAO perturbations. 5000 synthetic HI absorption line spectra are generated randomly piercing the box face. We calculate the cross-correlation function between widely separated pairs. We detect the BAO signal at z=3 where theory predicts to moderate statistical significance.

preprint2009arXiv

Cosmological AMR MHD with Enzo

In this work, we present MHDEnzo, the extension of the cosmological code Enzo to include the effects magnetic fields through the ideal MHD approximation. We use a higher order Godunov Riemann solver for the computation of interface fluxes. We use two constrained transport methods to compute the electric field from those interface fluxes, which simultaneously advances the induction equation and maintains the divergence of the magnetic field. A third order divergence free reconstruction technique is used to interpolate the magnetic fields in the block structured AMR framework already extant in Enzo. This reconstruction also preserves the divergence of the magnetic field to machine precision. We use operator splitting to include gravity and cosmological expansion. We then present a series of cosmological and non cosmological tests problems to demonstrate the quality of solution resulting from this combination of solvers.

preprint2009arXiv

Cosmological Radiation Hydrodynamics with ENZO

We describe an extension of the cosmological hydrodynamics code ENZO to include the self-consistent transport of ionizing radiation modeled in the flux-limited diffusion approximation. A novel feature of our algorithm is a coupled implicit solution of radiation transport, ionization kinetics, and gas photoheating, making the timestepping for this portion of the calculation resolution independent. The implicit system is coupled to the explicit cosmological hydrodynamics through operator splitting and solved with scalable multigrid methods. We summarize the numerical method, present a verification test on cosmological Stromgren spheres, and then apply it to the problem of cosmological hydrogen reionization.

preprint2009arXiv

Piecewise Parabolic Method on a Local Stencil for Magnetized Supersonic Turbulence Simulation

Stable, accurate, divergence-free simulation of magnetized supersonic turbulence is a severe test of numerical MHD schemes and has been surprisingly difficult to achieve due to the range of flow conditions present. Here we present a new, higher order-accurate, low dissipation numerical method which requires no additional dissipation or local "fixes" for stable execution. We describe PPML, a local stencil variant of the popular PPM algorithm for solving the equations of compressible ideal magnetohydrodynamics. The principal difference between PPML and PPM is that cell interface states are evolved rather that reconstructed at every timestep, resulting in a compact stencil. Interface states are evolved using Riemann invariants containing all transverse derivative information. The conservation laws are updated in an unsplit fashion, making the scheme fully multidimensional. Divergence-free evolution of the magnetic field is maintained using the higher order-accurate constrained transport technique of Gardiner and Stone. The accuracy and stability of the scheme is documented against a bank of standard test problems drawn from the literature. The method is applied to numerical simulation of supersonic MHD turbulence, which is important for many problems in astrophysics, including star formation in dark molecular clouds. PPML accurately reproduces in three-dimensions a transition to turbulence in highly compressible isothermal gas in a molecular cloud model. The low dissipation and wide spectral bandwidth of this method make it an ideal candidate for direct turbulence simulations.

preprint2009arXiv

Self-Consistent Solution of Cosmological Radiation-Hydrodynamics and Chemical Ionization

We consider a PDE system comprising compressible hydrodynamics, flux-limited diffusion radiation transport and chemical ionization kinetics in a cosmologically-expanding universe. Under an operator-split framework, the cosmological hydrodynamics equations are solved through the Piecewise Parabolic Method, as implemented in the Enzo community hydrodynamics code. The remainder of the model, including radiation transport, chemical ionization kinetics, and gas energy feedback, form a stiff coupled PDE system, which we solve using a fully-implicit inexact Newton approach, and which forms the crux of this paper. The inner linear Newton systems are solved using a Schur complement formulation, and employ a multigrid-preconditioned conjugate gradient solver for the inner Schur systems. We describe this approach and provide results on a suite of test problems, demonstrating its accuracy, robustness, and scalability to very large problems.

preprint2009arXiv

Self-organization in Turbulent Molecular Clouds: Compressional versus Solenoidal Modes

We use three-dimensional numerical simulations to study self-organization in supersonic turbulence in molecular clouds. Our numerical experiments describe decaying and driven turbulent flows with an isothermal equation of state, sonic Mach numbers from 2 to 10, and various degrees of magnetization. We focus on properties of the velocity field and, specifically, on the level of its potential (dilatational) component as a function of turbulent Mach number, magnetic field strength, and scale. We show how extreme choices of either purely solenoidal or purely potential forcing can reduce the extent of the inertial range in the context of periodic box models for molecular cloud turbulence. We suggest an optimized forcing to maximize the effective Reynolds number in numerical models.

preprint2009arXiv

The Power Spectrum of Turbulence in NGC 1333: Outflows or Large-Scale Driving?

Is the turbulence in cluster-forming regions internally driven by stellar outflows or the consequence of a large-scale turbulent cascade? We address this question by studying the turbulent energy spectrum in NGC 1333. Using synthetic 13CO maps computed with a snapshot of a supersonic turbulence simulation, we show that the VCS method of Lazarian and Pogosyan provides an accurate estimate of the turbulent energy spectrum. We then apply this method to the 13CO map of NGC 1333 from the COMPLETE database. We find the turbulent energy spectrum is a power law, E(k) k^-beta, in the range of scales 0.06 pc < ell < 1.5 pc, with slope beta=1.85\pm 0.04. The estimated energy injection scale of stellar outflows in NGC 1333 is ell_inj 0.3 pc, well resolved by the observations. There is no evidence of the flattening of the energy spectrum above the scale ell_inj predicted by outflow-driven simulations and analytical models. The power spectrum of integrated intensity is also a nearly perfect power law in the range of scales 0.16 pc < ell < 7.9 pc, with no feature above ell_inj. We conclude that the observed turbulence in NGC 1333 does not appear to be driven primarily by stellar outflows.

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.

preprint2009arXiv

Turbulence and Dynamo in Galaxy Cluster Medium: Implications on the Origin of Cluster Magnetic Fields

We present self-consistent cosmological magnetohydrodynamic (MHD) simulations that simultaneously follow the formation of a galaxy cluster and the magnetic field ejection by an active galactic nucleus (AGN). We find that the magnetic fields ejected by the AGNs, though initially distributed in relatively small volumes, can be transported throughout the cluster and be further amplified by the intra-cluster medium (ICM) turbulence during the cluster formation process. The ICM turbulence is shown to be generated and sustained by the frequent mergers of smaller halos. Furthermore, a cluster-wide dynamo process is shown to exist in the ICM and amplify the magnetic field energy and flux. The total magnetic energy in the cluster can reach $\sim$ $10^{61}$ ergs while micro Gauss ($μ$G) fields can distribute over $\sim$ Mpc scales throughout the whole cluster. This finding shows that magnetic fields from AGNs, being further amplified by the ICM turbulence through small-scale dynamo processes, can be the origin of cluster-wide magnetic fields.

preprint2008arXiv

The Lyman-alpha forest at redshifts 0.1 -- 1.6: good agreement between a large hydrodynamic simulation and HST spectra

We give a comprehensive statistical description of the Lyman-alpha absorption from the intergalactic medium in a hydrodynamic simulation at redshifts 0.1-1.6, the range of redshifts covered by HST spectra of QSOs. We use the ENZO code to make a 76 comoving Mpc cube simulation using 75 kpc cells, for a Hubble constant of 71 km/s/Mpc. The best prior work, by \citet{dave99},used an SPH simulation in a 15.6 Mpc box with an effective resolution of 245 kpc and slightly different cosmological parameters. At redshifts z=2 this simulation is different from data. \citet{tytler07b} found that the simulated spectra at z=2 have too little power on large scales, Lyman-alpha lines are too wide, there is a lack high column density lines, and there is a lack of pixels with low flux. Here we present statistics at z<1.6, including the flux distribution, the mean flux, the effective opacity, and the power and correlation of the flux. We also give statistics of the lyman alpha lines including the line width distribution, the column density distribution, the number of lines per unit equivalent width and redshift, and the correlation between the line width and column density. We find that the mean amount of absorption in the simulated spectra changes smoothly with redshift with DA(z)=0.01(1+z)^{2.25}. Both the trend and absolute values are close to measurements of HST spectra by \citet{kirkman07a}. The column density and line width distributions are also close to those measured from HST spectra by \citet{janknecht06a}, except for the mode of the line width distribution which is smaller in the HST spectra. Although some differences that we saw at z=2 are too subtle to be seen in existing HST spectra, overall, the simulation gives an good description of HST spectra at 0.1<z<1.6.

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)

preprint2007arXiv

Resonance Related Spiral Substructure in a Galactic Gaseous Disk

We use high resolution 2D hydrodynamic simulations to study the formation of spiral substructure in the gaseous disk of a galaxy. The obtained gaseous response is driven by a self-consistent non-axisymmetric potential obtained from an imposed spiral mass distribution. We highlight the importance of ultraharmonic resonances in generating these features. The temporal evolution of the system is followed with the parallel ZEUS-MP code, and we follow the steepening of perturbations induced by the spiral potential until large-scale shocks emerge. These shocks exhibit bifurcations that protrude from the gaseous arms and continue to steepen until new shocks are formed. When the contribution from the spiral potential relative to the axisymmetric background is increased from our default value, spurs protrude from the main arms after several revolutions of the gaseous disk. Such spurs overlap on top of the aforementioned shocks. These results support the hypothesis that a complicated gaseous response can coexist with an orderly spiral potential term, in the sense that the underlying background potential can be smooth yet drive a gaseous response that is far more spatially complex.

preprint2007arXiv

The Effect of Large-Scale Power on Simulated Spectra of the Lya forest

We study the effects of box size on ENZO simulations of the intergalactic medium (IGM) at z = 2. We follow statistics of the cold dark matter (CDM) and the Lya absorption. We find that the larger boxes have fewer pixels with significant absorption (flux < 0.96) and more pixels in longer stretches with little or no absorption, and they have wider Lya lines. We trace these effect back to the additional power in larger boxes from longer wavelength modes. The IGM in our larger boxes is hotter, from increased pressure heating due to faster hydrodynamical infall. When we increase the photoheating in smaller boxes to compensate, their Lya statistics change to mimic those of a box of twice the size. Statistics converge towards their value in the largest (76.8 Mpc) box, except for the most common value of the CDM density which continues to rise. When we compare to errors with data, we find that our 76.8 Mpc box is larger than we need for the mean flux, barely large enough for the column density distribution and the power spectrum of the flux, and too small for the line widths. This box with 75 kpc cells has approximately the same mean flux as QSO spectra, but the Lya lines are too wide by 2.6 km/s, there are too few lines with log H I column densities > 10^17 cm^-2, and the power of the flux is too low by 20 - 50%, from small to large scales. Four times smaller cell size does not resolve these differences, nor do simple changes to the ultraviolet background that drives the H and He II ionization. It is hard to see how simulations using popular cosmological and astrophysical parameters can match Lyman-alpha forest data at z=2.

preprint2007arXiv

Two regimes of Turbulent Fragmentation and the stellar IMF from Primordial to Present Day Star Formation

The Padoan and Nordlund model of the stellar initial mass function (IMF) is derived from low order statistics of supersonic turbulence, neglecting gravity (e.g. gravitational fragmentation, accretion and merging). In this work the predictions of that model are tested using the largest numerical experiments of supersonic hydrodynamic (HD) and magneto-hydrodynamic (MHD) turbulence to date (~1000^3 computational zones) and three different codes (Enzo, Zeus and the Stagger Code). The model predicts a power law distribution for large masses, related to the turbulence energy power spectrum slope, and the shock jump conditions. This power law mass distribution is confirmed by the numerical experiments. The model also predicts a sharp difference between the HD and MHD regimes, which is recovered in the experiments as well, implying that the magnetic field, even below energy equipartition on the large scale, is a crucial component of the process of turbulent fragmentation. These results suggest that the stellar IMF of primordial stars may differ from that in later epochs of star formation, due to differences in both gas temperature and magnetic field strength. In particular, we find that the IMF of primordial stars born in turbulent clouds may be narrowly peaked around a mass of order 10 solar masses, as long as the column density of such clouds is not much in excess of 10^22 cm^-2.

preprint2006arXiv

AMR Simulations of the Cosmological Light Cone: SZE Surveys of the Synthetic Universe

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

preprint2006arXiv

Beyond the Cool Core: The Formation of Cool Core Galaxy Clusters

Why do some clusters have cool cores while others do not? In this paper, cosmological simulations, including radiative cooling and heating, are used to examine the formation and evolution of cool core (CC) and non-cool core (NCC) clusters. Numerical CC clusters at z=0 accreted mass more slowly over time and grew enhanced cool cores via hierarchical mergers; when late major mergers occurred, the CCs survived the collisions. By contrast, NCC clusters of similar mass experienced major mergers early in their evolution that destroyed embryonic cool cores and produced conditions that prevent CC re-formation. We discuss observational consequences.

preprint2005arXiv

Adaptive Mesh Refinement for Supersonic Molecular Cloud Turbulence

We performed a series of three-dimensional numerical simulations of supersonic homogeneous Euler turbulence with adaptive mesh refinement (AMR) and effective grid resolution up to 1024^3 zones. Our experiments describe non-magnetized driven supersonic turbulent flows with an isothermal equation of state. Mesh refinement on shocks and shear is implemented to cover dynamically important structures with the highest resolution subgrids and calibrated to match the turbulence statistics obtained from the equivalent uniform grid simulations. We found that at a level of resolution slightly below 512^3, when a sufficient integral/dissipation scale separation is first achieved, the fraction of the box volume covered by the AMR subgrids first becomes smaller than unity. At the higher AMR levels subgrids start covering smaller and smaller fractions of the whole volume, which scale with the Reynolds number as Re^{-1/4}. We demonstrate the consistency of this scaling with a hypothesis that the most dynamically important structures in intermittent supersonic turbulence are strong shocks with a fractal dimension of two. We show that turbulence statistics derived from AMR simulations and simulations performed on uniform grids agree surprisingly well, even though only a fraction of the volume is covered by AMR subgrids. Based on these results, we discuss the signature of dissipative structures in the statistical properties of supersonic turbulence and their role in overall flow dynamics.

preprint2001arXiv

The Formation of the First Star in the Universe

We describe results from a fully self-consistent three dimensional hydrodynamical simulation of the formation of one of the first stars in the Universe. Dark matter dominated pre-galactic objects form because of gravitational instability from small initidal density perturbations. As they assemble via hierarchical merging, primordial gas cools through ro-vibrational lines of hydrogen molecules and sinks to the center of the dark matter potential well. The high redshift analog of a molecular cloud is formed. When the dense, central parts of the cold gas cloud become self-gravitating, a dense core of approximately 100 solar mass undergoes rapid contraction. At densities n>10^9 cm^-3 a one solar mass proto-stellar core becomes fully molecular due to three-body H_2 formation. Contrary to analytical expectations this process does not lead to renewed fragmentation and only one star is formed. The calculation is stopped when optical depth effects become important, leaving the final mass of the fully formed star somewhat uncertain. At this stage the protostar is acreting material very rapidly (~0.01 solar masses per year). Radiative feedback from the star will not only halt its growth but also inhibit the formation of other stars in the same pre-galactic object (at least until the first star ends its life, presumably as a supernova). We conclude that at most one massive (M >> 1 solar mass) metal free star forms per pre-galactic halo, consistent with recent abundance measurements of metal poor galactic halo stars.

preprint2000arXiv

Baryons in the Warm-Hot Intergalactic Medium

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

preprint1998arXiv

Cosmological Adaptive Mesh Refinement

We describe a grid-based numerical method for 3D hydrodynamic cosmological simulations which is adaptive in space and time and combines the best features of higher order--accurate Godunov schemes for Eulerian hydrodynamics with adaptive particle--mesh methods for collisionless particles. The basis for our method is the structured adaptive mesh refinement (AMR) algorithm of Berger & Collela (1989), which we have extended to cosmological hydro + N-body simulations. The resulting multiscale hybrid method is a powerful alternative to particle-based methods in current use. The choices we have made in constructing this algorithm are discussed, and its performance on the Zeldovich pancake test problem is given. We present a sample application of our method to the problem of first structure formation. We have achieved a spatial dynamic range $L_{box}/Δx > 250,000$ in a 3D multispecies gas + dark matter calculation, which is sufficient to resolve the formation of primordial protostellar cloud cores starting from linear matter fluctuations in an expanding FRW universe.

Michael L. Norman

What is connected

Connect this record

See the researcher in context

Building this map preview

64 published item(s)

Connecting Primordial Star Forming Regions and Second Generation Star Formation in the Phoenix Simulations

Massive Star Formation in Metal-Enriched Haloes at High Redshift

The Emergence of the First Star-free Atomic Cooling Haloes in the Universe

The Formation of Very Massive Stars in Early Galaxies and Implications for Intermediate Mass Black Holes

Galaxy Properties and UV Escape Fractions During Epoch of Reionization: Results from the Renaissance Simulations

Late Pop III Star Formation During the Epoch of Reionization: Results from the Renaissance Simulations

Length Scales and Turbulent Properties of Magnetic Fields in Simulated Galaxy Clusters

Probing The Ultraviolet Luminosity Function of the Earliest Galaxies with the Renaissance Simulations

Spatially Extended 21 cm Signal from Strongly Clustered UV and X-Ray Sources in the Early Universe

Fully-Coupled Simulation of Cosmic Reionization. I: Numerical Methods and Tests

Heating the IGM by X-rays from Population III Binaries in High Redshift Galaxies

Scaling Relations for Galaxies Prior to Reionization

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

A supersonic turbulence origin of Larson's laws

Direct Numerical Simulation of Reionization II: Recombinations, Clumping Factors, and the Photon Budget for Reionization

Energy cascade and scaling in supersonic isothermal turbulence

Enzo: An Adaptive Mesh Refinement Code for Astrophysics

Population III Stars and Remnants in High Redshift Galaxies

Comparisons of Cosmological MHD Galaxy Cluster Simulations to Radio Observations

Cosmological MHD Simulations of Galaxy Cluster Radio Relics: Insights and Warnings for Observations

Flux Correlations in Supersonic Isothermal Turbulence

Profiles of Dark Matter Velocity Anisotropy in Simulated Clusters

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

The Two States of Star Forming Clouds

What Shapes the Structure of Molecular Clouds: Turbulence or Gravity?

Comparing Numerical Methods for Isothermal Magnetized Supersonic Turbulence

Dispersal of Galactic Magnetic Fields into Intracluster Space

Evolution and Distribution of Magnetic Fields from AGNs in Galaxy Clusters II. The Effects of Cluster Size and Dynamical State

On the Density Distribution in Star-forming Interstellar Clouds

The Birth of a Galaxy: Primordial Metal Enrichment and Stellar Populations

Turbulent Clustering of Protoplanetary Dust and Planetesimal Formation

A Multi-Code Analysis Toolkit for Astrophysical Simulation Data

Evolution and Distribution of Magnetic Fields from AGNs in Galaxy Clusters. I. The Effect of Injection Energy and Redshift

High-Entropy Polar Regions Around the First Protostars

Historical perspective on astrophysical MHD simulations

Hydrodynamical Simulations of Galaxy Clusters with Galcons

Interstellar Turbulence and Star Formation

Mass and Magnetic distributions in Self Gravitating Super Alfvenic Turbulence with AMR

Parallel HOP: A Scalable Halo Finder for Massive Cosmological Data Sets

Pop III Stellar Masses and IMF

Radiation Hydrodynamical Instabilities in Cosmological and Galactic Ionization Fronts

Simulating Galaxy Clusters

The Effect of Projection on Derived Mass-Size and Linewidth-Size Relationships

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

Baryon Acoustic Oscillations in the Lyman Alpha Forest

Cosmological AMR MHD with Enzo

Cosmological Radiation Hydrodynamics with ENZO

Piecewise Parabolic Method on a Local Stencil for Magnetized Supersonic Turbulence Simulation

Self-Consistent Solution of Cosmological Radiation-Hydrodynamics and Chemical Ionization

Self-organization in Turbulent Molecular Clouds: Compressional versus Solenoidal Modes

The Power Spectrum of Turbulence in NGC 1333: Outflows or Large-Scale Driving?

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

Turbulence and Dynamo in Galaxy Cluster Medium: Implications on the Origin of Cluster Magnetic Fields

The Lyman-alpha forest at redshifts 0.1 -- 1.6: good agreement between a large hydrodynamic simulation and HST spectra

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

Resonance Related Spiral Substructure in a Galactic Gaseous Disk

The Effect of Large-Scale Power on Simulated Spectra of the Lya forest

Two regimes of Turbulent Fragmentation and the stellar IMF from Primordial to Present Day Star Formation

AMR Simulations of the Cosmological Light Cone: SZE Surveys of the Synthetic Universe

Beyond the Cool Core: The Formation of Cool Core Galaxy Clusters

Adaptive Mesh Refinement for Supersonic Molecular Cloud Turbulence

The Formation of the First Star in the Universe

Baryons in the Warm-Hot Intergalactic Medium

Cosmological Adaptive Mesh Refinement