Researcher profile

Thomas H. Greif

Thomas H. Greif contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2015arXiv

The numerical frontier of the high-redshift Universe

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

preprint2014arXiv

Formation of massive protostars in atomic cooling haloes

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

preprint2014arXiv

Multifrequency radiation hydrodynamics simulations of H2 line emission in primordial, star-forming clouds

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

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.

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

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

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.

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

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

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

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.