While magnetic fields are important in contemporary star formation, their role in primordial star formation is unknown. Magnetic fields of order 10^-16 G are produced by the Biermann battery due to the curved shocks and turbulence associated with the infall of gas into the dark matter minihalos that are the sites of formation of the first stars. These fields are rapidly amplified by a small-scale dynamo until they saturate at or near equipartition with the turbulence in the central region of the gas. Analytic results are given for the outcome of the dynamo, including the effect of compression in the collapsing gas. The mass-to-flux ratio in this gas is 2-3 times the critical value, comparable to that in contemporary star formation. Predictions of the outcomes of simulations using smooth particle hydrodynamics (SPH) and grid-based adaptive mesh refinement (AMR) are given. Because the numerical viscosity and resistivity for the standard resolution of 64 cells per Jeans length are several orders of magnitude greater than the physical values, dynamically significant magnetic fields affect a much smaller fraction of the mass in simulations than in reality. An appendix gives an analytic treatment of free-fall collapse, including that in a constant density background. Another appendix presents a new method of estimating the numerical viscosity; results are given for both SPH and grid-based codes.

Recent numerical simulations have suggested the probability of a fraction of the primordial stars to be ejected from the cluster of their origin. We explore the possibility that some of these can remain on the main sequence till the present epoch. We develop a semi-analytical model guided by results of cosmological simulations to study the mass accretion by these protostars as a function of the original stellar mass, and other parameters such as angular momentum and gravitational drag due to ambient gas. We also explore whether some of the protostars remain sufficiently low-mass and long-lived to survive to the present day. This requires that the protostars are ejected from the star forming region while their mass is less than $0.8\,M_{\odot}$. Assuming the protostars gain mass via the spherical Bondi-Hoyle accretion from the ambient medium, we show that Pop III protostars that initially form within a certain range of mass and are ejected with velocity larger than the escape velocity may survive to the present day on the main sequence. Thus, they may even be found in our Milky Way or its satellites. Our calculations also reveal that protostars that do not get ejected from the parent gas cloud accrete large amount of gas. We predict that these can become massive enough to be progenitors of black holes.

Recent numerical simulations have suggested the probability of a fraction of the primordial stars to be ejected from the cluster of their origin. We explore the possibility that some of these can remain on the main sequence till the present epoch. We develop a semi-analytical model guided by results of cosmological simulations to study the mass accretion by these protostars as a function of the original stellar mass, and other parameters such as angular momentum and gravitational drag due to ambient gas. We also explore whether some of the protostars remain sufficiently low-mass and long-lived to survive to the present day. This requires that the protostars are ejected from the star forming region while their mass is less than $0.8\,M_{\odot}$. Assuming the protostars gain mass via the spherical Bondi-Hoyle accretion from the ambient medium, we show that Pop III protostars that initially form within a certain range of mass and are ejected with velocity larger than the escape velocity may survive to the present day on the main sequence. Thus, they may even be found in our Milky Way or its satellites. Our calculations also reveal that protostars that do not get ejected from the parent gas cloud accrete large amount of gas. We predict that these can become massive enough to be progenitors of black holes.

We simulate the growth of a Population III stellar system, starting from cosmological initial conditions at z=100. We follow the formation of a minihalo and the subsequent collapse of its central gas to high densities, resolving scales as small as ~ 1 AU. Using sink particles to represent the growing protostars, we model the growth of the photodissociating and ionizing region around the first sink, continuing the simulation for ~ 5000 yr after initial protostar formation. Along with the first-forming sink, several tens of secondary sinks form before an ionization front develops around the most massive star. The resulting cluster has high rates of sink formation, ejections from the stellar disc, and sink mergers during the first ~ 2000 yr, before the onset of radiative feedback. By this time a warm ~ 5000 K phase of neutral gas has expanded to roughly the disc radius of 2000 AU, slowing mass flow onto the disc and sinks. By 5000 yr the most massive star grows to 20 M_sol, while the total stellar mass approaches 75 M_sol. Out of the ~ 40 sinks, approximately 30 are low-mass (M_* < 1 M_sol), and if the simulation had resolved smaller scales an even greater number of sinks might have formed. Thus, protostellar radiative feedback is insufficient to prevent rapid disc fragmentation and the formation of a high-member Pop III cluster before an ionization front emerges. Throughout the simulation, the majority of stellar mass is contained within the most massive stars, further implying that the Pop III initial mass function is top-heavy.

Gravitational waves (GWs) provide a revolutionary tool to investigate yet unobserved astrophysical objects. Especially the first stars, which are believed to be more massive than present-day stars, might be indirectly observable via the merger of their compact remnants. We develop a self-consistent, cosmologically representative, semi-analytical model to simulate the formation of the first stars. By extrapolating binary stellar-evolution models at 10% solar metallicity to metal-free stars, we track the individual systems until the coalescence of the compact remnants. We estimate the contribution of primordial stars to the merger rate density and to the detection rate of the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO). Owing to their higher masses, the remnants of primordial stars produce strong GW signals, even if their contribution in number is relatively small. We find a probability of $\gtrsim1\%$ that the current detection GW150914 is of primordial origin. We estimate that aLIGO will detect roughly 1 primordial BH-BH merger per year for the final design sensitivity, although this rate depends sensitively on the primordial initial mass function (IMF). Turning this around, the detection of black hole mergers with a total binary mass of $\sim 300\,\mathrm{M}_\odot$ would enable us to constrain the primordial IMF.

We explore the impact of a cosmic ray (CR) background generated by supernova explosions from the first stars on star-forming metal-free gas in a minihalo at $z\sim25$. Starting from cosmological initial conditions, we use the smoothed particle hydrodynamics code GADGET-2 to follow gas collapsing under the influence of a CR background up to densities of $n=10^{12}\,{\rm cm}^{-3}$, at which point we form sink particles. Using a suite of simulations with two sets of initial conditions and employing a range of CR background models, we follow each simulation for $5000\,$yr after the first sink forms. CRs both heat and ionise the gas, boosting ${\rm H}_2$ formation. Additional ${\rm H}_2$ enhances the cooling efficiency of the gas, allowing it to fulfil the Rees-Ostriker criterion sooner and expediting the collapse, such that each simulation reaches high densities at a different epoch. As it exits the loitering phase, the thermodynamic path of the collapsing gas converges to that seen in the absence of any CR background. By the time the gas approaches sink formation densities, the thermodynamic state of the gas is thus remarkably similar across all simulations. This leads to a robust characteristic mass that is largely independent of the CR background, of order $\sim$ a few $\times10\,{\rm M}_{\odot}$ even as the CR background strength varies by five orders of magnitude.

We investigate the impact of a cosmic X-ray background (CXB) on Population III stars forming in a minihalo at $z\simeq25$. Using the smoothed particle hydrodynamics code GADGET-2, we attain sufficient numerical resolution to follow gas collapsing into the centre of the minihalo from cosmological initial conditions up to densities of $10^{12}\,{\rm cm}^{-3}$, at which point we form sink particles. This allows us to study how the presence of a CXB affects the formation of H$_2$ and HD in the gas prior to becoming fully molecular. Using a suite of simulations for a range of possible CXB models, we follow each simulation for 5000\yr after the first sink particle forms. The CXB provides two competing effects, with X-rays both heating the gas and increasing the free electron fraction, allowing more H$_2$ to form. X-ray heating dominates below $n\sim1\,{\rm cm}^{-3}$, while the additional H$_2$ cooling becomes more important above $n\sim10^2\,{\rm cm}^{-3}$. The gas becomes optically thick to X-rays as it exits the quasi-hydrostatic `loitering phase,' such that the primary impact of the CXB is to cool the gas at intermediate densities, resulting in an earlier onset of baryonic collapse into the dark matter halo. At the highest densities, self-shielding results in similar thermodynamic behaviour across a wide range of CXB strengths. Consequently, we find that star formation is relatively insensitive to the presence of a CXB; both the number and the characteristic mass of the stars formed remains quite similar even as the strength of the CXB varies by several orders of magnitude.

We use cosmological simulations of high-redshift minihalos to investigate the effect of dark matter annihilation (DMA) on the collapse of primordial gas. We numerically investigate the evolution of the gas as it assembles in a Population III stellar disk. We find that when DMA effects are neglected, the disk undergoes multiple fragmentation events beginning at ~ 500 yr after the appearance of the first protostar. On the other hand, DMA heating and ionization of the gas speeds the initial collapse of gas to protostellar densities and also affects the stability of the developing disk against fragmentation, depending on the DM distribution. We compare the evolution when we model the DM density with an analytical DM profile which remains centrally peaked, and when we simulate the DM profile using N-body particles (the 'live' DM halo). When utilizing the analytical DM profile, DMA suppresses disk fragmentation for ~ 3500 yr after the first protostar forms, in agreement with earlier work. However, when using a 'live' DM halo, the central DM density peak is gradually flattened due to the mutual interaction between the DM and the rotating gaseous disk, reducing the effects of DMA on the gas, and enabling secondary protostars of mass ~ 1 M_sol to be formed within ~ 900 yr. These simulations demonstrate that DMA is ineffective in suppressing gas collapse and subsequent fragmentation, rendering the formation of long-lived dark stars unlikely. However, DMA effects may still be significant in the early collapse and disk formation phase of primordial gas evolution.

We perform numerical simulations of the growth of a Population III stellar system under photodissociating feedback. We start from cosmological initial conditions at z = 100, self-consistently following the formation of a minihalo at z = 15 and the subsequent collapse of its central gas to high densities. The simulations resolve scales as small as ~ 1 AU, corresponding to gas densities of 10^16 cm^-3. Using sink particles to represent the growing protostars, we evolve the stellar system for the next 5000 years. We find that this emerging stellar group accretes at an unusually low rate compared with minihalos which form at earlier times (z = 20 - 30), or with lower baryonic angular momentum. The stars in this unusual system will likely reach masses ranging from < 1 M_sun to 5 M_sun by the end of their main-sequence lifetimes, placing them in the mass range for which stars will undergo an asymptotic giant branch (AGB) phase. Based upon the simulation, we predict the rare existence of Population III stars that have survived to the present day and have been enriched by mass overflow from a previous AGB companion.

We perform a cosmological simulation in order to model the growth and evolution of Population III (Pop III) stellar systems in a range of host minihalo environments. A Pop III multiple system forms in each of the ten minihaloes, and the overall mass function is top-heavy compared to the currently observed initial mass function in the Milky Way. Using a sink particle to represent each growing protostar, we examine the binary characteristics of the multiple systems, resolving orbits on scales as small as 20 AU. We find a binary fraction of ~36%, with semi-major axes as large as 3000 AU. The distribution of orbital periods is slightly peaked at < 900 yr, while the distribution of mass ratios is relatively flat. Of all sink particles formed within the ten minihaloes, ~50% are lost to mergers with larger sinks, and ~50% of the remaining sinks are ejected from their star-forming disks. The large binary fraction may have important implications for Pop III evolution and nucleosynthesis, as well as the final fate of the first stars.

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

We numerically study the mutual interaction between dark matter (DM) and Population III (Pop III) stellar systems in order to explore the possibility of Pop III dark stars within this physical scenario. We perform a cosmological simulation, initialized at z ~ 100, which follows the evolution of gas and DM. We analyze the formation of the first minihalo at z ~ 20 and the subsequent collapse of the gas to densities of 10^12 cm^-3. We then use this simulation to initialize a set of smaller-scale `cut-out' simulations in which we further refine the DM to have spatial resolution similar to that of the gas. We test multiple DM density profiles, and we employ the sink particle method to represent the accreting star-forming region. We find that, for a range of DM configurations, the motion of the Pop III star-disk system serves to separate the positions of the protostars with respect to the DM density peak, such that there is insufficient DM to influence the formation and evolution of the protostars for more than ~ 5000 years. In addition, the star-disk system causes gravitational scattering of the central DM to lower densities, further decreasing the influence of DM over time. Any DM-powered phase of Pop III stars will thus be very short-lived for the typical multiple system, and DM will not serve to significantly prolong the life of Pop III stars.

We evaluate the effect of a supersonic relative velocity between the baryons and dark matter on the thermal and density evolution of the first gas clouds at z < 50. Through a series of cosmological simulations, initialized at z=100 with a range of relative streaming velocities and minihalo formation redshifts, we find that the typical streaming velocities will have little effect on the gas evolution. Once the collapse begins, the subsequent evolution of the gas will be nearly indistinguishable from the case of no streaming, and star formation will still proceed in the same way, with no change in the characteristic Pop III stellar masses. Reionization is expected to be dominated by halo masses of > 10^8 M_sun, for which the expected effect of streaming is negligible.

We estimate the rotation speed of Population III (Pop III) stars within a minihalo at z ~ 20 using a smoothed particle hydrodynamics (SPH) simulation, beginning from cosmological initial conditions. We follow the evolution of the primordial gas up to densities of 10^12 cm^-3. Representing the growing hydrostatic cores with accreting sink particles, we measure the velocities and angular momenta of all particles that fall onto these protostellar regions. This allows us to record the angular momentum of the sinks and estimate the rotational velocity of the Pop III stars expected to form within them. The rotation rate has important implications for the evolution of the star, the fate encountered at the end of its life, and the potential for triggering a gamma-ray burst (GRB). We find that there is sufficient angular momentum to yield rapidly rotating stars (> 1000 km s^-1, or near break-up speeds). This indicates that Pop III stars likely experienced strong rotational mixing, impacting their structure and nucleosynthetic yields. A subset of them was also likely to result in hypernova explosions, and possibly GRBs.

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

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.

We explore the implications of a possible cosmic ray (CR) background generated during the first supernova explosions that end the brief lives of massive Population III stars. We show that such a CR background could have significantly influenced the cooling and collapse of primordial gas clouds in minihaloes around redshifts of z ~ 15 - 20, provided the CR flux was sufficient to yield an ionization rate greater than about 10^-19 s^-1 near the center of the minihalo. The presence of CRs with energies less than approximately 10^7 eV would indirectly enhance the molecular cooling in these regions, and we estimate that the resulting lower temperatures in these minihaloes would yield a characteristic stellar mass as low as ~ 10 M_sun. CRs have a less pronounced effect on the cooling and collapse of primordial gas clouds inside more massive dark matter haloes with virial masses greater than approximately 10^8 M_sun at the later stages of cosmological structure formation around z ~ 10 - 15. In these clouds, even without CR flux the molecular abundance is already sufficient to allow cooling to the floor set by the temperature of the cosmic microwave background.

We discuss our radio (Australia Telescope Compact Array and Australian Long Baseline Array) and X-ray (XMM-Newton) monitoring observations of the unusual ultraluminous supernova SN 1978K in NGC 1313 at ~25 years after the explosion. SN 1978K is a rare example of a Type IIn supernova that has remained bright enough to have long-term X-ray and radio observations. The observations probe the dense medium that was ejected by the progenitor star prior to its explosion; the star might have been a Luminous Blue Variable. The radio imaging shows that the source remains compact, but it may be marginally resolved. The radio monitoring shows deviations from a smooth decay suggesting that inhomogeneities are present in the radio emitting region. It appears that a major change occurred in the mass-loss rate of the progenitor star shortly before the supernova event. The X-ray emission between 2000 and 2006 is consistent with the radiation coming from two shocks. All the X-ray data can be fit using the same model (with no systematic evolution or short-term variability) but this has a surprising requirement: the X-ray emitting regions have a very large abundance of helium. This would be consistent with the X-ray emitting shocks being located in a helium-rich layer that was ejected by the progenitor star, or helium-rich material was ejected in the supernova explosion. The unusual properties of the supernova motivated a search for an associated gamma-ray burst (GRB). We show that SN 1978K was inside the ~4 sigma error box of GRB 771029. If this association is correct, the GRB was exceptionally underluminous. However, the quality of the gamma-ray burst locations at that time was poor, and this is likely just a chance alignment.