Source author record

Volker Bromm

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

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph astro-ph.HE astro-ph.SR gr-qc

Catalog footprint

What is connected

77works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Little Red Dots and their Progenitors from Direct Collapse Black Holes

The James Webb Space Telescope (JWST) has discovered a new population of objects, the Little Red Dots (LRDs), characterized by V-shaped spectra indicative of strong breaks around the Balmer limit and compact morphology that gave them their name. A popular explanation is that they are a sub-population of active galactic nuclei/supermassive black holes (AGN/SMBHs) predominantly found in the high-redshift Universe ($z\gtrsim3$). Similarly, direct collapse black holes (DCBHs), theorized to form from collapsing massive, extremely metal-poor gas clouds, have been invoked to explain high-redshift quasars, the most massive AGN sub-population. Here, we employ the semi-analytical code A-SLOTH to produce a population of DCBHs and compare them against observed LRD demographics and properties. Specifically, we compare the DCBH-seeded SMBH population against the standard stellar-remnant seeds and find that DCBH models agree better with observed LRD population statistics and host halo properties. Furthermore, for the most extreme and earliest LRD detections, interpreted to be systems with an AGN but little stellar component, DCBHs are able to reproduce the observed spectral shape and properties under multiple scenarios - high dust attenuation or AGN surrounded by dense gas - that have been proposed to explain the unique shape of LRD spectra. Even when super-Eddington accretion, invoked previously to explain the nature of LRDs, is enforced on stellar remnant seeds, the spectral characteristics of extreme LRDs cannot be reproduced. We emphasize the importance of gas-metallicity observations as an additional dimension besides the widely used SMBH-stellar mass ratios to further constrain the progenitors of LRDs.

preprint2026arXiv

VENUS: Two Faint Little Red Dots Separated by $\sim70\,\mathrm{pc}$ Hidden in a Single Lensed Galaxy at $z\sim7$

We report the identification of a pair of faint little red dots (LRDs), dubbed Red Eyes, in a strongly-lensed galaxy at $z\sim7$ behind the PLCKG004.5-10.5 cluster, identified from the JWST Treasury program VENUS. Red Eyes are spatially resolved on the image plane with distinct colors, while the critical curve lies far north of Red Eyes, clearly requiring two different LRDs rather than a single LRD. Red Eyes is an extremely close pair of LRDs separated by $\sim70\,\mathrm{pc}$ in the source plane with a magnification of $μ\sim20$, which consistently explains another counter-image detected to the north-west. Red Eyes is hosted in a typical star-forming galaxy with $M_{\mathrm{UV,int}}\sim -19$, but its own UV emission is very faint ($M_{\mathrm{UV,int}} \gtrsim -16$). Moreover, Red Eyes does not reside at the galaxy center but lies at an offset position of approximately one effective radius $R_{\mathrm{e}}$ away from the galaxy center. If observed without lensing, Red Eyes would appear as a typical star-forming galaxy at $z\sim 7$ with $M_{\mathrm{UV}}\sim -19$, showing no apparent LRD signatures in either morphology or SED. These results suggest that multiple off-center LRDs, similar to Red Eyes, may be commonly hidden in a typical high-$z$ star-forming galaxy. In this case, various plausible scenarios may emerge, one of which is that intermediate-mass black holes (IMBHs) with $M_\mathrm{BH}\sim10^{4\text{--}6}\,M_\odot$ may form in star clusters on a stellar disk and contribute to the growth of the central supermassive black hole via mergers, with some IMBHs detectable as luminous LRDs in a sufficiently active and massive phase.

preprint2022arXiv

Effects of stellar-mass primordial black holes on first star formation

We use cosmological hydrodynamic zoom-in simulations and semi-analytical models to study the effects of primordial black holes (PBHs) on first star formation. Our models self-consistently combine two competing effects: initial (isocurvature) perturbations induced by PBHs and BH accretion feedback. Focusing on PBHs with masses $\sim 30\ \rm M_{\odot}$, we find that the standard picture of first star formation in molecular-cooling minihaloes is not changed by PBHs, as the simulated star-forming gas clouds in the central parsec are very similar to those in the $\rm ΛCDM$ case when PBHs make up $f_{\rm PBH}\sim 10^{-4}-0.1$ of dark matter. With a dynamical friction timescale of $\sim 2-10\ \rm Myr$ when the central gas density reaches $10^{5}\ \rm cm^{-3}$, it is also unlikely that PBHs can sink into star-forming discs and affect the evolution of protostars, although they may interact with the stars during the main-sequence stage. At larger scales, PBHs tend to shift star formation to more massive haloes, and accelerate structure formation. The latter effect is stronger in regions with higher initial overdensities. For $f_{\rm PBH}\sim 10^{-4}-0.01$ (allowed by observational constraints), the collapsed mass fraction of haloes hosting Population III stars is similar (within a factor of $\sim2$ at $z\lesssim 30$) to that in $\rm ΛCDM$, implying that the impact of stellar-mass PBHs on the cosmic star formation history at $z\gtrsim 10$ is small. We also find that the Lyman-Werner photons from PBH accretion in atomic-cooling haloes may facilitate the formation of direct-collapse BHs.

preprint2022arXiv

Maximal X-ray feedback in the pre-reionization universe

X-ray feedback in the pre-reionization Universe provided one of the major energy sources for reionization and the thermal evolution of the early intergalactic medium. However, X-ray sources at high redshift have remained largely inaccessible to observations. One alternative approach to study the overall effect of X-ray feedback in the early Universe is a full cosmological simulation. Toward this goal, in this paper we create an analytic model of X-ray feedback from accretion onto supermassive black holes (SMBHs), to be used as a sub-grid model in future cosmological simulations. Our analytic model provides a relation between the mass of a dark matter halo and the SMBH it hosts, where the efficiency is governed by an energy balance argument between thermal feedback and the confining gravitational potential of the halo. To calibrate the model, we couple the halo-level recipe with the Press-Schechter halo mass function and derive global mass and energy densities. We then compare our model to various observational constraints, such as the resulting soft X-ray and IR cosmic radiation backgrounds, to test our choice of model parameters. We in particular derive model parameters that do not violate any constraints, while providing maximal X-ray feedback prior to reionization. In addition, we consider the contribution of SMBH X-ray sources to reionization and the global 21 cm absorption signal.

preprint2022arXiv

On the probability of the extremely lensed $z$=6.2 Earendel source being a Population~III star

The recent discovery of the extremely lensed Earendel object at $z=6.2$ is remarkable in that it is likely a single star or stellar multiple, observed within the first billion years of cosmic history. Depending on its mass, which is still uncertain but will soon be more tightly constrained with the James Webb Space Telescope, the Earendel star might even be a member of the first generation of stars, the so-called Population~III (Pop~III). By combining results from detailed cosmological simulations of the assembly of the first galaxies, including the enrichment of the pristine gas with heavy chemical elements, with assumptions on key stellar parameters, we quantify the probability that Earendel has indeed a Pop~III origin. We find that this probability is non-negligible throughout the mass range inferred for Earendel, specifically ranging from a few percent at the lower-mass end to near unity for some Pop~III initial mass function (IMF) models towards the high-mass end of the allowed range. For models that extend the metal-enriched IMF to $500$\,M$_\odot$, the likelihood of Earendel being a Pop~III star stays at the few to ten percent level. We discuss the implications of such a discovery for the overall endeavor to probe the hitherto so elusive first stars in the universe.

preprint2022arXiv

Probing the initial mass function of the first stars with transients

The emergence of the first, so-called Population III (Pop III), stars shaped early cosmic history in ways that crucially depends on their initial mass function (IMF). However, because of the absence of direct observational constraints, the detailed IMF remains elusive. Nevertheless, numerical simulations agree in broad terms that the first stars were typically massive and should often end their lives in violent, explosive deaths. These fates include extremely luminous pair-instability supernovae (PISNe) and bright gamma-ray bursts (GRBs), the latter arising from the collapse of rapidly rotating progenitor stars into black holes. These high-redshift transients are expected to be within the detection limits of upcoming space telescope missions, allowing to place effective constraints on the shape of the primordial IMF that is not easily accessible with other probes. This paper presents a framework to probe the Pop III IMF, utilizing the cosmological source densities of high-redshift PISNe and GRBs. Considering these transients separately could provide useful constraints on the Pop III IMF, but tighter bounds are obtainable by combining PISN and GRB counts. This combined diagnostic is more robust as it is independent of the underlying Pop III star formation rate density, an unknown prior. Future surveys promise to capture most high-redshift GRBs across the entire sky, but high-redshift PISN searches with future telescopes, e.g. Roman Space Telescope, will likely be substantially incomplete. Nevertheless, we demonstrate that even such lower bounds on the PISN count will be able to provide key constraints on the primordial IMF, in particular, if it is top-heavy or not.

preprint2022arXiv

Public Release of A-SLOTH: Ancient Stars and Local Observables by Tracing Halos

The semi-analytical model A-SLOTH (Ancient Stars and Local Observables by Tracing Halos) is the first public code that connects the formation of the first stars and galaxies to observables. After several successful projects with this model, we publish the source code and describe the public version in this paper. The model is based on dark matter merger trees that can either be generated based on Extended Press-Schechter theory or that can be imported from dark matter simulations. On top of these merger trees, A-SLOTH applies analytical recipes for baryonic physics to model the formation of both metal-free and metal-poor stars and the transition between them with unprecedented precision and fidelity. A-SLOTH samples individual stars and includes radiative, chemical, and mechanical feedback. It is calibrated based on six observables, such as the optical depth to Thomson scattering, the stellar mass of the Milky Way and its satellite galaxies, the number of extremely-metal poor stars, and the cosmic star formation rate density at high redshift. A-SLOTH has versatile applications with moderate computational requirements. It can be used to constrain the properties of the first stars and high-z galaxies based on local observables, predicts properties of the oldest and most metal-poor stars in the Milky Way, can serve as a subgrid model for larger cosmological simulations, and predicts next-generation observables of the early Universe, such as supernova rates or gravitational wave events.

preprint2020arXiv

Dust emission from the first massive galaxies

We compare observed far-infrared (FIR) galaxy spectral energy distributions (SEDs) at $z\sim6$ and theoretical predictions for first massive population II (Pop II) galaxies. Observed FIR SEDs at $z\gtrsim5$ are broad and shifted to bluer wavelengths when compared to galaxies at $z\sim3$. By implementing an analytical model for dust emission from Pop II massive galaxies,we were able to reproduce the observed behaviour as a consequence of the high energy densities and silicate-rich dust composition of high-$z$ model galaxies. As noted in a previous work, the non-blackbody nature of galaxy SEDs at $z \sim 6$ should be taken into account when interpreting measurements of FIR luminosities to avoid underestimating star formation rates.

preprint2020arXiv

Formation of the first stars and black holes

We review the current status of knowledge concerning the early phases of star formation during cosmic dawn. This includes the first generations of stars forming in the lowest mass dark matter halos in which cooling and condensation of gas with primordial composition is possible at very high redshift ($z > 20$), namely metal-free Population III stars, and the first generation of massive black holes forming at such early epochs, the so-called black hole seeds. The formation of black hole seeds as end states of the collapse of Population III stars, or via direct collapse scenarios, is discussed. In particular, special emphasis is given to the physics of supermassive stars as potential precursors of direct collapse black holes, in light of recent results of stellar evolution models, and of numerical simulations of the early stages of galaxy formation. Furthermore, we discuss the role of the cosmic radiation produced by the early generation of stars and black holes at high redshift in the process of reionization.

preprint2020arXiv

Gravitational waves from Population III binary black holes formed by dynamical capture

We use cosmological hydrodynamic simulations to study the gravitational wave (GW) signals from high-redshift binary black holes (BBHs) formed by dynamical capture (ex-situ formation channel). We in particular focus on BHs originating from the first generation of massive, metal-poor, so-called Population~III (Pop~III) stars. An alternative (in-situ) formation pathway arises in Pop~III binary stars, whose GW signature has been intensively studied. In our optimistic model, we predict a local GW event rate density for ex-situ BBHs (formed at z>4) of ~0.04 /(yr Gpc^3). This is comparable to or even higher than the conservative predictions of the rate density for in-situ BBHs ~0.01-0.1 /(yr Gpc^3), indicating that the ex-situ formation channel may be as important as the in-situ one for producing GW events. We also evaluate the detectability of our simulated GW events for selected planned GW instruments, such as the Einstein Telescope (ET). For instance, we find the all-sky detection rate with signal-to-noise ratios above 10 to be <100 per year for the xylophone configuration of ET. However, our results are highly sensitive to the sub-grid models for BBH identification and evolution, such that the GW event efficiency (rate) is reduced by a factor of 4 (20) in the pessimistic case. The ex-situ channel of Pop III BBHs deserves further investigation with better modeling of the environments around Pop III-seeded BHs.

preprint2020arXiv

Lessons on early structure formation from a mature galaxy cluster observed at cosmic noon

We demonstrate a new approach of indirectly constraining both early star and structure formation via mature galaxy clusters at cosmic noon (z~2), using the cluster XLSSC 122 as an example. With the standard Press Schechter formalism, we infer a rapid evolution of the star formation efficiency (the ratio of stellar to halo mass) from 10^-4 to 0.01 during z~20-13, based on the age distribution of stars in post-starburst galaxies of XLSSC 122, measured by HST photometry assuming no dust extinction. Here, we consider all low-mass haloes, including minihaloes, that host the first stars and galaxies (5*10^5 Msun < M_halo < 10^10 Msun). We also place new constraints on fuzzy dark matter models of m_a < 5*10^-21 eV/c^2 for the ultra-light boson mass, from the abundance of galaxies with star formation at z > 13 in XLSSC 122. Our exploratory results are consistent with existing constraints. More comprehensive results will be obtained if our approach is extended to a large sample of clusters or field post-starburst galaxies at cosmic noon, with improved modelling of halo and stellar populations.

preprint2020arXiv

Redshift Horizon for the Origins Space Telescope from Primordial Dust Emission

We explore the possibility of detecting the first galaxies with the next generation Origins Space Telescope (OST) by applying an analytical model of primordial dust emission. By analysing source densities as a function of redshift (z), and considering deep-field exposures with the Origins Survey Spectrometer, we estimate that the redshift horizon for detecting one individual source would be above z~7 for systems with dust-to-metal ratios higher than those expected for typical primeval galaxies. On the other hand, if confusion limits could be overcome, the Far-infrared Imager and Polarimeter would enable the detection of more typical fainter systems at z>7. Given the dependence of the obtained results with the properties of primeval dust, we conclude that the OST could provide important clues to the nature of the interestellar medium in the early Universe.

preprint2020arXiv

Signatures of first galaxies at FIR/sub-mm wavelengths

We explore the possibility of detecting first galaxies with a generic far-infrared/sub-millimeter telescope by applying an analytical model of primordial dust emission. As shown in previous works, galaxies at redshifts $z>7$ experience a strong negative K-correction in such a way that systems of similar masses are brighter at higher $z$. In addition, at a given mass and $z$, our model predicts that luminosity fluxes increase proportionally to the dust-to-metal ratios ($D/M$) of primeval sources. We evaluate the observability of model sources at different $z>7$ as a function of the observed survey area ($ΔΩ$) and sensitivity ($S$) of a generic instrument. Assuming $ΔΩ\sim 10~{\rm deg}^2$ and a plausible $S \sim 1~μ{\rm Jy}$ for a near future survey, we could assure the detection of at least one typical source with $D/M \sim 5\times10^{-3}$ at $z>7$. For $S\gtrsim 1~μ{\rm Jy}$ and $ΔΩ\lesssim 10~{\rm deg}^2$, higher than typical $D/M$ are required to detect at least one individual source at $z>7$. The observability of model galaxies is also affected by the size distribution of dust grains, specially towards higher $z$.

preprint2020arXiv

The role of faint population III supernovae in forming CEMP stars in ultra-faint dwarf galaxies

CEMP-no stars, a subset of carbon enhanced metal poor (CEMP) stars ($\rm [C/Fe]\geq0.7$ and $\rm [Fe/H]\lesssim-1$) have been discovered in ultra-faint dwarf (UFD) galaxies, with $M_{\rm vir} \sim 10^8$ Msun and $M_{\ast}\sim10^3-10^4$ Msun at $z=0$, as well as in the halo of the Milky Way (MW). These CEMP-no stars are local fossils that may reflect the properties of the first (Pop~III) and second (Pop~II) generation of stars. However, cosmological simulations have struggled to reproduce the observed level of carbon enhancement of the known CEMP-no stars. Here we present new cosmological hydrodynamic zoom-in simulations of isolated UFDs that achieve a gas mass resolution of $m_{\rm gas}\sim60$ Msun. We include enrichment from Pop~III faint supernovae (SNe), with $ E_{\rm SN}=0.6\times10^{51}$ erg, to understand the origin of CEMP-no stars. We confirm that Pop~III and Pop~II stars are mainly responsible for the formation of CEMP and C-normal stars respectively. New to this study, we find that a majority of CEMP-no stars in the observed UFDs and the MW halo can be explained by Pop~III SNe with normal explosion energy ($ E_{\rm SN}=1.2\times10^{51}$~erg) and Pop~II enrichment, but faint SNe might also be needed to produce CEMP-no stars with $\rm [C/Fe]\gtrsim2$, corresponding to the absolute carbon abundance of $\rm A(C)\gtrsim6.0$. Furthermore, we find that while we create CEMP-no stars with high carbon ratio $\rm [C/Fe]\approx3-4$, by adopting faint SNe, it is still challenging to reproduce CEMP-no stars with extreme level of carbon abundance of $\rm A(C)\approx 7.0-7.5$, observed both in the MW halo and UFDs.

preprint2020arXiv

When did Population III star formation end?

We construct a theoretical framework to study Population III (Pop III) star formation in the post-reionization epoch ($z\lesssim 6$) by combining cosmological simulation data with semi-analytical models. We find that due to radiative feedback (i.e. Lyman-Werner and ionizing radiation) massive haloes ($M_{\rm halo}\gtrsim 10^{9}\ \rm M_{\odot}$) are the major ($\gtrsim 90$%) hosts for potential Pop III star formation at $z\lesssim 6$, where dense pockets of metal-poor gas may survive to form Pop III stars, under inefficient mixing of metals released by supernovae. Metal mixing is the key process that determines not only when Pop III star formation ends, but also the total mass, $M_{\rm PopIII}$, of \textit{active} Pop III stars per host halo, which is a crucial parameter for direct detection and identification of Pop III hosts. Both aspects are still uncertain due to our limited knowledge of metal mixing during structure formation. Current predictions range from early termination at the end of reionization ($z\sim 5$) to continuous Pop III star formation extended to $z=0$ at a non-negligible rate $\sim 10^{-7}\ \rm M_{\odot}\ yr^{-1}\ Mpc^{-3}$, with $M_{\rm PopIII}\sim 10^{3}-10^{6}\ \rm M_{\odot}$. This leads to a broad range of redshift limits for direct detection of Pop III hosts, $z_{\rm PopIII}\sim 0.5-12.5$, with detection rates $\lesssim 0.1-20\ \rm arcmin^{-2}$, for current and future space telescopes (e.g. HST, WFIRST and JWST). Our model also predicts that the majority ($\gtrsim 90$%) of the cosmic volume is occupied by metal-free gas. Measuring the volume filling fractions of this metal-free phase can constrain metal mixing parameters and Pop III star formation.

preprint2016arXiv

Building up the Population III initial mass function from cosmological initial conditions

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.

preprint2016arXiv

Contribution of the first galaxies to the cosmic far-infrared/sub-millimeter background - I. Mean background level

We study the contribution of the first galaxies to the far-infrared/sub-millimeter (FIR/sub-mm) extragalactic background light (EBL) by implementing an analytical model for dust emission. We explore different dust models, assuming different grain size distributions and chemical compositions. According to our findings, observed re-radiated emission from dust in dwarf-size galaxies at $z \sim 10$ would peak at a wavelength of $\sim 500 μ{\rm m}$ with observed fluxes of $\sim 10^{-3} - 10^{-2}$ nJy, which is below the capabilities of current observatories. In order to be detectable, model sources at these high redshifts should exhibit luminosities of $\gtrsim 10^{12} L_{\odot}$, comparable to that of local ultra-luminous systems. The FIR/sub-mm EBL generated by primeval galaxies peaks at $\sim 500 μ{\rm m}$, with an intensity ranging from $\sim 10^{-4}$ to $10^{-3} {\rm nW \ m^{-2} \ sr^{-1}}$, depending on dust properties. These values are $\sim 3 - 4$ orders of magnitude below the absolute measured cosmic background level, suggesting that the first galaxies would not contribute significantly to the observed FIR/sub-mm EBL. Our model EBL exhibits a strong correlation with the dust-to-metal ratio, where we assume a fiducial value of $D = 0.005$, increasing almost proportionally to it. Thus, measurements of the FIR/sub-mm EBL could provide constraints on the amount of dust in the early Universe. Even if the absolute signal from primeval dust emission may be undetectable, it might still be possible to obtain information about it by exploring angular fluctuations at $\sim 500 μ{\rm m}$, close to the peak of dust emission from the first galaxies.

preprint2016arXiv

Evidence for a direct collapse black hole in the Lyman-alpha source CR7

Throughout the epoch of reionization the most luminous Lyα emitters are capable of ionizing their own local bubbles. The CR7 galaxy at $z \approx 6.6$ stands out for its combination of exceptionally bright Lyα and HeII 1640 Angstrom line emission but absence of metal lines. As a result CR7 may be the first viable candidate host of a young primordial starburst or direct collapse black hole. High-resolution spectroscopy reveals a +160 km s$^{-1}$ velocity offset between the Lyα and HeII line peaks while the spatial extent of the Lyα emitting region is $\sim 16$ kpc. The observables are indicative of an outflow signature produced by a strong central source. We present one-dimensional radiation-hydrodynamics simulations incorporating accurate Lyα feedback and ionizing radiation to investigate the nature of the CR7 source. We find that a Population III star cluster with $10^5$ K blackbody emission ionizes its environment too efficiently to generate a significant velocity offset. However, a massive black hole with a nonthermal Compton-thick spectrum is able to reproduce the Lyα signatures as a result of higher photon trapping and longer potential lifetime. For both sources, Lyα radiation pressure turns out to be dynamically important.

preprint2016arXiv

Exploring the nature of the Lyman-$α$ emitter CR7

CR7 is the brightest Lyman-$α$ emitter observed at $z>6$, which shows very strong Lyman-$α$ and HeII 1640Å line luminosities, but no metal line emission. Previous studies suggest that CR7 hosts either young primordial stars with a total stellar mass of $\sim 10^7\,\mathrm{M}_\odot$ or a black hole of $\gtrsim 10^6\,\mathrm{M}_\odot$. Here, we explore different formation scenarios for CR7 with a semianalytical model, based on the random sampling of dark matter merger trees. We are unable to reproduce the observational constraints with a primordial stellar source, given our model assumptions, due to the short stellar lifetimes and the early metal enrichment. Black holes that are the remnants of the first stars are either not massive enough, or reside in metal-polluted haloes, ruling out this possible explanation of CR7. Our models instead suggest that direct collapse black holes, which form in metal-free haloes exposed to large Lyman-Werner fluxes, are more likely the origin of CR7. However, this result is derived under optimistic assumptions and future observations are necessary to further constrain the nature of CR7.

preprint2016arXiv

Gamma-ray bursts and Population III stars

Gamma-ray bursts (GRBs) are ideal probes of the epoch of the first stars and galaxies. We review the recent theoretical understanding of the formation and evolution of the first (so-called Population III) stars, in light of their viability of providing GRB progenitors. We proceed to discuss possible unique observational signatures of such bursts, based on the current formation scenario of long GRBs. These include signatures related to the prompt emission mechanism, as well as to the afterglow radiation, where the surrounding intergalactic medium might imprint a telltale absorption spectrum. We emphasize important remaining uncertainties in our emerging theoretical framework.

preprint2016arXiv

Gravitational Waves from the Remnants of the First Stars

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.

preprint2016arXiv

Lyman-alpha radiation hydrodynamics of galactic winds before cosmic reionization

The dynamical impact of Lyman-alpha (Lyα) radiation pressure on galaxy formation depends on the rate and duration of momentum transfer between Lyα photons and neutral hydrogen gas. Although photon trapping has the potential to multiply the effective force, ionizing radiation from stellar sources may relieve the Lyα pressure before appreciably affecting the kinematics of the host galaxy or efficiently coupling Lyα photons to the outflow. We present self-consistent Lyα radiation-hydrodynamics simulations of high-$z$ galaxy environments by coupling the Cosmic Lyα Transfer code (COLT) with spherically symmetric Lagrangian frame hydrodynamics. The accurate but computationally expensive Monte-Carlo radiative transfer calculations are feasible under the one-dimensional approximation. The initial starburst drives an expanding shell of gas from the centre and in certain cases Lyα feedback significantly enhances the shell velocity. Radiative feedback alone is capable of ejecting baryons into the intergalactic medium (IGM) for protogalaxies with a virial mass of $M_{\rm vir} \lesssim 10^8~{\rm M}_\odot$. We compare the Lyα signatures of Population III stars with $10^5$ K blackbody emission to that of direct collapse black holes with a nonthermal Compton-thick spectrum and find substantial differences if the Lyα spectra are shaped by gas pushed by Lyα radiation-driven winds. For both sources, the flux emerging from the galaxy is reprocessed by the IGM such that the observed Lyα luminosity is reduced significantly and the time-averaged velocity offset of the Lyα peak is shifted redward.

preprint2016arXiv

The First Stars: formation under cosmic ray feedback

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.

preprint2016arXiv

Toward $ab\,initio$ extremely metal poor stars

Extremely metal poor stars have been the focus of much recent attention owing to the expectation that their chemical abundances can shed light on the metal and dust yields of the earliest supernovae. We present our most realistic simulation to date of the astrophysical pathway to the first metal enriched stars. We simulate the radiative and supernova hydrodynamic feedback of a $60\,M_\odot$ Population III star starting from cosmological initial conditions realizing Gaussian density fluctuations. We follow the gravitational hydrodynamics of the supernova remnant at high spatial resolution through its freely-expanding, adiabatic, and radiative phases, until gas, now metal-enriched, has resumed runaway gravitational collapse. Our findings are surprising: while the Population III progenitor exploded with a low energy of $10^{51}\,\text{erg}$ and injected an ample metal mass of $6\,M_\odot$, the first cloud to collapse after the supernova explosion is a dense surviving primordial cloud on which the supernova blastwave deposited metals only superficially, in a thin, unresolved layer. The first metal-enriched stars can form at a very low metallicity, of only $2-5\times10^{-4}\,Z_\odot$, and can inherit the parent cloud's highly elliptical, radially extended orbit in the dark matter gravitational potential.

preprint2015arXiv

Abundance anomalies in metal-poor stars from Population III supernova ejecta hydrodynamics

We present a simulation of the long-term evolution of a Population III supernova remnant in a cosmological minihalo. Employing passive Lagrangian tracer particles, we investigate how chemical stratification and anisotropy in the explosion can affect the abundances of the first low-mass, metal-enriched stars. We find that reverse shock heating can leave the inner mass shells at entropies too high to cool, leading to carbon-enhancement in the re-collapsing gas. This hydrodynamic selection effect could explain the observed incidence of carbon-enhanced metal-poor (CEMP) stars at low metallicity. We further explore how anisotropic ejecta distributions, recently seen in direct numerical simulations of core-collapse explosions, may translate to abundances in metal-poor stars. We find that some of the observed scatter in the Population II abundance ratios can be explained by an incomplete mixing of supernova ejecta, even in the case of only one contributing enrichment event. We demonstrate that the customary hypothesis of fully-mixed ejecta clearly fails if post-explosion hydrodynamics prefers the recycling of some nucleosynthetic products over others. Furthermore, to fully exploit the stellar-archaeological program of constraining the Pop III initial mass function from the observed Pop II abundances, considering these hydrodynamical transport effects is crucial. We discuss applications to the rich chemical structure of ultra-faint dwarf satellite galaxies, to be probed in unprecedented detail with upcoming spectroscopic surveys.

preprint2015arXiv

Confined Population III Enrichment and the Prospects for Prompt Second-Generation Star Formation

It is widely recognized that nucleosynthetic output of the first, Population III supernovae was a catalyst defining the character of subsequent stellar generations. Most of the work on the earliest enrichment was carried out assuming that the first stars were extremely massive and that the associated supernovae were unusually energetic, enough to completely unbind the baryons in the host cosmic minihalo and disperse the synthesized metals into the intergalactic medium. Recent work, however, suggests that the first stars may in fact have been somewhat less massive, with a characteristic mass scale of a few tens of solar masses. We present a cosmological simulation following the transport of the metals synthesized in a Population III supernova assuming that it had an energy of 1e51 ergs, compatible with standard Type II supernovae. A young supernova remnant is inserted in the first star's relic HII region in the free expansion phase and is followed for 40 Myr employing adaptive mesh refinement and Lagrangian tracer particle techniques. The supernova remnant remains partially trapped within the minihalo and the thin snowplow shell develops pronounced instability and fingering. Roughly half of the ejecta turn around and fall back toward the center of the halo, with 1% of the ejecta reaching the center in 30 kyr and 10% in 10 Myr. The average metallicity of the combined returning ejecta and the pristine filaments feeding into the halo center from the cosmic web is 0.001 - 0.01 Z_sun, but the two remain unmixed until accreting onto the central hydrostatic core that is unresolved at the end of the simulation. We conclude that if Population III stars had less extreme masses, they promptly enriched the host minihalos with metals and triggered Population II star formation.

preprint2015arXiv

Constraining the primordial initial mass function with stellar archaeology

We present a new near-field cosmological probe of the initial mass function (IMF) of the first stars. Specifically, we constrain the lower-mass limit of the Population III (Pop III) IMF with the total number of stars in large, unbiased surveys of the Milky Way. We model the early star formation history in a Milky Way-like halo with a semi-analytic approach, based on Monte-Carlo sampling of dark matter merger trees, combined with a treatment of the most important feedback mechanisms. Assuming a logarithmically flat Pop III IMF and varying its low mass limit, we derive the number of expected survivors of these first stars, using them to estimate the probability to detect any such Pop III fossil in stellar archaeological surveys. Following our analysis, the most promising region to find possible Pop III survivors is the stellar halo of the Milky Way, which is the best target for future surveys. We find that if no genuine Pop III survivor is detected in a sample size of $4 \times 10^6$ ($2 \times 10^7$) halo stars with well-controlled selection effects, then we can exclude the hypothesis that the primordial IMF extended down below $0.8 M_\odot$ at a confidence level of 68% (99%). With the sample size of the Hamburg/ESO survey, we can tentatively exclude Pop III stars with masses below $0.65 M_\odot$ with a confidence level of 95%, although this is subject to significant uncertainties. To fully harness the potential of our approach, future large surveys are needed that employ uniform, unbiased selection strategies for high-resolution spectroscopic follow-up.

preprint2015arXiv

Cosmological Impact of Population III Binaries

We present the results of the stellar feedback from Pop III binaries by employing improved, more realistic Pop III evolutionary stellar models. To facilitate a meaningful comparison, we consider a fixed mass of 60 solar masses (Msun) incorporated in Pop III stars, either contained in a single star, or split up in binary stars of 30 Msun each or an asymmetric case of one 45 Msun and one 15 Msun star. Whereas the sizes of the resulting HII regions are comparable across all cases, the HeIII regions around binary stars are significantly smaller than that of the single star. Consequently, the He$^{+}$ 1640 angstrom recombination line is expected to become much weaker. Supernova feedback exhibits great variety due to the uncertainty in possible explosion pathways. If at least one of the component stars dies as a hypernova about ten times more energetic than conventional core-collapse supernovae, the gas inside the host minihalo is effectively blown out, chemically enriching the intergalactic medium (IGM) to an average metallicity of $10^{-4}-10^{-3}$ solar metallicity (Zsun), out to $\sim 2$ kpc. The single star, however, is more likely to collapse into a black hole, accompanied by at most very weak explosions. The effectiveness of early chemical enrichment would thus be significantly reduced, in difference from the lower mass binary stars, where at least one component is likely to contribute to heavy element production and dispersal. Important new feedback physics is also introduced if close binaries can form high-mass x-ray binaries, leading to the pre-heating and -ionization of the IGM beyond the extent of the stellar HII regions.

preprint2015arXiv

Metal Transport and Chemical Heterogeneity in Early Star Forming Systems

To constrain the properties of the first stars with the chemical abundance patterns observed in metal-poor stars, one must identify any non-trivial effects that the hydrodynamics of metal dispersal can imprint on the abundances. We use realistic cosmological hydrodynamic simulations to quantify the distribution of metals resulting from one Population III supernova and from a small number of such supernovae exploding in close succession. Overall, supernova ejecta are highly inhomogeneously dispersed throughout the simulations. When the supernova bubbles collapse, quasi-virialized metal-enriched clouds, fed by fallback from the bubbles and by streaming of metal-free gas from the cosmic web, grow in the centers of the dark matter halos. Partial turbulent homogenization on scales resolved in the simulation is observed only in the densest clouds where the vortical time scales are short enough to ensure true homogenization on subgrid scales. However, the abundances in the clouds differ from the gross yields of the supernovae. Continuing the simulations until the cloud have gone into gravitational collapse, we predict that the abundances in second-generation stars will be deficient in the innermost mass shells of the supernova (if only one has exploded) or in the ejecta of the latest supernovae (when multiple have exploded). This indicates that hydrodynamics gives rise to biases complicating the identification of nucleosynthetic sources in the chemical abundance spaces of the surviving stars.

preprint2015arXiv

On the Physical Requirements for a Pre-Reionization Origin of the Unresolved Near-Infrared Background

The study of the Cosmic Near-Infrared Background (CIB) light after subtraction of resolved sources can push the limits of current observations and infer the level of galaxy and black hole activity in the early universe. However, disentangling the relative contribution from low- and high-redshift sources is not trivial. Spatial fluctuations of the CIB exhibit a clustering excess at angular scales $\sim 1^\circ$ whose origin has not been conclusively identified. We explore the likelihood that this signal is dominated by emission from galaxies and accreting black holes in the early Universe. We find that, if the first small mass galaxies have a normal IMF, the light of their ageing stars (fossils) integrated over cosmic time contributes a comparable amount to the CIB as their pre-reionization progenitors. However, the measured fluctuation signal is too large to be produced by galaxies at redshifts $z>8$ unless their star formation efficiencies are much larger than those inferred from the observed Lyman-dropout population. In order to produce the observed level of CIB fluctuation without violating constraints from galaxy counts and the electron optical depth of the IGM, minihalos at $z>12$ must form stars with efficiency $f_\star \gtrsim 0.1$ and, although a top-heavy IMF is preferred, have a very low escape fraction of ionizing radiation, $f_{\rm esc}<0.01$. If instead the CIB fluctuations are produced by high-$z$ black holes, one requires vigorous accretion in the early universe reaching $ρ_{\rm acc} \gtrsim 10^5M_\odot{\rm Mpc^{-3}}$ by $z\simeq 10$. This growth must stop by $z \sim 6$ and be significantly obscured not to overproduce the soft cosmic X-ray background (CXB) and its observed coherence with the CIB. We therefore find the range of suitable possibilities at high-$z$ to be narrow, but could possibly be widened by including additional physics and evolution at those epochs.

preprint2015arXiv

Preserving chemical signatures of primordial star formation in the first low-mass stars

We model early star forming regions and their chemical enrichment by Population III (Pop III) supernovae with nucleosynthetic yields featuring high [C/Fe] ratios and pair-instability supernova (PISN) signatures. We aim to test how well these chemical abundance signatures are preserved in the gas prior to forming the first long-lived low-mass stars (or second-generation stars). Our results show that second-generation stars can retain the nucleosynthetic signature of their Pop III progenitors, even in the presence of nucleosynthetically normal Pop III core-collapse supernovae. We find that carbon-enhanced metal-poor stars are likely second-generation stars that form in minihaloes. Furthermore, it is likely that the majority of Pop III supernovae produce high [C/Fe] yields. In contrast, metals ejected by a PISN are not concentrated in the first star forming haloes, which may explain the absence of observed PISN signatures in metal-poor stars. We also find that unique Pop III abundance signatures in the gas are quickly wiped out by the emergence of Pop II supernovae. We caution that the observed fractions of stars with Pop III signatures cannot be directly interpreted as the fraction of Pop III stars producing that signature. Such interpretations require modelling the metal enrichment process prior to the second-generation stars' formation, including results from simulations of metal mixing. The full potential of stellar archaeology can likely be reached in ultra-faint dwarf galaxies, where the simple formation history may allow for straightforward identification of second-generation stars.

preprint2015arXiv

Star formation in the first galaxies - III. Formation, evolution, and characteristics of the first stellar cluster

We simulate the formation of a low metallicity (0.01 Zsun) stellar cluster in a dwarf galaxy at redshift z~14. Beginning with cosmological initial conditions, the simulation utilizes adaptive mesh refinement and sink particles to follow the collapse and evolution of gas past the opacity limit for fragmentation, thus resolving the formation of individual protostellar cores. A time- and location-dependent protostellar radiation field, which heats the gas by absorption on dust, is computed by integration of protostellar evolutionary tracks with the MESA code. The simulation also includes a robust non-equilibrium chemical network that self-consistently treats gas thermodynamics and dust-gas coupling. The system is evolved for 18 kyr after the first protostellar source has formed. In this time span, 30 sink particles representing protostellar cores form with a total mass of 81 Msun. Their masses range from ~0.1 Msun to 14.4 Msun with a median mass ~0.5-1 Msun. Massive protostars grow by competitive accretion while lower-mass protostars are stunted in growth by close encounters and many-body ejections. In the regime explored here, the characteristic mass scale is determined by the temperature floor set by the cosmic microwave background and by the onset of efficient dust-gas coupling. It seems unlikely that host galaxies of the first bursts of metal-enriched star formation will be detectable with the James Webb Space Telescope or other next-generation infrared observatories. Instead, the most promising access route to the dawn of cosmic star formation may lie in the scrutiny of metal-poor, ancient stellar populations in the Galactic neighborhood. The observable targets that correspond to the system simulated here are ultra-faint dwarf satellite galaxies such as Bootes II, Segue I and II, and Willman I.

preprint2015arXiv

The first galaxies: simulating their feedback-regulated assembly

We investigate the formation of a galaxy reaching a virial mass of $~ 10^8$ solar mass at $z=10$ by carrying out a zoomed radiation-hydrodynamical cosmological simulation. This simulation traces Population~III (Pop~III) star formation, characterized by a modestly top-heavy initial mass function (IMF), and considers stellar feedback such as photoionization heating from Pop III and Population~II (Pop~II) stars, mechanical and chemical feedback from supernovae (SNe), and X-ray feedback from accreting black holes (BHs) and high-mass X-ray binaries (HMXBs). We self-consistently impose a transition in star formation mode from top-heavy Pop III to low-mass Pop~II, and find that the star formation rate in the computational box is dominated by Pop~III until $z=13$, and by Pop~II thereafter. The simulated galaxy experiences bursty star formation, with a substantially reduced gas content due to photoionization heating from Pop~III and Pop~II stars, together with SN feedback. All the gas within the simulated galaxy is metal-enriched above $10^{-5}$ solar, such that there are no remaining pockets of primordial gas. The simulated galaxy has an estimated observed flux of $~10^{-3} nJy$, which is too low to be detected by the James Webb Space Telescope (JWST) without strong lensing amplification. We also show that our simulated galaxy is similar in terms of stellar mass to Segue 2, the least luminous dwarf known in the Local Group.

preprint2015arXiv

The First Stars: formation under X-ray feedback

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.

preprint2015arXiv

The Lyman-alpha signature of the first galaxies

We present the Cosmic Lyman-$α$ Transfer code (COLT), a massively parallel Monte-Carlo radiative transfer code, to simulate Lyman-$α$ (Ly$α$) resonant scattering through neutral hydrogen as a probe of the first galaxies. We explore the interaction of centrally produced Ly$α$ radiation with the host galactic environment. Ly$α$ photons emitted from the luminous starburst region escape with characteristic features in the line profile depending on the density distribution, ionization structure, and bulk velocity fields. For example, anisotropic ionization exhibits a tall peak close to line centre with a skewed tail that drops off gradually. Idealized models of first galaxies explore the effect of mass, anisotropic H II regions, and radiation pressure driven winds on Ly$α$ observables. We employ mesh refinement to resolve critical structures. We also post-process an ab initio cosmological simulation and examine images captured at various escape distances within the 1 Mpc$^3$ comoving volume. Finally, we discuss the emergent spectra and surface brightness profiles of these objects in the context of high-$z$ observations. The first galaxies will likely be observed through the red damping wing of the Ly$α$ line. Observations will be biased toward galaxies with an intrinsic red peak located far from line centre that reside in extensive H II super bubbles, which allows Hubble flow to sufficiently redshift photons away from line centre and facilitate transmission through the intergalactic medium (IGM). Even with gravitational lensing to boost the luminosity this preliminary work indicates that Ly$α$ emission from stellar clusters within haloes of $M_{\rm vir}<10^9~{\rm M}_\odot$ is generally too faint to be detected by the James Webb Space Telescope (JWST).

preprint2014arXiv

Dwarf Spheroidal Satellite Formation in a Reionized Local Group

Dwarf spheroidal satellite galaxies have emerged a powerful probe of small-scale dark matter clustering and of cosmic reionization. They exhibit structural and chemical continuity with dwarf irregular galaxies in the field and with spheroidal galaxies in high-density environments. By combining empirical constraints derived for star formation at low gas column densities and metallicities in the local universe with a model for dark matter and baryonic mass assembly, we provide an analytical description of how the dwarf spheroidals acquired their stellar content. Their progenitors formed stars until the gas content, initially reduced from the cosmic average by the thermal pressure of the reionized intergalactic medium, was finally ram pressure stripped during the progenitors' accretion on to the host galaxy. Dwarf spheroidal satellites of differing luminosities seem to share very similar most massive progenitor histories that reach thresholds for gas cooling by atomic line emission at epochs at which the Lagrangian volume of the Local Group should have been reionized. We hypothesize that dwarf spheroidals formed the bulk of their stars in partially rotationally supported HI disks in a reionized universe. This model provides an explanation for the "common mass scale" relation and reproduces the empirical luminosity-size and luminosity-metallicity relations. Explosive feedback phenomena, such as outflows driven by the concerted action of supernovae, need not have been significant in the dwarf spheroidals' formation. We further speculate that the true pre-reionization fossils should exhibit a structure distinct from that of the dwarf spheroidals, e.g., in the form of dense isolated or nuclear star clusters.

preprint2014arXiv

Formation of the First Low-Mass Stars from Cosmological Initial Conditions

We simulate the formation of a metal-poor (10^-2 Zsun) stellar cluster in one of the first galaxies to form in the early Universe, specifically a high-redshift atomic cooling halo (z~14). This is the first calculation that resolves the formation of individual metal-enriched stars in simulations starting from realistic cosmological initial conditions. We follow the evolution of a single dense clump among several in the parent halo. The clump forms a cluster of ~40 stars and sub-stellar objects within 7000 years and could continue forming stars ~5 times longer. Protostellar dust heating has a negligible effect on the star formation efficiency, at least during the early evolutionary stages, but it moderately suppresses gaseous fragmentation and brown dwarf formation. We observe fragmentation in thin gaseous filaments and sustained accretion in larger, rotating structures as well as ejections by binary interactions. The stellar initial mass function above 0.1 Msun, evaluated after ~10^4 years of fragmentation and accretion, seems in agreement with the recent measurement in ultra-faint dwarf spheroidal Galactic satellites of Geha et al. (2013).

preprint2014arXiv

Radiative Feedback from high mass X-ray binaries on the formation of the first galaxies and early reionization

Recent work suggests that the first generation of stars, the so-called Population III (Pop III), could have formed primarily in binaries or as members of small multiple systems. Here we investigate the impact of X-ray feedback from High-Mass X-ray Binaries (HMXBs) left behind in stellar binary systems after the primary forms a black hole (BH), accreting gas at a high rate from the companion, a process that is thought to be favored at the low metallicities characteristic of high-redshift gas. Thanks to their large mean free path, X-rays are capable of preionizing and preheating the gas in the intergalactic medium (IGM) and in haloes long before the reionization of the Universe is complete, and thus could have strongly affected the formation of subsequent generations of stars as well as reionization. We have carried out zoomed hydrodynamical cosmological simulations of minihaloes, accounting for the formation of Pop III stars and their collapse into BHs and HMXBs, and the associated radiation-hydrodynamic feedback from UV and X-ray photons. We find no strong net feedback from HMXBs on the simulated star formation history. On the other hand, the preheating of the IGM by HMXBs leads to a strong suppression of small-scale structures and significantly lowers the recombination rate in the IGM, thus yielding a net positive feedback on reionization. We further show that X-ray feedback from HMXBs can augment the ionizing feedback from the Pop III progenitor stars to suppress gas accretion onto the first BHs, limiting their growth into supermassive BHs. Finally, we show that X-ray ionization by HMXBs leaves distinct signatures in the properties of the high-redshift hydrogen that may be probed in upcoming observations of the redshifted 21cm spin-flip line.

preprint2014arXiv

Recovery from population III supernova explosions and the onset of second generation star formation

We use cosmological simulations to assess how the explosion of the first stars in supernovae (SNe) influences early cosmic history. Specifically, we investigate the impact by SNe on the host systems for Population~III (Pop~III) star formation and explore its dependence on halo environment and Pop~III progenitor mass. We then trace the evolution of the enriched gas until conditions are met to trigger second-generation star formation. To this extent, we quantify the recovery timescale, which measures the time delay between a Pop~III SN explosion and the appearance of cold, dense gas, out of which second-generation stars can form. We find that this timescale is highly sensitive to the Pop~III progenitor mass, and less so to the halo environment. For more massive progenitors, including those exploding in pair instability SNe, second-generation star formation is delayed significantly, for up to a Hubble time. The dependence of the recovery time on the mass of the SN progenitor is mainly due to the ionizing impact of the progenitor star. Photoionization heating increases the gas pressure and initiates a hydrodynamical response that reduces the central gas density, an effect that is stronger in more massive. The gas around lower mass Pop~III stars remains denser and hence the SN remnants cool more rapidly, facilitating the subsequent re-condensation of the gas and formation of a second generation of stars. In most cases, the second-generation stars are already metal-enriched to ~2-5 X 10^{-4}\zsun, thus belonging to Population~II. The recovery timescale is a key quantity governing the nature of the first galaxies, able to host low-mass, long-lived stellar systems. These in turn are the target of future deep-field campaigns with the James Webb Space Telescope.

preprint2014arXiv

The Chemical Imprint of Silicate Dust on the Most Metal-Poor Stars

We investigate the impact of dust-induced gas fragmentation on the formation of the first low-mass, metal-poor stars (< 1Msun) in the early universe. Previous work has shown the existence of a critical dust-to-gas ratio, below which dust thermal cooling cannot cause gas fragmentation. Assuming the first dust is silicon-based, we compute critical dust-to-gas ratios and associated critical silicon abundances ([Si/H]crit). At the density and temperature associated with protostellar disks, we find that a standard Milky Way grain size distribution gives [Si/H]crit = -4.5 +/- 0.1, while smaller grain sizes created in a supernova reverse shock give [Si/H]crit = -5.3 +/- 0.1. Other environments are not dense enough to be influenced by dust cooling. We test the silicate dust cooling theory by comparing to silicon abundances observed in the most iron-poor stars ([Fe/H] < -4.0). Several stars have silicon abundances low enough to rule out dust-induced gas fragmentation with a standard grain size distribution. Moreover, two of these stars have such low silicon abundances that even dust with a shocked grain size distribution cannot explain their formation. Adding small amounts of carbon dust does not significantly change these conclusions. Additionally, we find that these stars exhibit either high carbon with low silicon abundances or the reverse. A silicate dust scenario thus suggests that the earliest low-mass star formation in the most metal-poor regime may have proceeded through two distinct cooling pathways: fine structure line cooling and dust cooling. This naturally explains both the carbon-rich and carbon-normal stars at extremely low [Fe/H].

preprint2014arXiv

The First Stars: A Low-Mass Formation Mode

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.

preprint2014arXiv

The Mutual Interaction Between Population III Stars and Self-Annihilating Dark Matter

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.

preprint2013arXiv

Constraining the Statistics of Population III Binaries

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.

preprint2013arXiv

Formation of the First Stars

Understanding the formation of the first stars is one of the frontier topics in modern astrophysics and cosmology. Their emergence signaled the end of the cosmic dark ages, a few hundred million years after the Big Bang, leading to a fundamental transformation of the early Universe through the production of ionizing photons and the initial enrichment with heavy chemical elements. We here review the state of our knowledge, separating the well understood elements of our emerging picture from those where more work is required. Primordial star formation is unique in that its initial conditions can be directly inferred from the Lambda Cold Dark Matter (LCDM) model of cosmological structure formation. Combined with gas cooling that is mediated via molecular hydrogen, one can robustly identify the regions of primordial star formation, the so-called minihalos, having total masses of ~10^6 M_sun and collapsing at redshifts z~20-30. Within this framework, a number of studies have defined a preliminary standard model, with the main result that the first stars were predominantly massive. This model has recently been modified to include a ubiquitous mode of fragmentation in the protostellar disks, such that the typical outcome of primordial star formation may be the formation of a binary or small multiple stellar system. We will also discuss extensions to this standard picture due to the presence of dynamically significant magnetic fields, of heating from self-annihalating WIMP dark matter, or cosmic rays. We conclude by discussing possible strategies to empirically test our theoretical models.

preprint2013arXiv

On the operation of the chemothermal instability in primordial star-forming clouds

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.

preprint2013arXiv

Star Formation in the First Galaxies - II: Clustered Star Formation and the Influence of Metal Line Cooling

Population III stars are believed to have been more massive than typical stars today and to have formed in relative isolation. The thermodynamic impact of metals is expected to induce a transition leading to clustered, low-mass Population II star formation. In this work, we present results from three cosmological simulations, only differing in gas metallicity, that focus on the impact of metal fine-structure line cooling on the formation of stellar clusters in a high-redshift atomic cooling halo. Introduction of sink particles allows us to follow the process of gas hydrodynamics and accretion onto cluster stars for 4 Myr corresponding to multiple local free-fall times. At metallicities at least $10^{-3}\, Z_{\odot}$, gas is able to reach the CMB temperature floor and fragment pervasively resulting in a stellar cluster of size $\sim1$ pc and total mass $\sim1000\, M_{\odot}$. The masses of individual sink particles vary, but are typically $\sim100\, M_{\odot}$, consistent with the Jeans mass when gas cools to the CMB temperature, though some solar mass fragments are also produced. At the low metallicity of $10^{-4}\, Z_{\odot}$, fragmentation is completely suppressed on scales greater than 0.01 pc and total stellar mass is lower by a factor of 3 than in the higher metallicity simulations. The sink particle accretion rates, and thus their masses, are determined by the mass of the gravitationally unstable gas cloud and the prolonged gas accretion over many Myr. The simulations thus exhibit features of both monolithic collapse and competitive accretion. Even considering possible dust induced fragmentation that would occur at higher densities, the formation of a bona fide stellar cluster seems to require metal line cooling and metallicities of at least $10^{-3}\, Z_{\odot}$.

preprint2013arXiv

The First Galaxies: Assembly under Radiative Feedback from the First Stars

We investigate how radiative feedback from the first stars affects the assembly of the first dwarf galaxies. We perform cosmological zoomed SPH simulations of a dwarf galaxy assembling inside a halo of virial mass 10^9 solar at z = 10. The simulations follow the non-equilibrium chemistry/cooling of primordial gas and the conversion of the gas into metal-free stars. To quantify the radiative feedback, we compare a simulation in which stars emit both molecular hydrogen dissociating and hydrogen/helium ionizing radiation with a simulation in which stars emit only dissociating radiation, and with a simulation in which stars remain dark. Photodissociation and -ionization exert a strong negative feedback on the assembly of the simulated galaxy. Gas condensation is strongly impeded, and star formation is strongly suppressed in comparison with the simulation in which stars remain dark. The feedback on the gas implies a suppression of the central dark matter densities in the minihalo progenitor by factors of up to a few, which is a significant deviation from the singular isothermal density profile characterizing the dark matter distribution in the absence of radiative feedback. The evolution of gas densities, star formation rates, and the distribution of dark matter becomes insensitive to the inclusion of dissociating radiation in the late stages of the minihalo assembly, and it becomes insensitive to the inclusion of ionizing radiation once the minihalo turns into an atomically cooling galaxy. The formation of an extended disk inside the dwarf galaxy is a robust outcome not affected by the inclusion of radiation. We estimate that dwarf galaxies such as simulated here will be among the faintest galaxies the upcoming James Webb Space Telescope will detect. Our conclusions are subject to our neglect of feedback from supernovae and chemical enrichment as well as to cosmic variance. [abridged]

preprint2012arXiv

Chemical Signatures of the First Galaxies: Criteria for One-Shot Enrichment

We utilize metal-poor stars in the local, ultra-faint dwarf galaxies (UFDs; L_tot < 10^5 L_sun) to empirically constrain the formation process of the first galaxies. Since UFDs have much simpler star formation histories than the halo of the Milky Way, their stellar populations should preserve the fossil record of the first supernova (SN) explosions in their long-lived, low-mass stars. Guided by recent hydrodynamical simulations of first galaxy formation, we develop a set of stellar abundance signatures that characterize the nucleosynthetic history of such an early system if it was observed in the present-day universe. Specifically, we argue that the first galaxies are the product of chemical "one-shot" events, where only one (long-lived) stellar generation forms after the first, Population III, SN explosions. Our abundance criteria thus constrain the strength of negative feedback effects inside the first galaxies. We compare the stellar content of UFDs with these one-shot criteria. Several systems (Ursa Major II, and also Coma Berenices, Bootes I, Leo IV, Segue 1) largely fulfill the requirements, indicating that their high-redshift predecessors did experience strong feedback effects that shut off star formation. We term the study of the entire stellar population of a dwarf galaxy for the purpose of inferring details about the nature and origin of the first galaxies "dwarf galaxy archaeology". This will provide clues to the connection of the first galaxies, the surviving, metal-poor dwarf galaxies, and the building blocks of the Milky Way.

preprint2012arXiv

Effect of Population III Multiplicity on Dark Star Formation

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.

preprint2012arXiv

Formation and evolution of primordial protostellar systems

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.

preprint2012arXiv

Pre-galactic metal enrichment - The chemical signatures of the first stars

The emergence of the first sources of light at redshifts of z ~ 10-30 signaled the transition from the simple initial state of the Universe to one of increasing complexity. We review recent progress in our understanding of the formation of the first stars and galaxies, starting with cosmological initial conditions, primordial gas cooling, and subsequent collapse and fragmentation. We emphasize the important open question of how the pristine gas was enriched with heavy chemical elements in the wake of the first supernovae. We conclude by discussing how the chemical abundance patterns conceivably allow us to probe the properties of the first stars and subsequent stellar generations, and allow us to test models of early metal enrichment.

preprint2012arXiv

Probing Pre-galactic Metal Enrichment with High-Redshift 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.

preprint2012arXiv

Star Formation in the First Galaxies I: Collapse Delayed by Lyman-Werner Radiation

We investigate the process of metal-free star formation in the first galaxies with a high-resolution cosmological simulation. We consider the cosmologically motivated scenario in which a strong molecule-destroying Lyman-Werner (LW) background inhibits effective cooling in low-mass haloes, delaying star formation until the collapse or more massive haloes. Only when molecular hydrogen (H2) can self-shield from LW radiation, which requires a halo capable of cooling by atomic line emission, will star formation be possible. To follow the formation of multiple gravitationally bound objects, at high gas densities we introduce sink particles which accrete gas directly from the computational grid. We find that in a 1 Mpc^3 (comoving) box, runaway collapse first occurs in a 3x10^7 M_sun dark matter halo at z~12 assuming a background intensity of J21=100. Due to a runaway increase in the H2 abundance and cooling rate, a self-shielding, supersonically turbulent core develops abruptly with ~10^4 M_sun in cold gas available for star formation. We analyze the formation of this self-shielding core, the character of turbulence, and the prospects for star formation. Due to a lack of fragmentation on scales we resolve, we argue that LW-delayed metal-free star formation in atomic cooling haloes is very similar to star formation in primordial minihaloes, although in making this conclusion we ignore internal stellar feedback. Finally, we briefly discuss the detectability of metal-free stellar clusters with the James Webb Space Telescope.

preprint2012arXiv

The First Galaxies: Assembly with Black Hole Feedback

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.

preprint2012arXiv

The First Stars and Galaxies - Basic Principles

Understanding the formation of the first stars and galaxies is a key problem in modern cosmology. In these lecture notes, we will derive some of the basic physical principles underlying this emerging field. We will consider the basic cosmological context, the cooling and chemistry in primordial gas, the physics of gravitational instability, and the main properties of the first stars. We will conclude with a discussion of the observational signature of the first sources of light, to be probed with future telescopes, such as the James Webb Space Telescope.

preprint2012arXiv

The Source Density and Observability of Pair-Instability Supernovae from the First Stars

Theoretical models predict that some of the first stars ended their lives as extremely energetic pair-instability supernovae (PISNe). With energies approaching 10^53 ergs, these supernovae are expected to be within the detection limits of the upcoming James Webb Space telescope (JWST), allowing observational constraints to be placed on the properties of the first stars. We estimate the source density of PISNe using a semi-analytic halo mass function based approach, accounting for the effects of feedback from star formation on the PISN rates using cosmological simulations. We estimate an upper limit of ~0.2 PISNe per JWST field of view at any given time. Feedback can reduce this rate significantly, e.g., lowering it to as little as one PISN per 4000 JWST fields of view for the most pessimistic explosion models. We also find that the main obstacle to observing PISNe from the first stars is their scarcity, not their faintness; exposures longer than a few times 10^4s will do little to increase the number of PISNe found. Given this we suggest a mosaic style search strategy for detecting PISNe from the first stars. Even rather high redshift PISNe are unlikely to be missed by moderate exposures, and a large number of pointings will be required to ensure a detection.

preprint2011arXiv

Effect of Streaming Motion of Baryons Relative to Dark Matter on the Formation of the First 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.

preprint2011arXiv

Simulations on a Moving Mesh: The Clustered Formation of Population III Protostars

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.

preprint2011arXiv

The First Galaxies

We review our current understanding of how the first galaxies formed at the end of the cosmic dark ages, a few 100 million years after the Big Bang. Modern large telescopes discovered galaxies at redshifts greater than seven, whereas theoretical studies have just reached the degree of sophistication necessary to make meaningful predictions. A crucial ingredient is the feedback exerted by the first generation of stars, through UV radiation, supernova blast waves, and chemical enrichment. The key goal is to derive the signature of the first galaxies to be observed with upcoming or planned next-generation facilities, such as the James Webb Space Telescope or Atacama Large Millimeter Array. From the observational side, ongoing deep-field searches for very high-redshift galaxies begin to provide us with empirical constraints on the nature of the first galaxies.

preprint2011arXiv

The First Galaxies: Assembly of Disks and Prospects for Direct Detection

The James Webb Space Telescope (JWST) will enable observations of galaxies at redshifts z > 10 and hence allow to test our current understanding of structure formation at very early times. Previous work has shown that the very first galaxies inside halos with virial temperatures T < 10^4 K and masses M < 10^8 M_sun at z > 10 are probably too faint, by at least one order of magnitude, to be detected even in deep exposures with JWST. The light collected with JWST may therefore be dominated by radiation from galaxies inside ten times more massive halos. We use cosmological zoomed smoothed particle hydrodynamics simulations to investigate the assembly of such galaxies and assess their observability with JWST. We compare two simulations that are identical except for the inclusion of non-equilibrium H/D chemistry and radiative cooling by molecular hydrogen. In both simulations a large fraction of the halo gas settles in two nested, extended gas disks which surround a compact massive gas core. The presence of molecular hydrogen allows the disk gas to reach low temperatures and to develop marked spiral structure but does not qualitatively change its stability against fragmentation. We post-process the simulated galaxies by combining idealized models for star formation with stellar population synthesis models to estimate the luminosities in nebular recombination lines as well as in the ultraviolet continuum. We demonstrate that JWST will be able to constrain the nature of the stellar populations in galaxies such as simulated here based on the detection of the He1640 recombination line. Extrapolation of our results to halos with masses both lower and higher than those simulated shows that JWST may find up to a thousand star-bursting galaxies in future deep exposures of the z > 10 universe.

preprint2011arXiv

The Formation and Fragmentation of Disks around Primordial Protostars

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.

preprint2010arXiv

Effects of Varying the Three-Body Molecular Hydrogen Formation Rate in Primordial Star Formation

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.

preprint2010arXiv

Fragmentation in the First Galaxies

We study the gravitational fragmentation of cold accretion streams flowing into a typical first galaxy. We use a one-zone hydrodynamical model to examine the thermal evolution of the gas entering a 10^8 M_sun DM halo at z=10. The goal is to find the expected fragmentation mass scale and thus a characteristic mass of the first population of stars to form by shock fragmentation at high redshift. Our model accurately describes the chemical and thermal evolution of the gas as we are specifically concerned with how the cooling of the gas alters its fragmentation properties. We find there to be a sharp drop in the fragmentation mass at a metallicity of ~10^-4 Z_sun when a strong molecule destroying, LW background is present. However, If molecules can efficiently form, they dominate the cooling at T < 10^4 K, demonstrating no 'critical metallicity'. Dust grains are not included in our chemical model, but we argue their inclusion would not significantly the results. We also find that this physical scenario allows for the formation of a cluster of solar mass fragments, or a single 10^4 M_sun fragment, possibly the precursors to primeval clusters and SMBHs. Lastly, we conclude that the usual assumption of isobaricity for galactic shocks breaks down in gas of sufficiently high metallicity, suggesting that metal cooling may lead to thermal instabilities.

preprint2010arXiv

Gravitational fragmentation in turbulent primordial gas and the initial mass function of Population III stars

We report results from numerical simulations of star formation in the early universe that focus on the dynamical behavior of metal-free gas under different initial and environmental conditions. In particular we investigate the role of turbulence, which is thought to ubiquitously accompany the collapse of high-redshift halos. We distinguish between two main cases: the birth of Population III.1 stars - those which form in the pristine halos unaffected by prior star formation - and the formation of Population III.2 stars - those forming in halos where the gas is still metal free but has an increased ionization fraction. This latter case can arise either from exposure to the intense UV radiation of stellar sources in neighboring halos, or from the high virial temperatures associated with the formation of massive halos, that is, those with masses greater than 1e8 solar masses. We find that turbulent primordial gas is highly susceptible to fragmentation in both cases, even for turbulence in the subsonic regime, i.e. for rms velocity dispersions as low as 20 % of the sound speed. Contrary to our original expectations, fragmentation is more vigorous and more widespread in pristine halos compared to pre-ionized ones. We therefore predict Pop III.1 stars to be on average of somewhat lower mass, and form in larger groups, than Pop III.2 stars. We find that fragment masses cover over two orders of magnitude, indicating that the resulting Population III initial mass function was significantly extended in mass as well. This prompts the need for a large, high-resolution study of the formation of dark matter minihalos that is capable of resolving the turbulent flows in the gas at the moment when the baryons become self-gravitating. This would help determine which, if any, of the initial conditions presented in our study are realized in nature.

preprint2010arXiv

Rotation Speed of the First Stars

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.

preprint2010arXiv

The First Galaxies: Chemical Enrichment, Mixing, and Star Formation

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.

preprint2009arXiv

Accretion onto Seed Black Holes in the First Galaxies

The validity of the hypothesis that the massive black holes in high redshift quasars grew from stellar-sized "seeds" is contingent on a seed's ability to double its mass every few ten million years. This requires that the seed accrete at approximately the Eddington-limited rate. In the specific case of radiatively efficient quasiradial accretion in a metal-poor protogalactic medium, for which the Bondi accretion rate is often prescribed in cosmological simulations of massive black hole formation, we examine the effects of the radiation emitted near the black hole's event horizon on the structure of the surrounding gas flow. We find that the radiation pressure from photoionization significantly reduces the steady-state accretion rate and renders the quasiradial accretion flow unsteady and inefficient. The time-averaged accretion rates are a small fraction of the Eddington-limited accretion rate for Thomson scattering. The pressure of Ly-alpha photons trapped near the HII region surrounding the black hole may further attenuate the inflow. These results suggest that an alternative to quasiradial, radiatively efficient Bondi-like accretion should be sought to explain the rapid growth of quasar-progenitor seed black holes.

preprint2009arXiv

The First Galaxies: Signatures of the Initial Starburst

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.

preprint2009arXiv

The First Supernova Explosions: Energetics, Feedback, and Chemical Enrichment

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.

preprint2009arXiv

The formation of the first stars and galaxies

Observations made using large ground-based and space-borne telescopes have probed cosmic history all the way from the present-day to a time when the Universe was less than a tenth of its present age. Earlier on lies the remaining frontier, where the first stars, galaxies, and massive black holes formed. They fundamentally transformed the early Universe by endowing it with the first sources of light and chemical elements beyond the primordial hydrogen and helium produced in the Big Bang. The interplay of theory and upcoming observations promises to answer the key open questions in this emerging field.

preprint2009arXiv

The observational signature of the first H II regions

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.

preprint2007arXiv

GRB Cosmology

Current observations are about to open up a direct window into the final frontier of cosmology: the first billion years in cosmic history when the first stars and galaxies formed. Even before the launch of the James Webb Space Telescope, it might be possible to utilize Gamma-ray Bursts (GRBs) as unique probes of cosmic star formation and the state of the intergalactic medium (IGM) up to redshifts of several tens, when the first (Population III) stars had formed. The Swift mission, or future satellites such as EXIST, might be the first observatories to detect individual Population III stars, provided that massive metal-free stars were able to trigger GRBs. Spectroscopic follow-up observations of the GRB afterglow emission would allow to probe the ionization state and metal enrichment of the IGM as a function of redshift.

preprint2007arXiv

Local radiative feedback in the formation of the first protogalaxies

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.

preprint2006arXiv

GRB Cosmology and the First Stars

Gamma-Ray Bursts (GRBs) are unique probes of the cosmic star formation history and the state of the intergalactic medium up to the redshifts of the first stars. In particular, the ongoing {\it Swift} mission might be the first observatory to detect individual Population~III stars, provided that the massive, metal-free stars were able to trigger GRBs. {\it Swift} will empirically constrain the redshift at which Population~III star formation was terminated, thus providing crucial input to models of cosmic reionization and metal enrichment.

preprint2005arXiv

The cooling of shock-compressed primordial gas

We find that at redshifts z > 10, HD line cooling allows strongly-shocked primordial gas to cool to the temperature of the cosmic microwave background (CMB). This temperature is the minimum value attainable via radiative cooling. Provided that the abundance of HD, normalized to the total number density, exceeds a critical level of ~ 10^{-8}, the CMB temperature floor is reached in a time which is short compared to the Hubble time. We estimate the characteristic masses of stars formed out of shocked primordial gas in the wake of the first supernovae, and resulting from the mergers of dark matter haloes during hierarchical structure formation to be ~ 10 M_{solar}. In addition, we show that cooling by HD enables the primordial gas in relic H II regions to cool to temperatures considerably lower than those reached via H_2 cooling alone. We confirm that HD cooling is unimportant in cases where the primordial gas does not go through an ionized phase, as in the formation process of the very first stars in z ~ 20 minihaloes of mass ~ 10^{6} M_{solar}.

preprint2002arXiv

The Expected Redshift Distribution of Gamma-Ray Bursts

We predict the redshift distribution of Gamma-Ray Bursts (GRBs) assuming that they trace the cosmic star formation history. We find that a fraction >~50% of all GRBs on the sky originate at a redshift z >~5, even though the fraction of the total stellar mass formed by z~5 is only ~15%. These two fractions are significantly different because they involve different cosmological factors when integrating the star formation rate over redshift. Hence, deep observations of transient events, such as GRB afterglows or supernovae, provide an ideal strategy for probing the high-redshift universe. We also show that the intrinsic distribution of GRB durations is bimodal but significantly narrower and shifted towards shorter durations than the observed distribution.

preprint2002arXiv

The Formation of a Star Cluster: Predicting the Properties of Stars and Brown Dwarfs

We present results from the largest numerical simulation of star formation to resolve the fragmentation process down to the opacity limit. The simulation follows the collapse and fragmentation of a large-scale turbulent molecular cloud to form a stellar cluster and, simultaneously, the formation of circumstellar discs and binary stars. This large range of scales enables us to predict a wide variety of stellar properties for comparison with observations. The calculation clearly demonstrates that star formation is a highly-dynamic and chaotic process. Star-disc encounters form binaries and truncate discs. Stellar encounters disrupt bound multiple systems. The cloud produces roughly equal numbers of stars and brown dwarfs, with masses down to the opacity limit for fragmentation (~5 Jupiter masses). The initial mass function is consistent with a Salpeter slope (Gamma=-1.35) above 0.5 Msun, a roughly flat distribution (Gamma=0) in the range 0.006-0.5 Msun, and a sharp cutoff below ~0.005 Msun. This is consistent with recent observational surveys. The brown dwarfs form by the dynamical ejection of low-mass fragments from dynamically unstable multiple systems before the fragments have been able to accrete to stellar masses. Close binary systems (with separations <10 AU) are not formed by fragmentation in situ. Rather, they are produced by hardening of initially wider multiple systems through a combination of dynamical encounters, gas accretion, and/or the interaction with circumbinary and circumtriple discs. Finally, we find that the majority of circumstellar discs have radii less than 20 AU due to truncation in dynamical encounters. This is consistent with observations of the Orion Trapezium Cluster and implies that most stars and brown dwarfs do not form large planetary systems.

Volker Bromm

What is connected

Connect this record

See the researcher in context

Building this map preview

77 published item(s)

Little Red Dots and their Progenitors from Direct Collapse Black Holes

VENUS: Two Faint Little Red Dots Separated by $\sim70\,\mathrm{pc}$ Hidden in a Single Lensed Galaxy at $z\sim7$

Effects of stellar-mass primordial black holes on first star formation

Maximal X-ray feedback in the pre-reionization universe

On the probability of the extremely lensed $z$=6.2 Earendel source being a Population~III star

Probing the initial mass function of the first stars with transients

Public Release of A-SLOTH: Ancient Stars and Local Observables by Tracing Halos

Dust emission from the first massive galaxies

Formation of the first stars and black holes

Gravitational waves from Population III binary black holes formed by dynamical capture

Lessons on early structure formation from a mature galaxy cluster observed at cosmic noon

Redshift Horizon for the Origins Space Telescope from Primordial Dust Emission

Signatures of first galaxies at FIR/sub-mm wavelengths

The role of faint population III supernovae in forming CEMP stars in ultra-faint dwarf galaxies

When did Population III star formation end?

Building up the Population III initial mass function from cosmological initial conditions

Contribution of the first galaxies to the cosmic far-infrared/sub-millimeter background - I. Mean background level

Evidence for a direct collapse black hole in the Lyman-alpha source CR7

Exploring the nature of the Lyman-$α$ emitter CR7

Gamma-ray bursts and Population III stars

Gravitational Waves from the Remnants of the First Stars

Lyman-alpha radiation hydrodynamics of galactic winds before cosmic reionization

The First Stars: formation under cosmic ray feedback

Toward $ab\,initio$ extremely metal poor stars

Abundance anomalies in metal-poor stars from Population III supernova ejecta hydrodynamics

Confined Population III Enrichment and the Prospects for Prompt Second-Generation Star Formation

Constraining the primordial initial mass function with stellar archaeology

Cosmological Impact of Population III Binaries

Metal Transport and Chemical Heterogeneity in Early Star Forming Systems

On the Physical Requirements for a Pre-Reionization Origin of the Unresolved Near-Infrared Background

Preserving chemical signatures of primordial star formation in the first low-mass stars

Star formation in the first galaxies - III. Formation, evolution, and characteristics of the first stellar cluster

The first galaxies: simulating their feedback-regulated assembly

The First Stars: formation under X-ray feedback

The Lyman-alpha signature of the first galaxies

Dwarf Spheroidal Satellite Formation in a Reionized Local Group

Formation of the First Low-Mass Stars from Cosmological Initial Conditions

Radiative Feedback from high mass X-ray binaries on the formation of the first galaxies and early reionization

Recovery from population III supernova explosions and the onset of second generation star formation

The Chemical Imprint of Silicate Dust on the Most Metal-Poor Stars

The First Stars: A Low-Mass Formation Mode

The Mutual Interaction Between Population III Stars and Self-Annihilating Dark Matter

Constraining the Statistics of Population III Binaries

Formation of the First Stars

On the operation of the chemothermal instability in primordial star-forming clouds

Star Formation in the First Galaxies - II: Clustered Star Formation and the Influence of Metal Line Cooling

The First Galaxies: Assembly under Radiative Feedback from the First Stars

Chemical Signatures of the First Galaxies: Criteria for One-Shot Enrichment

Effect of Population III Multiplicity on Dark Star Formation

Formation and evolution of primordial protostellar systems

Pre-galactic metal enrichment - The chemical signatures of the first stars

Probing Pre-galactic Metal Enrichment with High-Redshift Gamma-Ray Bursts

Star Formation in the First Galaxies I: Collapse Delayed by Lyman-Werner Radiation

The First Galaxies: Assembly with Black Hole Feedback

The First Stars and Galaxies - Basic Principles

The Source Density and Observability of Pair-Instability Supernovae from the First Stars

Effect of Streaming Motion of Baryons Relative to Dark Matter on the Formation of the First Stars

Simulations on a Moving Mesh: The Clustered Formation of Population III Protostars

The First Galaxies

The First Galaxies: Assembly of Disks and Prospects for Direct Detection

The Formation and Fragmentation of Disks around Primordial Protostars

Effects of Varying the Three-Body Molecular Hydrogen Formation Rate in Primordial Star Formation

Fragmentation in the First Galaxies

Gravitational fragmentation in turbulent primordial gas and the initial mass function of Population III stars

Rotation Speed of the First Stars

The First Galaxies: Chemical Enrichment, Mixing, and Star Formation

Accretion onto Seed Black Holes in the First Galaxies

The First Galaxies: Signatures of the Initial Starburst

The First Supernova Explosions: Energetics, Feedback, and Chemical Enrichment

The formation of the first stars and galaxies

The observational signature of the first H II regions

GRB Cosmology

Local radiative feedback in the formation of the first protogalaxies

GRB Cosmology and the First Stars