We explore the use of field solvers as approximations of classical Vlasov-Poisson systems. This correspondence is investigated in both electrostatic and gravitational contexts. We demonstrate the ability of field solvers to be excellent approximations of problems with cold initial condition into the non linear regime. We also investigate extensions of the Schr\"odinger-Poisson system that employ multiple stacked cold streams, and the von Neumann-Poisson equation as methods that can successfully reproduce the classical evolution of warm initial conditions. We then discuss how appropriate simulation parameters need to be chosen to avoid interference terms, aliasing, and wave behavior in the field solver solutions. We present a series of criteria clarifying how parameters need to be chosen in order to effectively approximate classical solutions.

As computational resolution of modern cosmological simulations reach ever so close to resolving individual star-forming clumps in a galaxy, a need for "resolution-appropriate" physics for a galaxy-scale simulation has never been greater. To this end, we introduce a self-consistent numerical framework that includes explicit treatments of feedback from star-forming molecular clouds (SFMCs) and massive black holes (MBHs). In addition to the thermal supernovae feedback from SFMC particles, photoionizing radiation from both SFMCs and MBHs is tracked through full 3-dimensional ray tracing. A mechanical feedback channel from MBHs is also considered. Using our framework, we perform a state-of-the-art cosmological simulation of a quasar-host galaxy at z~7.5 for ~25 Myrs with all relevant galactic components such as dark matter, gas, SFMCs, and an embedded MBH seed of ~> 1e6 Ms. We find that feedback from SFMCs and an accreting MBH suppresses runaway star formation locally in the galactic core region. Newly included radiation feedback from SFMCs, combined with feedback from the MBH, helps the MBH grow faster by retaining gas that eventually accretes on to the MBH. Our experiment demonstrates that previously undiscussed types of interplay between gas, SFMCs, and a MBH may hold important clues about the growth and feedback of quasars and their host galaxies in the high-redshift Universe.

We present BAM: a novel Bias Assignment Method envisaged to generate mock catalogs. Combining the statistics of dark matter tracers from a high resolution cosmological $N$-body simulation and the dark matter density field calculated from down-sampled initial conditions using efficient structure formation solvers, we extract the halo-bias relation on a mesh of a $3\,h^{-1}$ Mpc cell side resolution as a function of properties of the dark matter density field (e.g. local density, cosmic web type), automatically including stochastic, deterministic, local and non-local components. We use this information to sample the halo density field, accounting for ignored dependencies through an iterative process. By construction, our approach reaches $\sim 1\%$ accuracy in the majority of the $k$-range up to the Nyquist frequency without systematic deviations for power spectra (about $k \sim 1\, h$ Mpc$^{-1}$) using either particle mesh or Lagrangian perturbation theory based solvers. When using phase-space mapping to compensate the low resolution of the approximate gravity solvers, our method reproduces the bispectra of the reference within $10\%$ precision studying configurations tracing the quasi-nonlinear regime. BAM has the potential to become a standard technique to produce mock halo and galaxy catalogs for future galaxy surveys and cosmological studies being highly accurate, efficient and parameter free.

In hierarchical structure formation, merging of galaxies is frequent and known to dramatically affect their properties. To comprehend these interactions high-resolution simulations are indispensable because of the nonlinear coupling between pc and Mpc scales. To this end, we present the first adaptive mesh refinement (AMR) simulation of two merging, low mass, initially gas-rich galaxies (1.8e10 Ms each), including star formation and feedback. With galaxies resolved by ~2e7 total computational elements, we achieve unprecedented resolution of the multiphase interstellar medium, finding a widespread starburst in the merging galaxies via shock-induced star formation. The high dynamic range of AMR also allows us to follow the interplay between the galaxies and their embedding medium depicting how galactic outflows and a hot metal-rich halo form. These results demonstrate that AMR provides a powerful tool in understanding interacting galaxies.

Using an isolated Milky Way-mass galaxy simulation, we compare results from 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical community. We utilize the infrastructure we have built for the AGORA High-resolution Galaxy Simulations Comparison Project. This includes the common disk initial conditions, common physics models (e.g., radiative cooling and UV background by the standardized package Grackle) and common analysis toolkit yt, all of which are publicly available. Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production are carefully constrained across code platforms. With numerical accuracy that resolves the disk scale height, we find that the codes overall agree well with one another in many dimensions including: gas and stellar surface densities, rotation curves, velocity dispersions, density and temperature distribution functions, disk vertical heights, stellar clumps, star formation rates, and Kennicutt-Schmidt relations. Quantities such as velocity dispersions are very robust (agreement within a few tens of percent at all radii) while measures like newly-formed stellar clump mass functions show more significant variation (difference by up to a factor of ~3). Systematic differences exist, for example, between mesh-based and particle-based codes in the low density region, and between more diffusive and less diffusive schemes in the high density tail of the density distribution. Yet intrinsic code differences are generally small compared to the variations in numerical implementations of the common subgrid physics such as supernova feedback. Our experiment reassures that, if adequately designed in accordance with our proposed common parameters, results of a modern high-resolution galaxy formation simulation are more sensitive to input physics than to intrinsic differences in numerical schemes.

We investigate the spatially-resolved star formation relation using a galactic disk formed in a comprehensive high-resolution (3.8 pc) simulation. Our new implementation of stellar feedback includes ionizing radiation as well as supernova explosions, and we handle ionizing radiation by solving the radiative transfer equation rather than by a subgrid model. Photoheating by stellar radiation stabilizes gas against Jeans fragmentation, reducing the star formation rate. Because we have self-consistently calculated the location of ionized gas, we are able to make spatially-resolved mock observations of star formation tracers, such as H-alpha emission. We can also observe how stellar feedback manifests itself in the correlation between ionized and molecular gas. Applying our techniques to the disk in a galactic halo of 2.3e11 Msun, we find that the correlation between star formation rate density (estimated from mock H-alpha emission) and molecular hydrogen density shows large scatter, especially at high resolutions of <~ 75 pc that are comparable to the size of giant molecular clouds (GMCs). This is because an aperture of GMC size captures only particular stages of GMC evolution, and because H-alpha traces hot gas around star-forming regions and is displaced from the molecular hydrogen peaks themselves. By examining the evolving environment around star clusters, we speculate that the breakdown of the traditional star formation laws of the Kennicutt-Schmidt type at small scales is further aided by a combination of stars drifting from their birthplaces, and molecular clouds being dispersed via stellar feedback.

We describe a new method for simulating ionizing radiation and supernova feedback in the analogues of low-redshift galactic disks. In this method, which we call star-forming molecular cloud (SFMC) particles, we use a ray-tracing technique to solve the radiative transfer equation for ultraviolet photons emitted by thousands of distinct particles on the fly. Joined with high numerical resolution of 3.8 pc, the realistic description of stellar feedback helps to self-regulate star formation. This new feedback scheme also enables us to study the escape of ionizing photons from star-forming clumps and from a galaxy, and to examine the evolving environment of star-forming gas clumps. By simulating a galactic disk in a halo of 2.3e11 Msun, we find that the average escape fraction from all radiating sources on the spiral arms (excluding the central 2.5 kpc) fluctuates between 0.08% and 5.9% during a ~20 Myr period with a mean value of 1.1%. The flux of escaped photons from these sources is not strongly beamed, but manifests a large opening angle of more than 60 degree from the galactic pole. Further, we investigate the escape fraction per SFMC particle, f_esc(i), and how it evolves as the particle ages. We discover that the average escape fraction f_esc is dominated by a small number of SFMC particles with high f_esc(i). On average, the escape fraction from a SFMC particle rises from 0.27% at its birth to 2.1% at the end of a particle lifetime, 6 Myrs. This is because SFMC particles drift away from the dense gas clumps in which they were born, and because the gas around the star-forming clumps is dispersed by ionizing radiation and supernova feedback. The framework established in this study brings deeper insight into the physics of photon escape fraction from an individual star-forming clump, and from a galactic disk.

There is mounting evidence for the coevolution of galaxies and their embedded massive black holes (MBHs) in a hierarchical structure formation paradigm. To tackle the nonlinear processes of galaxy-MBH interaction, we describe a self-consistent numerical framework which incorporates both galaxies and MBHs. The high-resolution adaptive mesh refinement (AMR) code Enzo is modified to model the formation and feedback of molecular clouds at their characteristic scale of 15.2 pc and the accretion of gas onto a MBH. Two major channels of MBH feedback, radiative feedback (X-ray photons followed through full 3D adaptive ray tracing) and mechanical feedback (bipolar jets resolved in high-resolution AMR), are employed. We investigate the coevolution of a 9.2e11 Msun galactic halo and its 1e5 Msun embedded MBH at redshift 3 in a cosmological LCDM simulation. The MBH feedback heats the surrounding ISM up to 1e6 K through photoionization and Compton heating and locally suppresses star formation in the galactic inner core. The feedback considerably changes the stellar distribution there. This new channel of feedback from a slowly growing MBH is particularly interesting because it is only locally dominant, and does not require the heating of gas globally on the disk. The MBH also self-regulates its growth by keeping the surrounding ISM hot for an extended period of time.

In the standard approach to studying cosmological structure formation, the overall expansion of the Universe is assumed to be homogeneous, with the gravitational effect of inhomogeneities encoded entirely in a Newtonian potential. A topic of ongoing debate is to what degree this fully captures the dynamics dictated by general relativity, especially in the era of precision cosmology. To quantitatively assess this, we directly compare standard N-body Newtonian calculations to full numerical solutions of the Einstein equations, for cold matter with various magnitude initial inhomogeneities on scales comparable to the Hubble horizon. We analyze the differences in the evolution of density, luminosity distance, and other quantities defined with respect to fiducial observers. This is carried out by reconstructing the effective spacetime and matter fields dictated by the Newtonian quantities, and by taking care to distinguish effects of numerical resolution. We find that the fully general relativistic and Newtonian calculations show excellent agreement, even well into the nonlinear regime. They only notably differ in regions where the weak gravity assumption breaks down, which arise when considering extreme cases with perturbations exceeding standard values.

We present a new method for generating initial conditions for numerical cosmological simulations in which massive neutrinos are treated as an extra set of N-body (collisionless) particles. It allows us to accurately follow the density field for both Cold Dark Matter (CDM) and neutrinos at both high and low redshifts. At high redshifts, the new method is able to reduce the shot noise in the neutrino power spectrum by a factor of more than $10^7$ compared to previous methods, where the power spectrum was dominated by shot noise at all scales. We find that our new approach also helps to reduce the noise on the total matter power spectrum on large scales, whereas on small scales the results agree with previous simulations. Our new method also allows for a systematic study of clustering of the low velocity tail of the distribution function of neutrinos. This method also allows for the study of the evolution of the overall velocity distribution as a function of the environment determined by the CDM field.

Primordial stars are likely to be very massive >30 Msun, form in isolation, and will likely leave black holes as remnants in the centers of their host dark matter halos. We expect primordial stars to form in halos in the mass range 10^6-10^10 Msun. Some of these early black holes, formed at redshifts z>10, could be the seed black hole for a significant fraction of the supermassive black holes found in galaxies in the local universe. If the black hole descendants of the primordial stars exist, their mergers with nearby supermassive black holes may be a prime candidate for long wavelength gravitational wave detectors. We simulate formation and evolution of dark matter halos in LambdaCDM universe. We seed high-redshift dark matter halos with early black holes, and explore the merger history of the host halos and the implications of black hole's kick velocities arising from their coalescence. The central concentration of low mass early black holes in present day galaxies is reduced if they experience even moderate kicks of tens of km/s. Even such modest kicks allow the black holes to leave their parent halo, which consequently leads to dynamical friction being less effective on the low mass black holes that were ejected, compared to those still embedded in their parent halos. Therefore, merger rates with central supermassive black holes in the largest halos may be reduced by more than an order of magnitude. Using analytical and illustrative cosmological N-body simulations, we quantify the role of kicks on the merger rates of black holes formed from massive metal free stars with supermassive black holes in present day galaxies.

Indirect detection of dark matter is a major avenue for discovery. However, baryonic backgrounds are diverse enough to mimic many possible signatures of dark matter. In this work, we study the newly proposed technique of dark matter velocity spectroscopy\,\cite{speckhard2016}. The non-rotating dark matter halo and the Solar motion produce a distinct longitudinal dependence of the signal which is opposite in direction to that produced by baryons. Using collisionless dark matter only simulations of Milky Way like halos, we show that this new signature is robust and holds great promise. We develop mock observations by high energy resolution X-ray spectrometer on a sounding rocket, the Micro-X experiment, to our test case, the 3.5 keV line. We show that by using six different pointings, Micro-X can exclude a constant line energy over various longitudes at $\geq$ 3$\sigma$. The halo triaxiality is an important effect and it will typically reduce the significance of this signal. We emphasize that this new {\it smoking gun in motion} signature of dark matter is general, and is applicable to any dark matter candidate which produces a sharp photon feature in annihilation or decay.

We study the magnetic field evolution of an isolated spiral galaxy, using isolated Milky Way-mass galaxy formation simulations and a novel prescription for magnetohydrodynamic (MHD) supernova feedback. Our main result is that a galactic dynamo can be seeded and driven by supernova explosions, resulting in magnetic fields whose strength and morphology is consistent with observations. In our model, supernovae supply thermal energy, and a low level magnetic field along with their ejecta. The thermal expansion drives turbulence, which serves a dual role by efficiently mixing the magnetic field into the interstellar medium, and amplifying it by means of turbulent dynamo. The computational prescription for MHD supernova feedback has been implemented within the publicly available ENZO code, and is fully described in this paper. This improves upon ENZO's existing modules for hydrodynamic feedback from stars and active galaxies. We find that the field attains $\mu G$-levels over Gyr-time scales throughout the disk. The field also develops large-scale structure, which appears to be correlated with the disk's spiral arm density structure. We find that seeding of the galactic dynamo by supernova ejecta predicts a persistent correlation between gas metallicity and magnetic field strength. We also generate all-sky maps of the Faraday rotation measure from the simulation-predicted magnetic field, and present a direct comparison with observations.

The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web -- depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper we bring twelve of these methods together and apply them to the same data set in order to understand how they compare. In general these cosmic web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore one would not {\it a priori} expect agreement between different techniques however, many of these methods do converge on the identification of specific features. In this paper we study the agreements and disparities of the different methods. For example, each method finds that knots inhabit higher density regions than filaments, etc. and that voids have the lowest densities. For a given web environment, we find substantial overlap in the density range assigned by each web classification scheme. We also compare classifications on a halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass haloes (i.e. $M_{\rm halo}\sim10^{13.5}h^{-1}M_{\odot}$) as being in filaments. Lastly, so that any future cosmic web classification scheme can be compared to the 12 methods used here, we have made all the data used in this paper public.

We present the Grackle chemistry and cooling library for astrophysical simulations and models. Grackle provides a treatment of non-equilibrium primordial chemistry and cooling for H, D, and He species, including H2 formation on dust grains; tabulated primordial and metal cooling; multiple UV background models; and support for radiation transfer and arbitrary heat sources. The library has an easily implementable interface for simulation codes written in C, C++, and Fortran as well as a Python interface with added convenience functions for semi-analytical models. As an open-source project, Grackle provides a community resource for accessing and disseminating astrochemical data and numerical methods. We present the full details of the core functionality, the simulation and Python interfaces, testing infrastructure, performance, and range of applicability. Grackle is a fully open-source project and new contributions are welcome.

We study evolution of voids in cosmological simulations using a new method for tracing voids over cosmic time. The method is based on tracking watershed basins (contiguous regions around density minima) of well developed voids at low redshift, on a regular grid of density field. It enables us to construct a robust and continuous mapping between voids at different redshifts, from initial conditions to the present time. We discuss how the new approach eliminates strong spurious effects of numerical origin when voids evolution is traced by matching voids between successive snapshots (by analogy to halo merger trees). We apply the new method to a cosmological simulation of a standard LambdaCDM cosmological model and study evolution of basic properties of typical voids (with effective radii between 6Mpc/h and 20Mpc/h at redshift z=0) such as volumes, shapes, matter density distributions and relative alignments. The final voids at low redshifts appear to retain a significant part of the configuration acquired in initial conditions. Shapes of voids evolve in a collective way which barely modifies the overall distribution of the axial ratios. The evolution appears to have a weak impact on mutual alignments of voids implying that the present state is in large part set up by the primordial density field. We present evolution of dark matter density profiles computed on iso-density surfaces which comply with the actual shapes of voids. Unlike spherical density profiles, this approach enables us to demonstrate development of theoretically predicted bucket-like shape of the final density profiles indicating a wide flat core and a sharp transition to high-density void walls.

The possibility of studying non-thermal electron energization in laser-driven plasma experiments of magnetic reconnection is studied using two- and three-dimensional particle-in-cell simulations. It is demonstrated that non-thermal electrons with energies more than an order of magnitude larger than the initial thermal energy can be produced in plasma conditions currently accessible in the laboratory. Electrons are accelerated by the reconnection electric field, being injected at varied distances from the X-points, and in some cases trapped in plasmoids, before escaping the finite-sized system. Trapped electrons can be further energized by the electric field arising from the motion of the plasmoid. This acceleration gives rise to a non-thermal electron component that resembles a power-law spectrum, containing up to ~ 8% of the initial energy of the interacting electrons and ~ 24 % of the initial magnetic energy. Estimates of the maximum electron energy and of the plasma conditions required to observe suprathermal electron acceleration are provided, paving the way for a new platform for the experimental study of particle acceleration induced by reconnection.

Understanding the velocity field is very important for modern cosmology: it gives insights to structure formation in general, and also its properties are crucial ingredients in modelling redshift-space distortions and in interpreting measurements of the kinetic Sunyaev-Zeldovich effect. Unfortunately, characterising the velocity field in cosmological N-body simulations is inherently complicated by two facts: i) The velocity field becomes manifestly multi-valued after shell-crossing and has discontinuities at caustics. This is due to the collisionless nature of dark matter. ii) N-body simulations sample the velocity field only at a set of discrete locations, with poor resolution in low-density regions. In this paper, we discuss how the associated problems can be circumvented by using a phase-space interpolation technique. This method provides extremely accurate estimates of the cosmic velocity fields and its derivatives, which can be properly defined without the need of the arbitrary "coarse-graining" procedure commonly used. We explore in detail the configuration-space properties of the cosmic velocity field on very large scales and in the highly nonlinear regime. In particular, we characterise the divergence and curl of the velocity field, present their one-point statistics, analyse the Fourier-space properties and provide fitting formulae for the velocity divergence bias relative to the non-linear matter power spectrum. We furthermore contrast some of the interesting differences in the velocity fields of warm and cold dark matter models. We anticipate that the high-precision measurements carried out here will help to understand in detail the dynamics of dark matter and the structures it forms.

We present an exact general remeshing scheme to compute analytic integrals of polynomial functions over the intersections between convex polyhedral cells of old and new meshes. In physics applications this allows one to ensure global mass, momentum, and energy conservation while applying higher-order polynomial interpolation. We elaborate on applications of our algorithm arising in the analysis of cosmological N-body data, computer graphics, and continuum mechanics problems. We focus on the particular case of remeshing tetrahedral cells onto a Cartesian grid such that the volume integral of the polynomial density function given on the input mesh is guaranteed to equal the corresponding integral over the output mesh. We refer to this as "physically conservative voxelization". At the core of our method is an algorithm for intersecting two convex polyhedra by successively clipping one against the faces of the other. This algorithm is an implementation of the ideas presented abstractly by Sugihara (1994), who suggests using the planar graph representations of convex polyhedra to ensure topological consistency of the output. This makes our implementation robust to geometric degeneracy in the input. We employ a simplicial decomposition to calculate moment integrals up to quadratic order over the resulting intersection domain. We also address practical issues arising in a software implementation, including numerical stability in geometric calculations, management of cancellation errors, and extension to two dimensions. In a comparison to recent work, we show substantial performance gains. We provide a C implementation intended to be a fast, accurate, and robust tool for geometric calculations on polyhedral mesh elements.

The Lyman-alpha (Lya) emission has played an important role in detecting high-redshift galaxies, including recently distant ones at redshift z > 7. It may also contain important information on the origin of these galaxies. Here, we investigate the formation of a typical L* galaxy and its observational signatures at the earliest stage, by combining a cosmological hydrodynamic simulation with three-dimensional radiative transfer calculations using the newly improved ART^2 code. Our cosmological simulation uses the Aquila initial condition which zooms in onto a Milky Way-like halo with high resolutions, and our radiative transfer couples multi-wavelength continuum, Lya line, and ionization of hydrogen. We find that the modeled galaxy starts to form at redshift z ~ 24 through efficient accretion of cold gas, which produces a strong Lya line with a luminosity of L(Lya) ~ 10^42 erg/s as early as z ~ 14. The Lya emission appears to trace the cold, dense gas. The lines exhibit asymmetric, single-peak profiles, and are shifted to the blue wing, a characteristic feature of gas inflow. Moreover, the contribution to the total Lya luminosity by excitation cooling increases with redshift, and it becomes dominant at z >~ 6. We predict that L* galaxies such as the modeled one may be detected at z <~ 8 by JWST and ALMA with a reasonable integration time. Beyond redshift 12, however, only Lya line may be observable by spectroscopic surveys. Our results suggest that Lya line is one of the most powerful tools to detect the first generation of galaxies, and to decipher their formation mechanism.

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.

A large number of high-redshift galaxies have been discovered via their narrow-band Lya line or broad-band continuum colors in recent years. The nature of the escaping process of photons from these early galaxies is crucial to understanding galaxy evolution and the cosmic reionization. Here, we investigate the escape of Lya, non-ionizing UV-continuum (l = 1300 - 1600 angstrom in rest frame), and ionizing photons (l < 912 angstrom) from galaxies by combining a cosmological hydrodynamic simulation with three-dimensional multi-wavelength radiative transfer calculations. The galaxies are simulated in a box of 5^3 h^-3 Mpc^3 with high resolutions using the Aquila initial condition which reproduces a Milky Way-like galaxy at redshift z=0. We find that the escape fraction (fesc) of these different photons shows a complex dependence on redshift and galaxy properties: fesc(Lya) and fesc(UV) appear to evolve with redshift, and they show similar, weak correlations with galaxy properties such as mass, star formation, metallicity, and dust content, while fesc(Ion) remains roughly constant at ~ 0.2 from z ~ 0 - 10, and it does not show clear dependence on galaxy properties. fesc(Lya) correlates more strongly with fesc(UV) than with fesc(Ion). In addition, we find a relation between the emergent Lya luminosity and the ionizing photon emissivity of Lyman Alpha Emitters (LAEs). By combining this relation with the observed luminosity functions of LAEs at different redshift, we estimate the contribution from LAEs to the reionization of intergalactic medium (IGM). Our result suggests that ionizing photons from LAEs alone are not sufficient to ionize IGM at z > 6, but they can maintain the ionization of IGM at z ~ 0 - 5.

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.

We investigate the large-scale clustering and gravitational interaction of baryons and dark matter (DM) over cosmic time using a set of collisionless N-body simulations. Both components, baryons and DM, are evolved from distinct primordial density and velocity power spectra as predicted by early-universe physics. We first demonstrate that such two-component simulations require an unconventional match between force and mass resolution (i.e. force softening on at least the mean particle separation scale). Otherwise, the growth on any scale is not correctly recovered because of a spurious coupling between the two species at the smallest scales. With these simulations, we then demonstrate how the primordial differences in the clustering of baryons and DM are progressively diminished over time. In particular, we explicitly show how the BAO signature is damped in the spatial distribution of baryons and imprinted in that of DM. This is a rapid process, yet it is still not fully completed at low redshifts. On large scales, the overall shape of the correlation function of baryons and DM differs by 2% at z = 9 and by 0.2% at z = 0. The differences in the amplitude of the BAO peak are approximately a factor of 5 larger: 10% at z = 9 and 1% at z = 0. These discrepancies are, however, smaller than effects expected to be introduced by galaxy formation physics in both the shape of the power spectrum and in the BAO peak, and are thus unlikely to be detected given the precision of the next generation of galaxy surveys. Hence, our results validate the standard practice of modelling the observed galaxy distribution using predictions for the total mass clustering in the Universe.

The microphysical properties of the DM particle can, in principle, be constrained by the properties and abundance of substructures in DM halos, as measured through strong gravitational lensing. Unfortunately, there is a lack of accurate theoretical predictions for the lensing signal of substructures, mainly because of the discreteness noise inherent to N-body simulations. Here we present Recursive-TCM, a method that is able to provide lensing predictions with an arbitrarily low discreteness noise, without any free parameters or smoothing scale. This solution is based on a novel way of interpreting the results of N-body simulations, where particles simply trace the evolution and distortion of Lagrangian phase-space volume elements. We discuss the advantages of this method over the widely used cloud-in-cells and adaptive-kernel smoothing density estimators. Applying the new method to a cluster-sized DM halo simulated in warm and cold DM scenarios, we show how the expected differences in their substructure population translate into differences in the convergence and magnification maps. We anticipate that our method will provide the high-precision theoretical predictions required to interpret and fully exploit strong gravitational lensing observations.

Warm Dark Matter (WDM) cosmologies are a viable alternative to the Cold Dark Matter (CDM) scenario. Unfortunately, an accurate scrutiny of the WDM predictions with N-body simulations has proven difficult due to numerical artefacts. Here, we report on cosmological simulations that, for the first time, are devoid of those problems, and thus, are able to accurately resolve the WDM halo mass function well below the cut-off. We discover a complex picture, with perturbations at different evolutionary stages populating different ranges in the halo mass function. On the smallest mass scales we can resolve, identified objects are typically centres of filaments that are starting to collapse. On intermediate mass scales, objects typically correspond to fluctuations that have collapsed and are in the process of relaxation, whereas the high mass end is dominated by objects similar to haloes identified in CDM simulations. We then explicitly show how the formation of low-mass haloes is suppressed, which translates into a strong cut-off in the halo mass function. This disfavours some analytic formulations that predict a halo mass function that would extend well below the free streaming mass. We argue for a more detailed exploration of the formation of the smallest structures expected to form in a given cosmology, which, we foresee, will advance our overall understanding of structure formation.

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.

Recently, we have shown how current cosmological N-body codes already follow the fine grained phase-space information of the dark matter fluid. Using a tetrahedral tesselation of the three-dimensional manifold that describes perfectly cold fluids in six-dimensional phase space, the phase-space distribution function can be followed throughout the simulation. This allows one to project the distribution function into configuration space to obtain highly accurate densities, velocities, and velocity dispersions. Here, we exploit this technique to show first steps on how to devise an improved particle-mesh technique. At its heart, the new method thus relies on a piecewise linear approximation of the phase space distribution function rather than the usual particle discretisation. We use pseudo-particles that approximate the masses of the tetrahedral cells up to quadrupolar order as the locations for cloud-in-cell (CIC) deposit instead of the particle locations themselves as in standard CIC deposit. We demonstrate that this modification already gives much improved stability and more accurate dynamics of the collisionless dark matter fluid at high force and low mass resolution. We demonstrate the validity and advantages of this method with various test problems as well as hot/warm-dark matter simulations which have been known to exhibit artificial fragmentation. This completely unphysical behaviour is much reduced in the new approach. The current limitations of our approach are discussed in detail and future improvements are outlined.

We investigate the connection between the epoch of reionization and the present day universe, by examining the extended mass reionization histories of dark matter halos identified at z=0. We combine an N-body dark matter simulation of a 600 Mpc volume with a three-dimensional, seminumerical reionization model. This provides reionization redshifts for each particle, which can then be connected with the properties of their halos at the present time. We find that the vast majority of present-day halos with masses larger than ~ few x 10^11 Msun reionize earlier than the rest of the universe. We also find significant halo-to-halo diversity in mass reionization histories, and find that in realistic inhomogenous models, the material within a given halo is not expected to reionize at the same time. In particular, the scatter in reionization times within individual halos is typically larger than the scatter among halos. From our fiducial reionization model, we find that the typical 68% scatter in reionization times within halos is ~ 115 Myr for 10^(12 \pm 0.25) Msun halos, decreasing slightly to ~ 95 Myr for 10^(15 \pm 0.25) Msun halos. We find a mild correlation between reionization history and environment: halos with shorter reionization histories are typically in more clustered environments, with the strongest trend on a scale of ~ 20 Mpc. Material in Milky Way mass halos with short reionization histories is preferentially reionized in relatively large HII regions, implying reionization mostly by sources external to the progenitors of the present-day halo. We investigate the impact on our results of varying the reionization model parameters, which span a range of reionization scenarios with varying timing and morphology.

Structured Adaptive Mesh Refinement (SAMR) is a popular numerical technique to study processes with high spatial and temporal dynamic range. It reduces computational requirements by adapting the lattice on which the underlying differential equations are solved to most efficiently represent the solution. Particularly in astrophysics and cosmology such simulations now can capture spatial scales ten orders of magnitude apart and more. The irregular locations and extensions of the refined regions in the SAMR scheme and the fact that different resolution levels partially overlap, poses a challenge for GPU-based direct volume rendering methods. kD-trees have proven to be advantageous to subdivide the data domain into non-overlapping blocks of equally sized cells, optimal for the texture units of current graphics hardware, but previous GPU-supported raycasting approaches for SAMR data using this data structure required a separate rendering pass for each node, preventing the application of many advanced lighting schemes that require simultaneous access to more than one block of cells. In this paper we present a single-pass GPU-raycasting algorithm for SAMR data that is based on a kD-tree. The tree is efficiently encoded by a set of 3D-textures, which allows to adaptively sample complete rays entirely on the GPU without any CPU interaction. We discuss two different data storage strategies to access the grid data on the GPU and apply them to several datasets to prove the benefits of the proposed method.

In the last decades cosmological N-body dark matter simulations have enabled ab initio studies of the formation of structure in the Universe. Gravity amplified small density fluctuations generated shortly after the Big Bang, leading to the formation of galaxies in the cosmic web. These calculations have led to a growing demand for methods to analyze time-dependent particle based simulations. Rendering methods for such N-body simulation data usually employ some kind of splatting approach via point based rendering primitives and approximate the spatial distributions of physical quantities using kernel interpolation techniques, common in SPH (Smoothed Particle Hydrodynamics)-codes. This paper proposes three GPU-assisted rendering approaches, based on a new, more accurate method to compute the physical densities of dark matter simulation data. It uses full phase-space information to generate a tetrahedral tessellation of the computational domain, with mesh vertices defined by the simulation's dark matter particle positions. Over time the mesh is deformed by gravitational forces, causing the tetrahedral cells to warp and overlap. The new methods are well suited to visualize the cosmic web. In particular they preserve caustics, regions of high density that emerge, when several streams of dark matter particles share the same location in space, indicating the formation of structures like sheets, filaments and halos. We demonstrate the superior image quality of the new approaches in a comparison with three standard rendering techniques for N-body simulation data.

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.

The primordial velocity dispersion of dark matter is small compared to the velocities attained during structure formation. The initial density distribution is close to uniform and it occupies an initial sheet in phase space that is single valued in velocity space. Because of gravitational forces this three dimensional manifold evolves in phase space without ever tearing, conserving phase-space volume and preserving the connectivity of nearby points. N-body simulations already follow the motion of this sheet in phase space. This fact can be used to extract full fine-grained phase-space-structure information from existing cosmological N-body simulations. Particles are considered as the vertices of an unstructured three dimensional mesh, moving in six dimensional phase-space. On this mesh, mass density and momentum are uniquely defined. We show how to obtain the space density of the fluid, detect caustics, and count the number of streams as well as their individual contributions to any point in configuration-space. We calculate the bulk velocity, local velocity dispersions, and densities from the sheet - all without averaging over control volumes. This gives a wealth of new information about dark matter fluid flow which had previously been thought of as inaccessible to N-body simulations. We outline how this mapping may be used to create new accurate collisionless fluid simulation codes that may be able to overcome the sparse sampling and unphysical two-body effects that plague current N-body techniques.

Narrow-band Lya line and broad-band continuum have played important roles in the discovery of high-redshift galaxies in recent years. Hence, it is crucial to study the radiative transfer of both Lya and continuum photons in the context of galaxy formation and evolution in order to understand the nature of distant galaxies. Here, we present a three-dimensional Monte Carlo radiative transfer code, All-wavelength Radiative Transfer with Adaptive Refinement Tree (ART^2), which couples Lya line and multi-wavelength continuum, for the study of panchromatic properties of galaxies and interstellar medium. This code is based on the original version of Li et al., and features three essential modules: continuum emission from X-ray to radio, Lya emission from both recombination and collisional excitation, and ionization of neutral hydrogen. The coupling of these three modules, together with an adaptive refinement grid, enables a self-consistent and accurate calculation of the Lya properties. As an example, we apply ART^2 to a cosmological simulation that includes both star formation and black hole growth, and study in detail a sample of massive galaxies at redshifts z=3.1 - 10.2. We find that these galaxies are Lya emitters (LAEs), whose Lya emission traces the dense gas region, and that their Lya lines show a shape characteristic of gas inflow. Furthermore, the Lya properties, including photon escape fraction, emergent luminosity, and equivalent width, change with time and environment. Our results suggest that LAEs evolve with redshift, and that early LAEs such as the most distant one detected at z ~ 8.6 may be dwarf galaxies with a high star formation rate fueled by infall of cold gas, and a low Lya escape fraction.

The Lya emission has been observed from galaxies over a redshift span z ~ 0 - 8.6. However, the evolution of high-redshift Lya emitters (LAEs), and the link between these populations and local galaxies, remain poorly understood. Here, we investigate the Lya properties of progenitors of a local L* galaxy by combining cosmological hydrodynamic simulations with three-dimensional radiative transfer calculations using the new ART^2 code. We find that the main progenitor (the most massive one) of a Milky Way-like galaxy has a number of Lya properties close to those of observed LAEs at z ~ 2 - 6, but most of the fainter ones appear to fall below the detection limits of current surveys. The Lya photon escape fraction depends sensitively on a number of physical properties of the galaxy, such as mass, star formation rate, and metallicity, as well as galaxy morphology and orientation. Moreover, we find that high-redshift LAEs show blue-shifted Lya line profiles characteristic of gas inflow, and that the Lya emission by excitation cooling increases with redshift, and becomes dominant at z > 6. Our results suggest that some observed LAEs at z ~ 2-6 with luminosity of L_Lya ~ 10^{42-43} ergs/s may be similar to the main progenitor of the Milky Way at high redshift, and that they may evolve into present-day L* galaxies.

We use large-scale simulations to investigate the morphology of reionization during the final, overlap phase. Our method uses an efficient three-dimensional smoothing technique which takes into account the finite mean free path due to absorption systems, lambda, by only smoothing over scales R_s<lambda. The large dynamic range of our calculations is necessary to resolve the neutral patches left at the end of reionization within a representative volume; we find that simulation volumes exceeding several hundred Mpc on a side are necessary in order to properly model reionization when the neutral fraction is ~0.01-0.3. Our results indicate a strong dependence of percolation morphology on a large and uncertain region of model parameter space. The single most important parameter is the mean free path to absorption systems, which serve as opaque barriers to ionizing radiation. If these absorption systems were as abundant as some realistic estimates indicate, the spatial structure of the overlap phase is considerably more complex than previously predicted. In view of the lack of constraints on the mean free path at the highest redshifts, current theories that do not include absorption by Lyman-limit systems, and in particular three-dimensional simulations, may underestimate the abundance of neutral clouds at the end of reionization. This affects predictions for the 21 cm signal associated with reionization, interpretation of absorption features in quasar spectra at z ~5-6, the connection between reionization and the local universe, and constraints on the patchiness and duration of reionization from temperature fluctuations measured in the cosmic microwave background arising from the kinetic Sunyaev-Zel'dovich effect.

We study the buildup of magnetic fields during the formation of Population III star-forming regions, by conducting cosmological simulations from realistic initial conditions and varying the Jeans resolution. To investigate this in detail, we start simulations from identical initial conditions, mandating 16, 32 and 64 zones per Jeans length, and studied the variation in their magnetic field amplification. We find that, while compression results in some amplification, turbulent velocity fluctuations driven by the collapse can further amplify an initially weak seed field via dynamo action, provided there is sufficient numerical resolution to capture vortical motions (we find this requirement to be 64 zones per Jeans length, slightly larger than, but consistent with previous work run with more idealized collapse scenarios). We explore saturation of amplification of the magnetic field, which could potentially become dynamically important in subsequent, fully-resolved calculations. We have also identified a relatively surprising phenomena that is purely hydrodynamic: the higher-resolved simulations possess substantially different characteristics, including higher infall-velocity, increased temperatures inside 1000 AU, and decreased molecular hydrogen content in the innermost region. Furthermore, we find that disk formation is suppressed in higher-resolution calculations, at least at the times that we can follow the calculation. We discuss the effect this may have on the buildup of disks over the accretion history of the first clump to form as well as the potential for gravitational instabilities to develop and induce fragmentation.

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.

It is shown that small-scale magnetic fields present before recombination induce baryonic density inhomogeneities of appreciable magnitude. The presence of such inhomogeneities changes the ionization history of the Universe, which in turn decreases the angular scale of the Doppler peaks and increases Silk damping by photon diffusion. This unique signature could be used to (dis)prove the existence of primordial magnetic fields of strength as small as B~10^(-11) Gauss by upcoming cosmic microwave background observations.

We discuss a new algorithm to generate multi-scale initial conditions with multiple levels of refinements for cosmological "zoom-in" simulations. The method uses an adaptive convolution of Gaussian white noise with a real space transfer function kernel together with an adaptive multi-grid Poisson solver to generate displacements and velocities following first (1LPT) or second order Lagrangian perturbation theory (2LPT). The new algorithm achieves RMS relative errors of order 10^(-4) for displacements and velocities in the refinement region and thus improves in terms of errors by about two orders of magnitude over previous approaches. In addition, errors are localized at coarse-fine boundaries and do not suffer from Fourier-space induced interference ringing. An optional hybrid multi-grid and Fast Fourier Transform (FFT) based scheme is introduced which has identical Fourier space behaviour as traditional approaches. Using a suite of re-simulations of a galaxy cluster halo our real space based approach is found to reproduce correlation functions, density profiles, key halo properties and subhalo abundances with per cent level accuracy. Finally, we generalize our approach for two-component baryon and dark-matter simulations and demonstrate that the power spectrum evolution is in excellent agreement with linear perturbation theory. For initial baryon density fields, it is suggested to use the local Lagrangian approximation in order to generate a density field for mesh based codes that is consistent with Lagrangian perturbation theory instead of the current practice of using the Eulerian linearly scaled densities.

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.

We describe a photon-conserving radiative transfer algorithm, using a spatially-adaptive ray tracing scheme, and its parallel implementation into the adaptive mesh refinement (AMR) cosmological hydrodynamics code, Enzo. By coupling the solver with the energy equation and non-equilibrium chemistry network, our radiation hydrodynamics framework can be utilised to study a broad range of astrophysical problems, such as stellar and black hole (BH) feedback. Inaccuracies can arise from large timesteps and poor sampling, therefore we devised an adaptive time-stepping scheme and a fast approximation of the optically-thin radiation field with multiple sources. We test the method with several radiative transfer and radiation hydrodynamics tests that are given in Iliev et al. (2006, 2009). We further test our method with more dynamical situations, for example, the propagation of an ionisation front through a Rayleigh-Taylor instability, time-varying luminosities, and collimated radiation. The test suite also includes an expanding H II region in a magnetised medium, utilising the newly implemented magnetohydrodynamics module in Enzo. This method linearly scales with the number of point sources and number of grid cells. Our implementation is scalable to 512 processors on distributed memory machines and can include radiation pressure and secondary ionisations from X-ray radiation. It is included in the newest public release of Enzo.

We present the results of a dedicated search for extremely metal-poor stars in the Fornax, Sculptor and Sextans dSphs. Five stars were selected from two earlier VLT/Giraffe and HET/HRS surveys and subsequently followed up at high spectroscopic resolution with VLT/UVES. All of them turned out to have [Fe/H] <= -3 and three stars are below [Fe/H] -3.5. This constitutes the first evidence that the classical dSphs Fornax and Sextans join Sculptor in containing extremely metal-poor stars and suggests that all of the classical dSphs contain extremely metal-poor stars. One giant in Sculptor at [Fe/H]=-3.96 +- 0.10 is the most metal-poor star ever observed in an external galaxy. We carried out a detailed analysis of the chemical abundances of the alpha, iron peak, and the heavy elements, and we performed a comparison with the Milky Way halo and the ultra faint dwarf stellar populations. Carbon, barium and strontium show distinct features characterized by the early stages of galaxy formation and can constrain the origin of their nucleosynthesis.

We suggest a novel discretisation of the momentum equation for Smoothed Particle Hydrodynamics (SPH) and show that it significantly improves the accuracy of the obtained solutions. Our new formulation which we refer to as relative pressure SPH, rpSPH, evaluates the pressure force in respect to the local pressure. It respects Newtons first law of motion and applies forces to particles only when there is a net force acting upon them. This is in contrast to standard SPH which explicitly uses Newtons third law of motion continuously applying equal but opposite forces between particles. rpSPH does not show the unphysical particle noise, the clumping or banding instability, unphysical surface tension, and unphysical scattering of different mass particles found for standard SPH. At the same time it uses fewer computational operations. and only changes a single line in existing SPH codes. We demonstrate its performance on isobaric uniform density distributions, uniform density shearing flows, the Kelvin-Helmholtz and Rayleigh-Taylor instabilities, the Sod shock tube, the Sedov-Taylor blast wave and a cosmological integration of the Santa Barbara galaxy cluster formation test. rpSPH is an improvement these cases. The improvements come at the cost of giving up exact momentum conservation of the scheme. Consequently one can also obtain unphysical solutions particularly at low resolutions.

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.

The transformation of atomic hydrogen to molecular hydrogen through three-body reactions is a crucial stage in the collapse of primordial, metal-free halos, where the first generation of stars (Population III stars) in the Universe are formed. However, in the published literature, the rate coefficient for this reaction is uncertain by nearly an order of magnitude. We report on the results of both adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of the collapse of metal-free halos as a function of the value of this rate coefficient. For each simulation method, we have simulated a single halo three times, using three different values of the rate coefficient. We find that while variation between halo realizations may be greater than that caused by the three-body rate coefficient being used, both the accretion physics onto Population III protostars as well as the long-term stability of the disk and any potential fragmentation may depend strongly on this rate coefficient.

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.

Magnetic fields are generally expected to increase the characteristic mass of stars formed in stellar clusters, because they tend to increase the effective Jeans mass. We test this expectation using adaptive mesh refinement (AMR) magnetohydrodynamic simulations of cluster formation in turbulent magnetized clumps of molecular clouds, treating stars as accreting sink particles. We find that, contrary to the common expectation, a magnetic field of strength in the observed range decreases, rather than increases, the characteristic stellar mass. It (1) reduces the number of intermediate-mass stars that are formed through direct turbulent compression, because sub-regions of the clump with masses comparable to those of stars are typically magnetically subcritical and cannot be compressed directly into collapse, and (2) increases the number of low-mass stars that are produced from the fragmentation of dense filaments. The filaments result from mass accumulation along the field lines. In order to become magnetically supercritical and fragment, the filament must accumulate a large enough column density (proportional to the field strength), which yields a high volume density (and thus a small thermal Jeans mass) that is conducive to forming low-mass stars. We find, in addition, that the characteristic stellar mass is reduced further by outflow feedback. The conclusion is that both magnetic fields and outflow feedback are important in shaping the stellar initial mass function (IMF).

We use the publicly available subhalo catalogs from the Via Lactea simulation along with a Gpc-scale N-body simulation to understand the impact of inhomogeneous reionization on the satellite galaxy population of the Milky Way. The large-volume simulation is combined with a model for reionization that allows us to predict the distribution of reionization times for Milky Way mass halos. Motivated by this distribution, we identify candidate satellite galaxies in the simulation by requiring that any subhalo must grow above a specified mass threshold before it is reionized; after this time the photoionizing background will suppress both the formation of stars and the accretion of gas. We show that varying the reionization time over the range expected for Milky Way mass halos can change the number of satellite galaxies by roughly two orders of magnitude. This conclusion is in contradiction with a number of studies in the literature, and we conclude that this is a result of inconsistent application of the results of Gnedin (2000). We compare our satellite galaxies to observations using both abundance matching and stellar population synthesis methods to assign luminosities to our subhalos and account for observational completeness effects. Additionally, if we assume that the mass threshold is set by the virial temperature Tvir = 8e3K we find that our model accurately matches the vmax distribution, radial distribution, and luminosity function of observed Milky Way satellites for a reionization time zreion = 9.6^{1.0}_{-2.1}, assuming that the Via Lacteasubhalo distribution is representative of the Milky Way. This results in the presence of 119^{+202}_{-50} satellite galaxies.

We describe an implementation of compressible inviscid fluid solvers with block-structured adaptive mesh refinement on Graphics Processing Units using NVIDIA's CUDA. We show that a class of high resolution shock capturing schemes can be mapped naturally on this architecture. Using the method of lines approach with the second order total variation diminishing Runge-Kutta time integration scheme, piecewise linear reconstruction, and a Harten-Lax-van Leer Riemann solver, we achieve an overall speedup of approximately 10 times faster execution on one graphics card as compared to a single core on the host computer. We attain this speedup in uniform grid runs as well as in problems with deep AMR hierarchies. Our framework can readily be applied to more general systems of conservation laws and extended to higher order shock capturing schemes. This is shown directly by an implementation of a magneto-hydrodynamic solver and comparing its performance to the pure hydrodynamic case. Finally, we also combined our CUDA parallel scheme with MPI to make the code run on GPU clusters. Close to ideal speedup is observed on up to four GPUs.

Eulerian hydrodynamical simulations are a powerful and popular tool for modeling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.

(Abridged) We investigate massive star formation in turbulent, magnetized, parsec-scale clumps of molecular clouds including protostellar outflow feedback using Enzo-based MHD simulations with accreting sink particles and effective resolution $2048^3$. We find that, in the absence of regulation by magnetic fields and outflow feedback, massive stars form readily in a turbulent, moderately condensed clump of $\sim 1,600$ solar masses, along with a cluster of hundreds of lower mass stars. The massive stars are fed at high rates by (1) transient dense filaments produced by large-scale turbulent compression at early times, and (2) by the clump-wide global collapse resulting from turbulence decay at late times. In both cases, the bulk of the massive star's mass is supplied from outside a 0.1 pc-sized "core" that surrounds the star. In our simulation, the massive star is clump-fed rather than core-fed. The need for large-scale feeding makes the massive star formation prone to regulation by outflow feedback, which directly opposes the feeding processes. The outflows reduce the mass accretion rates onto the massive stars by breaking up the dense filaments that feed the massive star formation at early times, and by collectively slowing down the global collapse that fuel the massive star formation at late times. The latter is aided by a moderate magnetic field of strength in the observed range. We conclude that the massive star formation in our simulated turbulent, magnetized, parsec-scale clump is outflow-regulated and clump-fed (ORCF for short). An important implication is that the formation of low-mass stars in a dense clump can affect the formation of massive stars in the same clump, through their outflow feedback on the clump dynamics.

Previous high resolution cosmological simulations predict the first stars to appear in the early universe to be very massive and to form in isolation. Here we discuss a cosmological simulation in which the central 50 solar mass clump breaks up into two cores, having a mass ratio of two to one, with one fragment collapsing to densities of 10^{-8} g/cc. The second fragment, at a distance of 800 astronomical units, is also optically thick to its own cooling radiation from molecular hydrogen lines, but is still able to cool via collision-induced emission. The two dense peaks will continue to accrete from the surrounding cold gas reservoir over a period of 10^5 years and will likely form a binary star system.

This White Paper to the National Academy of Sciences Astro2010 Decadal Review Committee outlines some of the outstanding questions regarding the assembly history of Massive Black Holes in the nuclei of galaxies and the revolutionary contributions anticipated in this field from low-frequency gravitational wave astronomy.

We present results of combined N-body and three-dimensional reionization calculations to determine the relationship between reionization history and local environment in a volume 1 Gpc/h across and a resolution of about 1 Mpc. We resolve the formation of about 2x10^6 halos of mass greater than ~10^12 Msun at z=0, allowing us to determine the relationship between halo mass and reionization epoch for galaxies and clusters. For our fiducial reionization model, in which reionization begins at z~15 and ends by z~6, we find a strong bias for cluster-size halos to be in the regions which reionized first, at redshifts 10<z<15. Consequently, material in clusters was reionized within relatively small regions, on the order of a few Mpc, implying that all clusters in our calculation were reionized by their own progenitors. Milky Way mass halos were on average reionized later and by larger regions, with a distribution similar to the global one, indicating that low mass halos are relatively uncorrelated with reionization when only their mass is taken as a prior. On average, we find that most halos with mass less than 10^13 Msun were reionized internally, while almost all halos with masses greater than 10^14 Msun were reionized by their own progenitors. We briefly discuss the implications of this work in light of the "missing satellites" problem and how this new approach may be extended further.

The first stars in the universe, forming at redshifts z>15 in minihalos with masses of order 10^6 Msun, may leave behind black holes as their remnants. These objects could conceivably serve as "seeds" for much larger black holes observed at redshifts z~6. We study the growth of the remnant black holes through accretion including for the first time the emitted accretion radiation with adaptive mesh refinement cosmological radiation-hydrodynamical simulations. The effects of photo-ionization and heating dramatically affect the accretion flow from large scales, resulting in negligible mass growth of the black hole. We compare cases with the accretion luminosity included and neglected to show that the accretion radiation drastically changes the environment within 100 pc of the black hole, where gas temperatures are increased by an order of magnitude. The gas densities are reduced and further star formation in the same minihalo prevented for the two hundred million years of evolution we followed. These calculations show that even without the radiative feedback included most seed black holes do not gain mass as efficiently as has been hoped for in previous theories, implying that black hole remnants of Pop III stars that formed in minihalos are not likely to be the origin of miniquasars. Most importantly, however, these calculations demonstrate that if early stellar mass black holes are indeed accreting close to the Bondi-Hoyle rate with ten percent efficiency they have a dramatic local effect in regulating star formation in the first galaxies.

(Abridged) Using ab initio cosmological Eulerian adaptive mesh refinement radiation hydrodynamical calculations, we discuss how very massive stars start the process of cosmological reionization. The models include non-equilibrium primordial gas chemistry and cooling processes and accurate radiation transport in the Case B approximation using adaptively ray traced photon packages, retaining the time derivative in the transport equation. Supernova feedback is modeled by thermal explosions triggered at parsec scales. All calculations resolve the local Jeans length by at least 16 grid cells at all times and as such cover a spatial dynamic range of ~10^6. These first sources of reionization are highly intermittent and anisotropic and first photoionize the small scales voids surrounding the halos they form in, rather than the dense filaments they are embedded in. As the merging objects form larger, dwarf sized galaxies, the escape fraction of UV radiation decreases and the HII regions only break out on some sides of the galaxies making them even more anisotropic. In three cases, SN blast waves induce star formation in overdense regions that were formed earlier from ionization front instabilities. These stars form tens of parsecs away from the center of their parent DM halo. Approximately 5 ionizing photons are needed per sustained ionization when star formation in 10^6 M_sun halos are dominant in the calculation. As the halos become larger than ~10^7 M_sun, the ionizing photon escape fraction decreases, which in turn increases the number of photons per ionization to 15-50, in calculations with stellar feedback only. Supernova feedback in these more massive halos creates a more diffuse medium, allowing the stellar radiation to escape more easily and maintaining the ratio of 5 ionizing photons per sustained ionization.

The first stars form in dark matter halos of masses ~10^6 M_sun as suggested by an increasing number of numerical simulations. Radiation feedback from these stars expels most of the gas from their shallow potential well of their surrounding dark matter halos. We use cosmological adaptive mesh refinement simulations that include self-consistent Population III star formation and feedback to examine the properties of assembling early dwarf galaxies. Accurate radiative transport is modeled with adaptive ray tracing. We include supernova explosions and follow the metal enrichment of the intergalactic medium. The calculations focus on the formation of several dwarf galaxies and their progenitors. In these halos, baryon fractions in 10^8 solar mass halos decrease by a factor of 2 with stellar feedback and by a factor of 3 with supernova explosions. We find that radiation feedback and supernova explosions increase gaseous spin parameters up to a factor of 4 and vary with time. Stellar feedback, supernova explosions, and H_2 cooling create a complex, multi-phase interstellar medium whose densities and temperatures can span up to 6 orders of magnitude at a given radius. The pair-instability supernovae of Population III stars alone enrich the halos with virial temperatures of 10^4 K to approximately 10^{-3} of solar metallicity. We find that 40% of the heavy elements resides in the intergalactic medium (IGM) at the end of our calculations. The highest metallicity gas exists in supernova remnants and very dilute regions of the IGM.

Numerous cosmological hydrodynamic studies have addressed the formation of galaxies. Here we choose to study the first stages of galaxy formation, including non-equilibrium atomic primordial gas cooling, gravity and hydrodynamics. Using initial conditions appropriate for the concordance cosmological model of structure formation, we perform two adaptive mesh refinement simulations of ~10^8 M_sun galaxies at high redshift. The calculations resolve the Jeans length at all times with more than 16 cells and capture over 14 orders of magnitude in length scales. In both cases, the dense, 10^5 solar mass, one parsec central regions are found to contract rapidly and have turbulent Mach numbers up to 4. Despite the ever decreasing Jeans length of the isothermal gas, we only find one site of fragmentation during the collapse. However, rotational secular bar instabilities transport angular momentum outwards in the central parsec as the gas continues to collapse and lead to multiple nested unstable fragments with decreasing masses down to sub-Jupiter mass scales. Although these numerical experiments neglect star formation and feedback, they clearly highlight the physics of turbulence in gravitationally collapsing gas. The angular momentum segregation seen in our calculations plays an important role in theories that form supermassive black holes from gaseous collapse.

At low temperatures, the main coolant in primordial gas is molecular hydrogen, H2. Recent work has shown that primordial gas that is not collapsing gravitationally but is cooling from an initially ionized state forms hydrogen deuteride, HD, in sufficient amounts to cool the gas to the temperature of the cosmic microwave background. This extra cooling can reduce the characteristic mass for gravitational fragmentation and may cause a shift in the characteristic masses of population III stars. Motivated by the importance of the atomic and molecular data for the cosmological question, we assess several chemical and radiative processes that have hitherto been neglected: the sensitivity of the low temperature H2 cooling rate to the ratio of ortho-H2 to para-H2, the uncertainty in the low temperature cooling rate of H2 excited by collisions with H, the effects of cooling from H2 excited by collisions with H+ and e-, and the large uncertainties in the rates of several of the reactions responsible for determining the H2 fraction in the gas. We show that the most important of the neglected processes is the excitation of H2 by collisions with protons and electrons. This cools the gas more rapidly at early times, and so it forms less H2 and HD at late times. This fact, as well as several of the chemical uncertainties presented here, significantly affects the thermal evolution of the gas. We anticipate that this may lead to clear differences in future detailed 3D studies of first structure formation. Finally, we show that although the thermal evolution of the gas is in principle sensitive to the ortho-para ratio, in practice the standard assumption of a 3:1 ratio produces results that are almost indistinguishable from those produced by a more detailed treatment. (abridged)

The understanding of the formation, life, and death of Population III stars, as well as the impact that these objects had on later generations of structure formation, is one of the foremost issues in modern cosmological research and has been an active area of research during the past several years. We summarize the results presented at "First Stars III," a conference sponsored by Los Alamos National Laboratory, the Kavli Institute for Particle Astrophysics and Cosmology, and the Joint Institute for Nuclear Astrophysics. This conference, the third in a series, took place in July 2007 at the La Fonda Hotel in Santa Fe, New Mexico, U.S.A.

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

Astrophysical relativistic flow problems require high resolution three-dimensional numerical simulations. In this paper, we describe a new parallel three-dimensional code for simulations of special relativistic hydrodynamics (SRHD) using both spatially and temporally structured adaptive mesh refinement (AMR). We used the method of lines to discretize the SRHD equations spatially and a total variation diminishing (TVD) Runge-Kutta scheme for time integration. For spatial reconstruction, we have implemented piecewise linear method (PLM), piecewise parabolic method (PPM), third order convex essentially non-oscillatory (CENO) and third and fifth order weighted essentially non-oscillatory (WENO) schemes. Flux is computed using either direct flux reconstruction or approximate Riemann solvers including HLL, modified Marquina flux, local Lax-Friedrichs flux formulas and HLLC. The AMR part of the code is built on top of the cosmological Eulerian AMR code {\sl enzo}. We discuss the coupling of the AMR framework with the relativistic solvers. Via various test problems, we emphasize the importance of resolution studies in relativistic flow simulations because extremely high resolution is required especially when shear flows are present in the problem. We also present the results of two 3d simulations of astrophysical jets: AGN jets and GRB jets. Resolution study of those two cases further highlights the need of high resolutions to calculate accurately relativistic flow problems.

Using magnetohydrodynamic (MHD) adaptive mesh refinement simulations, we study the formation and early evolution of disk galaxies with a magnetized interstellar medium. For a $10^{10}$ \msun halo with initial NFW dark matter and gas profiles, we impose a uniform $10^{-9}$ G magnetic field and follow its collapse, disk formation and evolution up to 1 Gyr. Comparing to a purely hydrodynamic simulation with the same initial condition, we find that a protogalactic field of this strength does not significantly influence the global disk properties. At the same time, the initial magnetic fields are quickly amplified by the differentially rotating turbulent disk. After the initial rapid amplification lasting $\sim500$ Myr, subsequent field amplification appears self-regulated. As a result, highly magnetized material begin to form above and below the disk. Interestingly, the field strengths in the self-regulated regime agrees well with the observed fields in the Milky Way galaxy both in the warm and the cold HI phase and do not change appreciably with time. Most of the cold phase shows a dispersion of order ten in the magnetic field strength. The global azimuthal magnetic fields reverse at different radii and the amplitude declines as a function of radius of the disk. By comparing the estimated star formation rate (SFR) in hydrodynamic and MHD simulations, we find that after the magnetic field strength saturates, magnetic forces provide further support in the cold gas and lead to a decline of the SFR.

Cosmic reionization progresses as HII regions form around sources of ionizing radiation. Their average size grows continuously until they percolate and complete reionization. We demonstrate how this typical growth can be calculated around the largest, biased sources of UV emission, such as quasars, by further developing an analytical model based on the excursion set formalism. This approach allows us to calculate the sizes and growth of the HII regions created by the progenitors of any dark matter halo of given mass and redshift with a minimum of free parameters. Statistical variations in the size of these pre-existing HII regions are an additional source of uncertainty in the determination of very high redshift quasar properties from their observed HII region sizes. We use this model to demonstrate that the transmission gaps seen in very high redshift quasars can be understood from the radiation of only their progenitors and associated clustered small galaxies. The fit sets a lower limit on the redshift of overlap at z = 5.8 +/- 0.1. This interpretation makes the transmission gaps independent of the age of the quasars observed. If this interpretation were correct it would raise the prospects of using radio interferometers currently under construction to detect the epoch of reionization.

In a hierarchical picture of galaxy formation virialization continually transforms gravitational potential energy into kinetic energies of the baryonic and dark matter. For the gaseous component the kinetic, turbulent energy is transformed eventually into internal thermal energy through shocks and viscous dissipation. Traditionally this virialization and shock heating has been assumed to occur instantaneously allowing an estimate of the gas temperature to be derived from the the virial temperature defined from the embedding dark matter halo velocity dispersion. As the mass grows the virial temperature of a halo grows. Mass accretion hence can be translated into a heating term. We derive this heating rate from the extended Press Schechter formalism and demonstrate its usefulness in semi-analytical models of galaxy formation. Our method explicitly conserves energy unlike the previous impulsive heating assumptions. Our formalism can trivially be applied in all current semi-analytical models as the heating term can be computed directly from the underlying merger trees. Our analytic results for the first cooling halos and the transition from cold to hot accretion are in agreement with numerical simulations.

The first luminous objects in the concordance cosmology form by molecular hydrogen cooling in dark matter dominated halos of masses ~10^6 M_sun. We use Eulerian adaptive mesh refinement simulations to demonstrate that in the presence of a large soft ultraviolet radiation background, molecular hydrogen is the dominant coolant. Even for very large radiation backgrounds, the halo masses that cool and collapse are up to two orders of magnitude smaller than the halos that cool via atomic hydrogen line cooling. The abundance of cooling halos and the cosmic mass fraction contained within them depends exponentially on this critical mass scale. Consequently, the majority of current models of cosmological reionization, chemical evolution, supermassive black hole formation, and galaxy formation underestimate the number of star forming progenitors of a given system by orders of magnitude. At the highest redshifts, this disagreement is largest. We also show that even in the absence of residual electrons, collisional ionization in central shocks create a sufficient amount of electrons to form molecular hydrogen and cool the gas in halos of virial temperatures far below the atomic cooling limit.

The nature of ultracompact H II regions (UCHRs) remains poorly determined. In particular, they are about an order of magnitude more common than would be expected if they formed around young massive stars and lasted for one dynamical time, around 10^4 yr. We here perform three-dimensional numerical simulations of the expansion of an H II region into self-gravitating, radiatively cooled gas, both with and without supersonic turbulent flows. In the laminar case, we find that H II region expansion in a collapsing core produces nearly spherical shells, even if the ionizing source is off-center in the core. This agrees with analytic models of blast waves in power-law media. In the turbulent case, we find that the H II region does not disrupt the central collapsing region, but rather sweeps up a shell of gas in which further collapse occurs. Although this does not constitute triggering, as the swept-up gas would eventually have collapsed anyway, it does expose the collapsing regions to ionizing radiation. We suggest that these regions of secondary collapse, which will not all themselves form massive stars, may form the bulk of observed UCHRs. As the larger shell will take over 10^5 years to complete its evolution, this could solve the timescale problem. Our suggestion is supported by the ubiquitous observation of more diffuse emission surrounding UCHRs.

(Abridged) Galaxies form in hierarchically assembling dark matter halos. With cosmological three dimensional adaptive mesh refinement simulations, we explore in detail the virialization of baryons in the concordance cosmology, including optically thin primordial gas cooling. We focus on early protogalaxies with virial temperatures of 10^4 K and their progenitors. Without cooling, virial heating occurs in shocks close to the virial radius for material falling in from voids. Material in dense filaments penetrates deeper to about half that radius. With cooling the virial shock position shrinks and also the filaments reach scales as small as a third the virial radius. The temperatures in protogalaxies found in adiabatic simulations decrease by a factor of two from the center and show flat entropy cores. In cooling halos the gas reaches virial equilibrium with the dark matter potential through its turbulent velocities. We observe turbulent Mach numbers ranging from one to three in the cooling cases. This turbulence is driven by the large scale merging and interestingly remains supersonic in the centers of these early galaxies even in the absence of any feedback processes. The virial theorem is shown to approximately hold over 3 orders of magnitude in length scale with the turbulent pressure prevailing over the thermal energy. The turbulent velocity distributions are Maxwellian and by far dominate the small rotation velocities associated with the total angular momentum of the galaxies. Decomposing the velocity field using the Cauchy-Stokes theorem, we show that ample amounts of vorticity are present around shocks even at the very centers of these objects.

We have performed a large set of high-resolution cosmological simulations using smoothed particle hydrodynamics to study the formation of the first luminous objects in the LCDM cosmology. We follow the collapse of primordial gas clouds in eight early structures and document the scatter in the properties of the first star-forming clouds. Our first objects span formation redshifts from z~10 to z~50 and cover an order of magnitude in halo mas. We find that the physical properties of the central star-forming clouds are very similar in all of the simulated objects despite significant differences in formation redshift and environment. The physical properties of the clouds have little correlation with spin, mass, or assembly history of the host halo. The collapse of protostellar objects at higher redshifts progresses much more rapidly due to the higher densities, which accelerates the formation of molecular hydrogen, enhances initial cooling and shortens the dynamical timescales. The mass of the star-forming clouds cover a broad range, from a few hundred to a few thousand solar masses, and exhibit various morphologies: some of have disk-like structures nearly rotational supported; others form flattened spheroids; still others form bars. All of them develop a single proto-stellar `seed' which does not fragment into multiple objects up to the moment that the central gas becomes optically thick to H2 cooling lines. At this time, the instantaneous mass accretion rate onto the centre varies significantly from object to object, with disk-like structuress have the smallest mass accretion rates. The formation epoch and properties of the star-forming clouds are sensitive to the values of cosmological parameters.

Understanding how seed black holes grow into intermediate and supermassive black holes (IMBHs and SMBHs, respectively) has important implications for the duty-cycle of active galactic nuclei (AGN), galaxy evolution, and gravitational wave astronomy. Most studies of the cosmological growth and merger history of black holes have used semianalytic models and have concentrated on SMBH growth in luminous galaxies. Using high resolution cosmological N-body simulations, we track the assembly of black holes over a large range of final masses -- from seed black holes to SMBHs -- over widely varying dynamical histories. We used the dynamics of dark matter halos to track the evolution of seed black holes in three different gas accretion scenarios. We have found that growth of Sagittarius A* - size SMBH reaches its maximum mass M_{SMBH}~10^6Msun at z~6 through early gaseous accretion episodes, after which it stays at near constant mass. At the same redshift, the duty-cycle of the host AGN ends, hence redshift z=6 marks the transition from an AGN to a starburst galaxy which eventually becomes the Milky Way. By tracking black hole growth as a function of time and mass, we estimate that the IMBH merger rate reaches a maximum of R_{max}=55 yr^-1 at z=11. From IMBH merger rates we calculate N_{ULX}=7 per Milky Way type galaxy per redshift in redshift range 2<z<6.

Recent observations show that the number of stars with very low metallicities in the dwarf spheroidal satellites of the Milky Way is low, despite the low average metallicities of stars in these systems. We undertake numerical simulations of star formation and metal enrichment of dwarf galaxies in order to verify whether this result can be reproduced with "standard" assumptions. The answer is likely to be negative, unless some selection bias against very low metallicity stars is present in the observations.

As part of the Dwarf galaxies Abundances and Radial-velocities Team (DART) program, we have measured the metallicities of a large sample of stars in four nearby dwarf spheroidal galaxies (dSph's): Sculptor, Sextans, Fornax, and Carina. The low mean metal abundances and the presence of very old stellar populations in these galaxies have supported the view that they are fossils from the early universe. However, contrary to naive expectations, we find a significant lack of stars with metallicities below [Fe/H] ~ -3 dex in all four systems. This suggests that the gas that made up the stars in these systems had been uniformly enriched prior to their formation. Furthermore, the metal-poor tail of the dSph metallicity distribution is significantly different from that of the Galactic halo. These findings show that the progenitors of nearby dSph's appear to have been fundamentally different from the building blocks of the Milky Way, even at the earliest epochs.

Primordial stars are likely to be very massive $\geq30\Msun$, form in isolation, and will likely leave black holes as remnants in the centers of their host dark matter halos in the mass range $10^{6}-10^{10}\Ms$. Such early black holes, at redshifts z$\gtsim10$, could be the seed black holes for the many supermassive black holes found in galaxies in the local universe. If they exist, their mergers with nearby supermassive black holes may be a prime signal for long wavelength gravitational wave detectors. We simulate formation of black holes in the center of high redshift dark matter halos and explore implications of initial natal kick velocities conjectured by some formation models. The central concentration of early black holes in present day galaxies is reduced if they are born even with moderate kicks of tens of km/s. The modest kicks allow the black holes to leave their parent halo, which consequently leads to dynamical friction being less effective on the lower mass black holes as compared to those still embedded in their parent halos. Therefore, merger rates may be reduced by more than an order of magnitude. Using analytical and illustrative cosmological N--body simulations we quantify the role of natal kicks of black holes formed from massive metal free stars on their merger rates with supermassive black holes in present day galaxies. Our results also apply to black holes ejected by the gravitational slingshot mechanism.

As part of the DART project we have used the ESO/2.2m Wide Field Imager in conjunction with the VLT/FLAMES GIRAFFE spectrograph to study the detailed properties of the resolved stellar population of the Fornax dwarf spheroidal galaxy out to and beyond its tidal radius. We re-derived the structural parameters of the Fornax dwarf spheroidal using our wide field imaging covering the galaxy out to its tidal radius, and analysed the spatial distribution of the Fornax stars of different ages as selected from Colour-Magnitude Diagram analysis. We have obtained accurate velocities and metallicities from spectra in the CaII triplet wavelength region for 562 Red Giant Branch stars which have velocities consistent with membership in Fornax dwarf spheroidal. We have found evidence for the presence of at least three distinct stellar components: a young population (few 100 Myr old) concentrated in the centre of the galaxy, visible as a Main Sequence in the Colour-Magnitude Diagram; an intermediate age population (2-8 Gyr old); and an ancient population (> 10Gyr),which are distinguishable from each other kinematically, from the metallicity distribution and in the spatial distribution of stars found in the Colour-Magnitude Diagram. From our spectroscopic analysis we find that the ``metal rich'' stars ([Fe/H]> -1.3) show a less extended and more concentrated spatial distribution, and display a colder kinematics than the ``metal poor'' stars ([Fe/H<-1.3). There is tentative evidence that the ancient stellar population in the centre of Fornax does not exhibit equilibrium kinematics. This could be a sign of a relatively recent accretion of external material.

We study the formation of the first generation of stars in the standard cold dark matter model, using a very high-resolution hydordynamic simulations. Our simulation achieves a dynamic range of 10^{10} in length scale. With accurate treatment of atomic and molecular physics, it allows us to study the chemo-thermal evolution of primordial gas clouds to densities up to n = 10^{16}/cc without assuming any a priori equation of state; a six orders of magnitudes improvement over previous three-dimensional calculations. All the relevant atomic and molecular cooling and heating processes, including cooling by collision-induced continuum emission, are implemented. For calculating optically thick H2 cooling at high densities, we use the Sobolev method. To examine possible gas fragmentation owing to thermal instability, we compute explicitly the growth rate of isobaric perturbations. We show that the cloud core does not fragment in either the low-density or high-density regimes. We also show that the core remains stable against gravitational deformation and fragmentation. We obtain an accurate gas mass accretion rate within a 10 Msun innermost region around the protostar. The protostar is accreting the surrounding hot gas at a rate of 0.001-0.01 Msun/yr. From these findings we conclude that primordial stars formed in early minihalos are massive. We carry out proto-stellar evolution calculations using the obtained accretion rate. The resulting mass of the first star is M_ZAMS = 60-100 Msun, with the exact mass dependent on the actual accretion rate.

The concordance model of cosmology and structure formation predicts the formation of isolated very massive stars at high redshifts in dark matter dominated halos of 10^5 to 10^6 Msun. These stars photo-ionize their host primordial molecular clouds, expelling all the baryons from their halos. When the stars die, a relic HII region is formed within which large amounts of molecular hydrogen form which will allow the gas to cool efficiently when gravity assembles it into larger dark matter halos. The filaments surrounding the first star hosting halo are largely shielded and provide the pathway for gas to stream into the halo when the star has died. We present the first fully three dimensional cosmological radiation hydrodynamical simulations that follow all these effects. A novel adaptive ray casting technique incorporates the time dependent radiative transfer around point sources. This approach is fast enough so that radiation transport, kinetic rate equations, and hydrodynamics are solved self-consistently. It retains the time derivative of the transfer equation and is explicitly photon conserving. This method is integrated with the cosmological adaptive mesh refinement code enzo, and runs on distributed and shared memory parallel architectures. Where applicable the three dimensional calculation not only confirm expectations from earlier one dimensional results but also illustrate the multi-fold hydrodynamic complexities of HII regions. In the absence of stellar winds the circumstellar environments of the first supernovae and putative early gamma-ray bursts will be of low density ~1 cm^-3. Albeit marginally resolved, ionization front instabilities lead to cometary and elephant trunk like small scale structures reminiscent of nearby star forming regions.

We discuss an analytic approach for modeling structure formation in sheets, filaments and knots. This is accomplished by combining models of triaxial collapse with the excursion set approach: sheets are defined as objects which have collapsed along only one axis, filaments have collapsed along two axes, and halos are objects in which triaxial collapse is complete. In the simplest version of this approach, which we develop here, large scale structure shows a clear hierarchy of morphologies: the mass in large-scale sheets is partitioned up among lower mass filaments, which themselves are made-up of still lower mass halos. Our approach provides analytic estimates of the mass fraction in sheets, filaments and halos, and its evolution, for any background cosmological model and any initial fluctuation spectrum. In the currently popular $\Lambda$CDM model, our analysis suggests that more than 99% of the cosmic mass is in sheets, and 72% in filaments, with mass larger than $10^{10} M_{\odot}$ at the present time. For halos, this number is only 46%. Our approach also provides analytic estimates of how halo abundances at any given time correlate with the morphology of the surrounding large-scale structure, and how halo evolution correlates with the morphology of large scale structure.

In these lecture notes we review the current knowledge about the formation of the first luminous objects. We start from the cosmological context of hierarchical models of structure formation, and discuss the main physical processes which are believed to lead to primordial star formation, i.e. the cooling processes and the chemistry of molecules (especially H2) in a metal-free gas. We then describe the techniques and results of numerical simulations, which indicate that the masses of the first luminous objects are likely to be much larger than that of present-day stars. Finally, we discuss the scenario presented above, exposing some of the most interesting problems which are currently being investigated, such as that of the feedback effects of these objects.

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.

Recent simulations of the formation of the first luminous objects in the universe predict isolated very massive stars to form in dark matter halos with virial temperatures large enough to allow significant amounts of molecular hydrogen to form. We construct a semi-analytic model based on the Press-Schechter formalism and calibrate the minimum halos mass that may form a primordial star with the results from extensive adaptive mesh refinement simulations. The model also includes star formation in objects with virial temperatures in excess of ten thousand Kelvin. The free parameters are tuned to match the optical depth measurements by the WMAP satellite. The models explicitly includes the negative feedback of the destruction of molecular hydrogen by a soft UV background which is computed self-consistently. We predict high redshift supernova rates as one of the most promising tools to test the current scenario of primordial star formation. The supernova rate from primordial stars peaks at redshifts ~20. Using an analytic model for the luminosities of pair-instability supernovae we predict observable magnitudes and discuss possible observational strategies. Such supernovae would release enough metals corresponding to a uniform enrichment to a few hundred thousands of solar metalicity. If some of these stars produce gamma ray bursts our rates will be directly applicable to understanding the anticipated results from the SWIFT satellite. This study highlights the great potential for the James Webb space telescope in probing cosmic structure at redshifts greater than 20.

We have found evidence for the presence of two distinct ancient stellar components (both geq 10 Gyr old) in the Sculptor dwarf spheroidal galaxy. We used the ESO Wide Field Imager (WFI) in conjunction with the VLT/FLAMES spectrograph to study the properties of the resolved stellar population of Sculptor out to and beyond the tidal radius. We find that two components are discernible in the spatial distribution of Horizontal Branch stars in our imaging, and in the [Fe/H] and v_hel distributions for our large sample of spectroscopic measurements. They can be generally described as a ``metal-poor'' component ([Fe/H]< -1.7) and a ``metal-rich'' component ([Fe/H]> -1.7). The metal-poor stars are more spatially extended than the metal-rich stars, and they also appear to be kinematically distinct. These results provide an important insight into the formation processes of small systems in the early universe and the conditions found there. Even this simplest of galaxies appears to have had a surprisingly complex early evolution.

We study the effects of the initial conditions of turbulent molecular clouds on the ionization structure in newly formed H_{ii} regions, using three-dimensional, photon-conserving radiative transfer in a pre-computed density field from three-dimensional compressible turbulence. Our results show that the initial density structure of the gas cloud can play an important role in the resulting structure of the H_{ii} region. The propagation of the ionization fronts, the shape of the resulting H_{ii} region, and the total mass ionized depend on the properties of the turbulent density field. Cuts through the ionized regions generally show ``butterfly'' shapes rather than spherical ones, while emission measure maps are more spherical if the turbulence is driven on scales small compared to the size of the H_{ii} region. The ionization structure can be described by an effective clumping factor $\zeta=< n > \cdot < n^2>/< n>^2$, where $n$ is number density of the gas. The larger the value of $\zeta$, the less mass is ionized, and the more irregular the H_{ii} region shapes. Because we do not follow dynamics, our results apply only to the early stage of ionization when the speed of the ionization fronts remains much larger than the sound speed of the ionized gas, or Alfv\'en speed in magnetized clouds if it is larger, so that the dynamical effects can be negligible.

We study the escape of Lyman continuum radiation from the disks of dwarf starburst galaxies, with and without supershells, by solving the radiation transfer problem of stellar radiation through them. We model disks with Md=10^8-10Msun, with exponential surface density profiles as a function of redshift, and model the of repeated supernova explosions driving supershells out of the disks, using the hydrodynamic simulation code, ZEUS-3D. The amount of star formation is assumed proportional to mass above some density threshold. We vary the threshold to explore the range of star formation efficiencies, f*=0.006, 0.06, and 0.6. We find that the interstellar gas swept up in dense supershells can effectively trap the ionizing photons, before the supershells blow out of the disks. The blow-outs then create galactic outflows, chimneys which allow the photons to escape directly to the intergalactic medium. Our results are consistent with escape fractions of less than 0.1 measured in local dwarf starburst galaxies, because they are likely observed while the starbursts are young, before blow-out. We suggest that high-redshift dwarf starburst galaxies may make a substantial contribution to the UV background radiation with total escape fractions >0.2, as expected if star formation efficiencies >0.06.

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

We combine fast radiative transfer calculations with high resolution hydrodynamical simulations to study an epoch of early hydrogen reionisation by primordial stellar sources at redshifts 15<z<30. With relatively conservative assumptions, population III star formation proceeds in a self-regulated manner both locally and globally and, for a conventional LCDM cosmology, can significantly reionise the intergalactic medium between 15<z<20 as long as a large fraction of ionising photons can escape from these earliest galaxies. We then combine these results with our earlier work focusing on the role of population II stars in galaxies with virial temperatures >10^4K at redshifts 5<z< 20. Hence, we construct a complete reionisation history of the Universe which matches the Thomson optical depths as measured by the WMAP satellite as well as the evolution of the Gunn Peterson optical depth as seen in the asborption spectra of the higest redshift quasars. We find that even with conservative estimates for the impact of negative feedback mechanisms, primordial stellar sources contribute significantly to early reionisation. Future observations of a Thomson optical depth of tau_e>~0.13 would bolster the claim for the existence of population III stars similar to the ones studied here.

Gamma-ray bursts are promising tools for tracing the formation of high redshift stars, including the first generation. At very high redshifts the reverse shock emission lasts longer in the observer frame, and its importance for detection and analysis purposes relative to the forward shock increases. We consider two different models for the GRB environment, based on current ideas about the redshift dependence of gas properties in galaxies and primordial star formation. We calculate the observed flux as a function of the redshift and observer time for typical GRB afterglows, taking into account intergalactic photoionization and Lyman-$\alpha$ absorption opacity as well as extinction by the Milky Way Galaxy. The fluxes in the X-ray and near IR bands are compared with the sensitivity of different detectors such as Chandra, XMM, Swift XRT and JWST. Using standard assumptions, we find that Chandra, XMM and Swift XRT can potentially detect GRBs in the X-ray band out to very high redshifts $z\gtrsim$ 30. In the K and M bands, the JWST and ground-based telescopes are potentially able to detect GRBs even one day after the trigger out to $z\sim$ 16 and 33, if present. While the X-ray band is insensitive to the external density and to reverse shocks, the near IR bands provide a sensitive tool for diagnosing both the environment and the reverse shock component.

The mechanisms which could lead to chemo-thermal instabilities and fragmentation during the formation of primordial protostars are investigated analytically. We introduce approximations for H2 cooling rates bridging the optically thin and thick regimes. These allow us to discuss instabilities up to densities when protostars become optically thick to continuum radiation (n~10^16 cm^-3). During the collapse, instability arises at two different stages: at low density (n~10^8-10^11 cm^-3), it is due to fast 3-body reactions converting H into H2; at high density (n>10^13 cm^-3), it is due to Collisional Induced Emission (CIE). In agreement with the 3D simulations, we find that the instability at low densities cannot lead to fragmentation, because fluctuations do not survive turbulent mixing, and because their growth is slow. The situation at high density is similar. The CIE-induced instability is as weak as the low density one, with similar ratios of growth and dynamical time scales. Fluctuation growth time is longer than free fall time, and fragmentation seems unlikely. One then expects the first stars to be massive, not to form binaries nor harbour planets. Nevertheless, full 3D simulations are required. They could become possible using simplified estimates of radiative transfer effects, which we show to work very well in the 1D case. This indicates that the effects of radiative transfer during the initial stages of formation of primordial protostars can be treated as local corrections to cooling. (Abridged)

We examine the effect of preheating of the intergalactic medium on galaxy formation using cosmological hydrodynamical simulations. By performing simulations both with and without a simple model for preheating, we analyse and compare the angular momentum distributions of the dark matter and the baryons. Preheating unbinds baryons from their dark matter haloes, yielding a baryonic mass fraction that declines with decreasing halo mass. In addition, the spin parameter of the gas is reduced with respect to the case without preheating, while the misalignment between the angular momentum directions of the gas and dark matter increases strongly. The angular momentum distributions of individual haloes reveal that preheating decreases (increases) the mass fraction with low (negative) specific angular momentum. We discuss the implications of these findings for the formation of disk galaxies in a preheated intergalactic medium, and compare our results to the predictions of Maller & Dekel (2002), who propose an alternative interpretation for the origin of the angular momentum of (proto)-galaxies.

(abridged) We use large cosmological simulations to study the origin of primordial star-forming clouds in a Lambda CDM universe, by following the formation of dark matter halos and the cooling of gas within them. To model the physics of chemically pristine gas, we employ a non-equilibrium treatment of the chemistry of 9 species and include cooling by molecular hydrogen. We explore the hierarchical growth of bound structures forming at redshifts z = 25 - 30 with total masses in the range 10^5 - 10^6 Msun. The complex interplay between the gravitational formation of dark halos and the thermodynamic and chemical evolution of the gas clouds compromises analytic estimates of the critical H2 fraction. Dynamical heating from mass accretion and mergers opposes relatively inefficient cooling by molecular hydrogen, delaying the production of star-forming clouds in rapidly growing halos. We also investigate the impact of photo-dissociating ultra-violet (UV) radiation on the formation of primordial gas clouds. We consider two extreme cases by first including a uniform radiation field in the optically thin limit and secondly by accounting for the maximum effect of gas self-shielding in virialized regions. In both the cases we consider, the overall impact can be described by computing an equilibrium H2 abundance for the radiation flux and defining an effective shielding factor. Based on our numerical results, we develop a semi-analytic model of the formation of the first stars, and demonstrate how it can be coupled with large N-body simulations to predict the star formation rate in the early universe.

We study the reionisation of the Universe by stellar sources using a numerical approach that combines fast 3D radiative transfer calculations with high resolution hydrodynamical simulations. Ionising fluxes for the sources are derived from intrinsic star formation rates computed in the underlying hydrodynamical simulations. Our mass resolution limit for sources is M~ 4.0 x 10^7 h^-1 M_sol, which is roughly an order of magnitude smaller than in previous studies of this kind. Our calculations reveal that the reionisation process is sensitive to the inclusion of dim sources with masses below ~10^9 h^-1 M_sol. We present the results of our reionisation simulation assuming a range of escape fractions for ionising photons and make statistical comparisons with observational constraints on the neutral fraction of hydrogen at z~6 derived from the z=6.28 SDSS quasar of Becker and coworkers. Our best fitting model has an escape fraction of ~20% and causes reionisation to occur by z~8, although the IGM remains fairly opaque until z~6. In order to simultaneously match the observations from the z=6.28 SDSS quasar and the optical depth measurement from WMAP with the sources modeled here, we require an evolving escape fraction that rises from f_esc=0.20 near z~6 to f_esc>~10 at z~18.

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 simulations to investigate the effects of a soft X-ray background, that could have been produced by an early generation of mini-quasars, on the subsequent cooling and collapse of high redshift pregalactic clouds. The simulations use an Eulerian adaptive mesh refinement technique with initial conditions drawn from a flat Lambda-dominated cold dark matter model cosmology to follow the nonequilibrium chemistry of nine chemical species in the presence of both a soft ultraviolet Lyman-Werner H_2 photodissociating flux and a soft X-ray background extending to 7.2 keV. Although we vary the normalization of the X-ray background by two orders of magnitude, the positive feedback effect of the X-rays on cooling and collapse of the pregalactic cloud expected due to the increased electron fraction is quite mild, only weakly affecting the mass threshold for collapse and the fraction of gas within the cloud that is able to cool. Inside most of the cloud we find that H_2 is in photodissociation equilibrium with the soft UV flux. The net buildup of the electron density needed to enhance H_2 formation occurs too slowly compared to the H_2 photodissociation and dynamical timescales within the cloud to overcome the negative impact of the soft UV photodissociating flux on cloud collapse. However, we find that even in the most extreme cases the first objects to form do rely on molecular hydrogen as coolant and stress that our results do not justify the neglect of these objects in models of galaxy formation.

Using radiative transfer calculations and cosmological simulations of structure formation, we study constraints that can be placed on the nature of the cosmic ultraviolet (UV) background in the redshift interval 2.5< z <5. Our approach makes use of observational estimates of the opacities of hydrogen and singly ionised helium in the intergalactic medium during this epoch. In particular, we model the reionisation of He II by sources of hard ultraviolet radiation, i.e. quasars, and infer values for our parameterisation of this population from observational estimates of the opacity of the He II Lyman-alpha forest. Next, we estimate the photoionisation rate of H I from these sources and find that their contribution to the ionising background is insufficient to account for the measured opacity of the H I Lyman-alpha forest at a redshift z 3. This motivates us to include a soft, stellar component to the ionising background to boost the hydrogen photoionisation rate, but which has a negligible impact on the He II opacity. In order to simultaneously match observational estimates of the H I and He II opacities, we find that galaxies and quasars must contribute about equally to the ionising background in H I at z~ 3.

We use non-radiative N-body/SPH simulations of structure formation in a LCDM cosmology to compare the angular momentum distributions of dark matter and gas in a large sample of halos. We show that the two components have identical spin parameter distributions and that their angular momentum distributions within individual halos are very similar, all in excellent agreement with standard assumptions. Despite these similarities, however, we find that the angular momentum vectors of the gas and dark matter are poorly aligned, with a median misalignment angle of about 30 degrees, which might have important implications for spin correlation statistics used in weak lensing studies. We present distributions for the component of the angular momentum that is aligned with the total angular momentum of each halo, and find that for between 5 and 50 percent of the mass this component is negative. This disagrees with the `Universal' angular momentum distribution suggested by Bullock et al. (2001), for which the mass fraction with negative specific angular momentum is zero. We discuss the implications of our results for the formation of disk galaxies. Since galactic disks generally do not contain counter-rotating stars or gas, disk formation cannot occur under detailed conservation of specific angular momentum. We suggest that the material with negative specific angular momentum combines with positive angular momentum material to build a bulge component, and show that in such a scenario the remaining material can form a disk with a density distribution that is very close to exponential.

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 study the reionization of Helium II by quasars using a numerical approach that combines 3D radiative transfer calculations with cosmological hydrodynamical simulations. Sources producing the ionizing radiation are selected according to an empirical quasar luminosity function and are assigned luminosities according to their intrinsic masses. We present models in which these parameters are varied and examine characteristics of the resultant reionization process that distinguish the various cases. In addition, we extract artificial spectra from the simulations and quantify statistical properties of the spectral features in each model. We find that the most important factor affecting the evolution of He II reionization is the cumulative number of ionizing photons that are produced by the sources. Comparisons between He II opacities measured observationally and those obtained by our analysis reveal that the available ranges in plausible values for the parameters provide enough leeway to provide a satisfactory match. However, one property common to all our calculations is that the epoch of Helium II reionization must have occurred at a redshift between 3 < z < 4. If so, future observational programs will be able to directly trace the details of the ionization history of helium and probe the low density phase of the intergalactic medium during this phase of the evolution of the Universe.

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

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.

While the modern stellar IMF shows a rapid decline with increasing mass, theoretical investigations suggest that very massive stars (>100 solar masses) may have been abundant in the early universe. Other calculations also indicate that, lacking metals, these same stars reach their late evolutionary stages without appreciable mass loss. After central helium burning, they encounter the electron-positron pair instability, collapse, and burn oxygen and silicon explosively. If sufficient energy is released by the burning, these stars explode as brilliant supernovae with energies up to 100 times that of an ordinary core collapse supernova. They also eject up to 50 solar masses of radioactive Ni56. Stars less massive than 140 solar masses or more massive than 260 solar masses should collapse into black holes instead of exploding, thus bounding the pair-creation supernovae with regions of stellar mass that are nucleosynthetically sterile. Pair-instability supernovae might be detectable in the near infrared out to redshifts of 20 or more and their ashes should leave a distinctive nucleosynthetic pattern.

Spherical density profiles and specific angular momentum profiles of Dark Matter halos found in cosmological N-body simulations have been measured extensively. The distribution of the total angular momentum of dark matter halos is also used routinely in semi-analytic modeling of the formation of disk galaxies. However, it is unclear whether the initial (i.e. at the time the halo is assembled) angular momentum distributions of baryons is related to the dark matter at all. Theoretical models for ellipticities in weak lensing studies often rely on an assumed correlation of the angular momentum vectors of dark matter and gas in galaxies. Both of these assumptions are shown to be in reasonable agreement with high resolution cosmological smoothed particle hydrodynamical simulations that follow the dark matter as long as only adiabatic gas physics are included. However, we argue that in more realistic models of galaxy formation one expects pressure forces to play a significant role at turn--around. Consequently the torquing force on DM and baryons will be uncorrelated and their respective angular momenta are not expected to align. An SPH simulation with ad-hoc feedback is presented that illustrates these effects. Massive low redshift elliptical galaxies may be a notable exception where "light may trace mass".

We describe a novel adaptive ray tracing scheme to solve the equation of radiative transfer around point sources in hydrodynamical simulations. The angular resolution adapts to the local hydrodynamical resolution and hence is of use for adaptive meshes as well as adaptive smooth particle hydrodynamical simulations. Recursive creation of rays ensures ease of implementation. The multiple radial integrations needed to solve the time dependent radiative transfer are sped up significantly using a quad-tree once the rays are cast. Simplifications advantageous for methods with one radiation source are briefly discussed. The suggested method is easily generalized to speed up Monte Carlo radiative transfer techniques. In summary a nearly optimal use of long characteristics is presented and aspects of its implementation and comparison to other methods are given.

Absorption or emission against the cosmic microwave background radiation (CMB) may be observed in the redshifted 21cm line if the spin temperature of the neutral intergalactic medium prior to reionization differs from the CMB temperature. This so-called 21cm tomography should reveal important information on the physical state of the intergalactic medium at high redshifts. The fluctuations in the redshifted 21 cm, due to gas density inhomogeneities at early times, should be observed at meter wavelengths by the next generation radio telescopes such as the proposed {\it Square Kilometer Array (SKA)}. Here we show that the extra-galactic radio sources provide a serious contamination to the brightness temperature fluctuations expected in the redshifted 21 cm emission from the IGM at high redshifts. Unless the radio source population cuts off at flux levels above the planned sensitivity of SKA, its clustering noise component will dominate the angular fluctuations in the 21 cm signal. The integrated foreground signal is smooth in frequency space and it should nonetheless be possible to identify the sharp spectral feature arising from the non-uniformities in the neutral hydrogen density during the epoch when the first UV sources reionize the intergalactic medium.

We present a new approach to numerically model continuum radiative transfer based on the Optically Thin Variable Eddington Tensor (OTVET) approximation. Our method insures the exact conservation of the photon number and flux (in the explicit formulation) and automatically switches from the optically thick to the optically thin regime. It scales as N logN with the number of hydrodynamic resolution elements and is independent of the number of sources of ionizing radiation (i.e. works equally fast for an arbitrary source function). We also describe an implementation of the algorithm in a Soften Lagrangian Hydrodynamic code (SLH) and a multi--frequency approach appropriate for hydrogen and helium continuum opacities. We present extensive tests of our method for single and multiple sources in homogeneous and inhomogeneous density distributions, as well as a realistic simulation of cosmological reionization.

The incorporation of radiative transfer effects into cosmological hydrodynamical simulations is essential for understanding how the intergalactic medium (IGM) makes the transition from a neutral medium to one that is almost fully ionized. Here, we present an approximate numerical method designed to study in a statistical sense how a cosmological density field is ionized by a set of discrete point sources. A diffuse background radiation field is also computed self-consistently in our procedure. The method requires relatively few time steps and can be employed with simulations having high resolution. We describe the details of the algorithm and provide a description of how the method can be applied to the output from a pre-existing cosmological simulation to study the systematic reionization of a particular ionic species. As a first application, we compute the reionization of He II by quasars in the redshift range 3 to 6.

At the earliest epochs of structure formation in cold dark matter (CDM) cosmologies, the smallest nonlinear objects are the numerous small halos that condense with virial temperatures below 10,000 K. Such ``minihalos'' are not yet resolved in large-scale three-dimensional cosmological simulations. Here we employ a semi-analytic method, combined with three-dimensional simulations of individual minihalos, to examine their importance during cosmological reionization. We show that, depending on when reionization takes place, they potentially play an important role as sinks of ionizing radiation. If reionization occurs at sufficiently high redshifts (z_r > 20), the intergalactic medium is heated to 10,000 K and most minihalos never form. On the other hand, if z_r < 20, a significant fraction (>10 percent) of all baryons have already collapsed into minihalos, and are subsequently removed from the halos by photoevaporation as the ionizing background flux builds up. We show that this process can require a significant budget of ionizing photons; exceeding the production by a straightforward extrapolation back in time of known quasar and galaxy populations by a factor of up to 10 and 3, respectively.

We present results from three-dimensional hydrodynamic simulations of the high redshift collapse of pregalactic clouds including feedback effects from a soft H2 photodissociating UV radiation field. The simulations use an Eulerian adaptive mesh refinement technique to follow the nonequilibrium chemistry of nine chemical species with cosmological initial conditions drawn from a popular Lambda-dominated cold dark matter model. The results confirm that the soft UV background can delay the cooling and collapse of small halos (~10^6 Msun). For reasonable values of the photo-dissociating flux, the H2 fraction is in equilibrium throughout most of the objects we simulate. We determine the mass threshold for collapse for a range of soft-UV fluxes and also derive a simple analytic expression. Continuing the simulations beyond the point of initial collapse demonstrates that the fraction of gas which can cool depends mostly on the virial mass of the halo and the amount of soft-UV flux, with remarkably little scatter. We parameterize this relation, for use in semi-analytic models.

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.

Radiative Transfer (RT) effects play a crucial role in the thermal history of the intergalactic medium. Here I discuss recent advances in the development of numerical methods that introduce RT to cosmological hydrodynamics. These methods can also readily be applied to time dependent problems on interstellar and galactic scales.

We review the relevance of H2 molecules for structure formation in cosmology. Molecules are important at high redshifts, when the first collapsed structures appear with typical temperatures of a few hundred Kelvin. In these chemically pristine clouds, radiative cooling is dominated H2 molecules. As a result, H2 ``astro-chemistry'' is likely to determine the epoch when the first astrophysical objects appear. We summarize results of recent three-dimensional simulations. A discussion of the effects of feedback, and implications for the reionization of the universe is also given.

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.

We study the effects of the ionizing and dissociating photons produced by PopIII objects on the surrounding intergalactic medium. We find that the typical size of a H_2 photodissociated region, R_d ~ 1-5 kpc, is smaller than the mean distance between sources at z ~ 20-30, but larger than the ionized region by a factor depending on the detailed properties of the emission spectrum. This implies that clearing of intergalactic H_2 occurs before reionization of the universe is complete. In the same redshift range, the soft-UV background in the Lyman-Werner bands, when the intergalactic H and H_2 opacity is included, is found to be J_LW ~ 1d-30 - 1d-27 erg cm^{-2} s^{-1} Hz^{-1}. This value is well below the threshold required for the negative feedback of PopIII objects on the subsequent galaxy formation to be effective in that redshift range.

In hierarchical models of structure formation, an early cosmic UV background (UVB) is produced by the small (T_vir < 10^4 K) halos that collapse before reionization. The UVB at energies below 13.6eV suppresses the formation of stars or black holes inside small halos, by photo-dissociating their only cooling agent, molecular H2. We self-consistently compute the buildup of the early UVB in Press-Schechter models, coupled with H2 photo-dissociation both in the intergalactic medium (IGM), and inside virialized halos. We find that the intergalactic H2 has a negligible effect on the UVB, both because its initial optical depth is small (tau<0.1), and because it is photo-dissociated at an early stage. If the UV sources in the first collapsed halos are stars, then their UV flux suppresses further star-formation inside small halos. This results in a pause in the buildup of the UVB, and reionization is delayed until larger halos (T_vir> 10^4 K) collapse. If the small halos host mini-quasars with hard spectra extending to approximately 1 keV, then their X-rays balance the effects of the UVB, the negative feedback does not occur, and reionization can be caused by the small halos.

The thermal history of the intergalactic medium (IGM) after reionization is to a large extent determined by photoheating. Here we demonstrate that calculations of the photoheating rate which neglect radiative transfer effects substantially underestimate the energy input during and after reionization. The neglect of radiative transfer effects results in temperatures of the IGM which are too low by a factor of two after HeII reionization. We briefly discuss implications for the absorption properties of the IGM and the distribution of baryons in shallow potential wells.

Recent observations and theoretical calculations suggest that some QSO absorption line systems may be due to gas in small dark matter halos with circular velocities on the order of 30 km/s. Additional observational evidence suggests that, in general, many absorption line systems may also be multi-phase in nature. Thus, computing the absorption lines of mini-halos, in addition to providing signatures of small halos, is a natural way to explore multi-phase behavior. The state of gas in mini-halos is strongly affected by the background UV radiation field. To address this issue a code was developed that includes many of the chemical and radiative processes found in CLOUDY and also incorporates spherically symmetric multi-wavelength radiative transfer of an isotropic field, non-equilibrium chemistry, heating, cooling and self-consistent quasi hydro-static equilibrium gas dynamics. With this code detailed simulations were conducted of gas in mini-halos using different types of background spectra. From these simulations the absorption line signatures of the gas were computed and compared with a variety of observations: high redshift metal lines, He lines and low redshift metal line systems. Based on these results the mini-halo model absorption line signatures appear to be consistent with many current observations given a sufficiently soft spectrum.

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.

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 propose that a large fraction of the QSO Lyman limit absorption systems (LLS) observed at high redshift (z > 3) originate from gas trapped in small objects, such as minihalos, that form prior to reionization. In the absence of a strong UV flux, the gas is predominantly neutral and may form clouds with HI column density NHI > 10^18 cm^-2. Due to their high densities and high HI column densities, these clouds are not destroyed by the onset of the UV background at a later time. Thus, if not disrupted by other processes, such as mergers into larger systems or `blow away' by supernovae, they will produce LLS. We show that the observed number density of LLS at high redshifts can well be reproduced by the survived `minihalos' in hierarchical clustering models such as the standard cold dark matter model. The number density of LLS in such a population increases with z even beyond the redshifts accessible to current observations and dies off quickly at z < 2. This population is distinct from other populations because the absorbing systems have small velocity widths and a close to primordial chemical composition. The existence of such a population requires that the reionization of the universe occurs late, at z < 20.

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

The minimum mass that a virialized gas cloud must have in order to be able to cool in a Hubble time is computed, using a detailed treatment of the chemistry of molecular hydrogen. With a simple model for halo profiles, we reduce the problem to that of numerically integrating a system of chemical equations. The results agree well with numerically expensive 3D simulations, and our approach has the advantage of rapidly being able to explore large regions of parameter space. The minimum baryonic mass M_b is found to be strongly redshift dependent, dropping from 10^6 solar masses at z=15 to 5000 solar masses at z=100 as molecular cooling becomes effective. For z>>100, M_b rises again, as CMB photons inhibit H_2-formation through the H^- channel. Finally, for z>>200, the H_2^+ channel for H_2-formation becomes effective, driving M_b down towards 10^3 solar masses. With a standard CDM power spectrum with sigma_8=0.7, this implies that a fraction 10^{-3} of all baryons may have formed luminous objects by z=30, which could be sufficient to reheat the universe.

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.