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.

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.

The origin of supermassive black holes (SMBHs) that inhabit the centers of massive galaxies is largely unconstrained. Remnants from supermassive stars (SMSs) with masses around 10,000 solar masses provide the ideal seed candidates, known as direct collapse black holes. However, their very existence and formation environment in the early Universe are still under debate, with their supposed rarity further exacerbating the problem of modeling their ab-initio formation. SMS models have shown that rapid collapse, with an infall rate above a critical value, in metal-free haloes is a requirement for the formation of a proto-stellar core which will then form an SMS. Using a radiation hydrodynamics simulation of early galaxy formation, we show the natural emergence of metal-free haloes both massive enough, and with sufficiently high infall rates, to form an SMS. We find that haloes that are exposed to both a Lyman-Werner intensity of J_LW ~ 3 J_21 and that undergo at least one period of rapid growth early in their evolution are ideal cradles for SMS formation. This rapid growth induces substantial dynamical heating, amplifying the existing Lyman-Werner suppression originating from a group of young galaxies 20 kiloparsecs away. Our results strongly indicate that structure formation dynamics, rather than a critical Lyman-Werner (LW) flux, may be the main driver of massive black hole formation in the early Universe. We find that massive black hole seeds may be much more common in overdense regions of the early Universe than previously considered with a comoving number density up to 10^-3 Mpc^-3.

We review two decades of progress using the Enzo hydrodynamic cosmology code to simulate the Cosmic Dawn, a period of roughly 1 billion years beginning with the formation of the first stars in the universe, and ending with cosmic reionization. Using simulations of increasing size and complexity, working up in length and mass scale and to lower redshifts, a connected narrative is built up covering the entire epoch. In the first part of the paper, we draw on results we and our collaborators have achieved using the Enzo cosmological adaptive mesh refinement code. Topics include the formation of Population III stars, the transition to Population II star formation, chemical enrichment, the assembly of the first galaxies, their high redshift galaxy statistics, and their role in reionization. In the second part of the paper we highlight physical difficulties that will require new, more physically complex simulations to address, drawing from a broader literature survey. We discuss the healthy interplay between self-consistent numerical simulations and analytic and semi-analytic approaches. Finally, we discuss technical advances in hardware and software that will enable a new class of more realistic simulations to be carried out on exascale supercomputers in the future.

Astrophysical and cosmological phenomena involve a large variety of physical processes, and can encompass an enormous range of scales. To effectively investigate these phenomena computationally, applications must similarly support modeling these phenomena on enormous ranges of scales; furthermore, they must do so efficiently on high-performance computing platforms of ever-increasing parallelism and complexity. We describe Enzo-P, a Petascale redesign of the ENZO adaptive mesh refinement astrophysics and cosmology application, along with Cello, a reusable and scalable adaptive mesh refinement software framework, on which Enzo-P is based. Cello's scalability is enabled by the Charm++ Parallel Programming System, whose data-driven asynchronous execution model is ideal for taking advantage of the available but irregular parallelism in adaptive mesh refinement-based applications. We present scaling results on the NSF Blue Waters supercomputer, and outline our future plans to bring Enzo-P to the Exascale Era by targeting highly-heterogeneous accelerator-based platforms.

Previously we identified a new class of early galaxy that we estimate contributes up to 30\% of the ionizing photons responsible for reionization. These are low mass halos in the range $M_h =10^{6.5}-10^{8} M_{\odot}$ that have been chemically enriched by supernova ejecta from prior Pop III star formation. Despite their low star formation rates, these Metal Cooling halos (MCs) are significant sources of ionizing radiation, especially at the onset of reionization, due to their high number density and ionizing escape fractions. Here we present a fully-coupled radiation hydrodynamic simulation of reionization that includes these MCs as well the more massive hydrogen atomic line cooling halos. Our method is novel: we perform halo finding inline with the radiation hydrodynamical simulation, and assign escaping ionizing fluxes to halos using a probability distribution function (PDF) measured from the galaxy-resolving Renaissance Simulations. The PDF captures the mass dependence of the ionizing escape fraction as well as the probability that a halo is actively forming stars. With MCs, reionization starts earlier than if only halos of $10^8 M_{\odot}$ and above are included, however the redshift when reionization completes is only marginally affected as this is driven by more massive galaxies. Because star formation is intermittent in MCs, the earliest phase of reionization exhibits a stochastic nature, with small H II regions forming and recombining. Only later, once halos of mass $\sim 10^9 M_{\odot}$ and above begin to dominate the ionizing emissivity, does reionization proceed smoothly in the usual manner deduced from previous studies. This occurs at $z\approx 10$ in our simulation.

The formation of stellar mass black holes from the remnants of Population III stars provides a source of initial black hole seeds with the potential to grow into intermediate or, in rare cases, possibly supermassive black holes. We use the Renaissance simulation suite to follow the growth of over 15,000 black holes born into mini-haloes in the early Universe. We compute the evolution of the black holes by post-processing individual remnant Population III star particles in the Renaissance simulation snapshots. The black holes populate haloes from 10$^{6}$ M$_{\odot}$ up to 10$^{9}$ M$_{\odot}$. We find that all of the black holes display very inefficient growth. On average the black holes increase their initial mass by a factor 10$^{-5}$, with the most active black holes increasing their mass by approximately 10%. Only a single black hole experiences any period of super-Eddington accretion, but the duration is very short and not repeated. Furthermore, we find no correlation of black hole accretion with halo mass in the mass range sampled. Within most haloes, we identify clumps of cool, dense gas for which accretion rates would be high, but instances of black holes encountering these clumps are rare and short-lived. Star formation competes with black hole growth by consuming available gas and driving down accretion rates through feedback. We conclude that the black holes born from Population III remnants do not form a significant population of intermediate mass black holes in the early Universe and will need to wait until later times to undergo significant accretion, if at all.

We produce synthetic spectra and observations for metal-free stellar populations and high mass X-ray binaries in the Renaissance Simulations at a redshift of 15. We extend our methodology from the first paper in the series by modelling the production and extinction of emission lines throughout a dusty and metal-enriched interstellar and circum-galactic media extracted from the simulation, using a Monte Carlo calculation. To capture the impact of high-energy photons, we include all frequencies from hard X-ray to far infrared with enough frequency resolution to discern line emission and absorption profiles. The most common lines in our sample in order of their rate of occurrence are Ly$\alpha$, the C IV $\lambda\lambda1548,1551$ doublet, H-$\alpha$, and the Ca II $\lambda\lambda\lambda8498,8542,8662$ triplet. The best scenario for a direct observation of a metal-free stellar population is a merger between two Population III galaxies. In mergers between metal-enriched and metal-free stellar populations, some characteristics may be inferred indirectly. Single Population III galaxies are too dim to be observed photometrically at $z = 15$. Ly$\alpha$ emission is discernible by $JWST$ as an increase in $\rm{J_{200w} - J_{277w}}$ colour off the intrinsic stellar tracks. Observations of metal-free stars will be difficult, though not impossible, with the next generation of space telescopes.

We present synthetic observations for the first generations of galaxies in the Universe and make predictions for future deep field observations for redshifts greater than 6. Due to the strong impact of nebular emission lines and the relatively compact scale of HII regions, high resolution cosmological simulations and a robust suite of analysis tools are required to properly simulate spectra. We created a software pipeline consisting of FSPS, Hyperion, Cloudy and our own tools to generate synthetic IR observations from a fully three-dimensional arrangement of gas, dust, and stars. Our prescription allows us to include emission lines for a complete chemical network and tackle the effect of dust extinction and scattering in the various lines of sight. We provide spectra, 2-D binned photon imagery for both HST and JWST IR filters, luminosity relationships, and emission line strengths for a large sample of high redshift galaxies in the Renaissance Simulations. Our resulting synthetic spectra show high variability between galactic halos with a strong dependence on stellar mass, metallicity, gas mass fraction, and formation history. Halos with the lowest stellar mass have the greatest variability in [OIII]/H$\beta$, [OIII] and CIII] while halos with higher masses are seen to show consistency in their spectra and [OIII] equivalent widths (EW) between 1\AA\ and 10\AA. Viewing angle accounted for three-fold difference in flux due to the presence of ionized gas channels in a halo. Furthermore, JWST color plots show a discernible relationship between redshift, color, and mean stellar age.

We explore the structure and statistics of multiphase, magnetized ISM turbulence in the local Milky Way by means of driven periodic box numerical MHD simulations. Using the higher order-accurate piecewise-parabolic method on a local stencil (PPML), we carry out a small parameter survey varying the mean magnetic field strength and density while fixing the rms velocity to observed values. We quantify numerous characteristics of the transient and steady-state turbulence, including its thermodynamics and phase structure, kinetic and magnetic energy power spectra, structure functions, and distribution functions of density, column density, pressure, and magnetic field strength. The simulations reproduce many observables of the local ISM, including molecular clouds, such as the ratio of turbulent to mean magnetic field at 100 pc scale, the mass and volume fractions of thermally stable HI, the lognormal distribution of column densities, the mass-weighted distribution of thermal pressure, and the linewidth-size relationship for molecular clouds. Our models predict the shape of magnetic field probability density functions (PDFs), which are strongly non-Gaussian, and the relative alignment of magnetic field and density structures. Finally, our models show how the observed low rates of star formation per free-fall time are controlled by the multiphase thermodynamics and large-scale turbulence.

Due to their long mean free paths, X-rays are expected to have many significant impacts globally on the properties of the intergalactic medium (IGM) by their heating and ionizing processes on large scales. At high redshifts, X-rays from Population (Pop) III binaries might have important effects on cosmic reionization and the Lyman alpha forest. As a continuation of our previous work on Pop III binary X-rays (Xu et al. 2014), we use the Pop III distribution and evolution from the Renaissance Simulations, a suite of self-consistent cosmological radiation hydrodynamics simulations of the formation of the first galaxies, to calculate the X-ray luminosity density and background over the redshift range 20 > z > 7.6. As we find that Pop III star formation continues at a low, nearly constant rate to the end of reionization, X-rays are being continuously produced at significant rates compared to other possible X-ray sources, such as AGNs and normal X-ray binaries during the same period of time. We estimate that Pop III binaries produce approximately 6 eV of energy in the X-rays per hydrogen atom. We calculate the X-ray background for different monochromatic photon energies. KeV X-rays redshift and accumulate to produce a strong X-ray background spectrum extending to roughly 500 eV. The X-ray background is strong enough to heat the IGM to ~ 1000 K and to ionize a few percent of the neutral hydrogen. These effects are important for an understanding of the neutral hydrogen hyperfine transition 21-cm line signatures, the Ly alpha forest, and optical depth of the CMB to Thomson scattering.

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.

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.

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 $\alpha$ 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.

We present a simulation of the formation of the earliest Population II stars, starting from cosmological initial conditions and ending when metals created in the first supernovae are incorporated into a collapsing gas-cloud. This occurs after a supernova blast-wave collides with a nearby mini-halo, inducing further turbulence that efficiently mixes metals into the dense gas in the center of the halo. The gas that first collapses has been enriched to a metallicity of Z ~ 2e-5 Zsun. Due to the extremely low metallicity, collapse proceeds similarly to metal-free gas until dust cooling becomes efficient at high densities, causing the cloud to fragment into a large number of low mass objects. This external enrichment mechanism provides a plausible origin for the most metal-poor stars observed, such as SMSS J031300.36-670839.3, that appear to have formed out of gas enriched by a single supernova. This mechanism operates on shorter timescales than the time for low-mass mini-halos (M < 5e5 Msun) to recover their gas after experiencing a supernova. As such, metal-enriched stars will likely form first via this channel if the conditions are right for it to occur. We identify a number of other externally enriched halos that may form stars in this manner. These halos have metallicities as high as 0.01 Zsun, suggesting that some members of the first generation of metal-enriched stars may be hiding in plain sight in current stellar surveys.

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.

We present our prediction for the local 21 cm differential brightness temperature ($\delta 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 $\delta 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 $\delta T_{b}$. A combination of emission ($\delta T_{b}>0$) and absorption ($\delta T_{b}<0$) regions appears in varying amplitudes and angular scales. The boost of the signal by the high-density environment ($\delta\sim0.64$) and on a relatively large scale combine to make Rarepeak a discernible, spatially-extended ($\theta\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.

Galaxy clusters are unique laboratories to investigate turbulent fluid motions and large scale magnetic fields. Synchrotron radio halos at the center of merging galaxy clusters provide the most spectacular and direct evidence of the presence of relativistic particles and magnetic fields associated with the intracluster medium. The study of polarized emission from radio halos is extremely important to constrain the properties of intracluster magnetic fields and the physics of the acceleration and transport of the relativistic particles. However, detecting this polarized signal is a very hard task with the current radio facilities.We use cosmological magneto-hydrodynamical simulations to predict the expected polarized surface brightness of radio halos at 1.4 GHz. We compare these expectations with the sensitivity and the resolution reachable with the SKA1. This allows us to evaluate the potential for studying intracluster magnetic fields in the surveys planned for SKA1.

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.

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.

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.

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

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

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.

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.

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.

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.

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.

Synchrotron radio halos at the center of merging galaxy clusters provide the most spectacular and direct evidence of the presence of relativistic particles and magnetic fields associated with the intracluster medium. The study of polarized emission from radio halos has been shown to be extremely important to constrain the properties of intracluster magnetic fields. However, detecting this polarized signal is a very hard task with the current radio facilities.We investigate whether future radio observatories, such as the Square Kilometer Array (SKA) and its precursors and pathfinders, will be able to detect the polarized emission of radio halos in galaxy clusters.On the basis of cosmological magnetohydrodynamical simulations with initial magnetic fields injected by active galactic nuclei, we predict the expected radio halo polarized signal at 1.4 GHz. We compare these expectations with the limits of current radio facilities and explore the potential of the forthcoming radio interferometers to investigate intracluster magnetic fields through the detection of polarized emission from radio halos.The resolution and sensitivity values that are expected to be obtained in future sky surveys performed at 1.4 GHz using the SKA precursors and pathfinders (like APERTIF and ASKAP) are very promising for the detection of the polarized emission of the most powerful (L1.4GHz>10e25 Watt/Hz) radio halos. Furthermore, the JVLA have the potential to already detect polarized emission from strong radio halos, at a relatively low resolution.However, the possibility of detecting the polarized signal in fainter radio halos (L1.4GHz~10e24 Watt/Hz) at high resolution requires a sensitivity reachable only with SKA.

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, \rho~10^(-28)-10^(-27) g/cm^3, with magnetic field strengths of 0.1-1.0 \mu 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.

We report statistical results for dark matter (DM) velocity anisotropy, \beta, from a sample of some 6000 cluster-size halos (at redshift zero) identified in a \Lambda CDM hydrodynamical adaptive mesh refinement simulation performed with the Enzo code. These include profiles of \beta\ 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, \gamma, and \beta\ at large radii, $0.3 r_{vir} < r < r_{vir}$.

Radio observations of galaxy clusters show that there are $\mu$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 ($\sigma_{RM}$), are consistent at a first-order with the radial distribution from observations. The correlations between the $\sigma_{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.

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.

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.

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.

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\'enic and one trans-Alfv\'enic. 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({\rho},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\'en 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.

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.

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.

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, $\tau_{\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 $\tau_{\rm p}$ in the inertial range. Clustering of these particles occurs primarily around a scale where the eddy turnover time is $\sim\tau_{\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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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(\rho) \propto \rho^{-1.67}$ and $P(B)\propto B^{-2.74}$. A power law is also found between magnetic field strength and density, with $B\propto \rho^{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.

We use extensive measurements of the cluster A1689 to assess the expected similarity in the dynamics of galaxies and dark matter (DM) in their motion as collisionless `particles' in the cluster gravitational potential. To do so we derive the radial profile of the specific kinetic energy of the cluster galaxies from the Jeans equation and observational data. Assuming that the specific kinetic energies of galaxies and DM are roughly equal, we obtain the mean value of the DM velocity anisotropy parameter, and the DM density profile. Since this deduced profile has a scale radius that is higher than inferred from lensing observations, we tested the validity of the assumption by repeating the analysis using results of simulations for the profile of the DM velocity anisotropy. Results of both analyses indicate a significant difference between the kinematics of galaxies and DM within $r \lesssim 0.3r_{\rm vir}$. This finding is reflected also in the shape of the galaxy number density profile, which flattens markedly with respect to the steadily rising DM profile at small radii. Thus, $r \sim 0.3r_{\rm vir}$ seems to be a transition region interior to which collisional effects significantly modify the dynamical properties of the galaxy population with respect to those of DM in A1689

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.

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.

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.

Gravitational and electronic forces produce a correlation between the mass and shape of objects in the universe. For example, at an average radius of ~ 200 km - 300 km, the icy moons and rocky asteroids of our Solar System transition from a rounded potato shape to a sphere. We derive this potato-to-sphere transition radius -- or "potato radius" -- from first principles. Using the empirical potato radii of asteroids and icy moons, we derive a constraint on the yield strength of these bodies during their formative years when their shapes were determined. Our proposed ~ 200 km potato radius for icy moons would substantially increase the number of trans-Neptunian objects classified as dwarf planets.

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.

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.

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.

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.

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.

We present results of large-scale three-dimensional simulations of weakly magnetized supersonic turbulence at grid resolutions up to 1024^3 cells. Our numerical experiments are carried out with the Piecewise Parabolic Method on a Local Stencil and assume an isothermal equation of state. The turbulence is driven by a large-scale isotropic solenoidal force in a periodic computational domain and fully develops in a few flow crossing times. We then evolve the flow for a number of flow crossing times and analyze various statistical properties of the saturated turbulent state. We show that the energy transfer rate in the inertial range of scales is surprisingly close to a constant, indicating that Kolmogorov's phenomenology for incompressible turbulence can be extended to magnetized supersonic flows. We also discuss numerical dissipation effects and convergence of different turbulence diagnostics as grid resolution refines from 256^3 to 1024^3 cells.

We have implemented an Adaptive Mesh Refinement criterion explicitly designed to increase spatial resolution around discontinuities in the velocity field in ENZO cosmological simulations. With this technique, shocks and turbulent eddies developed during the hierarchical assembly of galaxy clusters are followed with unprecedented spatial resolution, even at large distances from the clusters center. By measuring the spectral properties of the gas velocity field, its time evolution and the properties of shocks for a reference galaxy cluster, we investigate the connection between accretion processes and the onset of chaotic motions in the simulated Inter Galactic Medium over a wide range of scales

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.

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

We use three-dimensional simulations to study the statistics of supersonic turbulence in molecular clouds. Our numerical experiments describe driven turbulent flows with an isothermal equation of state, Mach numbers around 10, and various degrees of magnetization. We first support the so-called 1/3-rule of Kritsuk et al. 2007 with our new data from a larger 2048^3 simulation. We then attempt to extend the 1/3-rule to supersonic MHD turbulence and get encouraging preliminary results based on a set of 512^3 simulations. Our results suggest an interesting new approach to tackle universal scaling relations and intermittency in supersonic MHD turbulence.

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.

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 ($\mu$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.

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.

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.

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)

We present new results characterizing cosmological shocks within adaptive mesh refinement N-Body/hydrodynamic simulations that are used to predict non-thermal components of large-scale structure. This represents the first study of shocks using adaptive mesh refinement. We propose a modified algorithm for finding shocks from those used on unigrid simulations that reduces the shock frequency of low Mach number shocks by a factor of ~3. We then apply our new technique to a large, (512 Mpc/h)^3, cosmological volume and study the shock Mach number (M) distribution as a function of pre-shock temperature, density, and redshift. Because of the large volume of the simulation, we have superb statistics that results from having thousands of galaxy clusters. We find that the Mach number evolution can be interpreted as a method to visualize large-scale structure formation. Shocks with Mach<5 typically trace mergers and complex flows, while 5<Mach<20 and Mach>20 generally follow accretion onto filaments and galaxy clusters, respectively. By applying results from nonlinear diffusive shock acceleration models using the first-order Fermi process, we calculate the amount of kinetic energy that is converted into cosmic ray protons. The acceleration of cosmic ray protons is large enough that in order to use galaxy clusters as cosmological probes, the dynamic response of the gas to the cosmic rays must be included in future numerical simulations.

Current numerical studies suggest that the first galaxies formed a few stars at a time and were enriched only gradually by the first heavy elements. However, the large box sizes in these models cannot resolve primordial supernova explosions or the mixing of their metals with ambient gas, which could result in intervening, prompt generations of low-mass stars. We present multiscale 1D models of Population III supernovae in cosmological halos that evolve the blast from its earliest stages as a free expansion. We find that if the star ionizes the halo, the ejecta strongly interacts with the dense shell swept up by the H II region, potentially cooling and fragmenting it into clumps that are unstable to gravitational collapse. If the star fails to ionize the halo, the explosion propagates metals out to 20 - 40 pc and then collapses, enriching tens of thousands of solar masses of primordial gas, in contrast to previous models that suggest such explosions 'fizzle'. Rapid formation of low-mass stars trapped in the gravitational potential well of the halo appears to be inevitable in such circumstances. Consequently, it is possible that far more stars were swept up into the first galaxies, at earlier times and with distinct chemical signatures, than in present models. Upcoming measurements by JWST and ALMA may discriminate between these two paradigms.

We present cosmological magnetohydrodynamic simulations of the formation of a galaxy cluster with magnetic energy feedback from an active galactic nuclei (AGN). We demonstrate that X-ray cavities can be produced by the magnetically dominated jet-lobe system that is supported by a central axial current. The cavities are magnetically dominated and their morphology is determined jointedly by the magnetic fields and the background cluster pressure profile. The expansion and motion of the cavities are driven initially by the Lorentz force of the magnetic fields, and the cavities only become buoyant at late stages ($> 500$ Myr). We find that up to $80%-90%$ of the injected magnetic energy goes into doing work against the hot cluster medium, heating it, and lifting it in the cluster potential.

We report the results of the first self-consistent three-dimensional adaptive mesh refinement magnetohydrodynamical simulations of Population III star formation including the Biermann Battery effect. We find that the Population III stars formed including this effect are both qualitatively and quantitatively similar to those from hydrodynamics-only (non-MHD) cosmological simulations. We observe peak magnetic fields of ~10^-9 G in the center of our star-forming halo at z ~ 17.55. The magnetic fields created by the Biermann Battery effect are predominantly formed early in the evolution of the primordial halo at low density and large spatial scales, and then grow through compression and by shear flows. The fields seen in this calculation are never large enough to be dynamically important (with beta >= 10^{15} at all times), and should be considered the minimum possible fields in existence during Population III star formation, and may be seed fields for the stellar dynamo or the magnetorotational instability at higher densities and smaller spatial scales.

We describe a subgrid model for including galaxies into hydrodynamical cosmological simulations of galaxy cluster evolution. Each galaxy construct- or galcon- is modeled as a physically extended object within which star formation, galactic winds, and ram pressure stripping of gas are modeled analytically. Galcons are initialized at high redshift (z~3) after galaxy dark matter halos have formed but before the cluster has virialized. Each galcon moves self-consistently within the evolving cluster potential and injects mass, metals, and energy into intracluster (IC) gas through a well-resolved spherical interface layer. We have implemented galcons into the Enzo adaptive mesh refinement code and carried out a simulation of cluster formation in a LambdaCDM universe. With our approach, we are able to economically follow the impact of a large number of galaxies on IC gas. We compare the results of the galcon simulation with a second, more standard simulation where star formation and feedback are treated using a popular heuristic prescription. One advantage of the galcon approach is explicit control over the star formation history of cluster galaxies. Using a galactic SFR derived from the cosmic star formation density, we find the galcon simulation produces a lower stellar fraction, a larger gas core radius, a more isothermal temperature profile, and a flatter metallicity gradient than the standard simulation, in better agreement with observations.

We present numerical simulations of primordial supernovae in cosmological minihalos at $z \sim$ 20. We consider Type II supernovae, hypernovae, and pair instability supernovae (PISN) in halos from 6.9 $\times$ 10$^5$ - 1.2 $\times$ 10$^7$ $\Ms$, those in which Population III stars are expected to form via H$_2$ cooling. The supernovae evolve along two evolutionary paths according to whether they explode in \ion{H}{2} regions or neutral halos. Those in \ion{H}{2} regions first expand adiabatically and then radiate strongly upon collision with baryons ejected from the halo during its photoevaporation by the progenitor. Explosions in neutral halos promptly emit most of their kinetic energy as x-rays, but retain enough momentum to seriously disrupt the halo. We find that the least energetic of the supernovae are capable of destroying halos $\lesssim$ 10$^7$ $\Ms$, while a single PISN can destroy even more massive halos. Blasts in \ion{H}{2} regions disperse heavy elements into the IGM, but neutral halos confine the explosion and its metals. In \ion{H}{2} regions, a prompt second generation of stars may form in the remnant at radii of 100 - 200 pc in the halo. Explosions confined by large halos instead recollapse, with infall rates in excess of 10$^{-2}$ $\Ms$ yr$^{-1}$ that heavily contaminate their interior. This fallback may either fuel massive black hole growth at very high redshifts or create the first globular cluster with a radius of 10 - 20 pc at the center of the halo. Our findings allow the possibility that the first primitive galaxies formed sooner, with greater numbers of stars and distinct chemical abundance patterns, than in current models.

Magnetic fields play an important role in almost all astrophysical phenomena including star formation. But due to the difficulty in analytic modeling and observation, magnetic fields are still poorly studied and numerical simulation has become a major tool. We have implemented a cosmological magnetohydrodynamics package for Enzo which is an AMR hydrodynamics code designed to simulate structure formation. We use the TVD solver developed by S. Li as the base solver. In addition, we employ the constrained transport(CT) algorithm as described by D. Balsara. For interpolation magnetic fields to fine grids we used a divergence free quadratic reconstruction, also described by Balsara. We present results from several test problems including MHD caustics, MHD pancake and galaxy cluster formation with magnetic fields. We also discuss possible applications of our AMR MHD code to first star research.

The median observed velocity width v_90 of low-ionization species in damped Ly-alpha systems is close to 90 km/s, with approximately 10% of all systems showing v_90 > 210 km/s at z=3. We show that a relative shortage of such high-velocity neutral gas absorbers in state-of-the-art galaxy formation models is a fundamental problem, present both in grid-based and particle-based numerical simulations. Using a series of numerical simulations of varying resolution and box size to cover a wide range of halo masses, we demonstrate that energy from gravitational infall alone is insufficient to produce the velocity dispersion observed in damped Ly-alpha systems, nor does this dispersion arise from an implementation of star formation and feedback in our highest resolution (~ 45 pc) models, if we do not put any galactic winds into our models by hand. We argue that these numerical experiments highlight the need to separate dynamics of different components of the multiphase interstellar medium at z=3.

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.

We present numerical simulations of how a 120 M$_\odot$ primordial star regulates star formation in nearby cosmological halos at $z \sim$ 20 by photoevaporation. Our models include nine-species primordial chemistry and self-consistent multifrequency conservative transfer of UV photons with all relevant radiative processes. Whether or not new stars form in halos clustered around a Population III star ultimately depends on their core densities and proximity to the star. Diffuse halos with central densities below 2 - 3 cm$^{-3}$ are completely ionized and evaporated anywhere in the cluster. Evolved halos with core densities above 2000 cm$^{-3}$ are impervious to both ionizing and Lyman-Werner flux at most distances from the star and collapse as quickly as they would in its absence. Star formation in halos of intermediate density can be either promoted or suppressed depending on how the I-front remnant shock compresses, deforms and enriches the core with H$_2$. We find that the 120 M$_\odot$ star photodissociates H$_2$ in most halos in the cluster but that catalysis by H- restores it a few hundred kyr after the death of the star, with little effect on star formation. Our models exhibit significant departures from previous one-dimensional spherically-symmetric simulations, which are prone to serious errors due to unphysical geometric focusing effects.

We review recent 3D cosmological hydrodynamic simulations of primordial star formation from cosmological initial conditions (Pop III.1) and from initial conditions that have been altered by radiative feedback from stellar sources (Pop III.2). We concentrate on simulations that resolve the formation of the gravitationally unstable cloud cores in mini-halos over the mass range $10^5 < M/\Msun < 10^7 $ and follow their evolution to densities of at least $10^{10} \cmm3$ and length scales of $<10^{-2}$ pc such that accretion rates can be estimated. The advent of ensembles of such simulations exploring a variety of conditions permits us to assess the robustness of the standard model for Pop III.1 star formation and investigate scatter in their formation redshifts and accretion rates, thereby providing much needed information about the Pop III IMF. The simulations confirm the prediction that Pop III.1 stars were massive ($\sim 100 \Msun$), and form in isolation in primordial mini-halos. Simulations of Pop III.2 star forming in relic HII regions suggest somewhat lower masses ($\sim 30 \Msun$) which may help explain the chemical abundances of extremely metal poor stars. We note that no 3D simulation at present has achieved stellar density let alone followed the entire accretion history of the star in any scenario, and thus the IMF of Pop III stars remains poorly determined theoretically.

Radiative cooling by metals in shocked gas mediates the formation of ionization front instabilities in the galaxy today that are responsible for a variety of phenomena in the interstellar medium, from the morphologies of nebulae to triggered star formation in molecular clouds. An important question in early reionization and chemical enrichment of the intergalactic medium is whether such instabilities arose in the H II regions of the first stars and primeval galaxies, which were devoid of metals. We present three-dimensional numerical simulations that reveal both shadow and thin-shell instabilities readily formed in primordial gas. We find that the hard UV spectra of Population III stars broadened primordial ionization fronts, causing H2 formation capable of inciting violent thin- shell instabilities in D-type fronts, even in the presence of intense Lyman-Werner flux. The high post- front gas temperatures associated with He ionization sustained and exacerbated shadow instabilities, unaided by molecular hydrogen cooling. Our models indicate that metals eclipsed H2 cooling in I-front instabilities at modest concentrations, from 0.001- 0.01 solar. We conclude that ionization front instabilities were prominent in the H II regions of the first stars and galaxies, influencing the escape of ionizing radiation and metals into the early universe.

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.

Flux-limited X-ray samples indicate that about half of rich galaxy clusters have cool cores. Why do only some clusters have cool cores while others do not? In this paper, cosmological N-body + Eulerian hydrodynamic simulations, including radiative cooling and heating, are used to address this question as we examine the formation and evolution of cool core (CC) and non-cool core (NCC) clusters. These adaptive mesh refinement simulations produce both CC and NCC clusters in the same volume. They have a peak resolution of 15.6 h^{-1} kpc within a (256 h^{-1} Mpc)^3 box. Our simulations suggest that there are important evolutionary differences between CC clusters and their NCC counterparts. Many of the numerical CC clusters accreted mass more slowly over time and grew enhanced cool cores via hierarchical mergers; when late major mergers occurred, the CC's survived the collisions. By contrast, NCC clusters experienced major mergers early in their evolution that destroyed embryonic cool cores and produced conditions that prevented CC re-formation. As a result, our simulations predict observationally testable distinctions in the properties of CC and NCC beyond the core regions in clusters. In particular, we find differences between CC versus NCC clusters in the shapes of X-ray surface brightness profiles, between the temperatures and hardness ratios beyond the cores, between the distribution of masses, and between their supercluster environs. It also appears that CC clusters are no closer to hydrostatic equilibrium than NCC clusters, an issue important for precision cosmology measurements.

We study the ionization and thermal evolution of the intergalactic medium during the epoch of \heii reionization by means of radiation hydrodynamical cosmological simulations. We post-process baryonic density fields from a standard optically-thin IGM simulation with a homogeneous galaxy-dominated UV background (UVB) which reionizes \hi and \hei at z=6.5 but does not have any contribution to the ionization of \heii. Quasars with luminosities proportional to the mass of the host halos are then introduced as point sources throughout the 100 Mpc simulation volume consistent with the Pei luminosity function. We evolve the spatial distribution of the \heii ionizing radiation field using a time-implicit variable tensor Eddington factor radiative transfer scheme. Simultaneously, we also solve for the local ionization of \heii to \heii and the associated photoheating of the gas. We find that the percolation of the \heiii regions is essentially complete by z=2.5. When comparing to a self-consistent optically thin simulation we find that in optically thick calculation the gas temperature is higher by a factor of approximately 1.7 at the mean gas density level. We use 300 random lines of sight to compute at $\bar{z} = 2.5 \pm 0.1$ a mean \heii \lya line transmission of $\bar{F} = 0.304 \pm 0.002$. We compare the broadening width of the \hi and \heii \lya lines to the results from the self-consistent optically thin simulation and find a shift by approximately 1.25 km/s of the b-parameter distribution. Estimating the relative broadening width between the two forests shows that the \heii median b-parameter is about 0.8 times the median \hi broadening width. This implies that the \heii absorbers are physically extended consistent with conclusions from observed lines of sight.

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.

We examine aspects of primordial star formation in the presence of a molecular hydrogen-dissociating ultraviolet background. We compare a set of AMR hydrodynamic cosmological simulations using a single cosmological realization but with a range of ultraviolet background strengths in the Lyman-Werner band. This allows us to study the effects of Lyman-Werner radiation on suppressing H2 cooling at low densities as well as the high-density evolution of the collapsing core in a self-consistent cosmological framework. We find that the addition of a photodissociating background results in a delay of the collapse of high density gas at the center of the most massive halo in the simulation and, as a result, an increase in the virial mass of this halo at the onset of baryon collapse. We find that, contrary to previous results, Population III star formation is not suppressed for J$_{21} \geq 0.1$, but occurs even with backgrounds as high as J$_{21} = 1$. We find that H2 cooling leads to collapse despite the depressed core molecular hydrogen fractions due to the elevated H2 cooling rates at $T=2-5 \times 10^3$ K. We observe a relationship between the strength of the photodissociating background and the rate of accretion onto the evolving protostellar cloud core, with higher LW background fluxes resulting in higher accretion rates. Finally, we find that the collapsing halo cores in our simulations do not fragment at densities below $n \sim 10^{10}$ cm$^{-3}$ regardless of the strength of the LW background, suggesting that Population III stars forming in halos with T$_{vir} \sim 10^4$ K may still form in isolation.

We present the first results from a new generation of 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 simulated a very large (512^3h^{-3}Mpc^3) volume with unprecedented dynamic range. We have generated simulated light cones to match the resolution and sensitivity of current and future SZE instruments. Unlike many previous studies of this type, our simulation includes unbound gas, where an appreciable fraction of the baryons in the universe reside. We have found that cluster line-of-sight overlap may be a significant issue in upcoming single-dish SZE surveys. Smaller beam surveys (~1 arcmin) have more than one massive cluster within a beam diameter 5-10% of the time, and a larger beam experiment like Planck has multiple clusters per beam 60% of the time. We explore the contribution of unresolved halos and unbound gas to the SZE signature at the maximum decrement. We find that there is a contribution from gas outside clusters of ~16% per object on average for upcoming surveys. This adds both bias and scatter to the deduced value of the integrated SZE, adding difficulty in accurately calibrating a cluster Y-M relationship. Finally, we find that in images where objects with M > 5x10^{13} M_{\odot} have had their SZE signatures removed, roughly a third of the total SZE flux still remains. This gas exists at least partially in the Warm Hot Intergalactic Medium (WHIM), and will possibly be detectable with the upcoming generation of SZE surveys.

We present numerical simulations of the photoevaporation of cosmological halos clustered around a 120 M$_\odot$ primordial star, confining our study to structures capable of hosting Population III star formation. The calculations include self-consistent multifrequency conservative transfer of UV photons together with nine-species primordial chemistry and all relevant radiative processes. The ultimate fates of these halos varies with central density and proximity to the central source but generally fall into one of four categories. Diffuse halos with central densities below 2 - 3 cm$^{-3}$ are completely ionized and evaporated by the central star anywhere in the cluster. More evolved halo cores at densities above 2000 cm$^{-3}$ are impervious to both ionizing and Lyman-Werner flux at most distances from the star and collapse of their cores proceeds without delay. Radiative feedback in halos of intermediate density can be either positive or negative, depending on how the I-front remnant shock both compresses and deforms the core and enriches it with H$_2$. We find that the 120 M$_\odot$ star photodissociates H$_2$ in most halos within the cluster but that catalysis by H- rapidly restores molecular hydrogen within a few hundred Kyr after the death of the star, with little delay in star formation. Our models exhibit significant departures from previous one-dimensional spherically-symmetric simulations, which are prone to serious errors due to unphysical geometric focusing effects.

Ionization front instabilities have long been of interest for their suspected role in a variety of phenomena in the galaxy, from the formation of bright rims and 'elephant trunks' in nebulae to triggered star formation in molecular clouds. Numerical treatments of these instabilities have historically been limited in both dimensionality and input physics, leaving important questions about their true evolution unanswered. We present the first three-dimensional radiation hydrodynamical calculations of both R-type and D-type ionization front instabilities in galactic environments (i.e., solar metallicity gas). Consistent with linear stability analyses of planar D-type fronts, our models exhibit many short-wavelength perturbations growing at early times that later evolve into fewer large-wavelength structures. The simulations demonstrate that both self-consistent radiative transfer and three-dimensional flow introduce significant morphological differences to unstable modes when compared to earlier two-dimensional approximate models. We find that the amplitude of the instabilities in the nonlinear regime is primarily determined by the efficiency of cooling within the shocked neutral shell. Strong radiative cooling leads to long, extended structures with pronounced clumping while weaker cooling leads to saturated modes that devolve into turbulent flows. These results suggest that expanding H II regions may either promote or provide turbulent support against the formation of later generations of stars, with potential consequences for star formation rates in the galaxy today.

We present results of large-scale three-dimensional simulations of supersonic Euler turbulence with the piecewise parabolic method and multiple grid resolutions up to 2048^3 points. Our numerical experiments describe non-magnetized driven turbulent flows with an isothermal equation of state and an rms Mach number of 6. We discuss numerical resolution issues and demonstrate convergence, in a statistical sense, of the inertial range dynamics in simulations on grids larger than 512^3 points. The simulations allowed us to measure the absolute velocity scaling exponents for the first time. The inertial range velocity scaling in this strongly compressible regime deviates substantially from the incompressible Kolmogorov laws. The slope of the velocity power spectrum, for instance, is -1.95 compared to -5/3 in the incompressible case. The exponent of the third-order velocity structure function is 1.28, while in incompressible turbulence it is known to be unity. We propose a natural extension of Kolmogorov's phenomenology that takes into account compressibility by mixing the velocity and density statistics and preserves the Kolmogorov scaling of the power spectrum and structure functions of the density-weighted velocity v=\rho^{1/3}u. The low-order statistics of v appear to be invariant with respect to changes in the Mach number. For instance, at Mach 6 the slope of the power spectrum of v is -1.69, and the exponent of the third-order structure function of v is unity. We also directly measure the mass dimension of the "fractal" density distribution in the inertial subrange, D_m = 2.4, which is similar to the observed fractal dimension of molecular clouds and agrees well with the cascade phenomenology.

We use large-scale three-dimensional simulations of supersonic Euler turbulence to study the physics of a highly compressible cascade. Our numerical experiments describe non-magnetized driven turbulent flows with an isothermal equation of state and an rms Mach number of 6. We find that the inertial range velocity scaling deviates strongly from the incompressible Kolmogorov laws. We propose an extension of Kolmogorov's K41 phenomenology that takes into account compressibility by mixing the velocity and density statistics and preserves the K41 scaling of the density-weighted velocity v=rho^{1/3}u. We show that low-order statistics of 'v' are invariant with respect to changes in the Mach number. For instance, at Mach 6 the slope of the power spectrum of 'v' is -1.69 and the third-order structure function of 'v' scales linearly with separation. We directly measure the mass dimension of the "fractal" density distribution in the inertial subrange, D_m=2.4, which is similar to the observed fractal dimension of molecular clouds and agrees well with the cascade phenomenology.

In this paper we describe our massively parallel version of Enzo, a multiphysics, parallel, AMR application for simulating cosmological structure formation developed at UCSD and Columbia. We describe its physics, numerical algorithms, implementation, and performance on current terascale platforms. We also discuss our future plans and some of the challenges we face as we move to the petascale.

We have analyzed a large sample of numerically simulated clusters to demonstrate the adverse effects resulting from use of X-ray fitted beta-model parameters with Sunyaev-Zeldovich effect (SZE) data. There is a fundamental incompatibility between model fits to X-ray surface brightness profiles and those done with SZE profiles. Since observational SZE radial profiles are in short supply, the X-ray parameters are often used in SZE analysis. We show that this leads to biased estimates of the integrated Compton y-parameter inside r_{500} calculated from clusters. We suggest a simple correction of the method, using a non-isothermal beta-model modified by a universal temperature profile, which brings these calculated quantities into closer agreement with the true values.

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.

We test a method of estimating the power spectrum of turbulence in molecular clouds based on the comparison of power spectra of integrated intensity maps and single-velocity-channel maps, suggested by Lazarian and Pogosyan. We use synthetic 13CO data from non-LTE radiative transfer calculations based on density and velocity fields of a simulation of supersonic hydrodynamic turbulence. We find that the method yields the correct power spectrum with good accuracy. We then apply the method to the Five College Radio Astronomy Observatory 13CO map of the Perseus region, from the COMPLETE website. We find a power law power spectrum with slope beta=1.81+-0.10. The values of beta as a function of velocity resolution are also confirmed using the lower resolution map of the same region obtained with the AT&T Bell Laboratories antenna. Because of its small uncertainty, this result provides a useful constraint for numerical codes used to simulate molecular cloud turbulence.

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.

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.

(abridged) We perform 12 extremely high resolution adaptive mesh refinement cosmological hydrodynamic simulations of Population III star formation in a Lambda CDM universe, varying the box size and large-scale structure, to understand systematic effects in the formation of primordial protostellar cores. We find results that are qualitatively similar to those observed previously. We observe that the threshold halo mass for formation of a Population III protostar does not evolve significantly with time in the redshift range studied (33 > z > 19) but exhibits substantial scatter due to different halo assembly histories: Halos which assembled more slowly develop cooling cores at lower mass than those that assemble more rapidly, in agreement with Yoshida et al. (2003). We do, however, observe significant evolution in the accretion rates of Population III protostars with redshift, with objects that form later having higher maximum accretion rates, with a variation of two orders of magnitude (10^-4 - 10^-2 Msolar/year). This can be explained by considering the evolving virial properties of the halos with redshift and the physics of molecular hydrogen formation at low densities. Our result implies that the mass distribution of Population III stars inferred from their accretion rates may be broader than previously thought, and may evolve with redshift. Finally, we observe that our collapsing protostellar cloud cores do not fragment, consistent with previous results, which suggests that Population III stars which form in halos of mass 10^5 - 10^6 Msun always form in isolation.

We present the results of three-dimensional simulations of supersonic Euler turbulence with grid resolutions up to 1024^3 points. Our numerical experiments describe nonmagnetized driven turbulent flows with an isothermal equation of state and an rms Mach number of 6. We demonstrate that the inertial range scaling properties of turbulence in this strongly compressible regime deviate substantially from a Kolmogorov-like behavior previously recovered for mildly compressible transonic flows.

Withdrawn at the current time due to errors.

Maps of estimated dust column density in molecular clouds are usually assumed to reliably trace the total gas column density structure. In this work we present results showing a clear discrepancy between the dust and the gas distribution in the Taurus molecular cloud complex. We compute the power spectrum of a 2MASS extinction map of the Taurus region and find it is much shallower than the power spectrum of a 13CO map of the same region previously analyzed. This discrepancy may be explained as the effect of grain growth on the grain extinction efficiency. However, this would require a wide range of maximum grain sizes, which is ruled out based on constraints from the extinction curve and the available grain models. We show that major effects due to CO formation and depletion are also ruled out. Our result may therefore suggest the existence of intrinsic spatial fluctuations of the dust to gas ratio, with amplitude increasing toward smaller scales. Preliminary results of numerical simulations of trajectories of inertial particles in turbulent flows illustrate how the process of clustering of dust grains by the cloud turbulence may lead to observable effects. However, these results cannot be directly applied to large scale supersonic and magnetized turbulence at present.

In this paper we examine aspects of primordial star formation in a gravitino warm dark matter universe with a cosmological constant. We compare a set of simulations using a single cosmological realization but with a wide range of warm dark matter particle masses which have not yet been conclusively ruled out by observations. The addition of a warm dark matter component to the initial power spectrum results in a delay in the collapse of high density gas at the center of the most massive halo in the simulation and, as a result, an increase in the virial mass of this halo at the onset of baryon collapse. Both of these effects become more pronounced as the warm dark matter particle mass becomes smaller. A cosmology using a gravitino warm dark matter power spectrum assuming a particle mass of m_{WDM} ~ 40keV is effectively indistinguishable from the cold dark matter case, whereas the m_{WDM} ~ 15 keV case delays star formation by approx. 10^8 years. There is remarkably little scatter between simulations in the final properties of the primordial protostar which forms at the center of the halo, possibly due to the overall low rate of halo mergers which is a result of the WDM power spectrum. The detailed evolution of the collapsing halo core in two representative WDM cosmologies is described. At low densities (n_{b} <= 10^5 cm^{-3}), the evolution of the two calculations is qualitatively similar, but occurs on significantly different timescales, with the halo in the lower particle mass calculation taking much longer to evolve over the same density range and reach runaway collapse. Once the gas in the center of the halo reaches relatively high densities (n_{b} >= 10^5 cm^{-3}) the overall evolution is essentially identical in the two calculations.

We critically analyze the measurement of galaxy cluster gas masses, which is central to cosmological studies that rely on the galaxy cluster gas mass fraction. Using synthetic observations of numerically simulated clusters viewed through their X-ray emission and thermal Sunyaev-Zeldovich effect (SZE), we reduce the observations to obtain measurements of the cluster gas mass. We are thus able to quantify the possible sources of uncertainty and systematic bias associated with the common simplifying assumptions used in reducing real cluster observations including isothermality and hydrostatic equilibrium. We find that intrinsic variations in clusters limit the precision of observational gas mass estimation to ~10% to 1 sigma confidence excluding instrumental effects. Gas mass estimates performed via all methods surprisingly show little or no trending in the scatter as a function of cluster redshift. For the full cluster sample, methods that use SZE profiles out to roughly the virial radius are the simplest, most accurate, and unbiased way to estimate cluster mass. X-ray methods are systematically more precise mass estimators than are SZE methods if merger and cool core systems are removed, but X-ray methods slightly overestimate (5-10%) the cluster gas mass on average. We find that cool core clusters in our samples are particularly poor candidates for observational mass estimation, even when excluding emission from the core region. The effects of cooling in the cluster gas alter the radial profile of the X-ray and SZE surface brightness outside the cool core region, leading to poor gas mass estimates in cool core clusters. Finally, we find that methods using a universal temperature profile estimate cluster masses to higher precision than those assuming isothermality.

This paper describes ZEUS-MP, a multi-physics, massively parallel, message- passing implementation of the ZEUS code. ZEUS-MP differs significantly from the ZEUS-2D code, the ZEUS-3D code, and an early "version 1" of ZEUS-MP distributed publicly in 1999. ZEUS-MP offers an MHD algorithm better suited for multidimensional flows than the ZEUS-2D module by virtue of modifications to the Method of Characteristics scheme first suggested by Hawley and Stone (1995), and is shown to compare quite favorably to the TVD scheme described by Ryu et. al (1998). ZEUS-MP is the first publicly-available ZEUS code to allow the advection of multiple chemical (or nuclear) species. Radiation hydrodynamic simulations are enabled via an implicit flux-limited radiation diffusion (FLD) module. The hydrodynamic, MHD, and FLD modules may be used in one, two, or three space dimensions. Self gravity may be included either through the assumption of a GM/r potential or a solution of Poisson's equation using one of three linear solver packages (conjugate-gradient, multigrid, and FFT) provided for that purpose. Point-mass potentials are also supported. Because ZEUS-MP is designed for simulations on parallel computing platforms, considerable attention is paid to the parallel performance characteristics of each module. Strong-scaling tests involving pure hydrodynamics (with and without self-gravity), MHD, and RHD are performed in which large problems (256^3 zones) are distributed among as many as 1024 processors of an IBM SP3. Parallel efficiency is a strong function of the amount of communication required between processors in a given algorithm, but all modules are shown to scale well on up to 1024 processors for the chosen fixed problem size.

We use high resolution Eulerian hydrodynamics simulations to study kinematic properties of the low ionization species in damped Ly-alpha systems at redshift z=3. Our adaptive mesh refinement simulations include most key ingredients relevant for modeling neutral gas in high-column density absorbers: hydrodynamics, gravitational collapse, continuum radiative transfer and gas chemistry, but no star formation. We model high-resolution Keck spectra with unsaturated low ion transitions in two Si II lines (1526 and 1808 A), and compare simulated line profiles to the data from the SDSS DLA survey. We find that with increasing grid resolution the models show a trend in convergence towards the observed distribution of HI column densities. While in our highest resolution model we recover the cumulative number of DLAs per unit absorption distance, none of our models predicts DLA velocity widths as high as indicated by the data, suggesting that feedback from star formation might be important. At z=3 a non-negligible fraction of DLAs with column densities below 10^21 cm^-2 is caused by tidal tails due to galaxy-galaxy interactions in more massive halo environments. Lower column density absorbers with N_HI < 10^21.4 cm^-2 are sensitive to changes in the UV background resulting in a 10% reduction of the cumulative number of DLAs for twice the quasar background relative to the fiducial value. We find that the mass cut-off below which a large fraction of dwarf galaxies cannot retain gas after reionization is 7*10^7 msun, lower than the previous estimates. Finally, we show that models with self-shielding commonly used in the literature produce significantly lower DLA velocity widths than the full radiative transfer runs.

Radiative transfer simulations are now at the forefront of numerical astrophysics. They are becoming crucial for an increasing number of astrophysical and cosmological problems; at the same time their computational cost has come to the reach of currently available computational power. Further progress is retarded by the considerable number of different algorithms (including various flavours of ray-tracing and moment schemes) developed, which makes the selection of the most suitable technique for a given problem a non-trivial task. Assessing the validity ranges, accuracy and performances of these schemes is the main aim of this paper, for which we have compared 11 independent RT codes on 5 test problems: (0) basic physics, (1) isothermal H II region expansion and (2) H II region expansion with evolving temperature, (3) I-front trapping and shadowing by a dense clump, (4) multiple sources in a cosmological density field. The outputs of these tests have been compared and differences analyzed. The agreement between the various codes is satisfactory although not perfect. The main source of discrepancy appears to reside in the multi-frequency treatment approach, resulting in different thicknesses of the ionized-neutral transition regions and different temperature structure. The present results and tests represent the most complete benchmark available for the development of new codes and improvement of existing ones. To this aim all test inputs and outputs are made publicly available in digital form.

Scattering of the cosmic microwave background (CMB) in clusters of galaxies polarizes the radiation. We explore several polarization components which have their origin in the kinematic quadrupole moments induced by the motion of the scattering electrons, either directed or random. Polarization levels and patterns are determined in a cluster simulated by the hydrodynamical Enzo code. We find that polarization signals can be as high as $\sim 1 \mu$K, a level that may be detectable by upcoming CMB experiments.

We review the theory of the formation of galaxy clusters and discuss their role as cosmological probes. We begin with the standard cosmological framework where we discuss the origin of the CDM matter power spectrum and the growth of density fluctuations in the linear regime. We then summarize the spherical top-hat model for the nonlinear growth of fluctuations from which scaling relations and halo statistics are derived. Numerical methods for simulating gas in galaxy clusters are then overviewed with an emphasis on multiscale hydrodynamic simulations of cluster ensembles. Results of hydrodynamic AMR simulations are described which compare cluster internal and statistical properties as a function of their assumed baryonic processes. Finally, we compare various methods of measuring cluster masses using X-ray and the thermal Sunyaev-Zeldovich effect (SZE). We find that SZE offers great promise for precision measurements in raw samples of high-z clusters.

Radiation hydrodynamical transport of ionization fronts in the next generation of cosmological reionization simulations holds the promise of predicting UV escape fractions from first principles as well as investigating the role of photoionization in feedback processes and structure formation. We present a multistep integration scheme for radiative transfer and hydrodynamics for accurate propagation of I-fronts and ionized flows from a point source in cosmological simulations. The algorithm is a photon-conserving method which correctly tracks the position of I-fronts at much lower resolutions than non-conservative techniques. The method applies direct hierarchical updates to the ionic species, bypassing the need for the costly matrix solutions required by implicit methods while retaining sufficient accuracy to capture the true evolution of the fronts. We review the physics of ionization fronts in power-law density gradients, whose analytical solutions provide excellent validation tests for radiation coupling schemes. The advantages and potential drawbacks of direct and implicit schemes are also considered, with particular focus on problem timestepping which if not properly implemented can lead to morphologically plausible I-front behavior that nonetheless departs from theory. We also examine the effect of radiation pressure from very luminous central sources on the evolution of I-fronts and flows.

There has been considerable theoretical debate over whether photoionization and supernova feedback from the first Population III stars facilitate or suppress the formation of the next generation of stars. We present results from an Eulerian adaptive mesh refinement simulation demonstrating the formation of a primordial star within a region ionized by an earlier nearby star. Despite the higher temperatures of the ionized gas and its flow out of the dark matter potential wells, this second star formed within 23 million years of its neighbor's death. The enhanced electron fraction within the HII region catalyzes rapid molecular hydrogen formation that leads to faster cooling in the subsequent star forming halos than in the first halos. This "second generation" primordial protostar has a much lower accretion rate because, unlike the first protostar, it forms in a rotationally supported disk of approx. 10-100 solar masses. This is primarily due to the much higher angular momentum of the halo in which the second star forms. In contrast to previously published scenarios, such configurations may allow binaries or multiple systems of lower mass stars to form. These first high resolution calculations offer insight into the impact of feedback upon subsequent populations of stars and clearly demonstrate how primordial chemistry promotes the formation of subsequent generations of stars even in the presence of the entropy injected by the first stars into the IGM.

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.

We compare two cosmological hydrodynamic simulation codes in the context of hierarchical galaxy formation: The SPH code GADGET, and the Eulerian AMR code ENZO. Both codes represent dark matter with the N-body method, but use different gravity solvers and fundamentally different approaches to hydrodynamics. We compare the GADGET `entropy conserving' SPH formulation with two ENZO methods: The piecewise parabolic method (PPM), and the artificial viscosity-based scheme used in the ZEUS code. In this paper we focus on a comparison of cosmological simulations that follow either only dark matter, or also adiabatic baryonic gas. The dark matter-only runs agree generally quite well, provided ENZO is run with a comparatively fine root grid and a low overdensity threshold for mesh refinement, otherwise the abundance of low-mass halos is suppressed. This is due to the hierarchical particle-mesh method used to compute gravitational forces in ENZO, which tends to deliver lower force resolution than the tree algorithm of GADGET. At comparable force resolution, we find that the latter offers substantially better performance and lower memory consumption than the present gravity solver in ENZO. In simulations that include adiabatic gas dynamics, we find general agreement in the distribution functions of temperature, entropy, and density for gas of moderate to high overdensity, as found inside dark matter halos. However, there are some significant differences at lower overdensities. We argue that these discrepancies are presumably owing to differences in the shock-capturing abilities of the different methods. In particular, ZEUS hydro leads to some unphysical heating at early times in preshock regions. Overall, the GADGET hydro results are bracketed by those for ENZO/ZEUS and ENZO/PPM. (abridged)

We use high quality echelle spectra of 24 QSOs to provide a calibrated measurement of the total amount of Lyman-alpha forest absorption (DA) over the redshift range 2.2 < z < 3.2. Our measurement of DA excludes absorption from metal lines or the Lyman-alpha lines of Lyman limit systems and damped Lyman-alpha systems. We use artificial spectra with realistic flux calibration errors to show that we are able to place continuum levels that are accurate to better than 1%. When we combine our results with our previous results between 1.6 < z < 2.2, we find that the redshift evolution of DA is well described over 1.6 < z < 3.2 as A*(1+z)^gamma, where A = 0.0062 and gamma = 2.75. We detect no significant deviations from a smooth power law evolution over the redshift range studied. We find less H I absorption than expected at z=3, implying that the UV background is about 40% higher than expected. Our data appears to be consistent with an H I ionization rate of Gamma ~ 1.4E-12 1/s.

We investigate empirical scaling relations between the thermal Sunyaev-Zeldovich effect (SZE) and cluster mass in simulated clusters of galaxies. The simulated clusters have been compiled from four different samples that differ only in their assumed baryonic physics. We show that the strength of the thermal SZE integrated over a significant fraction of the virialized region of the clusters is relatively insensitive to the detailed heating and cooling processes in the cores of clusters by demonstrating that the derived scaling relations are nearly identical between the four cluster samples considered. For our synthetic images, the central Comptonization parameter shows significant boosting during transient merging events, but the integrated SZE appears to be relatively insensitive to these events. Most importantly, the integrated SZE closely tracks the underlying cluster mass. Observations through the thermal SZE allow a strikingly accurate mass estimation from relatively simple measurements that do not require either parametric modeling or geometric deprojection and thus avoid assumptions regarding the physics of the ICM or the symmetry of the cluster. This result offers significant promise for precision cosmology using clusters of galaxies.

Accretion rates of order 10^-8 M_\odot/yr are observed in young protostars of approximately a solar mass with evidence of circumstellar disks. The accretion rate is significantly lower for protostars of smaller mass, approximately proportional to the second power of the stellar mass, \dot{M}_accr\propto M^2. The traditional view is that the observed accretion is the consequence of the angular momentum transport in isolated protostellar disks, controlled by disk turbulence or self--gravity. However, these processes are not well understood and the observed protostellar accretion, a fundamental aspect of star formation, remains an unsolved problem. In this letter we propose the protostellar accretion rate is controlled by accretion from the large scale gas distribution in the parent cloud, not by the isolated disk evolution. Describing this process as Bondi--Hoyle accretion, we obtain accretion rates comparable to the observed ones. We also reproduce the observed dependence of the accretion rate on the protostellar mass. These results are based on realistic values of the ambient gas density and velocity, as inferred from numerical simulations of star formation in self--gravitating turbulent clouds.

We present 40 fully hydrodynamical numerical simulations of the intergalactic gas that gives rise to the Ly-alpha forest. We make artificial spectra from each and measure three output parameters: the mean flux, a measure of the most common Ly-alpha line width (b-value), and the 1D power spectrum of the flux on scales from 0.01 - 0.1 s/km. We investigate how these output parameters are correlated with the input parameters that describe the simulations. These input parameters include the main cosmological parameters; and two astrophysical parameters gammah and gammahe. The parameter gammah controls the rate of ionization of H I, He I and He II while gammahe controls the rate of heating. We compare the output from the simulations to data from HIRES and UVES spectra and we give a new measurement of the flux power from these spectra at z=1.95. We find that simulations with a wide variety of sigma_8 values, from at least 0.8 - 1.1, can fit the small scale flux power and b-values when we adjust gammahe to compensate for the sigma_8 change. We can also use gammah to adjust the H I ionization rate to simultaneously match the mean flux. We derive scaling relations that give the output parameter values expected for a variety of input parameters with errors of 2% for optical depth and 1.4% for b-balues. We can break the degeneracy between gammahe and sigma_8 using large scale power or other data to fix sigma_8.

We use hydrodynamic cosmological simulations in a 9.6 Mpc box to explore the evolution of the intergalactic medium (IGM) transmissivity from z=2 through the epoch of reionization. Reionization is achieved through an ultraviolet background (UVB) that includes evolving stellar and QSO source populations. We construct and analyze noiseless synthetic HI Lya absorption spectra along lines of sight through our continuously evolving box and find a smooth evolution of the effective optical depth under a power law up to the epoch of reionization. Crossing into the epoch of reionization, both the mean transmitted flux (MTF) and variance to the mean transmitted flux sharply deviate from a smooth evolution and hence such observables have large margins of error. Despite the statistical uncertainty in inferring the reionization profile from spectra, the end of an opacity phase transition of the IGM correlates well with the redshift when both the mean and variance of the transmitted flux rapidly deviate from the post-reionization profile. However, an unobtainable number of lines of sight is needed to allow an estimate of the MTF with less than 10% relative margin of error. In addition to optical transmission, we compare the predicted dark gap length distribution with observations. We show that this statistic is sensitive to spectral resolution at reionization redshifts, but overall in agreement with results by Songaila and Cowie (2002). Finally, we derive a positive correlation between the mean optical depth within a gap and the size of the gap, in attempt to relate "transmission statistics" to "dark gap statistics" in high redshift studies of the IGM.

The origin of brown dwarfs (BDs) is an important component of the theory of star formation, because BDs are approximately as numerous as solar mass stars. It has been suggested that BDs originate from the gravitational fragmentation of protostellar disks, a very different mechanism from the formation of hydrogen burning stars. We propose that BDs are instead formed by the process of turbulent fragmentation, like more massive stars. In numerical simulations of turbulence and star formation we find that gravitationally unstable density peaks of BD mass are commonly formed by the turbulent flow. These density peaks collapse into BD mass objects with circumstellar disks, like more massive protostars. We rely on numerical experiments with very large resolution, achieved with adaptive mesh refinement (AMR). The turbulence simulation presented here is the first AMR turbulence experiment ever attempted and achieves an effective resolution of 1024^3 computational zones. The star formation simulation achieves an effective resolution of (10^6)^3 computational zones, from a cloud size of 5pc to protostellar disks resolved down to 1AU.

Using H I absorption alone, we attempt to separate H I absorption lines in quasar spectra into two categories; HDLs (Higher Density Lines) and LDLs (Lower Density Lines), and we discuss the difference in their physical properties. We deblend and fit all H I lines with Voigt profiles, and make an unbiased sample of H I lines covering a wide column density range (12 < log N_HI < 19 cm^-2). To reduce the influence of line blending, we simultaneously fit several Lyman series lines. As a result of a two-point correlation analysis, we found that higher column density H I lines are clustering at $\Delta$ v < 200 km/s, while lower ones at $\Delta$ v < 100 km/s. We define HDLs as H I lines with 15 < log N_HI < 19 cm^-2 and all H I lines within +/- 200 km/s of a line with log N_HI > 15 cm^-2, and LDLs as others with 12 < log N_HI < 15 cm^-2. We found that the HDLs have smaller minimum b-values for a given column density than the LDLs. This difference is successfully reproduced by our Hydrodynamic simulation. The LDLs seem to be cool or shock-heated diffuse IGM gas, while the HDLs are likely to be cooler dense gas near to galaxies.

We identify a concordant model for the intergalactic medium (IGM) at redshift z=1.9 that uses popular values for cosmological and astrophysical parameters and accounts for all baryons with an uncertainty of 6%. We have measured the amount of absorption, DA, in the Ly-alpha forest at redshift 1.9 in spectra of 77 QSO from the Kast spectrograph. We calibrated the continuum fits with realistic artificial spectra, and we found that averaged over all 77 QSOs the mean continuum level is within 1-2% of the correct value. Absorption from all lines in the Ly-alpha forest at z=1.9 removes DA=15.1 +/- 0.7% of the flux between 1070 and 1170 (rest) Angstroms. This is the first measurement using many QSOs at this z, and the first calibrated measurement at any redshift. Metal lines absorb 2.3 +/- 0.5%, and LLS absorb 1.0 +/- 0.4% leaving 11.8 +/- 1.0% from the lower density bulk of the IGM. Averaging over Delta z=0.1 or 154 Mpc, the dispersion is 6.1 +/- 0.3% including LLS and metal lines, or 3.9 (+0.5, -0.7)% for the lower density IGM alone, consistent with the usual description of large scale structure. LLS and metal lines are major contributors to the variation in the mean flux, and they make the flux field significantly non-Gaussian. We find that a hydrodynamic simulation on a 1024 cubed grid in a 75.7 Mpc box reproduces the observed DA from the low density IGM with parameters values H_o=71 km/s/Mpc, Omega_Lambda=0.73, Omega_m=0.27, Omega_b=0.044, sigma_8=0.9 and a UV background that has an ionization rate that is 1.08 +/- 0.08 times the prediction by Madau, Haardt & Rees (1999).

We employ a generalization of the She & Leveque model to study velocity scaling relations based on our simulations of thermal instability--induced turbulence. Being a by-product of the interstellar phase transition, such multiphase turbulence tends to be more intermittent than compressible isothermal turbulence. Due to radiative cooling, which promotes nonlinear instabilities in supersonic flows, the Hausdorff dimension of the most singular dissipative structures, D, can be as high as 2.3, while in supersonic isothermal turbulence D is limited by shock dissipation to D < 2. We also show that single-phase velocity statistics carry only incomplete information on the turbulent cascade in a multiphase medium. We briefly discuss the possible implications of these results on the hierarchical structure of molecular clouds and on star formation.

I review cosmological simulations of X-ray clusters. Simulations have increased in resolution dramatically and the effects of radiative cooling, star formation feedback, and chemical enrichment on the ICM are being simulated. The structure and evolution of non-radiative X-ray clusters is now well characterized. Such models fail to reproduce the observed L_x-T relation, implying the need for additional physics. Simulations adding radiative cooling produce too much cool gas and unreasonably high X-ray luminosities. Simulations including star formation and feedback appear more promising, but need further refinement. New observations should help in this regard.

We simulate the ionization environment of z ~ 20 luminous objects formed within the framework of the current CDM cosmology and compute their UV escape fraction. These objects are likely single very massive stars that are copious UV emitters. We present analytical estimates as well as one--dimensional radiation hydrodynamical calculations of the evolution of these first HII regions in the universe. The initially D--type ionization front evolves to become R--type within $\lesssim 10^5$ yrs at a distance $\sim1$ pc. This ionization front then completely overruns the halo, accelerating an expanding shell of gas outward to velocities in excess of 30 km s$^{-1}$, about ten times the escape velocity of the confining dark matter halo. We find that the evolution of the HII region depends only weakly on the assumed stellar ionizing luminosities. Consequently, most of the gas surrounding the first stars will leave the dark halo whether or not the stars produce supernovae. If they form the first massive seed black holes these are unlikely to accrete within a Hubble time after they formed until they are incorporated into larger dark matter halos that contain more gas. Because these I--fronts exit the halo on timescales much shorter than the stars' main sequence lifetimes their host halos have UV escape fractions of $\gtrsim 0.95$, fixing an important parameter for theoretical studies of cosmological hydrogen reionization.

In this paper we introduce Enzo, a 3D MPI-parallel Eulerian block-structured adaptive mesh refinement cosmology code. Enzo is designed to simulate cosmological structure formation, but can also be used to simulate a wide range of astrophysical situations. Enzo solves dark matter N-body dynamics using the particle-mesh technique. The Poisson equation is solved using a combination of fast fourier transform (on a periodic root grid) and multigrid techniques (on non-periodic subgrids). Euler's equations of hydrodynamics are solved using a modified version of the piecewise parabolic method. Several additional physics packages are implemented in the code, including several varieties of radiative cooling, a metagalactic ultraviolet background, and prescriptions for star formation and feedback. We also show results illustrating properties of the adaptive mesh portion of the code. Enzo is publicly available and can be downloaded at http://cosmos.ucsd.edu/enzo/ .

I survey the use and impact of adaptive mesh refinement (AMR) simulations in numerical astrophysics and cosmology. Two basic techniques are in use to extend the dynamic range of Eulerian grid simulations in multi-dimensions: cell refinement, and patch refinement, otherwise known as block-structured adaptive mesh refinement (SAMR). In this survey, no attempt is made to assess the relative merits of these two approaches. Rather, the discussion focuses on how AMR is being used and how AMR is making a scientific impact in a diverse set of fields from space physics to the cosmology of the early universe. The increased adoption of AMR techniques in the past decade is driven in part by the public availability of AMR codes and frameworks. I provide a partial list of resources for those interested in learning more about AMR simulations.

We examine the dynamical evolution and statistical properties of the supernova ejecta of massive primordial stars in a cosmological framework to determine whether this first population of stars could have enriched the universe to the levels and dispersions seen by the most recent observations of the Lyman-Alpha forest. We evolve a lambda CDM model in a 1 Mpc^3 volume to a redshift of z = 15 and add ``bubbles'' of metal corresponding to the supernova ejecta of the first generation of massive stars in all dark matter halos with masses greater than 5 times 10^5 solar masses. These initial conditions are then evolved to z = 3 and the distribution and levels of metals are compared to observations. In the absence of further star formation the primordial metal is initially contained in halos and filaments. Photoevaporation of metal-enriched gas due to the metagalactic ultraviolet background radiation at the epoch of reionization (z ~ 6) causes a sharp increase of the metal volume filling factor. At z = 3, ~ 2.5% of the simulation volume (approx. 20% of the total gas mass) is filled with gas enriched above a metallicity of 10^-4 Z_solar, and less than 0.6% of the volume is enriched above a metallicity of 10^-3 Z_solar. This suggests that, even with the most optimistic prescription for placement of primordial supernova and the amount of metals produced by each supernova, this population of stars cannot entirely be responsible for the enrichment of the Lyman-$\alpha$ forest to the levels and dispersions seen by current observations unless we have severely underestimated the duration of the Pop III epoch. However, comparison to observations show that Pop III supernovae can be significant contributors to the very low overdensity Lyman-Alpha forest.

We use ZEUS-MP to perform high resolution, three-dimensional, super-Alfvenic turbulent simulations in order to investigate the role of magnetic fields in self-gravitating core formation within turbulent molecular clouds. Statistical properties of our super-Alfvenic model without gravity agree with previous similar studies. Including self-gravity, our models give the following results. They are consistent with the turbulent fragmentation prediction of the core mass distribution of Padoan & Nordlund. They also confirm that local gravitational collapse is not prevented by magnetohydrodynamic waves driven by turbulent flows, even when the turbulent Jeans mass exceeds the mass in the simulation volume. Comparison of results between 256^3 and 512^3 zone simulations reveals convergence in the collapse rate. Analysis of self-gravitating cores formed in the simulation shows that: (1) All cores formed are magnetically supercritical by at least an order of magnitude. (2) A power law relation between central magnetic field strength and density B_c propto rho_c^{1/2} is observed despite the cores being strongly supercritical. (3) Specific angular momentum j propto R^{3/2} for cores with radius R. (4) Most cores are prolate and triaxial in shape, in agreement with the results of Gammie et al.

Cooling flows, cluster mergers and motions of galaxies through the cluster gas with supersonic and sonic velocities must lead to large scale motions of the intracluster medium (ICM). A high-resolution numerical simulation of X-ray cluster formation by Norman and Bryan predicts cluster-wide turbulence with v_turb ~ 300-600 ~km/s and eddy scales l_outer ~ 100-500 kpc, the larger numbers being characteristic of turbulence near the virial radius, while the smaller numbers pertain to the core. The simulation also predicts the existence of ordered bulk flows in the core with v ~ 400 km/s on scales of several hundred kpc. In this paper we consider the observability of such fluid motions via the distortions it induces in the CMB via the kinematic SZ effect, as well as via Doppler broadening and shifting of metal lines in the X-ray spectrum. We estimate |\Delta T/T|_kinematic < few X 10^{-6} -- at or below current limits of detectability. However, we find that an energy resolution of a few eV is sufficient to detect several Doppler shifted components in the 6.7 keV Fe line in the core of the cluster.

We present a new model for the creation of cool cores in rich galaxy clusters within a LambdaCDM cosmological framework using the results from high spatial dynamic range, adaptive mesh hydro/N-body simulations. It is proposed that cores of cool gas first form in subclusters and these subclusters merge to create rich clusters with cool, central X-Ray excesses. The rich cool clusters do not possess ``cooling flows'' due to the presence of bulk velocities in the intracluster medium in excess of 1000 km/sec produced by on-going accretion of gas from supercluster filaments. This new model has several attractive features including the presence of substantial core substructure within the cool cores, and it predicts the appearance of cool bullets, cool fronts, and cool filaments all of which have been recently observed with X-Ray satellites. This hierarchical formation model is also consistent with the observation that cool cores in Abell clusters occur preferentially in dense supercluster environments. On the other hand, our simulations overproduce cool cores in virtually all of our numerical clusters, the central densities are high, and physical core temperatures are often below 1 keV (in contrast to recent observations). We will discuss additional preliminary simulations to ``soften'' the cool cores involving star formation and supernova feedback.

We present results from recent simulations of the formation and evolution of clusters of galaxies in a LambdaCDM cosmology. These simulations contain our most physically complete input physics to date including radiative cooling, star formation that transforms rapidly cooling material into aggregate star particles and we also model the thermal feedback from resulting supernovae in the star particles. We use an adaptive mesh refinement (AMR) Eulerian hydrodynamics scheme to obtain very high spatial resolution (~ 2 kpc) in a computational volume 256 Mpc on a side with mass resolution for dark matter and star particles of ~ 10^8 M_solar. We examine in detail the appearance and evolution of the core region of our simulated clusters.

Measurements of the Hubble constant to distant galaxy clusters using the Sunyaev-Zeldovich effect are systematically low in comparison to values obtained by other means. These measurements usually assume a spherical isothermal $\beta$ model for the intracluster medium (ICM). We present the results of a statistical analysis of {\em temperature bias} in $H_0$ determinations in a sample of 27 numerically simulated X-ray clusters drawn from a $\Lambda$CDM model at z=0.5. We employ adaptive mesh refinement to provide high resolution (15.6 h$^{-1} kpc$) in cluster cores which dominate the X-ray and radio signals. Fitting synthetic X-ray and y-parameter maps to the standard isothermal beta model, we find a broad, skewed distribution in $f\equiv H_0(SZ)/H_0(actual)$ with a mean, median, and standard deviation of 0.89, 0.83 and 0.32 respectively, where $H_0(SZ)$ is the value of $H_0$ derived by using Sunyaev-Zeldovich effect method and $H_0(actual)$ is the value used in the cosmological simulation. We find that the clusters' declining temperature profiles systematically lower estimates of $H_0$ by 10% to 20%. The declining temperature profile of our adiabatic system is consistent with the result including radiative cooling and supernovae feedback (Loken et al. 2002). We thereby introduce a non-isothermal $\beta$ model as an improvement. Applying the non-isothermal $\beta$ model to the refined sample with well-fitted temperature profile, the value of f improves 9% relative to the actual value. The study of the morphology and the clumping effects conclude that these two factors combine to induce scatter in f of $\pm$ 30 %.

In recent years, the possibility of measuring the cosmological constant $\Omega_\Lambda$ through the application of the Alcock-Paczynski test to the Lyman Alpha (Ly$\alpha$) forest has been suggested (McDonald et al. 1999; Hui et al. 1999). Despite the theoretical uncertainties due to a few other cosmological parameters, some of the greatest difficulties we encounter concern the huge uncertainties due to cosmic variance and noise. In this paper, we propose a maximum likelihood estimation (MLE) method to deal with cosmic variance and noise using synthetic spectra of quasistellar objects (QSOs) from our cosmological hydrodynamic simulations. We demonstrate that the MLE method can overcome the cosmic variance problem. Applying the MLE method, we find that we have more than 90% probability to determine $\Omega_\Lambda$ within 20% error and approximately of 66% probability to determine $\Omega_\Lambda$ within 10% error by using 30 pairs QSO spectra when other cosmological parameters are assumed. Another important source of error is from noise in the flux spectra, and we have modeled the corresponding effect by studying artificial spectra with different kinds of noise added. We discover that the noise distribution does not have significant effect on the final cross-correlation functions as long as the signal-to-noise ratio (S/N) is fixed. Finally, a preliminary test and discussion about the sensitivities to other cosmological parameters are included in this paper as well.

We present a new scenario for the formation of cool cores in rich galaxy clusters based on results from recent high spatial dynamic range, adaptive mesh Eulerian hydrodynamic simulations of large-scale structure formation. We find that cores of cool gas, material that would be identified as a classical cooling flow based on its X-ray luminosity excess and temperature profile, are built from the accretion of discrete, stable subclusters. Any ``cooling flow'' present is overwhelmed by the velocity field within the cluster - the bulk flow of gas through the cluster typically has speeds up to about 2,000 km s^-1 and significant rotation is frequently present in the cluster core. The inclusion of consistent initial cosmological conditions for the cluster within its surrounding supercluster environment is crucial when simulating the evolution of cool cores in rich galaxy clusters. This new model for the hierarchical assembly of cool gas naturally explains the high frequency of cool cores in rich galaxy clusters despite the fact that a majority of these clusters show evidence of substructure which is believed to arise from recent merger activity. Furthermore, our simulations generate complex cluster cores in concordance with recent X-ray observations of cool fronts, cool ``bullets'', and filaments in a number of galaxy clusters. Our simulations were computed with a coupled N-body, Eulerian, adaptive mesh refinement, hydrodynamics cosmology code that properly treats the effects of shocks and radiative cooling by the gas. We employ up to seven levels of refinement to attain a peak resolution of 15.6 h^-1 kpc within a volume 256 h-1 Mpc on a side and assume a standard LambdaCDM cosmology.

We present first results from three-dimensional hydrodynamic simulations of the high redshift formation of dwarf galaxies. The simulations use an Eulerian adaptive mesh refinement technique to follow the non-equilibrium chemistry of hydrogen and helium with cosmological initial conditions drawn from a popular Lambda-dominated CDM model. We include the effects of reionization using a uniform radiation field, a phenomenological description of the effect of star formation and, in a separate simulation, the effects of stellar feedback. The results highlight the effects of stellar feedback and photoionization on the baryon content and star formation of galaxies with virial temperatures of approximately 10^4K. Dwarf sized dark matter halos that assemble prior to reionization are able to form stars. Most halos of similar mass that assemble after reionization do not form stars by redshift of three. Dwarf galaxies that form stars show large variations in their gas content because of stellar feedback and photoionization effects. Baryon-to-dark matter mass ratios are found to lie below the cosmic mean as a result of stellar feedback. The supposed substructure problem of CDM is critically assessed on the basis of these results. The star formation histories modulated by radiative and stellar feedbacks are discussed. In addition, metallicities of individual objects are shown to be naturally correlated with their mass-to-light ratios as is also evident in the properties of local dwarf galaxies.

We use three-dimensional hydrodynamic numerical simulations to study phase transformations occurring in a clumpy interstellar gas exposed to time-dependent volumetric heating. To mimic conditions in the Galactic interstellar medium, we take a numerical model of a turbulent multiphase medium from Kritsuk & Norman (2002) computed in a periodic box with mean density n_0=0.25 cm^-3 and mean pressure P_0/k~10^3.4 K cm^-3. A second model with n_0=1 cm^-3 is also considered. Variations of the heating rate on a timescale of 1-10 Myr applied thereafter cause pressure variations in the gas and shifting of the thermal equilibrium curve in the phase plane. This stimulates mass transfer between the gas phases via thermal instability, converting 5-10% of the thermal energy into kinetic energy of gas motions. The experiments demonstrate that recurrent substantial heating episodes can maintain turbulence at this level. Possible applications to the interstellar gas heated by variable far-ultraviolet background radiation produced by short-living massive stars are discussed.

This paper presents a new code for performing multidimensional radiation hydrodynamic (RHD) simulations on parallel computers involving anisotropic radiation fields and nonequilibrium effects. The radiation evolution modules described here encapsulate the physics provided by the serial algorithm of Stone et. al (ApJSupp, vol 80, pp. 819-845), but add new functionality with regard to physics and numerics. Physics enhancments include the addition of time dependence to the computation of the variable tensor Eddington factor (VTEF) closure term, and a matter-radiation coupling scheme which is particularly robust for nonequilibrium problems. Numerical highlights include a discussion of how our code is implemented for parallel execution and a description of our scalable linear solver module. We present a suite of numerical tests from which the virtues and vices of our method may be gleaned; these include nonequilibrium Marshak waves, 2-D "shadow" tests showing the one-sided illumination of an opaque cloud, and full RHD+VTEF simulations of radiating shocks. We conclude that radiation moment solutions closed with variable tensor Eddington factors show a dramatic qualitative improvement over results obtained with flux-limited diffusion, and further that this approach has a bright future in the context of parallel RHD simulations in astrophysics.

We investigate the predicted present-day temperature profiles of the hot, X-ray emitting gas in galaxy clusters for two cosmological models - a current best-guess LCDM model and standard cold dark matter (SCDM). Our numerically-simulated "catalogs" of clusters are derived from high-resolution (15/h kpc) simulations which make use of a sophisticated, Eulerian-based, Adaptive Mesh-Refinement (AMR) code that faithfully captures the shocks which are essential for correctly modelling cluster temperatures. We show that the temperature structure on Mpc-scales is highly complex and non-isothermal. However, the temperature profiles of the simulated LCDM and SCDM clusters are remarkably similar and drop-off as $T +AFw-propto (1+-r/a_x)^{-+AFw-delta}$ where $a_x +AFw-sim r_{vir}/1.5$ and $+AFw-delta +AFw-sim 1.6$. This decrease is in good agreement with the observational results of Markevitch et al.(1998) but diverges, primarily in the innermost regions, from their fit which assumes a polytropic equation of state. Our result is also in good agreement with a recent sample of clusters observed by BeppoSAX though there is some indication of missing physics at small radii ($r<0.2 r_{vir}$). We discuss the interpretation of our results and make predictions for new x-ray observations that will extend to larger radii than previously possible. Finally, we show that, for $r>0.2 r_{vir}$, our universal temperature profile is consistent with our most recent simulations which include both radiative cooling and supernovae feedback.

We study the dynamics of phase transitions in the interstellar medium by means of three-dimensional hydrodynamic numerical simulations. We use a realistic cooling function and generic nonequilibrium initial conditions to follow the formation history of a multiphase medium in detail in the absence of gravity. We outline a number of qualitatively distinct stages of this process, including a linear isobaric evolution, transition to an isochoric regime, formation of filaments and voids (also known as "thermal" pancakes), the development and decay of supersonic turbulence, an approach to pressure equilibrium, and final relaxation of the multiphase medium. We find that 1%-2% of the initial thermal energy is converted into gas motions in one cooling time. The velocity field then randomizes into turbulence that decays on a dynamical timescale E_k ~ t^-n, 1 < n < 2. While not all initial conditions yield a stable two-phase medium, we examine such a case in detail. We find that the two phases are well mixed with the cold clouds possessing a fine-grained structure near our numerical resolution limit. The amount of gas in the intermediate unstable phase roughly tracks the rms turbulent Mach number, peaking at 25% when M_rms ~ 8, decreasing to 11% when M_rms ~ 0.4.

We present new calculations of the inhomogeneous process of cosmological reionization by carefully following the radiative transfer in pre-computed hydrodynamical simulations of galaxy formation. These new computations represent an important step on the way towards fully self-consistent and adaptive calculations which will eventually cover the enormous range of scales from the sizes of individual mini-halos to the mean free path of ionizing photons in the post-overlap Universe. The goal of such simulations is to include enough realistic physics to accurately model the formation of early structures and the end of the `dark ages'. Our new calculations demonstrate that the process by which the ionized regions percolate the Universe is complex, and that the idea of voids being ionized before overdense regions is perhaps too simplistic. It seems that observational information pertaining to the reionization epoch may now be in our grasp, through the detection of Gunn-Peterson troughs at z~6. If so, then the comparison of information from many lines of sight with simulations such as ours may allow us to disentangle details of the ionization history and trace the early formation of structure.

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.

As an entry for the 2001 Gordon Bell Award in the "special" category, we describe our 3-d, hybrid, adaptive mesh refinement (AMR) code, Enzo, designed for high-resolution, multiphysics, cosmological structure formation simulations. Our parallel implementation places no limit on the depth or complexity of the adaptive grid hierarchy, allowing us to achieve unprecedented spatial and temporal dynamic range. We report on a simulation of primordial star formation which develops over 8000 subgrids at 34 levels of refinement to achieve a local refinement of a factor of 10^12 in space and time. This allows us to resolve the properties of the first stars which form in the universe assuming standard physics and a standard cosmological model. Achieving extreme resolution requires the use of 128-bit extended precision arithmetic (EPA) to accurately specify the subgrid positions. We describe our EPA AMR implementation on the IBM SP2 Blue Horizon system at the San Diego Supercomputer Center.

The dynamical state of star-forming molecular clouds cannot be understood without determining the structure and strength of their magnetic fields. Measurements of polarized far-infrared radiation from thermally aligned dust grains are used to map the orientation of the field and estimate its strength, but the accuracy of the results has remained in doubt. In order to assess the reliability of this method, we apply it to simulated far-infrared polarization maps derived from three-dimensional simulations of supersonic magnetohydrodynamical turbulence, and compare the estimated values to the known magnetic field strengths in the simulations. We investigate the effects of limited telescope resolution and self-gravity on the structure of the maps. Limited observational resolution affects the field structure such that small scale variations can be completely suppressed, thus giving the impression of a very homogeneous field. The Chandrasekhar-Fermi method of estimating the mean magnetic field in a turbulent medium is tested, and we suggest an extension to measure the rms field. Both methods yield results within a factor of 2 for field strengths typical of molecular clouds, with the modified version returning more reliable estimates for slightly weaker fields. However, neither method alone works well for very weak fields, missing them by a factor of up to 150. Taking the geometric mean of both methods estimates even the weakest fields accurately within a factor of 2.5. Limited telescope resolution leads to a systematic overestimation of the field strengths for all methods. We discuss the effects responsible for this overestimation and show how to extract information on the underlying (turbulent) power spectrum.

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.

We discuss the results of recent 3D simulations of first structure formation in relationship to the formation of the first stars. On the basis of a new, high-resolution AMR simulation (spatial dynamic range = 30,000,000), we conclude that the first stars are likely to be massive.

We have simulated the evolution of a large sample of X-ray clusters in a Lambda-CDM universe at high spatial resolution using adaptive mesh refinement and placed the results in an online archive for public access. The Simulated Cluster Archive website http://sca.ncsa.uiuc.edu provides tools for interactive 2D and 3D analysis of gas and dark matter fields, X-ray and SZ imaging, and data export. We encourage community use and solicit their feedback.

We describe ZEUS-MP: a Multi-Physics, Massively-Parallel, Message-Passing code for astrophysical fluid dynamics simulations in 3 dimensions. ZEUS-MP is a follow-on to the sequential ZEUS-2D and ZEUS-3D codes developed and disseminated by the Laboratory for Computational Astrophysics (lca.ncsa.uiuc.edu) at NCSA. V1.0 released 1/1/2000 includes the following physics modules: ideal hydrodynamics, ideal MHD, and self-gravity. Future releases will include flux-limited radiation diffusion, thermal heat conduction, two-temperature plasma, and heating and cooling functions. The covariant equations are cast on a moving Eulerian grid with Cartesian, cylindrical, and spherical polar coordinates currently supported. Parallelization is done by domain decomposition and implemented in F77 and MPI. The code is portable across a wide range of platforms from networks of workstations to massively parallel processors. Some parallel performance results are presented as well as an application to turbulent star formation.

We use cosmological hydrodynamic simulations to investigate formation of galactic bulges within the framework of hierarchical clustering in a representative CDM cosmological model. We show that largest objects forming at cosmological redshifts z~4 resemble observed bulges of spiral galaxies or moderate size ellipticals in their general properties like sizes, shapes, and density profiles. This is consistent with observational data indicating the existence of ``old'' bulges and ellipticals at more moderate redshifts. These bulges are gas dominated at redshift z=3, with high rates of star formation and would appear to be good candidates for small blue galaxies seen in the Hubble Deep Field.

Many questions in physical cosmology regarding the thermal history of the intergalactic medium, chemical enrichment, reionization, etc. are thought to be intimately related to the nature and evolution of pregalactic structure. In particular the efficiency of primordial star formation and the primordial IMF are of special interest. We present results from high resolution three--dimensional adaptive mesh refinement simulations that follow the collapse of primordial molecular clouds and their subsequent fragmentation within a cosmologically representative volume. Comoving scales from 128 kpc down to 1 pc are followed accurately. Dark matter dynamics, hydrodynamics and all relevant chemical and radiative processes (cooling) are followed self-consistently for a cluster normalized CDM structure formation model. Primordial molecular clouds with ~10^5 solar masses are assembled by mergers of multiple objects that have formed hydrogen molecules in the gas phase with a fractional abundance of ~10^-4. As the subclumps merge cooling lowers the temperature to ~200 Kelvin in a `cold pocket' at the center of the halo. Within this cold pocket, a quasi-hydrostatically contracting core with ~200 solar mass and number densities > 10^5 cm^-3 is found. We find that less than 1% of the primordial gas in such small scale structures cools and collapses to sufficiently high densities to be available for primordial star formation. Furthermore, it is worthwhile to note that this study achieved the highest dynamic range covered by structured adaptive mesh techniques in cosmological hydrodynamics to date.

New adaptive mesh refinement N-body + hydrodynamics numerical simulations are used to illustrate the complex and changing cluster environments in which many radio galaxies live and evolve. Groups and clusters of galaxies form at the intersections of filaments where they continue to accrete gas and dark matter to the present day. The accretion process produces shocks, turbulence, and transonic bulk flows forming a kind of stormy weather within the intracluster medium (ICM). Radio sources embedded within the stormy ICM form distorted, complex morphologies as observed in recent VLA cluster surveys. We show that the bending of wide-angle tailed radio sources can be understood as the result of recent cluster-subcluster mergers. We use new MHD simulations to illustrate how cluster radio halos can be formed by the shocks and turbulence produced during cluster mergers. Finally, we discuss new observations of distant Abell clusters that reveal a class of weak radio sources, probably starbursts, likely produced during the formation of the clusters as they accrete material from the supercluster environment.

We investigate the properties of the Lyman alpha forest as predicted by numerical simulations for a range of currently viable cosmological models. This is done in order to understand the dependencies of the forest on cosmological parameters. Focusing on the redshift range from two to four, we show that: (1) most of the evolution in the distributions of optical depth, flux and column density can be understood by simple scaling relations, (2) the shape of optical depth distribution is a sensitive probe of the amplitude of density fluctuations on scales of a few hundred kpc, (3) the mean of the b distribution (a measure of the width of the absorption lines) is also very sensitive to fluctuations on these scales, and decreases as they increase. We perform a preliminary comparison to observations, where available. A number of other properties are also examined, including the evolution in the number of lines, the two-point flux distribution and the HeII opacity.

We have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including SPH and grid methods with fixed, deformable or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gas dynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of two. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to about 10%, and the ratio of specific kinetic to thermal energies of the gas agree to about 5%. Various factors contribute to the spread in calculated cluster properties, including differences in the internal timing of the simulations. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modelled.

Hierarchical theories of structure formation predict that clusters of galaxies should be embedded in a web like structure, with filaments emanating from them to large distances. The amount of mass contained within such filaments near a cluster can be comparable to the collapsed mass of the cluster itself. Diffuse infalling material also contains a large amount of mass. Both these components can contribute to the cluster weak lensing signal. This ``projection bias'' is maximized if a filament lies close to the line-of-sight to a cluster. Using large--scale numerical simulations of structure formation in a cosmological constant dominated cold dark matter model, we show that the projected mass typically exceeds the actual mass by several tens of percent. This effect is significant for attempts to estimate cluster masses through weak lensing observations, and will affect weak lensing surveys aimed at constructing the cluster mass function.

The comparison of quasar absorption spectra with numerically simulated spectra from hierarchical cosmological models of structure formation promises to be a valuable tool to discriminate among these models. We present simulation results for the column density, Doppler b parameter, and optical depth probability distributions for five popular cosmological models.

In this paper we critically examine predictions of the Lya forest within the standard cold dark matter (SCDM) model, paying particular attention to the low end of the column-density distribution. We show in particular that the width of these lines, typically measured by the b-parameter of a Voigt profile, is sensitive to spatial resolution in numerical simulations and has previously been overestimated. The new result, which predicts a distribution with a median b of around 20-22 km/s at z=3, is substantially below that observed. We examine a number of possible causes of this discrepancy and argue that it is unlikely to be rectified by an increase in the thermal broadening of the absorbing gas, but is instead telling us something about the distribution of matter on these scales. Although the median differs, the shape of the b-parameter distribution agrees quite well with that observed, and the high-end tail is naturally produced by the filamentary nature of gravitational collapse in these models. In particular, we demonstrate that lines of sight which obliquely intersect a filament or sheet tend to produce absorption lines with larger b parameters. We also examine the physical nature of the gas which is responsible for the forest, showing that for lines with neutral column densities below N_HI ~ 10^14 cm^-2 (for this model at z=3), the peculiar infall velocity is actually slower than the Hubble flow, while larger lines have, on average, turned around and are collapsing.

Many questions in physical cosmology regarding the thermal and ionization history of the intergalactic medium are now successfully studied with the help of cosmological hydrodynamical simulations. Here we present a numerical method that solves the radiative transfer around point sources within a three dimensional cartesian grid. The method is energy conserving independently of resolution: this ensures the correct propagation speeds of ionization fronts. We describe the details of the algorithm, and compute as first numerical application the ionized region surrounding a mini-quasar in a cosmological density field at z=7.

Many questions in physical cosmology regarding the thermal history of the intergalactic medium, chemical enrichment, reionization, etc. are thought to be intimately related to the nature and evolution of pregalactic structure. In particular the efficiency of primordial star formation and the primordial IMF are of special interest. We present results from high resolution three-dimensional adaptive mesh refinement simulations that follow the collapse of primordial molecular clouds and their subsequent fragmentation within a cosmologically representative volume. Comoving scales from 128 kpc down to 0.5 pc are followed accurately. Dark matter dynamics, hydrodynamics and all relevant chemical and radiative processes (cooling) are followed self-consistently for a cluster normalized CDM structure formation model. Primordial molecular clouds with ~100,000 solar masses are assembled by mergers of multiple objects that have formed hydrogen molecules in the gas phase with a fractional abundance of <~ 10^-4. As the subclumps merge cooling decreases the temperature to ~200 Kelvin in multiple "cold pockets" of the merger product. Within these cold pockets, quasi-hydrostatically contracting cores with masses \~ 100 M_sun and number densities > 100,000/cm^3 are found. We find that less than 1% of the primordial gas in such small scale structures cools and collapses to sufficiently high densities to be available for primordial star formation. Furthermore, our results indicate that the formation of very massive objects, massive black holes, fragmentation of a large fraction of baryons into brown dwars or Jupiter size fragments seems, in contrast to various claims in the literature, very unlikely. The expected escape fraction of UV photons with (h nu) > 11eV is very small.

In this paper we review the importance of the Lyman alpha forest as a probe of structure formation in the universe. We first discuss the statistics used to describe the Lyman alpha forest and the numerical techniques used to produce simulated spectra of the forest from a given cosmological model. We then discuss the physical picture of the absorbing structures that emerges from these numerical simulations. Finally, we comment on how two of the statistics, the slope of the column density distribution and the b parameter distribution, may be used to constrain competing cosmologies.

We describe an approach for incorporating radiative transfer into 3D hydrodynamic cosmological simulations. The method, while approximate, allows for a self-consistent treatment of self-shielding and shadowing, diffuse and point sources of radiation, and frequency dependent transfer. Applications include photodissociation, photoheating, and photoionization of the IGM.

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}/\Delta 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.

We report on results of recent, high resolution hydrodynamic simulations of the formation and evolution of X-ray clusters of galaxies carried out within a cosmological framework. We employ the highly accurate piecewise parabolic method (PPM) on fixed and adaptive meshes which allow us to resolve the flow field in the intracluster gas. The excellent shock capturing and low numerical viscosity of PPM represent a substantial advance over previous studies using SPH. We find that in flat, hierarchical cosmological models, the ICM is in a turbulent state long after turbulence generated by the last major merger should have decayed away. Turbulent velocites are found to vary slowly with cluster radius, being $\sim 25%$ of $\sigma_{vir}$ in the core, increasing to $\sim 60%$ at the virial radius. We argue that more frequent minor mergers maintain the high level of turbulence found in the core where dynamical times are short. Turbulent pressure support is thus significant throughout the cluster, and results in a somewhat cooler cluster ($T/T_{vir} \sim .8$) for its mass. Some implications of cluster turbulence are discussed.

We describe an application of Berger and Colella's structured adaptive mesh refinement algorithm to the field of cosmological astrophysics. Simulations in this area must include not just a gaseous component which follows the hyperbolic equations of compressible gas dynamics, but also a collisionless component (such as dark matter or stars) described by the Newtonian dynamical equations. The two fluids interact via gravity which requires an elliptic solver. The challenge for AMR is twofold. First, the collisionless material is most easily modeled by following trajectories of individual objects, a method often referred to as an N-body scheme. The introduction of particles poses a number of difficulties, both physical (how do they interact accurately with the fluid in the mesh) and algorithmic (how to efficiently add a new data structure). The second challenge is to incorporate the gravitational interaction between the two components. We discuss our solutions to both of these issues and briefly present very encouraging results.

Gravitational instabilities naturally give rise to multi-scale structure, which is difficult for traditional Eulerian hydrodynamic methods to accurately evolve. This can be circumvented by adaptively adding resolution (in the form of multiple levels of finer meshes) to relatively small volumes as required. We describe an application of this adaptive mesh refinement (AMR) technique to cosmology, focusing on the formation and evolution of X-ray clusters. A set of simulations are performed on a single cluster, varying the initial resolution and refinement criteria. We find that although new, small scale structure continues to appear as the resolution is increased, bulk properties and radial profiles appear to converge at an effective resolution of $8192^3$.

We compare the results of Eulerian hydrodynamic simulations of cluster formation against virial scaling relations between four bulk quantities: the cluster mass, the dark matter velocity dispersion, the gas temperature and the cluster luminosity. The comparison is made for a large number of clusters at a range of redshifts in three different cosmological models (CHDM, CDM and OCDM). We find that the analytic formulae provide a good description of the relations between three of the four numerical quantities. The fourth (luminosity) also agrees once we introduce a procedure to correct for the fixed numerical resolution. We also compute the normalizations for the virial relations and compare extensively to the existing literature, finding remarkably good agreement. The Press-Schechter prescription is calibrated with the simulations, again finding results consistent with other authors. We also examine related issues such as the size of the scatter in the virial relations, the effect of metallicity with a fixed pass-band, and the structure of the halos. All of this is done in order to establish a firm groundwork for the use of clusters as cosmological probes. Implications for the models are briefly discussed.

We examine the structure of baryonic wakes in the cosmological fluid which would form behind GUT-scale cosmic strings at early times (redshifts z > 100) in a neutrino-dominated universe. We show, using simple analytical arguments as well as 1- and 2-dimensional hydrodynamical simulations, that these wakes will NOT be able to form interesting cosmological objects before the neutrino component collapses. The width of the baryonic wakes (< 10 kpc comoving) is smaller than the scale of wiggles on the strings and are probably not enhanced by the wiggliness of the string network.

We discuss the origin and physical nature of the Lya forest absorption systems as found in hydrodynamical simulations of the Intergalactic Medium in a standard Cold Dark Matter cosmology (Omega=1, H0=50 km/s/Mpc, sigma_8=0.7). The structures of the systems that give rise to the Lya forest span a wide range in morphologies, depending on the density contrast. The absorption properties of the forest clouds correlate strongly with those of the underlying physical systems from which they arise. The highest column density systems (log N_HI > 15), correspond to highly overdense spheroidal structures, moderate column density systems (13 < log N_HI < 14), correspond to filaments, and the lowest density absorption systems originate from discrete fluctuations within underdense regions a few megaparsecs across, cosmic minivoids. Most of the intergalactic HeII opacity arises from these underdense regions. We are able to account for the distribution of optical depths in our synthesized spectra entirely by absorption due to discrete systems. We find that virtually all the baryons in the simulation fragment into structures that we can identify with discrete absorption lines, with at most 5% remaining in a smoothly distributed component (the Gunn-Peterson effect). We compare our results with the cloud ionization parameters inferred from Keck HIRES measurements of carbon and silicon in the Lya forest. Combining with constraints imposed by measurements of the mean intergalactic HI opacity, we find 0.03 < Omega_b < 0.08 and 0.3 < (Gamma_HI/ 10^-12 s^-1) < 1 at z=3-3.5.

We investigate the possibility of very early formation of primordial star clusters from high-\sigma perturbations in cold dark matter dominated structure formation scenarios. For this we have developed a powerful 2-level hierarchical cosmological code with a realistic and robust treatment of multi-species primordial gas chemistry, paying special attention to the formation and destruction of hydrogen molecules, non-equilibrium ionization, and cooling processes. We performed 3-D simulations at small scales and at high redshifts and find that, analogous to simulations of large scale structure, a complex system of filaments, sheets, and spherical knots at the intersections of filaments form. On the mass scales covered by our simulations (5x10^5 - 1x10^9\Ms) that collapse at redshifts z>25, we find that only at the spherical knots can enough H2 be formed (n_{H_2}/n_H > 5x10^-4) to cool the gas appreciably. Quantities such as the time dependence of the formation of H2 molecules, the final H2 fraction, and central densities from the simulations are compared to the theoretical predictions of Abel (1995) and Tegmark et al. (1997) and found to agree remarkably well. Comparing the 3-D results to an isobaric collapse model we further discuss the possible implications of the extensive merging of small structure that is inherent in hierarchical models. Typically only 5-8% percent of the total baryonic mass in the collapsing structures is found to cool significanlty. Assuming the Padoan (1995) model for star formation our results would predict the first stellar systems to be as small as ~30\Ms. Some implications for primordial globular cluster formation scenarios are also discussed.

We simulate the Lya forest in a standard CDM universe using a 2-level hierarchical grid code to evolve the dark and baryonic matter components self-consistently. We solve the time-dependent ionization equations for hydrogen and helium, adopting the Haardt & Madau (1996) estimate for the metagalactic UV radiation background. We compare our simulation results with the measured properties of the Lya forest by constructing synthetic spectra which we analyze in a manner similar to the method used by observers. The HI column density and Doppler parameter distributions we obtain agree closely with those measured by the Keck HIRES and earlier high spectral resolution observations over the column density range 10^{12} < NHI < 10^{16}. In particular, we find a power-law column density distribution persists to the lowest values reported from the HIRES measurements, in agreement with the incompleteness-corrected observations. We find evolution in the cloud number density and opacity comparable to the observed evolution. We also compare with measured values of the intergalactic HeII opacity. Our results require a HeII ionizing background lower than the Haardt & Madau estimate by a factor of 4, corresponding to a soft intrinsic QSO spectrum of alpha_Q = 1.8 - 2.

We developed a three-dimensional 2-level hierarchical cosmological code with a realistic and robust treatment of multi-species non-equilibrium ionization and radiative cooling processes, and use it to investigate primordial star forming regions that originate from high-\sigma perturbations in a standard CDM dominated cosmology. We find it is possible to produce gravitationally bound and cooled structures at very high redshift (z ~ 40) with baryonic masses as small as ~1000Ms. The molecular hydrogen formation in these small scale structures follows very well the analytical predictions of Abel (1995) and Tegmark et al. (1996). We also discuss the minimum mass that cosmological structures must have in order to be able to cool and collapse.

We have reviewed the chemistry and cooling behaviour of low-density (n<10^4 cm^-3) primordial gas and devised a cooling model wich involves 19 collisional and 9 radiative processes and is applicable for temperatures in the range (1 K < T < 10^8 K). We derived new fits of rate coefficients for the photo-attachment of neutral hydrogen, the formation of molecular hydrogen via H-, charge exchange between H2 and H+, electron detachment of H- by neutral hydrogen, dissociative recombination of H2 with slow electrons, photodissociation of H2+, and photodissociation of H2. Further it was found that the molecular hydrogen produced through the gas-phase processes, H2+ + H -> H2 + H+, and H- + H -> H2 + e-, is likely to be converted into its para configuration on a faster time scale than the formation time scale. We have tested the model extensively and shown it to agree well with former studies. We further studied the chemical kinetics in great detail and devised a minimal model which is substantially simpler than the full reaction network but predicts correct abundances. This minimal model shows convincingly that 12 collisional processes are sufficient to model the H, He, H+, H-, He+, He++, and H2 abundances in low density primordial gas for applications with no radiation fields.

We have developed a method of solving for multi-species chemical reaction flows in non--equilibrium and self--consistently with the hydrodynamic equations in an expanding FLRW universe. The method is based on a backward differencing scheme for the required stability when solving stiff sets of equations and is designed to be efficient for three-dimensional calculations without sacrificing accuracy. In all, 28 kinetic reactions are solved including both collisional and radiative processes for the following nine separate species: H, H+, He, He+, He++, H-, H2+, H2, and e-. The method identifies those reactions (involving H- and H2+) ocurring on the shortest time scales, decoupling them from the rest of the network and imposing equilibrium concentrations to good accuracy over typical cosmological dynamical times. Several tests of our code are presented, including radiative shock waves, cosmological sheets, conservation constraints, and fully three-dimensional simulations of CDM cosmological evolutions in which we compare our method to results obtained when the packaged routine LSODAR is substituted for our algorithms.

We present the first X-ray Luminosity Function (XLF) for an optically-selected sample of 49 nearby poor clusters of galaxies and a sample of 67 Abell clusters with z < 0.15. We have extended the measured cluster XLF by more than a factor of 10 in X-ray luminosity. Our poor cluster sample was drawn from an optical catalog of groups with 0.01<z<0.03 composed of Zwicky galaxies. The X-ray emission was measured from the ROSAT all-sky survey. About 45% of the poor clusters were detected with 0.5-2.0 keV luminosities from (1.7-65)x 10^{41} h^{-2} ergs/sec. These are among the X-ray brightest, optically-selected poor clusters in the northern hemisphere. For this sample, the poor cluster XLF was found to be smooth extrapolation of the rich cluster XLF. A new Hydro/N-body simulation of a Hot +Cold dark matter model with \Omega_{total}=1, \Omega_{\nu}=0.2, and a baryon fraction of 7.5% was used to model and understand our observational selection effects. We found that the observed cluster Gas Mass Function was consistent with our model.

We discover a new dynamical mechanism that significantly enhances the growth of Rayleigh-Taylor fingers developed near the contact interface between the supernova ejecta and swept-up ambient gas in young supernova remnants if the supernova remnant expands into a clumpy (cloudy) circumstellar medium. Our numerical simulation demonstrates that large Rayleigh-Taylor fingers can obtain a sufficient terminal velocity to protrude through the forward shock front by taking extra kinetic energy from vorticies generated by shock-cloud interactions. We suggest this mechanism as a means to generate the aspherical expansion of the supernova ejecta. Ambient magnetic fields are stretched and amplified as the Rayleigh-Taylor fingers protrude, possibly leading to strongly enhanced radio emission. The material in the protrusions originates from the ejected stellar material with greatly enhanced heavy elements. Therefore, it can be a strong X-ray emitter. The timescale for the Rayleigh-Taylor fingers to reach the forward shock depends on the size, mass density and distribution of clouds being engulfed by the supernova shock, although the details will require further numerical investigation.

We study the radio emission from young supernova remnants by means of 3D numerical MHD simulations of the Rayleigh-Taylor instability in the shell of the remnant. The computation is carried out in spherical polar coordinates ($r, \theta, \phi$) by using a moving grid technique which allows us to finely resolve the shell. Three-dimensional result shows more turbulent (complex) structures in the mixing region than the two-dimensional result, and the instability is found to deform the reverse shock front. Stokes parameters (I,Q, and U) are computed to study the radio properties of the remnant. The total intensity map shows two distinctive regions (inner and outer shells). The inner shell appears to be complex and turbulent exhibiting loop structures and plumes as a result of the Rayleigh-Taylor instability, while the outer shell is faint and laminar due to the shocked uniform ambient magnetic fields. The inner shell resembles the observed radio structure in the main shell of young SNRs, which is evidence that the Rayleigh-Taylor instability is an ongoing process in young SNRs. When only the peculiar components of the magnetic fields generated by the instability are considered, the polarization B-vector in the inner radio shell is preferentially radial with about $20 \sim 50\%$ of fractional polarization which is higher than the observed value. The fractional polarization is lowest in the turbulent inner shell and increases outward, which is attributed to the geometric effect. The polarized intensity is found to be correlated with the total intensity. We demonstrate that the polarized intensity from the turbulent region can dominate over the polarized intensity from the shocked uniform fields if the amplified field is sufficiently strong. Therefore, we conclude that the Rayleigh-Taylor instability can explain the dominant radial magnetic

We perform a double line of sight (DLOS) analysis of the Lyman-alpha forest structures that form and evolve in cosmological N-body/hydrodynamic simulations. Pairs of simulated spectra, extracted from lines of sight separated by distances from D=12.5kpc up to 800kpc, and a ``control sample'' of unrelated lines of sight, are analyzed at redshifts 3, 2, and 1. Coincident line samples are defined for HI column density thresholds of Nco = 10^{12.5}, 10^{13}, and 10^{14} per square cm. We find that: 1) Under the assumption of a single structure size, a Bayesian analysis yields sizes that are larger for smaller Nco, and at fixed Nco the size decreases with decreasing redshift. However, these derived sizes are found to increase with increasing D indicating that the assumption of a single structure size is invalid. 2) The column densities of coincident pairs are highly correlated for small D, with increasing scatter as D is increased, consistent with structures that have a centrally peaked N(HI) that decreases gradually with radius. 3) The velocity difference distribution for coincident lines is very narrow for small D, and widens as D is increased to meet the expectation for chance coincidences in unrelated lines of sight. This behavior is indicative of organized motion within the structures. 4) For small D, the distribution of anticoincident line column densities, Nac, falls steeply as Nac increases from the cutoff value, but has a significant tail at large values which is inconsistent with a population of spherical absorbers with sharp edges, and consistent with a flattened geometry. The conclusions reached on the basis of the DLOS analysis are validated by an examination of the three-dimensional structures and velocity flows in the simulation data.

Soon after the first QSO was identified, Gunn & Peterson searched for the expected characteristic absorption trough on the blueward side of Lya in the spectrum of the QSO due to an Intergalactic Medium (IGM). They failed to find it, placing a constraint on the density of neutral hydrogen in the IGM that was less than 1 part in 10^5 of that residing in galaxies, and concluded that the IGM must be highly ionized. Soon afterwards absorption by Lya was detected in the IGM; not by a diffuse component, but by a clumpy component of intergalactic gas clouds, the Lya forest. The goal then became to search for `excess' absorption beyond that expected from the forest. No clear absorption by a diffuse HI component, however, was ever detected. The results of recent numerical hydrodynamics computations of the formation of the Lya forest appear to indicate that the division between a clumpy component and a diffuse one may be inappropriate. We find that in a CDM-dominated cosmology, no substantial diffuse medium should be expected. Instead, the medium condenses into a network of complex structures that reveal themselves as discrete absorption systems in the spectra of QSOs. The lowest column density lines arise from the fine structure in minivoids -- small regions with densities below the cosmic mean. These results suggest that the long sought for diffuse HI Gunn-Peterson effect may not exist.

Young supernova remnants such as Tycho generally exhibit a bright circular clumpy shell in both radio and X-ray emission. For several young remnants, various arguments suggest that the magnetic field is larger than can be explained by compression of a few $\mu G$ ambient magnetic field by the shock wave. Radio polarization studies reveal a net radial orientation of magnetic fields in the shell which cannot be explained by the simple compression either. We model Rayleigh-Taylor instability at the interface of the ejecta and the shocked ambient medium to explain these observations. We have performed multidimensional MHD simulations of the instability in the shell of a Type-I supernova remnant for the first time utilizing a moving grid technique which allows us to follow the growth of the instability and its effect on the local magnetic field in detail. We find that the evolution of the instability is very sensitive to the deceleration of the ejecta and the evolutionary stage of the remnant. As the reverse shock enters the inner uniform density region, Chevalier's self-similar stage ends and the thickness of radio shell increases and the instability weakens. Our simulation shows that Rayleigh-Taylor and Kelvin-Helmholtz instabilities amplify ambient magnetic fields locally by as much as a factor of 60 around dense fingers due to stretching, winding, and compression. Globally, the amount of magnetic-field amplification is low and the magnetic energy density reaches only about 0.3% of the turbulent energy density at the end of simulation. Strong magnetic field lines draped around the fingers produce the radial B-vector polarization, whereas thermal bremsstrahlung from the dense fingers themselves produce the clumpy X-ray emission. As a result, the X-ray emission peaks inside of the radio emission.

We have performed a multi-species hydrodynamical simulation of the formation and evolution of Lyman alpha clouds in a flat CDM dominated universe with an external flux of ionizing radiation. We solve the fully coupled non-equilibrium rate equations for the following species: H, H^+, H^-, H_2, H_2^+, He, He^+, He^{++}, and e^-. The statistical properties of the distribution and evolution of both hydrogen and helium absorption lines are extracted and compared to observed data. We find excellent agreement for the following neutral hydrogen data: the distribution of column densities is fit well by a power law with exponent beta=1.55 with a possible deficiency of lines above column density 10^{15} cm^{-2}; the integrated distribution matches observed data over a broad range of column densities 10^{13} to 10^{17} cm^{-2}; a Gaussian statistical fit to the Doppler parameter distribution yields a median of 35.6 km s^{-1}; the evolution of the number of clouds with column densities larger than 10^{14} cm^{-2} follows a power law with exponent gamma=2.22. Analogous calculations are presented for HeII absorption lines and we find the ratio of Doppler parameters b_{HeII}/b_{HI} = 0.87. Our data also suggests that Ly$\alpha$ clouds may belong to two morphologically different groups: small column density clouds which tend to reside in sheets or filamental structures and are very elongated and/or flattened, and the large column density clouds which are typically found at the intersections of these caustic structures and are quasi-spherical.

We present the first self-consistently computed results on the evolution of X-ray properties of galaxy clusters in a Cold + Hot Dark Matter (CHDM) model. We have performed a hydrodynamic plus N-body simulation for the COBE-compatible CHDM model with standard mass components: Omega(hot) = 0.3, Omega(cold) = 0.6 and Omega(baryon) = 0.1 (h = 0.5). In contrast with the CDM model, which fails to reproduce the observed temperature distribution function dN/dT (Bryan et al. 1994b), the CHDM model fits the observational dN/dT quite well. Our results on X-ray luminosity are less firm but even more intriguing. We find that the resulting X-ray luminosity functions at redshifts z = 0.0, 0.2, 0.4, 0.7 are well fit by observations, where they overlap. The fact that both temperatures and luminosities provide a reasonable fit to the available observational data indicates that, unless we are missing some essential physics, there is neither room nor need for a large fraction of gas in rich clusters: 10% (or less) in baryons is sufficient to explain their X-ray properties. We also see a tight correlation between X-ray luminosity and gas temperature.

We present a detailed comparison of the simulation results of various cosmological hydrodynamic codes. Starting with identical initial conditions based on the Cold Dark Matter scenario for the growth of structure, we integrate from redshift $z=20$ to $z=0$ to determine the physical state within a representative volume of size $L^3$ where $L=64 h^{-1} {\rm Mpc}$. Five independent codes are compared: three of them Eulerian mesh based and two variants of the Smooth Particle Hydrodynamics "SPH" Lagrangian approach. The Eulerian codes were run at $N^3=(32^3,~64^3,~128^3,~{\rm and},~256^3)$ cells, the SPH codes at $N^3= 32^3$ and $64^3$ particles. Results were then rebinned to a $16^3$ grid with the expectation that the rebinned data should converge, by all techniques, to a common and correct result as $N \rightarrow \infty$. We find that global averages of various physical quantities do, as expected, tend to converge in the rebinned model, but that uncertainties in even primitive quantities such as $\langle T \rangle$, $\langle \rho^2\rangle^{1/2}$ persists at the 3\%-17\% level after completion of very large simulations. The two SPH codes and the two shock capturing Eulerian codes achieve comparable and satisfactory accuracy for comparable computer time in their treatment of the high density, high temperature regions as measured in the rebinned data; the variance among the five codes (at highest resolution) for the mean temperature (as weighted by $\rho^2$) is only 4.5\%. Overall the comparison allows us to better estimate errors, it points to ways of improving this current generation of hydrodynamic codes and of suiting their use to problems which exploit their individually best features.