Andreas H. Pawlik

On the evolution of the H i column density distribution in cosmological simulations

We use a set of cosmological simulations combined with radiative transfer calculations to investigate the distribution of neutral hydrogen in the post-reionization Universe. We assess the contributions from the metagalactic ionizing background, collisional ionization and diffuse recombination radiation to the total ionization rate at redshifts z = 0 5. We find that the densities above which hydrogen self-shielding becomes important are consistent with analytic calculations and previous work. However, because of diffuse recombination radiation, whose intensity peaks at the same density, the transition between highly ionized and self-shielded regions is smoother than what is usually assumed. We provide fitting functions to the simulated photoionization rate as a function of density and show that post-processing simulations with the fitted rates yields results that are in excellent agreement with the original radiative transfer calculations. The predicted neutral hydrogen column density distributions agree very well with the observations. In particular, the simulations reproduce the remarkable lack of evolution in the column density distribution of Lyman limit and weak damped Lyα systems below z = 3. The evolution of the low column density end is affected by the increasing importance of collisional ionization with decreasing redshift. On the other hand, the simulations predict the abundance of strong damped Lyα systems to broadly track the cosmic star formation rate density.

CANDELS: The Contribution of the Observed Galaxy Population to Cosmic Reionization

We present measurements of the specific ultraviolet luminosity density from a sample of 483 galaxies at 6 z 8. These galaxies were selected from new deep near-infrared Hubble Space Telescope imaging from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey, Hubble UltraDeep Field 2009, and Wide Field Camera?3 Early Release Science programs. We investigate the contribution to reionization from galaxies that we observe directly, thus sidestepping the uncertainties inherent in complementary studies that have invoked assumptions regarding the intrinsic shape or the faint-end cutoff of the galaxy ultraviolet (UV) luminosity function. Due to our larger survey volume, wider wavelength coverage, and updated assumptions about the clumping of gas in the intergalactic medium (IGM), we find that the observable population of galaxies can sustain a fully reionized IGM at z?= 6, if the average ionizing photon escape fraction (f esc) is ~30%. Our result contrasts with a number of previous studies that have measured UV luminosity densities at these redshifts that vary by a factor of five, with many concluding that galaxies could not complete reionization by z?= 6 unless a large population of galaxies fainter than the detection limit were invoked, or extremely high values of f esc were present. The specific UV luminosity density from our observed galaxy samples at z?= 7 and 8 is not sufficient to maintain a fully reionized IGM unless f esc > 50%. We examine the contribution from galaxies in different luminosity ranges and find that the sub-L* galaxies we detect are stronger contributors to the ionizing photon budget than the L > L* population, unless f esc is luminosity dependent. Combining our observations with constraints on the emission rate of ionizing photons from Ly? forest observations at z?= 6, we find that we can constrain f esc < 34% (2?) if the observed galaxies are the only contributors to reionization, or <13% (2?) if the luminosity function extends to a limiting magnitude of M UV = ?13. These escape fractions are sufficient to sustain an ionized IGM by z?= 6. Current constraints on the high-redshift galaxy population imply that the volume ionized fraction of the IGM, while consistent with unity at z ? 6, appears to drop at redshifts not much higher than 7, consistent with a number of complementary reionization probes. If faint galaxies dominated the ionizing photon budget at z?= 6-7, future extremely deep observations with the James Webb Space Telescope will probe deep enough to directly observe them, providing an indirect constraint on the global ionizing photon escape fraction.

Power spectrum extraction for redshifted 21-cm Epoch of Reionization experiments: the LOFAR case

One of the aims of the Low Frequency Array (LOFAR) Epoch of Reionization (EoR) project is to measure the power spectrum of variations in the intensity of redshifted 21-cm radiation from the EoR. The sensitivity with which this power spectrum can be estimated depends on the level of thermal noise and sample variance, and also on the systematic errors arising from the extraction process, in particular from the subtraction of foreground contamination. We model the extraction process using realistic simulations of the cosmological signal, the foregrounds and noise, and so estimate the sensitivity of the LOFAR EoR experiment to the redshifted 21-cm power spectrum. Detection of emission from the EoR should be possible within 360 h of observation with a single station beam. Integrating for longer, and synthesizing multiple station beams within the primary (tile) beam, then enables us to extract progressively more accurate estimates of the power at a greater range of scales and redshifts. We discuss different observational strategies which compromise between depth of observation, sky coverage and frequency coverage. A plan in which lower frequencies receive a larger fraction of the time appears to be promising. We also study the nature of the bias which foreground fitting errors induce on the inferred power spectrum and discuss how to reduce and correct for this bias. The angular and line-of-sight power spectra have different merits in this respect, and we suggest considering them separately in the analysis of LOFAR data.

The impact of local stellar radiation on the H i column density distribution

It is often assumed that local sources of ionizing radiation have little impact on the distribution of neutral hydrogen in the post-reionization Universe. While this is a good assumption for the intergalactic medium, analytic arguments suggest that local sources may typically be more important than the meta-galactic background radiation for high column density absorbers (NHI > 10 17 cm 2 ). We post-process cosmological, hydrodynamical simulations with accurate radiation transport to investigate the impact of local stellar sources on the column density distribution function of neutral hydrogen. We demonstrate that the limited numerical resolution and the simplified treatment of the interstellar medium (ISM) that are typical of the current generation of cosmological simulations provide significant challenges, but that many of the problems can be overcome by taking two steps. First, using ISM particles rather than stellar particles as sources results in a much better sampling of the source distribution, effectively mimicking higher-resolution simulations. Second, by rescaling the source luminosities so that the amount of radiation escaping into the intergalactic medium agrees with that required to produce the observed background radiation, many of the results become insensitive to errors in the predicted fraction of the radiation that escapes the immediate vicinity of the sources. By adopting this strategy and by drastically varying the assumptions about the structure of the unresolved ISM, we conclude that we can robustly estimate the effect of local sources for column densities NHI ≪ 10 21 cm 2 . However, neither the escape fraction of ionizing radiation nor the effect of local sources on the abundance of NHI & 10 21 cm 2 systems can be predicted with confidence. We find that local stellar radiation is unimportant for NHI ≪ 10 17 cm 2 , but that it can affect Lyman Limit and Damped Lyα systems. For 10 18 < NHI < 10 21 cm 2 the impact of local sources increases with redshift. At z = 5 the abundance of absorbers is substantially reduced for NHI ≫ 10 17 cm 2 , but at z = 0 the effect only becomes significant for NHI & 10 21 cm 2 .

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. To this end, we perform cosmological zoomed smoothed particle hydrodynamics simulations of a dwarf galaxy assembling inside a halo reaching a virial mass ~109 M ☉ at z = 10. The simulations follow the non-equilibrium chemistry and cooling of primordial gas and the subsequent conversion of the cool dense gas into massive 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 molecular hydrogen dissociating radiation, and further with a simulation in which stars remain dark. Photodissociation and photoionization exert a strong negative feedback on the assembly of the galaxy inside the main minihalo progenitor. 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 from either dissociating or ionizing radiation 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 inside the virial radius 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 a rotationally supported extended disk inside the dwarf galaxy is a robust outcome of our simulations not affected by the inclusion of radiation. Low-mass galaxies in the neighborhood of the dwarf galaxy show a large scatter in the baryon fraction which is driven by radiative feedback from sources both internal and external to these galaxies. Our estimates of the observability of the first galaxies show that dwarf galaxies such as simulated here will be among the faintest galaxies the upcoming James Webb Space Telescope will detect. Our conclusions regarding the structure and observability of the first galaxies are subject to our neglect of feedback from supernovae and chemical enrichment as well as to statistical uncertainties implied by the limited number of galaxies in our simulations.

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.

Multifrequency, thermally coupled radiative transfer with traphic: method and tests

We present an extension of traphic, the method for radiative transfer of ionizing radiation in smoothed particle hydrodynamics simulations that we introduced in Pawlik & Schaye. The new version keeps all advantages of the original implementation: photons are transported at the speed of light, in a photon-conserving manner, directly on the spatially adaptive, unstructured grid traced out by the particles, in a computation time that is independent of the number of radiation sources, and in parallel on distributed memory machines. We extend the method to include multiple frequencies, both hydrogen and helium, and to model the coupled evolution of the temperature and ionization balance. We test our methods by performing a set of simulations of increasing complexity and including a small cosmological reionization run. The results are in excellent agreement with exact solutions, where available, and also with results obtained with other codes if we make similar assumptions and account for differences in the atomic rates used. We use the new implementation to illustrate the differences between simulations that compute photoheating in the grey approximation and those that use multiple frequency bins. We show that close to ionizing sources the grey approximation asymptotes to the multifrequency result if photoheating rates are computed in the optically thin limit, but that the grey approximation breaks down everywhere if, as is often done, the optically thick limit is assumed.

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 Pop III progenitor masses M . 40M , recovery is prompt, 10 Myr. 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 and hence more luminous progenitors. The gas around lower mass Pop III stars remains therefore 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, 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.

Spatially adaptive radiation-hydrodynamical simulations of galaxy formation during cosmological reionization

We present a suite of cosmological radiation-hydrodynamical simulations of the assembly of galaxies driving the reionization of the intergalactic medium (IGM) at z >~ 6. The simulations account for the hydrodynamical feedback from photoionization heating and the explosion of massive stars as supernovae (SNe). Our reference simulation, which was carried out in a box of size 25 comoving Mpc/h using 2 x 512^3 particles, produces a reasonable reionization history and matches the observed UV luminosity function of galaxies. Simulations with different box sizes and resolutions are used to investigate numerical convergence, and simulations in which either SNe or photoionization heating or both are turned off, are used to investigate the role of feedback from star formation. Ionizing radiation is treated using accurate radiative transfer at the high spatially adaptive resolution at which the hydrodynamics is carried out. SN feedback strongly reduces the star formation rates (SFRs) over nearly the full mass range of simulated galaxies and is required to yield SFRs in agreement with observations. Photoheating helps to suppress star formation in low-mass galaxies, but its impact on the cosmic SFR is small. Because the effect of photoheating is masked by the strong SN feedback, it does not imprint a signature on the UV galaxy luminosity function, although we note that our resolution is insufficient to model star-forming minihaloes cooling through molecular hydrogen transitions. Photoheating does provide a strong positive feedback on reionization because it smooths density fluctuations in the IGM, which lowers the IGM recombination rate substantially. Our simulations demonstrate a tight non-linear coupling of galaxy formation and reionization, motivating the need for the accurate and simultaneous inclusion of photoheating and SN feedback in models of the early Universe.

The first galaxies: simulating their feedback-regulated assembly

We investigate the formation of a galaxy reaching a virial mass of