In this paper I describe the visualization toolkit Sator, which is designed to read, analyze and visualize simulation data of the moving-mesh code Arepo. It is written in Python and employs a graphical user interface based on the Tkinter module, providing an interactive and intuitive user experience. All three snapshot formats employed by Arepo are supported, including the HDF5 format. Individual snapshot fields, for example cell coordinates, are read interactively only when they are needed. Refined fields used for plots, for example the temperature, are constructed from snapshot fields and can be introduced with only a few lines of code. Sator currently supports the generation of image slices and projections, which can be zoomed, moved and rotated with minimal computational effort. Furthermore, various phase space and line plots can be created. Due to its modular nature, additional analysis and plotting routines can be easily implemented and are under construction. The code is well documented and is intended to be free and open source, such that any member of the astrophysical community may contribute.

The first stars are believed to have formed a few hundred million years after the big bang in so-called dark matter minihalos with masses ~10^6 M_sun. Their radiation lit up the Universe for the first time, and the supernova explosions that ended their brief lives enriched the intergalactic medium with the first heavy elements. Influenced by their feedback, the first galaxies assembled in halos with masses ~10^8 M_sun, and hosted the first metal-enriched stellar populations. In this review, I summarize the theoretical progress made in the field of high-redshift star and galaxy formation since the turn of the millennium, with an emphasis on numerical simulations. These have become the method of choice to understand the multi-scale, multi-physics problem posed by structure formation in the early Universe. In the first part of the review, I focus on the formation of the first stars in minihalos - in particular the post-collapse phase, where disk fragmentation, protostellar evolution, and radiative feedback become important. I also discuss the influence of additional physical processes, such as magnetic fields and streaming velocities. In the second part of the review, I summarize the various feedback mechanisms exerted by the first stars, followed by a discussion of the first galaxies and the various physical processes that operate in them.

We present the highest-resolution three-dimensional simulation to date of the collapse of an atomic cooling halo in the early Universe. We use the moving-mesh code arepo with the primordial chemistry module introduced in Greif (2014), which evolves the chemical and thermal rate equations for over more than 20 orders of magnitude in density. Molecular hydrogen cooling is suppressed by a strong Lyman-Werner background, which facilitates the near-isothermal collapse of the gas at a temperature of about $10^4\,$K. Once the central gas cloud becomes optically thick to continuum emission, it settles into a Keplerian disc around the primary protostar. The initial mass of the protostar is about $0.1\,{\rm M}_\odot$, which is an order of magnitude higher than in minihaloes that cool via molecular hydrogen. The high accretion rate and efficient cooling of the gas catalyse the fragmentation of the disc into a small protostellar system with 5-10 members. After about 12 yr, strong gravitational interactions disrupt the disc and temporarily eject the primary protostar from the centre of the cloud. By the end of the simulation, a secondary clump has collapsed at a distance of $\simeq 150\,$au from the primary clump. If this clump undergoes a similar evolution as the first, the central gas cloud may evolve into a wide binary system. High accretion rates of both the primary and secondary clumps suggest that fragmentation is not a significant barrier for forming at least one massive black hole seed.

We investigate the collapse of primordial gas in a minihalo with three-dimensional radiation hydrodynamics simulations that accurately model the transfer of H2 line emission. For this purpose, we have implemented a multiline, multifrequency ray-tracing scheme in the moving-mesh code Arepo that is capable of adaptively refining rays based on the Healpix algorithm, as well as a hybrid equilibrium/non-equilibrium primordial chemistry solver. We find that a multifrequency treatment of the individual H2 lines is essential, since for high optical depths the smaller cross-section in the wings of the lines greatly increases the amount of energy that can escape. The influence of Doppler shifts due to bulk velocities is comparatively small, since systematic velocity differences in the cloud are typically smaller than the sound speed. During the initial collapse phase, the radially averaged escape fraction agrees relatively well with the fit of Ripamonti & Abel. However, in general it is not advisable to use a simple density-dependent fitting function, since the escape fraction depends on many factors and does not capture the suppression of density perturbations due to the diffusion of radiation. The Sobolev method overestimates the escape fraction by more than an order of magnitude, since the properties of the gas change on scales smaller than the Sobolev length.

We investigate the operation of the chemothermal instability in primordial star-forming clouds with a suite of three-dimensional, moving-mesh simulations. In line with previous studies, we find that the gas at the centre of high-redshift minihaloes becomes chemothermally unstable as three-body reactions convert the atomic hydrogen into a fully molecular gas. The competition between the increasing rate at which the gas cools and the increasing optical depth to H2 line emission creates a characteristic dip in the cooling time over the free-fall time on a scale of 100 au. As a result, the free-fall time decreases to below the sound-crossing time, and the cloud may become gravitationally unstable and fragment on a scale of a few tens of au during the initial free-fall phase. In three of the nine haloes investigated, secondary clumps condense out of the parent cloud, which will likely collapse in their own right before they are accreted by the primary clump. In the other haloes, fragmentation at such an early stage is less likely. However, given that previous simulations have shown that the infall velocity decreases substantially once the gas becomes rotationally supported, the amount of time available for perturbations to develop may be much greater than is evident from the limited period of time simulated here.

We analyze the cosmological simulations performed in the recent work of Greif et al. (2012), which followed the early growth and merger history of Pop III stars while resolving scales as small as 0.05 R_sol. This is the first set of cosmological simulations to self-consistently resolve the rotation and internal structure of Pop III protostars. We find that Pop III stars form under significant rotational support which is maintained for the duration of the simulations. The protostellar surfaces spin from ~50% to nearly 100% of Keplerian rotational velocity. These rotation rates persist after experiencing multiple stellar merger events. In the brief time period simulated (~ 10 yr), the protostars show little indication of convective instability, and their properties furthermore show little correlation with the properties of their host minihaloes. If Pop III protostars within this range of environments generally form with high degrees of rotational support, and if this rotational support is maintained for a sufficient amount of time, this has a number of crucial implications for Pop III evolution and nucleosynthesis, as well as the possibility for Pop III pair-instability supernovae, and the question of whether the first stars produced gamma-ray bursts.

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. We utilize the bright afterglow of a Pop III GRB exploding in a primordial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an existing highly-resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial temperature >10^4 K. We explore the unusual circumburst environment inside the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1+z) with values of about a few cm^{-3}. In more massive halos, corresponding to the first galaxies, the density may be larger, n>100 cm^{-3}. The resulting afterglow fluxes may be detectable with the JWST and VLA in the near-IR and radio wavebands, respectively, out to redshift z>20. We predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak fluxes from mJy to Jy at different observed times. GRBs are ideal tools for probing the metal enrichment in the early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages. We show that if the afterglow can be observed sufficiently early, analysis of the metal lines can distinguish whether the first heavy elements were produced in a PISN, or a core-collapse (Type II) SN, thus constraining the initial mass function of the first stars.

We investigate the formation of the first stars at the end of the cosmic dark ages with a suite of three-dimensional, moving mesh simulations that directly resolve the collapse of the gas beyond the formation of the first protostar at the centre of a dark matter minihalo. The simulations cover more than 25 orders of magnitude in density and have a maximum spatial resolution of 0.05 R_sun, which extends well below the radius of individual protostars and captures their interaction with the surrounding gas. In analogy to previous studies that employed sink particles, we find that the Keplerian disc around the primary protostar fragments into a number of secondary protostars, which is facilitated by H2 collisional dissociation cooling and collision-induced emission. The further evolution of the protostellar system is characterized by strong gravitational torques that transfer angular momentum between the secondary protostars formed in the disc and the surrounding gas. This leads to the migration of about half of the secondary protostars to the centre of the cloud in a free-fall time, where they merge with the primary protostar and enhance its growth to about five times the mass of the second most massive protostar. By the same token, a fraction of the protostars obtain angular momentum from other protostars via N-body interactions and migrate to higher orbits. On average, only every third protostar survives until the end of the simulation. However, the number of protostars present at any given time increases monotonically, suggesting that the system will continue to grow beyond the limited period of time simulated here.

We study how the first galaxies were assembled under feedback from the accretion onto a central black hole (BH) that is left behind by the first generation of metal-free stars through self-consistent, cosmological simulations. X-ray radiation from the accretion of gas onto BH remnants of Population III (Pop III) stars, or from high-mass X-ray binaries (HMXBs), again involving Pop III stars, influences the mode of second generation star formation. We track the evolution of the black hole accretion rate and the associated X-ray feedback starting with the death of the Pop III progenitor star inside a minihalo and following the subsequent evolution of the black hole as the minihalo grows to become an atomically cooling galaxy. We find that X-ray photoionization heating from a stellar-mass BH is able to quench further star formation in the host halo at all times before the halo enters the atomic cooling phase. X-ray radiation from a HMXB, assuming a luminosity close to the Eddington value, exerts an even stronger, and more diverse, feedback on star formation. It photoheats the gas inside the host halo, but also promotes the formation of molecular hydrogen and cooling of gas in the intergalactic medium and in nearby minihalos, leading to a net increase in the number of stars formed at early times. Our simulations further show that the radiative feedback from the first BHs may strongly suppress early BH growth, thus constraining models for the formation of supermassive BHs.

We perform three-dimensional cosmological simulations to examine the growth of metal-free, Population III (Pop III) stars under radiative feedback. We begin our simulation at z=100 and trace the evolution of gas and dark matter until the formation of the first minihalo. We then follow the collapse of the gas within the minihalo up to densities of n = 10^12 cm^-3, at which point we replace the high-density particles with a sink particle to represent the growing protostar. We model the effect of Lyman-Werner (LW) radiation emitted by the protostar, and employ a ray-tracing scheme to follow the growth of the surrounding H II region over the next 5000 yr. We find that a disk assembles around the first protostar, and that radiative feedback will not prevent further fragmentation of the disk to form multiple Pop III stars. Ionization of neutral hydrogen and photodissociation of H_2 by LW radiation leads to heating of the dense gas to several thousand Kelvin, and this warm region expands outward at the gas sound speed. Once the extent of this warm region becomes equivalent to the size of the disk, the disk mass declines while the accretion rate onto the protostars is reduced by an order of magnitude. This occurs when the largest sink has grown to ~ 20 M_sol while the second sink has grown to 7 M_sol, and we estimate the main sink will approach an asymptotic value of ~ 30 M_sol by the time it reaches the main sequence. Our simulation thus indicates that the most likely outcome is a massive Pop III binary. However, we simulate only one minihalo, and the statistical variation between minihaloes may be substantial. If Pop III stars were typically unable to grow to more than a few tens of solar masses, this would have important consequences for the occurence of pair-instability supernovae in the early Universe as well as the Pop III chemical signature in the oldest stars observable today.

It has recently been demonstrated that coherent relative streaming velocities of order 30 km / s between dark matter and gas permeated the universe on scales below a few Mpc directly after recombination. We here use a series of high-resolution moving-mesh calculations to show that these supersonic motions significantly influence the virialization of the gas in minihalos, and delay the formation of the first stars. As the gas streams into minihalos with bulk velocities around 1 km / s at z ~ 20, the additional momentum and energy input reduces the gas fractions and central densities of the halos, increasing the typical virial mass required for efficient cooling by a factor of three, and delaying Population III star formation by dz ~ 4. Since the distribution of the magnitude of the streaming velocities is narrowly peaked around a non-negligible value, this effect is important in most regions of the universe. As a consequence, the increased minimum halo mass implies a reduction of the absolute number of minihalos that can be expected to cool and form Population III stars by up to an order of magnitude. We further find that the streaming velocities increase the turbulent velocity dispersion of the minihalo gas, which could affect its ability to fragment and hence alter the mass function of the first stars.

The cosmic dark ages ended a few hundred million years after the Big Bang, when the first stars began to fill the universe with new light. It has generally been argued that these stars formed in isolation and were extremely massive - perhaps 100 times as massive as the Sun. In a recent study, Clark and collaborators showed that this picture requires revision. They demonstrated that the accretion disks that build up around Population III stars are strongly susceptible to fragmentation and that the first stars should therefore form in clusters rather than in isolation. We here use a series of high-resolution hydrodynamical simulations performed with the moving mesh code AREPO to follow up on this proposal and to study the influence of environmental parameters on the level of fragmentation. We model the collapse of five independent minihalos from cosmological initial conditions, through the runaway condensation of their central gas clouds, to the formation of the first protostar, and beyond for a further 1000 years. During this latter accretion phase, we represent the optically thick regions of protostars by sink particles. Gas accumulates rapidly in the circumstellar disk around the first protostar, fragmenting vigorously to produce a small group of protostars. After an initial burst, gravitational instability recurs periodically, forming additional protostars with masses ranging from ~ 0.1 to 10 M_sun. Although the shape, multiplicity, and normalization of the protostellar mass function depend on the details of the sink-particle algorithm, fragmentation into protostars with diverse masses occurs in all cases, confirming earlier reports of Population III stars forming in clusters. Depending on the efficiency of later accretion and merging, Population III stars may enter the main sequence in clusters and with much more diverse masses than are commonly assumed.

We introduce a prescription for the luminosity from accreting protostars into smoothed particle hydrodynamics simulation, and apply the method to simulations of five primordial minihalos generated from cosmological initial conditions. We find that accretion luminosity delays fragmentation within the halos, but does not prevent it. In halos that slowly form a low number of protostars, the accretion luminosity can reduce the number of fragments that are formed before the protostars start ionising their surroundings. However, halos that rapidly form many protostars become dominated by dynamical processes, and the effect of accretion luminosity becomes negligible. Generally the fragmentation found in the halos is highly dependent on the initial conditions. Accretion luminosity does not substantially affect the accretion rates experienced by the protostars, and is far less important than dynamical interactions, which can lead to ejections that effectively terminate the accretion. We find that the accretion rates onto the inner regions of the disks (20 AU) around the protostars are highly variable, in contrast to the constant or smoothly decreasing accretion rates currently used in models of the pre-main sequence evolution of Population III stars.

The very first stars to form in the Universe heralded an end to the cosmic dark ages and introduced new physical processes that shaped early cosmic evolution. Until now, it was thought that these stars lived short, solitary lives, with only one extremely massive star, or possibly a very wide binary system, forming in each dark matter minihalo. Here we describe numerical simulations that show that these stars were, to the contrary, often members of tight multiple systems. Our results show that the disks that formed around the first young stars were unstable to gravitational fragmentation, possibly producing small binary and higher-order systems that had separations as small as the distance between the Earth and the Sun.

Although the influence of magnetic fields is regarded as vital in the star formation process, only a few magnetohydrodynamics (MHD) simulations have been performed on this subject within the smoothed particle hydrodynamics (SPH) method. This is largely due to the unsatisfactory treatment of non-vanishing divergence of the magnetic field. Recently smoothed particle magnetohydrodynamics (SPMHD) simulations based on Euler potentials have proven to be successful in treating MHD collapse and fragmentation problems, however these methods are known to have some intrinsical difficulties. We have performed SPMHD simulations based on a traditional approach evolving the magnetic field itself using the induction equation. To account for the numerical divergence, we have chosen an approach that subtracts the effects of numerical divergence from the force equation, and additionally we employ artificial magnetic dissipation as a regularization scheme. We apply this realization of SPMHD to a widely known setup, a variation of the 'Boss & Bodenheimer standard isothermal test case', to study the impact of the magnetic fields on collapse and fragmentation. In our simulations, we concentrate on setups, where the initial magnetic field is parallel to the rotation axis. We examine different field strengths and compare our results to other findings reported in the literature. We are able to confirm specific results found elsewhere, namely the delayed onset of star formation for strong fields, accompanied by the tendency to form only single stars. We also find that the 'magnetic cushioning effect', where the magnetic field is wound up to form a 'cushion' between the binary, aids binary fragmentation in a case, where previously only formation of a single protostar was expected.

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.

One possible scenario for the formation of massive black holes (BHs) in the early Universe is from the direct collapse of primordial gas in atomic-cooling dark matter haloes in which the gas is unable to cool efficiently via molecular transitions. We study the formation of such BHs, as well as the accretion of gas onto these objects and the high energy radiation emitted in the accretion process, by carrying out cosmological radiation hydrodynamics simulations. In the absence of radiative feedback, we find an upper limit to the accretion rate onto the central object which forms from the initial collapse of hot (~ 10^4 K) gas of the order of 0.1 MSun per year. This is high enough for the formation of a supermassive star, the immediate precursor of a BH, with a mass of the order of 10^5 MSun. Assuming that a fraction of this mass goes into a BH, we track the subsequent accretion of gas onto the BH self-consistently with the high energy radiation emitted from the accretion disk. Using a ray-tracing algorithm to follow the propagation of ionizing radiation, we model in detail the evolution of the photoionized region which forms around the accreting BH. We find that BHs with masses of the order of 10^4 MSun initially accrete at close to the Eddington limit, but that the accretion rate drops to of order 1 percent of the Eddington limit after ~ 10^6 yr, due to the expansion of the gas near the BH in response to strong photoheating and radiation pressure. One signature of the accretion of gas onto BHs formed by direct collapse, as opposed to massive Pop III star formation, is an extremely high ratio of the luminosity emitted in He II 1640 to that emitted in H_alpha (or Ly_alpha); this could be detected by the James Webb Space Telescope. Finally, we briefly discuss implications for the coevolution of BHs and their host galaxies.

Using three-dimensional cosmological simulations, we study the assembly process of one of the first galaxies, with a total mass of 10^8 M_sun, collapsing at z = 10. Our main goal is to trace the transport of the heavy chemical elements produced and dispersed by a pair-instability supernova exploding in one of the minihalo progenitors. To this extent, we incorporate an efficient algorithm into our smoothed particle hydrodynamics code which approximately models turbulent mixing as a diffusion process. We study this mixing with and without the radiative feedback from Population III stars that subsequently form in neighboring minihalos. Our simulations allow us to constrain the initial conditions for second-generation star formation, within the first galaxy itself, and inside of minihalos that virialize after the supernova explosion. We find that most minihalos remain unscathed by ionizing radiation or the supernova remnant, while some are substantially photoheated and enriched to supercritical levels, likely resulting in the formation of low-mass Population III or even Population II stars. At the center of the newly formed galaxy, 10^5 M_sun of cold, dense gas uniformly enriched to 10^-3 Z_sun are in a state of collapse, suggesting that a cluster of Population II stars will form. The first galaxies, as may be detected by the James Webb Space Telescope, would therefore already contain stellar populations familiar from lower redshifts.

We investigate the formation of metal-free, Population III (Pop III), stars within a minihalo at z ~ 20 with a smoothed particle hydrodynamics (SPH) simulation, starting from cosmological initial conditions. Employing a hierarchical, zoom-in procedure, we achieve sufficient numerical resolution to follow the collapsing gas in the center of the minihalo up to number densities of 10^12 cm^-3. This allows us to study the protostellar accretion onto the initial hydrostatic core, which we represent as a growing sink particle, in improved physical detail. The accretion process, and in particular its termination, governs the final masses that were reached by the first stars. The primordial initial mass function (IMF), in turn, played an important role in determining to what extent the first stars drove early cosmic evolution. We continue our simulation for 5000 yr after the first sink particle has formed. During this time period, a disk-like configuration is assembled around the first protostar. The disk is gravitationally unstable, develops a pronounced spiral structure, and fragments into several other protostellar seeds. At the end of the simulation, a small multiple system has formed, dominated by a binary with masses ~ 40 M_Sun and ~ 10 M_Sun. If Pop III stars were to form typically in binaries or small multiples, the standard model of primordial star formation, where single, isolated stars are predicted to form in minihaloes, would have to be modified. This would have crucial consequences for the observational signature of the first stars, such as their nucleosynthetic pattern, and the gravitational-wave emission from possible Pop III black-hole binaries.

The formation of the first galaxies at redshifts z~10-15 signaled the transition from the simple initial state of the universe to one of ever increasing complexity. We here review recent progress in understanding their assembly process with numerical simulations, starting with cosmological initial conditions and modelling the detailed physics of star formation. In particular, we study the role of HD cooling in ionized primordial gas, the impact of UV radiation produced by the first stars, and the propagation of the supernova blast waves triggered at the end of their brief lives. We conclude by discussing promising observational diagnostics that will allow us to probe the properties of the first galaxies, such as their contribution to reionization and the chemical abundance pattern observed in extremely low-metallicity stars.

We construct star formation histories at redshifts z > 5 for two physically distinct populations of primordial, metal-free stars, motivated by theoretical and observational arguments that have hinted towards the existence of an intermediate stellar generation between Population III and Population I/II. Taking into account the cosmological parameters as recently revised by WMAP after three years of operation, we determine self-consistent reionization histories and discuss the resulting chemical enrichment from these early stellar generations. We find that the bulk of ionizing photons and heavy elements produced at high redshifts must have originated in Population II.5 stars, which formed out of primordial gas in haloes with virial temperatures >= 10^4 K, and had typical masses >= 10 M_sun. Classical Population III stars, formed in minihaloes and having masses >= 100 M_sun, on the other hand, had only a minor impact on reionization and early metal enrichment. Specifically, we find that only 10 per cent by mass of metal-free star formation went into Population III.

We perform three-dimensional smoothed particle hydrodynamics simulations in a realistic cosmological setting to investigate the expansion, feedback, and chemical enrichment properties of a 200 M_sun pair-instability supernova in the high-redshift universe. We find that the SN remnant propagates for a Hubble time at z = 20 to a final mass-weighted mean shock radius of 2.5 kpc (proper), roughly half the size of the HII region, and in this process sweeps up a total gas mass of 2.5*10^5 M_sun. The morphology of the shock becomes highly anisotropic once it leaves the host halo and encounters filaments and neighboring minihalos, while the bulk of the shock propagates into the voids of the intergalactic medium. The SN entirely disrupts the host halo and terminates further star formation for at least 200 Myr, while in our specific case it exerts positive mechanical feedback on neighboring minihalos by shock-compressing their cores. In contrast, we do not observe secondary star formation in the dense shell via gravitational fragmentation, due to the previous photoheating by the progenitor star. We find that cooling by metal lines is unimportant for the entire evolution of the SN remnant, while the metal-enriched, interior bubble expands adiabatically into the cavities created by the shock, and ultimately into the voids with a maximum extent similar to the final mass-weighted mean shock radius. Finally, we conclude that dark matter halos of at least M_vir > 10^8 M_sun must be assembled to recollect all components of the swept-up gas.

We investigate the properties of the first galaxies at z > 10 with highly resolved numerical simulations, starting from cosmological initial conditions and taking into account all relevant primordial chemistry and cooling. A first galaxy is characterized by the onset of atomic hydrogen cooling, once the virial temperature exceeds 10^4 K, and its ability to retain photoheated gas. We follow the complex accretion and star formation history of a 5*10^7 M_sun system by means of a detailed merger tree and derive an upper limit on the number of Population III (Pop III) stars formed prior to its assembly. We investigate the thermal and chemical evolution of infalling gas and find that partial ionization at temperatures > 10^4 K catalyses the formation of H2 and hydrogen deuteride, allowing the gas to cool to the temperature of the cosmic microwave background. Depending on the strength of radiative and chemical feedback, primordial star formation might be dominated by intermediate-mass Pop III stars formed during the assembly of the first galaxies. Accretion on to the nascent galaxy begins with hot accretion, where gas is accreted directly from the intergalactic medium and shock-heated to the virial temperature, but is quickly accompanied by a phase of cold accretion, where the gas cools in filaments before flowing into the parent halo with high velocities. The latter drives supersonic turbulence at the centre of the galaxy and could lead to very efficient chemical mixing. The onset of turbulence in the first galaxies thus likely marks the transition to Pop II star formation.

We use three-dimensional smoothed particle hydrodynamics simulations together with a dynamical ray-tracing scheme to investigate the build-up of the first H ii regions around massive Population III stars in minihaloes. We trace the highly anisotropic breakout of the ionising radiation into the intergalactic medium, allowing us to predict the resulting recombination radiation with greatly increased realism. Our simulations, together with Press-Schechter type arguments, allow us to predict the Population III contribution to the radio background at 100 MHz via bremsstrahlung and 21 cm emission. We find a global bremsstrahlung signal of around 1 mK, and a combined 21 cm signature which is an order of magnitude larger. Both might be within reach of the planned Square Kilometer Array experiment, although detection of the free-free emission is only marginal. The imprint of the first stars on the cosmic radio background might provide us with one of the few diagnostics to test the otherwise elusive minihalo star formation site.

The formation of the first galaxies at redshifts z ~ 10-15 signaled the transition from the simple initial state of the universe to one of ever increasing complexity. We here review recent progress in understanding their assembly process with numerical simulations, starting with cosmological initial conditions and modelling the detailed physics of star formation. In this context we emphasize the importance and influence of selecting appropriate initial conditions for the star formation process. We revisit the notion of a critical metallicity resulting in the transition from primordial to present-day initial mass functions and highlight its dependence on additional cooling mechanisms and the exact initial conditions. We also review recent work on the ability of dust cooling to provide the transition to present-day low-mass star formation. In particular, we highlight the extreme conditions under which this transition mechanism occurs, with violent fragmentation in dense gas resulting in tightly packed clusters.

Detection of the radiation emitted from the first galaxies at z > 10 will be made possible in the next decade, with the launch of the James Webb Space Telescope (JWST). We carry out cosmological radiation hydrodynamics simulations of Population III (Pop III) starbursts in a 10^8 M_Sun dwarf galaxy at z = 12.5. For different star formation efficiencies and stellar initial mass functions (IMFs), we calculate the luminosities and equivalent widths (EWs) of the recombination lines H_alpha, Ly_alpha, and He II 1640, under the simple assumption that the stellar population does not evolve over the first ~3 Myr of the starburst. Although only < 40 percent of the gas in the central 100 pc of the galaxy is photoionized, we find that photoheating by massive stars causes a strong dynamical response, which results in a weak correlation between luminosity emitted in hydrogen recombination lines and the total mass in stars. However, owing to the low escape fraction of He II-ionizing photons, the luminosity emitted in He II 1640 is much more strongly correlated with the total stellar mass. The ratio of the luminosity in He II 1640 to that in Ly_alpha or H_alpha is found to be a good indicator of the IMF in many cases. The ratio of observable fluxes is F_1640/F_Halpha ~ 1 for clusters of 100 M_Sun Pop III stars and F_1640/F_Halpha ~ 0.1 for clusters of 25 M_Sun Pop III stars. The EW of the He II 1640 emission line is the most reliable IMF indicator, its value varying between ~ 20 and ~ 200 angstrom for a massive and very massive Pop III IMF, respectively. Even the bright, initial stages of Pop III starbursts in the first dwarf galaxies will likely not be directly detectable by the JWST. Instead, the JWST may discover only more massive, and hence more chemically evolved, galaxies which host normal, Pop I/II, star formation.

The primary concern of this thesis is to understand the formation and properties of the first galaxies, as well as the influence of the first stars in terms of radiative, mechanical and chemical feedback. In particular, we elucidate the role of turbulence, ionizing radiation by massive Population III stars, mechanical feedback by highly energetic supernovae, and chemical enrichment. In light of the next generation of ground- and space based telescopes, we derive their observational signature in terms of recombination radiation, bremsstrahlung and 21 cm emission. We find that the cumulative 21 cm signal of the first H II regions will likely be observable by the planned SKA, while the recombination radiation from the first starbursts might be observable by JWST. These probes are essential to test the theoretical framework of the first stars and galaxies and shed some light on this elusive period of cosmic history.

We introduce a simple and efficient algorithm for diffusion in smoothed particle hydrodynamics (SPH) simulations and apply it to the problem of chemical mixing. Based on the concept of turbulent diffusion, we link the diffusivity of a pollutant to the local physical conditions and can thus resolve mixing in space and time. We apply our prescription to the evolution of an idealized supernova remnant and find that we can model the distribution of heavy elements without having to explicitly resolve hydrodynamic instabilities in the post-shock gas. Instead, the dispersal of the pollutant is implicitly modelled through its dependence on the local velocity dispersion. Our method can thus be used in any SPH simulation that investigates chemical mixing but lacks the necessary resolution on small scales. Potential applications include the enrichment of the interstellar medium in present-day galaxies, as well as the intergalactic medium at high redshifts.

The first stars were key drivers of early cosmic evolution. We review the main physical elements of the current consensus view, positing that the first stars were predominantly very massive. We continue with a discussion of important open questions that confront the standard model. Among them are uncertainties in the atomic and molecular physics of the hydrogen and helium gas, the multiplicity of stars that form in minihalos, and the possible existence of two separate modes of metal-free star formation.

The character of the first galaxies at redshifts z > 10 strongly depends on their level of pre-enrichment, which is in turn determined by the rate of primordial star formation prior to their assembly. In order for the first galaxies to remain metal-free, star formation in minihaloes must be highly suppressed, most likely by H2-dissociating Lyman-Werner (LW) radiation. We show that the build-up of such a strong LW background is hindered by two effects. Firstly, the level of the LW background is self-regulated, being produced by the Population III (Pop III) star formation which it, in turn, suppresses. Secondly, the high opacity to LW photons which is built up in the relic H II regions left by the first stars acts to diminish the global LW background. Accounting for a self-regulated LW background, we estimate a lower limit for the rate of Pop III star formation in minihaloes at z > 15. Further, we simulate the formation of a 'first galaxy' with virial temperature T > 10^4 K and total mass > 10^8 M_Sun at z > 10, and find that complete suppression of previous Pop III star formation is unlikely, with stars of > 100 M_Sun (Pop III.1) and > 10 M_Sun (Pop III.2) likely forming. Finally, we discuss the implications of these results for the nature of the first galaxies, which may be observed by future missions such as the James Webb Space Telescope.

The formation of the first generations of stars at redshifts z > 15-20 signaled the transition from the simple initial state of the universe to one of increasing complexity. We here review recent progress in understanding the assembly process of the first galaxies, starting with cosmological initial conditions and modelling the detailed physics of star formation. In particular, we study the role of HD cooling in ionized primordial gas, the impact of UV radiation produced by the first stars, and the propagation of the supernova blast waves triggered at the end of their brief lives. We conclude by discussing how the chemical abundance patterns observed in extremely low-metallicity stars allow us to probe the properties of the first stars.

The formation of the first galaxies is influenced by the radiative feedback from the first generations of stars. This feedback is manisfested by the heating and ionization of the gas which lies within the H II regions surrounding the first stars, as well as by the photodissociation of hydrogen molecules within the larger Lyman-Werner (LW) bubbles that surround these sources. Using a ray-tracing method in three-dimensional cosmological simulations, we self-consistently track the formation of, and radiative feedback from, individual stars in the course of the formation of a protogalaxy. We compute in detail the H II regions of each of these sources, as well as the regions affected by their molecule-dissociating radiation. We follow the thermal, chemical, and dynamical evolution of the primordial gas, as it becomes incorporated into the protogalaxy. While the IGM is, in general, optically thick to LW photons only over physical distances of > 30 kpc at redshifts z < 20, the high molecule fraction that is built up in relic H II regions and their increasing volume-filling fraction renders even the local IGM optically thick to LW photons over physical distances of the order of a few kiloparsecs. We find that efficient accretion onto Population III relic black holes may occur after ~ 60 Myr from the time of their formation, by which time the photo-heated relic H II region gas can cool and re-collapse into the 10^6 M_solar minihalo which hosts the black hole. Also, Pop II.5 stars, postulated to have masses of the order of 10 M_solar, can likely form from this re-collapsing relic H II region gas. Overall, we find that the local radiative feedback from the first generations of stars suppresses the star formation rate by only a factor of, at most, a few.