Source author record

Martin G. Haehnelt

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

55works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

55 published item(s)

preprint2022arXiv

Hydrogen reionisation ends by $z=5.3$: Lyman-$α$ optical depth measured by the XQR-30 sample

The presence of excess scatter in the Ly-$α$ forest at $z\sim 5.5$, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionisation. However, low data quality and systematic uncertainties complicate the use of Ly-$α$ transmission as a precision probe of reionisation's end stages. In this paper, we assemble a sample of 67 quasar sightlines at $z>5.5$ with high signal-to-noise ratios of $>10$ per $\leq 15$ km s$^{-1}$ spectral pixel, relying largely on the new XQR-30 quasar sample. XQR-30 is a large program on VLT/X-Shooter which obtained deep (SNR $>20$ per pixel) spectra of 30 quasars at $z>5.7$. We carefully account for systematics in continuum reconstruction, instrumentation, and contamination by damped Ly-$α$ systems. We present improved measurements of the mean Ly-$α$ transmission over $4.9<z<6.1$. Using all known systematics in a forward modelling analysis, we find excellent agreement between the observed Ly-$α$ transmission distributions and the homogeneous-UVB simulations Sherwood and Nyx up to $z\leq5.2$ ($<1 σ$), and mild tension ($\sim 2.5 σ$) at $z=5.3$. Homogeneous UVB models are ruled out by excess Ly-$α$ transmission scatter at $z\geq5.4$ with high confidence ($>3.5 σ$). Our results indicate that reionisation-related fluctuations, whether in the UVB, residual neutral hydrogen fraction, and/or IGM temperature, persist in the intergalactic medium until at least $z=5.3$ ($t=1.1$ Gyr after the Big Bang). This is further evidence for a late end to reionisation.

preprint2022arXiv

Limits on non-canonical heating and turbulence in the intergalactic medium from the low redshift Lyman-alpha forest

We examine the column density distribution function (CDDF) and Doppler parameter distribution from hydrodynamical simulations and Cosmic Origins Spectrograph (COS) observations of the Lyman-alpha forest at redshift $0\leq z\leq 0.2$. Allowing for a factor of two uncertainty in the metagalactic HI photoionisation rate, our hydrodynamical simulations are in good agreement ($1$-$1.5σ$) with the shape and amplitude of the observed CDDF at HI column densities $10^{13.3}\rm\,cm^{-2}\leq N_{\rm HI}\leq 10^{14.5}\rm\,cm^{-2}$. However, the Doppler widths of the simulated lines remain too narrow with respect to the COS data. We argue that invoking AGN feedback does not resolve this discrepancy. We also disfavour enhanced photoheating rates as a potential solution, as this requires an unphysically hard UV background spectrum. If instead appealing to a non-canonical source of heating, an additional specific heat injection of $u \lesssim 6.9\rm\,eV\,m_{\rm p}^{-1}$ is required at $z\lesssim 2.5$ for gas that has $N_{\rm HI}\simeq 10^{13.5}\rm\,cm^{-2}$ by $z=0.1$. Alternatively, there may be an unresolved line of sight turbulent velocity component of $v_{\rm turb}\lesssim 8.5\rm\,km\,s^{-1}(N_{\rm HI}/10^{13.5}\rm\,cm^{-2})^{0.21}$ for the coldest gas in the diffuse IGM.

preprint2022arXiv

Long Dark Gaps in the Ly$β$ Forest at $z<6$: Evidence of Ultra Late Reionization from XQR-30 Spectra

We present a new investigation of the intergalactic medium (IGM) near reionization using dark gaps in the Lyman-$β$ (Ly$β$) forest. With its lower optical depth, Ly$β$ offers a potentially more sensitive probe to any remaining neutral gas compared to commonly used Ly$α$ line. We identify dark gaps in the Ly$β$ forest using spectra of 42 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30 VLT Large Programme. Approximately $40\%$ of these QSO spectra exhibit dark gaps longer than $10h^{-1}{\rm Mpc}$ at $z\simeq5.8$. By comparing the results to predictions from simulations, we find that the data are broadly consistent both with models where fluctuations in the Ly$α$ forest are caused solely by ionizing ultraviolet background (UVB) fluctuations and with models that include large neutral hydrogen patches at $z<6$ due to a late end to reionization. Of particular interest is a very long ($L=28h^{-1}{\rm Mpc}$) and dark ($τ_{\rm eff} \gtrsim 6$) gap persisting down to $z\simeq 5.5$ in the Ly$β$ forest of the $z_{\rm}=5.85$ QSO PSO J025$-$11. This gap may support late reionization models with a volume-weighted average neutral hydrogen fraction of $ \langle x_{\rm HI}\rangle \gtrsim 5\%$ by $z=5.6$. Finally, we infer constraints on $\langle x_{\rm HI}\rangle$ over $5.5 \lesssim z \lesssim 6.0$ based on the observed Ly$β$ dark gap length distribution and a conservative relationship between gap length and neutral fraction derived from simulations. We find $\langle x_{\rm HI}\rangle \leq 0.05$, 0.17, and 0.29 at $z\simeq 5.55$, 5.75, and 5.95, respectively. These constraints are consistent with models where reionization ends significantly later than $z = 6$.

preprint2022arXiv

The effect of inhomogeneous reionisation on the Lyman-$α$ forest power spectrum at redshift $z>4$: implications for thermal parameter recovery

We use the Sherwood-Relics suite of hybrid hydrodynamical and radiative transfer simulations to model the effect of inhomogeneous reionisation on the 1D power spectrum of the \Lya forest transmitted flux at redshifts $4.2\leq z \leq 5$. Relative to models that assume a homogeneous UV background, reionisation suppresses the power spectrum at small scales, $k \sim 0.1\rm\,km^{-1}\,s$ by $\sim 10$ per cent because of spatial variations in the thermal broadening kernel and the divergent peculiar velocity field associated with over-pressurised intergalactic gas. On larger scales, $k<0.03\rm\,km^{-1}\,s$, the power spectrum is instead enhanced by $10$--$50$ per cent by large scale spatial variations in the neutral hydrogen fraction. The effect of inhomogeneous reionisation must therefore be accounted for in analyses of forthcoming high precision measurements. We provide a correction for the Ly$α$ forest power spectrum at $4.1\leq z \leq 5.4$ that can be easily applied within other parameter inference frameworks using similar reionisation models. We perform a Bayesian analysis of mock data to assess the extent of systematic biases that may arise in measurements of the intergalactic medium if ignoring this correction. At the scales probed by current high resolution Ly$α$ forest data at $z>4$, $0.006 \rm \,km^{-1}\,s\leq k \leq 0.2 \rm\, km^{-1}\,s$, we find inhomogeneous reionisation does not introduce any significant bias in thermal parameter recovery for the current measurement uncertainties of $\sim 10$ per cent. However, for $5$ per cent uncertainties, $\sim 1σ$ shifts between the estimated and true parameters occur.

preprint2022arXiv

The Nature of High ${\rm [OIII]}_{\rm 88μm}$/${\rm [CII]}_{\rm 158μm}$ Galaxies in the Epoch of Reionization: Low Carbon Abundance and a Top-Heavy IMF?

ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88μm}$/[CII]$_{\rm 158μm}$ ratios and [CII]$_{\rm 158μm}$ deficits compared to local galaxies. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties. In order to elucidate the underlying physics that drives this high-redshift phenomenon, we employ SPHINX$^{20}$, a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization. We find that the observed $z>6$ [OIII]-SFR and [CII]-SFR relations can only be reproduced when the C/O abundance ratio is $\sim8\times$ lower than Solar and the total metal production is $\sim4\times$ higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [CII]-SFR and [CII] luminosity functions are predicted to converge to observed values at $z\sim4.5$. While we demonstrate that ionisation parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [OIII]/[CII], observed values at $z>6$ can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [CII] and [OIII] can be used to predict the values of ionisation parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed $z>6$ galaxies. Finally, we demonstrate that [OI]$_{\rm 63μm}$ can be used as a replacement for [CII] in high-redshift galaxies where [CII] is unobserved and argue that more observation time should be used to target [OI] at $z>6$. (Abridged)

preprint2021arXiv

Revised estimates of CMB $B$-mode polarization induced by patchy reionization

The search for primordial gravitational waves through the $B$-mode polarization pattern in the CMB is one of the major goals of current and future CMB experiments. Besides foregrounds, a potential hurdle in this search is the anisotropic secondary $B$-mode polarization generated by the scattering of CMB photons off free electrons produced during patchy cosmological reionization. Robust predictions of these secondary anisotropies are challenging because of uncertainties in the reionization history. In this paper, we revise estimates of the reionization-induced $B$-mode signal by incorporating recent advances in the understanding of reionization through observations of the Lyman-$α$ forest. To derive these $B$-mode estimates, we use high-dynamic-range radiative transfer simulations of reionization that are calibrated to the Ly$α$ data. These simulations are also consistent with a variety of other high-redshift observations. We find that around multipoles $\ell\approx 100$, reionization induces $B$-mode power with $\ell(\ell+1)C_\ell^{BB}/2π\approx 4\times 10^{-6}\,μ$K$^2$. This secondary signal is thus at the level of the primordial signal with the tensor-to-scalar ratio $r<10^{-4}$, and can increase by a factor of $\sim 50$ if reionization is sourced by highly clustered sources residing in haloes with mass of $\sim 10^{11}$ M$_\odot$. Our findings suggest that the contribution of patchy reionization to the search for primordial gravitational waves is unlikely to be a concern for currently planned CMB experiments.

preprint2020arXiv

Constraining the second half of reionization with the Lyman-$β$ forest

We present an analysis of the evolution of the Lyman-series forest into the epoch of reionization using cosmological radiative transfer simulations in a scenario where reionization ends late. We explore models with different midpoints of reionization and gas temperatures. We find that once the simulations have been calibrated to match the mean flux of the observed Lyman-$α$ forest at $4 < z < 6$, they also naturally reproduce the distribution of effective optical depths of the Lyman-$β$ forest in this redshift range. We note that the tail of the largest optical depths that is most challenging to match corresponds to the long absorption trough of ULAS J0148+0600, which we have previously shown to be rare in our simulations. We consider the evolution of the Lyman-series forest out to higher redshifts, and show that future observations of the Lyman-$β$ forest at $z>6$ will discriminate between different reionization histories. The evolution of the Lyman-$α$ and Lyman-$γ$ forests are less promising as a tool for pushing studies of reionization to higher redshifts due to the stronger saturation and foreground contamination, respectively.

preprint2020arXiv

How to Quench a Dwarf Galaxy: The Impact of Inhomogeneous Reionization on Dwarf Galaxies and Cosmic Filaments

We use the SPHINX suite of high-resolution cosmological radiation hydrodynamics simulations to study how spatially and temporally inhomogeneous reionization impacts the baryonic content of dwarf galaxies and cosmic filaments. The SPHINX simulations simultaneously model an inhomogeneous reionization, follow the escape of ionising radiation from thousands of galaxies, and resolve haloes well below the atomic cooling threshold. This makes them an ideal tool for examining how reionization impacts star formation and the gas content of dwarf galaxies. We compare simulations with and without stellar radiation to isolate the effects of radiation feedback from that of supernova, cosmic expansion, and numerical resolution. We find that the gas content of cosmic filaments can be reduced by more than 80% following reionization. The gas inflow rates into haloes with $M_{vir}<10^8M_{\odot}$ are strongly affected and are reduced by more than an order of magnitude compared to the simulation without reionization. A significant increase in gas outflow rates is found for halo masses $M_{vir}<7\times10^7M_{\odot}$. Our simulations show that inflow suppression, rather than photoevaporation, is the dominant mechanism by which the baryonic content of high-redshift dwarf galaxies is regulated. At fixed redshift and halo mass, there is a large scatter in the halo baryon fractions that is entirely dictated by the timing of reionization in the local region surrounding a halo which can change by $Δz>3$ at fixed mass. Finally, although the gas content of high-redshift dwarf galaxies is significantly impacted by reionization, we find that most haloes with $M_{vir}<10^8M_{\odot}$ can remain self-shielded and form stars long after reionization, until their local gas reservoir is depleted, suggesting that local group dwarf galaxies do not necessarily exhibit star formation histories that peak prior to $z=6$...

preprint2020arXiv

New Methods for Identifying Lyman Continuum Leakers and Reionization-Epoch Analogues

Identifying low-redshift galaxies that emit Lyman Continuum radiation (LyC leakers) is one of the primary, indirect methods of studying galaxy formation in the epoch of reionization. However, not only has it proved challenging to identify such systems, it also remains uncertain whether the low-redshift LyC leakers are truly "analogues" of the sources that reionized the Universe. Here, we use high-resolution cosmological radiation hydrodynamics simulations to examine whether simulated galaxies in the epoch of reionization share similar emission line properties to observed LyC leakers at $z\sim3$ and $z\sim0$. We find that the simulated galaxies with high LyC escape fractions ($f_{\rm esc}$) often exhibit high O32 and populate the same regions of the R23-O32 plane as $z\sim3$ LyC leakers. However, we show that viewing angle, metallicity, and ionisation parameter can all impact where a galaxy resides on the O32-$f_{\rm esc}$ plane. Based on emission line diagnostics and how they correlate with $f_{\rm esc}$, lower-metallicity LyC leakers at $z\sim3$ appear to be good analogues of reionization-era galaxies. In contrast, local [SII]-deficient galaxies do not overlap with the simulated high-redshift LyC leakers on the SII-BPT diagram; however, this diagnostic may still be useful for identifying leakers. We use our simulated galaxies to develop multiple new diagnostics to identify LyC leakers using IR and nebular emission lines. We show that our model using only [CII]$_{\rm 158μm}$ and [OIII]$_{\rm 88μm}$ can identify potential leakers from non-leakers from the local Dwarf Galaxy Survey. Finally, we apply this diagnostic to known high-redshift galaxies and find that MACS1149_JD1 at $z=9.1$ is the most likely galaxy to be actively contributing to the reionization of the Universe.

preprint2020arXiv

Probing delayed-end reionization histories with the 21cm-LAE cross-power spectrum

We model the 21-cm signal and LAE population evolution during the epoch of reionization in order to predict the 21cm-LAE cross-power spectrum. We employ high-dynamic-range simulations of the IGM to create models that are consistent with constraints from the CMB, Lyman-$α$ forest and LAE population statistics. Using these models we consider the evolution of the cross-power spectrum for a selection of realistic reionization histories and predict the sensitivity of current and upcoming surveys to measuring this signal. We find that the imprint of a delayed-end to reionization can be observed by future surveys, and that strong constraints can be placed on the progression of reionization as late as $z=5.7$ using a Subaru-SKA survey. We make predictions for the signal-to-noise ratios achievable by combinations of Subaru/PFS with the MWA, LOFAR, HERA and SKA interferometers for an integration time of 1000 hours. We find that a Subaru-SKA survey could measure the cross-power spectrum for a late reionization at $z=6.6$ with a total signal-to-noise greater than 5, making it possible to constrain both the timing and bubble size at the end of reionization. Furthermore, we find that expanding the current Subaru/PFS survey area and depth by a factor of three would double the total signal-to-noise.

preprint2020arXiv

Probing the thermal state of the intergalactic medium at $z>5$ with the transmission spikes in high-resolution Ly$α$ forest spectra

We compare a sample of five high-resolution, high S/N Ly$α$ forest spectra of bright $6<z \lesssim 6.5$ QSOs aimed at spectrally resolving the last remaining transmission spikes at $z>5$ with those obtained from mock absorption spectra from the Sherwood and Sherwood-Relics suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile fitting procedure for the inverted transmitted flux, $1-F$, similar to the widely used Voigt profile fitting of the transmitted flux $F$ at lower redshifts, to characterise the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominant in low temeperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution ($\leq $ 8 km/s) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of $T_0= (11000 \pm 1600,10500\pm 2100,12000 \pm 2200)$ K at $z= (5.4,5.6,5.8)$. This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to $z\sim5.3$, the width distribution of the observed transmission spikes is consistent with the range of $T_0$ caused by spatial fluctuations in the temperature-density relation.

preprint2019arXiv

AEDGE: Atomic Experiment for Dark Matter and Gravity Exploration in Space

We propose in this White Paper a concept for a space experiment using cold atoms to search for ultra-light dark matter, and to detect gravitational waves in the frequency range between the most sensitive ranges of LISA and the terrestrial LIGO/Virgo/KAGRA/INDIGO experiments. This interdisciplinary experiment, called Atomic Experiment for Dark Matter and Gravity Exploration (AEDGE), will also complement other planned searches for dark matter, and exploit synergies with other gravitational wave detectors. We give examples of the extended range of sensitivity to ultra-light dark matter offered by AEDGE, and how its gravitational-wave measurements could explore the assembly of super-massive black holes, first-order phase transitions in the early universe and cosmic strings. AEDGE will be based upon technologies now being developed for terrestrial experiments using cold atoms, and will benefit from the space experience obtained with, e.g., LISA and cold atom experiments in microgravity. This paper is based on a submission (v1) in response to the Call for White Papers for the Voyage 2050 long-term plan in the ESA Science Programme. ESA limited the number of White Paper authors to 30. However, in this version (v2) we have welcomed as supporting authors participants in the Workshop on Atomic Experiments for Dark Matter and Gravity Exploration held at CERN: ({\tt https://indico.cern.ch/event/830432/}), as well as other interested scientists, and have incorporated additional material.

preprint2019arXiv

Long troughs in the Lyman-$α$ forest below redshift 6 due to islands of neutral hydrogen

A long (110 cMpc/$h$) and deep absorption trough in the Ly$α$ forest has been observed extending down to redshift 5.5 in the spectrum of ULAS J0148+0600. Although no Ly$α$ transmission is detected, Ly$β$ spikes are present which has led to claims that the gas along this trough must be ionized. Using high resolution cosmological radiative transfer simulations in large volumes, we show that in a scenario where reionization ends late ($z \sim 5.2$), our simulations can reproduce troughs as long as observed. In this model, we find that the troughs are caused by islands of neutral hydrogen. Small ionized holes within the neutral islands allow for the transmission of Ly$β$. We have also modelled the Ly$α$ emitter population around the simulated troughs, and show that there is a deficit of Ly$α$ emitters close to the trough as is observed.

preprint2016arXiv

Models of the Cosmological 21 cm Signal from the Epoch of Reionization Calibrated with Lyman-alpha and CMB Data

We present here 21 cm predictions from high dynamic range simulations for a range of reionization histories that have been tested against available Lyman-alpha and CMB data. We assess the observability of the predicted spatial 21 cm fluctuations by ongoing and upcoming experiments in the late stages of reionization in the limit in which the hydrogen spin temperature is significantly larger than the CMB temperature. Models consistent with the available Lyman-alpha data and CMB measurement of the Thomson optical depth predict typical values of 10--20 mK^2 for the variance of the 21 cm brightness temperature at redshifts z=7--10 at scales accessible to ongoing and upcoming experiments (k < 1 h/cMpc). This is only a factor of a few below the sensitivity claimed to have been already reached by ongoing experiments in the signal rms value. Our different models for the reionization history make markedly different predictions for the redshift evolution and thus frequency dependence of the 21 cm power spectrum and should be easily discernible by LOFAR (and later HERA and SKA1) at their design sensitivity. Our simulations have sufficient resolution to assess the effect of high-density Lyman limit systems that can self-shield against ionizing radiation and stay 21 cm bright even if the hydrogen in their surroundings is highly ionized. Our simulations predict that including the effect of the self-shielded gas in highly ionized regions reduces the large scale 21 cm power by about 30%.

preprint2016arXiv

Probing the end of reionization with the near-zones of z > 6 QSOs

QSO near-zones are an important probe of the the ionization state of the IGM at z ~ 6-7, at the end of reionization. We present here high-resolution cosmological 3D radiative transfer simulations of QSO environments for a wide range of host halo masses, 10$^{10-12.5}$ M_sun. Our simulated near-zones reproduce both the overall decrease of observed near-zone sizes at 6 < z < 7 and their scatter. The observable near-zone properties in our simulations depend only very weakly on the mass of the host halo. The size of the H II region expanding into the IGM is generally limited by (super-)Lyman Limit systems loosely associated with (low-mass) dark matter haloes. This leads to a strong dependence of near-zone size on direction and drives the large observed scatter. In the simulation centred on our most massive host halo, many sightlines show strong red damping wings even for initial volume averaged neutral hydrogen fractions as low as ~ 10$^{-3}$. For QSO lifetimes long enough to allow growth of the central supermassive black hole while optically bright, we can reproduce the observed near-zone of ULAS J1120+0641 only with an IGM that is initially neutral. Our results suggest that larger samples of z > 7 QSOs will provide important constraints on the evolution of the neutral hydrogen fraction and thus on how late reionization ends.

preprint2016arXiv

Testing the effect of galactic feedback on the IGM at z ~ 6 with metal-line absorbers

We present models of low- and high-ionization metal-line absorbers (O I, C II, C IV and Mg II) during the end of the reionization epoch, at z ~ 6. Using four cosmological hydrodynamical simulations with different feedback schemes (including the Illustris and Sherwood simulations) and two different choices of hydro-solver, we investigate how the overall incidence rate and equivalent width distribution of metal-line absorbers varies with the galactic wind prescription. We find that the O I and C II absorbers are reasonably insensitive to the feedback scheme. All models, however, struggle to reproduce the observations of C IV and Mg II, which are probing down to lower overdensities than O I and C II at z ~ 6, suggesting that the metals in the simulations are not being transported out into the IGM efficiently enough. The situation is improved but not resolved if we choose a harder (but still reasonable) and/or (locally) increased UV background at z ~ 6.

preprint2016arXiv

The Sherwood simulation suite: overview and data comparisons with the Lyman-alpha forest at redshifts 2 < z < 5

We introduce a new set of large scale, high resolution hydrodynamical simulations of the intergalactic medium: the Sherwood simulation suite. These are performed in volumes 10^3-160^3 h^-3 comoving Mpc^{3}, span almost four orders of magnitude in mass resolution with up to 17.2 billion particles, and employ a variety of physics variations including warm dark matter and galactic outflows. We undertake a detailed comparison of the simulations to high resolution, high signal-to-noise observations of the Lyman-alpha forest over the redshift range 2 < z < 5. The simulations are in very good agreement with the observational data, lending further support to the paradigm that the Lyman-alpha forest is a natural consequence of the web-like distribution of matter arising in LCDM cosmological models. Only a small number of minor discrepancies remain with respect to the observational data. Saturated Lyman-alpha absorption lines with column densities N_HI > 10^14.5 cm^-2 at 2 < z < 2.5 are underpredicted in the models. An uncertain correction for continuum placement bias is required to match the distribution and power spectrum of the transmitted flux, particularly at z > 4. Finally, the temperature of intergalactic gas in the simulations may be slightly too low at z=2.7 and a flatter temperature-density relation is required at z=2.4, consistent with the expected effects of non-equilibrium ionisation during He-II reionisation.

preprint2015arXiv

Fast cold gas in hot AGN outflows

Observations of the emission from spatially extended cold gas around bright high-redshift QSOs reveal surprisingly large velocity widths exceeding 2000 km s^(-1), out to projected distances as large as 30 kpc. The high velocity widths have been interpreted as the signature of powerful AGN-driven outflows. Naively, these findings appear in tension with hydrodynamic models in which AGN-driven outflows are energy-driven and thus very hot with typical temperatures T = 10^6-7 K. Using the moving-mesh code Arepo, we perform 'zoom-in' cosmological simulations of a z = 6 QSO and its environment, following black hole growth and feedback via energy-driven outflows. In the simulations, the QSO host galaxy is surrounded by a clumpy circum-galactic medium pre-enriched with metals due to supernovae-driven galactic outflows. As a result, part of the AGN-driven hot outflowing gas can cool radiatively, leading to large amounts (> 10^9 M_sun) of cold gas comoving with the hot bipolar outflow. This results in velocity widths of spatially extended cold gas similar to those observed. We caution, however, that gas inflows, random motions in the deep potential well of the QSO host galaxy and cooling of supernovae-driven winds contribute significantly to the large velocity width of the cold gas in the simulations, complicating the interpretation of observational data.

preprint2015arXiv

Feedback from Active Galactic Nuclei: Energy- versus momentum-driving

We employ hydrodynamical simulations using the moving-mesh code AREPO to investigate the role of energy and momentum input from Active Galactic Nuclei (AGN) in driving large-scale galactic outflows. We start by reproducing analytic solutions for both energy- and momentum-driven outflowing shells in simulations of a spherical isolated dark matter potential with gas in hydrostatic equilibrium and with no radiative cooling. We confirm that for this simplified setup, galactic outflows driven by a momentum input rate of order L_Edd/c can establish an M_BH - sigma relation with slope and normalisation similar to that observed. We show that momentum input at a rate of L_Edd/c is however insufficient to drive efficient outflows once cooling and gas inflows as predicted by cosmological simulations at resolved scales are taken into account. We argue that observed large-scale AGN-driven outflows are instead likely to be energy-driven and show that such outflows can reach momentum fluxes exceeding 10 L_Edd/c within the innermost 10 kpc of the galaxy. The outflows are highly anisotropic, with outflow rates and a velocity structure found to be inadequately described by spherical outflow models. We verify that the hot energy-driven outflowing gas is expected to be strongly affected by metal-line cooling, leading to significant amounts (>10^9 M_sun) of entrained cold gas.

preprint2015arXiv

Lyman-alpha emitters gone missing: evidence for late reionization?

We combine high resolution hydrodynamical simulations with an intermediate resolution, dark matter only simulation and an analytical model for the growth of ionized regions to estimate the large scale distribution and redshift evolution of the visibility of Lyman-alpha emission in 6<=z<=8 galaxies. The inhomogeneous distribution of neutral hydrogen during the reionization process results in significant fluctuations in the Lyman-alpha transmissivity on large scales. The transmissivity depends not only on the ionized fraction of the intergalactic medium by volume and the amplitude of the local ionizing background, but is also rather sensitive to the evolution of the relative velocity shift of the Lyman-alpha emission line due to resonant scattering. We reproduce a decline in the space density of Lyman-alpha emitting galaxies as rapid as observed with a rather rapidly evolving neutral fraction between z=6-8, and a typical Lyman-alpha line velocity offset of 100 km/s redward of systemic at z=6 which decreases toward higher redshift. The new (02/2015) Planck results indicate such a recent end to reionization is no longer disfavoured by constraints from the cosmic microwave background.

preprint2015arXiv

Seeding High Redshift QSOs by Collisional Runaway in Primordial Star Clusters

We study how runaway stellar collisions in high-redshift, metal-poor star clusters form very massive stars (VMSs) that can directly collapse to intermediate-mass black holes (IMBHs). We follow the evolution of a pair of neighbouring high-redshift mini-haloes with high-resolution, cosmological hydrodynamical zoom-in simulations using the adaptive mesh refinement code RAMSES combined with the non-equilibrium chemistry package KROME. The first collapsing mini-halo is assumed to enrich the central nuclear star cluster (NSC) of the other to a critical metallicity, sufficient for Population II (Pop. II) star formation at redshift $z\approx27$. Using the spatial configuration of the flattened, asymmetrical gas cloud forming in the core of the metal enriched halo, we set the initial conditions for simulations of an initially non-spherical star cluster with the direct summation code NBODY6 which are compared to about 2000 NBODY6 simulations of spherical star clusters for a wide range of star cluster parameters. The final mass of the VMS that forms depends strongly on the initial mass and initial central density of the NSC. For the initial central densities suggested by our RAMSES simulations, VMSs with mass > 400 M$_{\odot}$ can form in clusters with stellar masses of $\approx10^4$ M$_{\odot}$, and this can increase to well over 1000 M$_{\odot}$ for more massive and denser clusters. The high probability we find for forming a VMS in these mini-haloes at such an early cosmic time makes collisional runaway of Pop. II star clusters a promising channel for producing large numbers of high-redshift IMBHs that may act as the seeds of supermassive black holes.

preprint2015arXiv

The bias of DLAs at z ~ 2.3: contraining stellar feedback in shallow potential wells

We discuss the recent Baryon Oscillation Spectroscopic Survey measurement of a rather high bias factor for the host galaxies/haloes of Damped Lyman-alpha Absorbers (DLAs), in the context of our previous modelling of the physical properties of DLAs within the $Λ$ cold dark matter paradigm. Joint modelling of the column density distribution, the velocity width distribution of associated low ionization metal absorption, and the bias parameter suggests that DLAs are hosted by galaxies with dark matter halo masses in the range $10 < \log M_v < 12$, with a rather sharp cutoff at the lower mass end, corresponding to virial velocities of 35 km/sec. The observed properties of DLAs appear to suggest efficient (stellar) feedback in haloes with masses/virial velocities below the cutoff and a large retained baryon fraction (> 35 %) in haloes above the cutoff.

preprint2015arXiv

The environment of bright QSOs at z ~ 6: Star forming galaxies and X-ray emission

We employ cosmological hydrodynamical simulations to investigate models in which the supermassive black holes (BHs) powering luminous z ~ 6 QSOs grow from massive seeds. We simulate at high resolution 18 fields sampling regions with densities ranging from the mean cosmic density all the way to the highest sigma peaks in the Millennium simulation volume. Only in the most massive halos, BHs can grow to masses up to ~ 10^9 Msun by z ~ 6 without invoking super-Eddington accretion. Accretion onto the most massive BHs becomes limited by thermal AGN feedback by z ~ 9-8 with further BH growth proceeding in short Eddington limited bursts. Our modelling suggests that current flux-limited surveys of QSOs at high redshift preferentially detect objects at their peak luminosity and therefore miss a substantial population of QSOs powered by similarly massive BHs but with low accretion rates. To test whether the required host halo masses are consistent with the observed galaxy environment of z ~ 6 QSOs, we produce realistic rest-frame UV images of our simulated galaxies. Without strong stellar feedback, our simulations predict numbers of bright galaxies larger than observed by a factor ten or more. Supernova-driven galactic winds reduce the predicted numbers to a level consistent with observations indicating that stellar feedback was already very efficient at high redshifts. We have further investigated the effect of thermal AGN feedback on the surrounding gas. Our adopted AGN feedback prescription drives mostly energy-driven highly anisotropic outflows with gas speeds of >= 1000 km/s to distances of >= 10 kpc consistent with observations. The spatially extended thermal X-ray emission around bright QSOs powered by these outflows can exceed by large factors the emission expected without AGN feedback and is an important diagnostic of the mechanism whereby AGN feedback energy couples to surrounding gas.

preprint2015arXiv

The photoheating of the intergalactic medium in synthesis models of the UV background

We compare cosmological hydrodynamical simulations combined with the homogeneous metagalactic UV background (UVB) of Haardt & Madau (2012) (HM2012) to observations of the Lyman-alpha forest that are sensitive to the thermal and ionization state of the intergalactic medium (IGM). The transition from optically thick to thin photoheating predicted by the simple one-zone, radiative transfer model implemented by HM2012 predicts a thermal history that is in remarkably good agreement with the observed rise of the IGM temperature at z~3 if we account for the expected evolution of the volume filling factor of HeIII. Our simulations indicate that there may be, however, some tension between the observed peak in the temperature evolution and the rather slow evolution of the HeII opacities suggested by recent Hubble Space Telescope/COS measurements. The HM2012 UVB also underpredicts the metagalactic hydrogen photoionization rate required by our simulations to match the observed opacity of the forest at z>4 and z<2.

preprint2014arXiv

Massive Neutrinos and the Non-linear Matter Power Spectrum

We perform an extensive suite of N-body simulations of the matter power spectrum, incorporating massive neutrinos in the range M = 0.15-0.6 eV, probing the non-linear regime at scales k < 10 hMpc-1 at z < 3. We extend the widely used HALOFIT approximation to account for the effect of massive neutrinos on the power spectrum. In the strongly non-linear regime HALOFIT systematically over-predicts the suppression due to the free-streaming of the neutrinos. The maximal discrepancy occurs at k ~ 1 hMpc-1, and is at the level of 10% of the total suppression. Most published constraints on neutrino masses based on HALOFIT are not affected, as they rely on data probing the matter power spectrum in the linear or mildly non-linear regime. However, predictions for future galaxy, Lyman-alpha forest and weak lensing surveys extending to more non-linear scales will benefit from the improved approximation to the non-linear matter power spectrum we provide. Our approximation reproduces the induced neutrino suppression over the targeted scales and redshifts significantly better. We test its robustness with regard to changing cosmological parameters and a variety of modelling effects.

preprint2014arXiv

Numerical resolution effects on simulations of massive black hole seeds

We have performed high-resolution numerical simulations with the hydrodynamical AMR code Enzo to investigate the formation of massive seed black holes in a sample of six dark matter haloes above the atomic cooling threshold. The aim of this study is to illustrate the effects of varying the maximum refinement level on the final object formed. The virial temperatures of the simulated haloes range from $\rm{T} \sim 10000\ \rm{K} - 16000\ \rm{K}$ and they have virial masses in the range $\rm{M} \sim 2 \times 10^7 \rm{M_{\odot}}$ to $\rm{M} \sim 7 \times 10^7 \rm{M_{\odot}}$ at $z \sim 15$. The outcome of our six fiducial simulations is both generic and robust. A rotationally supported, marginally gravitationally stable, disk forms with an exponential profile. The mass and scale length of this disk depends strongly on the maximum refinement level used. Varying the maximum refinement level by factors between 1 / 64 to 256 times the fiducial level illustrates the care that must be taken in interpreting the results. The lower resolution simulations show tentative evidence that the gas may become rotationally supported out to 20 pc while the highest resolution simulations show only weak evidence of rotational support due to the shorter dynamical times for which the simulation runs. The higher resolution simulations do, however, point to fragmentation at small scales of the order of $\sim 100$ AU. In the highest resolution simulations a central object of a few times $10^2\ \rm{M_{\odot}}$ forms with multiple strongly bound, Jeans unstable, clumps of $\sim 10\ \rm{M_{\odot}}$ and radii of 10 - 20 AU suggesting the formation of dense star clusters in these haloes.

preprint2013arXiv

A consistent determination of the temperature of the intergalactic medium at redshift z=2.4

We present new measurements of the thermal state of the intergalactic medium (IGM) at $z\sim2.4$ derived from absorption line profiles in the Ly$α$ forest. We use a large set of high-resolution hydrodynamical simulations to calibrate the relationship between the temperature-density ($T$--$Δ$) relation in the IGM and the distribution of HI column densities, $N_{\rm HI}$, and velocity widths, $b_{\rm HI}$, of discrete Ly$α$ forest absorbers. This calibration is then applied to the measurement of the lower cut-off of the $b_{\rm HI}$--$N_{\rm HI}$ distribution recently presented by Rudie et al. (2012). We infer a power-law $T$--$Δ$ relation, $T=T_{0}Δ^{γ-1}$, with a temperature at mean density, $T_{0}=[1.00^{+0.32}_{-0.21}]\times10^{4}\rm\,K$ and slope $(γ-1)=0.54\pm0.11$. The slope is fully consistent with that advocated by the analysis of Rudie et al (2012); however, the temperature at mean density is lower by almost a factor of two, primarily due to an adjustment in the relationship between column density and physical density assumed by these authors. These new results are in excellent agreement with the recent temperature measurements of Becker et al. (2011), based on the curvature of the transmitted flux in the Ly$α$ forest. This suggests that the thermal state of the IGM at this redshift is reasonably well characterised over the range of densities probed by these methods.

preprint2013arXiv

A z=3.045 Lyman alpha emitting halo hosting a QSO and a possible candidate for AGN-triggered star-formation

In this third paper in a series on the nature of extended, asymmetric Lyman alpha emitters at z ~ 3 we report the discovery, in an ultra-deep, blind, spectroscopic long-slit survey, of a Lyman alpha emitting halo around a QSO at redshift 3.045. The QSO is a previously known, obscured AGN. The halo appears extended along the direction of the slit and exhibits two faint patches separated by 17 proper kpc in projection from the QSO. Comparison of the 2-dimensional spectrum with archival HST ACS images shows that these patches coincide spatially with emission from a peculiar, dumbbell-shaped, faint galaxy. The assumptions that the Lyman alpha emission patches are originating in the galaxy and that the galaxy is physically related to the QSO are at variance with photometric estimates of the galaxy redshift. We show, however, that a population of very young stars at the redshift of the QSO may fit the existing rest frame broad band UV photometry of the galaxy. If this scenario is correct then the symmetry of the galaxy in continuum and Lyman alpha emission, the extension of the QSO's Lyman alpha emission in its direction, and the likely presence of a young stellar population in close proximity to a (short-lived) AGN suggest that this may be an example of AGN feedback triggering external star formation in high redshift galaxies.

preprint2013arXiv

Probing the metallicity and ionization state of the circumgalactic medium at z~6 and beyond with OI absorption

Low ionization metal absorption due to OI has been identified as an important probe of the physical state of the inter-/circumgalactic medium at the tail-end of reionization. We use here high-resolution hydrodynamic simulations to interpret the incidence rate of OI absorbers at z~6 as observed by Becker et al. (2011). We infer weak OI absorbers (EW > 0.1 A) to have typical HI column densities in the range of sub-DLAs, densities of 80 times the mean baryonic density and metallicities of about 1/500 th solar. This is similar to the metallicity inferred at similar overdensities at z~3, suggesting that the metal enrichment of the circumgalactic medium around low-mass galaxies has already progressed considerably by z~6. The apparently rapid evolution of the incidence rates for OI absorption over the redshift range 5<z<6 mirrors that of self-shielded Lyman-Limit systems at lower redshift and is mainly due to the rapid decrease of the meta-galactic photo-ionization rate at z>5. We predict the incidence rate of OI absorbers to continue to rise rapidly with increasing redshift as the IGM becomes more neutral. If the distribution of metals extends to lower density regions, OI absorbers will allow the metal enrichment of the increasingly neutral filamentary structures of the cosmic web to be probed.

preprint2013arXiv

Star-forming Galactic Contrails at z=3.2 as a Source of Metal Enrichment and Ionizing Radiation

A spectroscopically detected Lyman alpha emitting halo at redshift 3.216 in the GOODS-N field is found to reside at the convergence of several Lyman alpha filaments. HST images show that some of the filaments are inhabited by galaxies. Several of the galaxies in the field have pronounced head-tail structures, which are partly aligned with each other. The blue colors of most tails suggest the presence of young stars, with the emission from at least one of the galaxies apparently dominated by high equivalent width Lyman alpha. Faint, more diffuse, and similarly elongated, apparently stellar features, can be seen over an area with a linear extent of at least 90 kpc. The region within several arcseconds of the brightest galaxy exhibits spatially extended emission by HeII, NV and various lower ionization metal lines. The gas-dynamical features present are strongly reminiscent of ram-pressure stripped galaxies, including evidence for recent star formation in the stripped contrails. Spatial gradients in the appearance of several galaxies may represent a stream of galaxies passing from a colder to a hotter intergalactic medium. The stripping of gas from the in-falling galaxies, in conjunction with the occurrence of star formation and stellar feedback in the galactic contrails suggests a mechanism for the metal enrichment of the high redshift intergalactic medium that does not depend on long-range galactic winds, at the same time opening a path for the escape of ionizing radiation from galaxies.

preprint2012arXiv

Constraints on Massive Neutrinos from the CFHTLS Angular Power Spectrum

We use the galaxy angular power spectrum at $z\sim0.5-1.2$ from the Canada-France-Hawaii-Telescope Legacy Survey Wide fields (CFHTLS-Wide) to constrain separately the total neutrino mass $\sum{m_ν}$ and the effective number of neutrino species $N_{\rm{eff}}$. This survey has recently benefited from an accurate calibration of the redshift distribution, allowing new measurements of the (non-linear) matter power spectrum in a unique range of scales and redshifts sensitive to neutrino free streaming. Our analysis makes use of a recent model for the effect of neutrinos on the weakly non-linear matter power spectrum derived from accurate N-body simulations. We show that CFHTLS, combined with WMAP7 and a prior on the Hubble constant provides an upper limit of $\sum{m_ν}<0.29\,$eV and $N_{\rm{eff}} =4.17^{+1.62}_{-1.26}$ (2$\,σ$ confidence levels). If we omit smaller scales which may be affected by non-linearities, these constraints become $\sum{m_ν}<0.41\,$eV and $N_{\rm{eff}} =3.98^{+2.02}_{-1.20}$ (2$\,σ$ confidence levels). Finally we show that the addition of other large scale structures probes can further improve these constraints, demonstrating that high redshift large volumes surveys such as CFHTLS are complementary to other cosmological probes of the neutrino mass.

preprint2012arXiv

Detection of dark galaxies and circum-galactic filaments fluorescently illuminated by a quasar at z=2.4

A deep narrow-band survey for Ly-alpha emission carried out on the VLT-FORS2 has revealed 98 Ly-alpha candidates down to a flux limit of 4.e-18 erg/s/cm^2 in a volume of 5500 comoving Mpc^3 at z=2.4 centered on the hyperluminous quasar HE0109-3518. The properties of the detected sources in terms of their i) equivalent width distribution, ii) luminosity function, and iii) the average luminosity versus projected distance from the quasar, all suggest that a large fraction of these objects have been fluorescently "illuminated" by HE0109-3518. This conclusion is supported by comparison with detailed radiative transfer simulations of the effects of the quasar illumination. 18 objects have a rest-frame Equivalent Width (EW0) larger than 240A, the expected limit for Ly-alpha emission powered by Population II star formation and 12 sources among these do not have any continuum counterpart in a deep V-band imaging of the same field. For these, a stacking analysis indicates EW0>800A, effectively ruling out Ly-alpha powered by internal star formation. These sources are thus the best candidates so far for proto-galactic clouds or "dark" galaxies at high-redshift, whose existence has recently been suggested by several theoretical studies. Assuming they are mostly ionized by the quasar radiation, we estimate that their gas masses would be about 10^9 Msun implying that their star formation efficiencies (SFE) are less than 1.e-11 yr^-1 several times below the SFE of the most gas-rich dwarf galaxies locally, and five hundred times lower than typical massive star-forming galaxies at z~2. We have also discovered extended, filamentary gas, also likely illuminated by the quasar, around some of the brightest continuum-detected sources with EW0>240A. This emission is compatible with the expectations for circum-galactic cold streams but other origins, including tidal stripping, are also possible.

preprint2012arXiv

Extended and Filamentary Lyman Alpha Emission from the Formation of a Protogalactic Halo at z=2.63

We report the observation of a further asymmetric, extended Lyman alpha emitting halo at z=2.63, from our ultra-deep, long-slit spectroscopic survey of faint high redshift emitters, undertaken with Magellan LDSS3 in the GOODS-S field. The Lya emission, detected over more than 30 kpc, is spatially coincident with a concentration of galaxies visible in deep broad-band imaging. While these faint galaxies without spectroscopic redshifts cannot with certainty be associated with one another or with the Lya emission, there are a number of compelling reasons why they very probably form a Milky Way halo-mass group at the Lya redshift. A filamentary structure, possibly consisting of Lya emission at very high equivalent width, and evidence for disturbed stellar populations, suggest that the properties of the emitting region reflect ongoing galaxy assembly, with recent galaxy mergers and star formation occurring in the group. Hence, the Lya provides unique insights into what is probably a key mode of galaxy formation at high redshifts. The Lya emission may be powered by cooling radiation or spatially extended star formation in the halo, but is unlikely to be fluorescence driven by either an AGN or one of the galaxies. The spatial profile of the emission is conspicuously different from that of typical Lya emitters or Lyman break galaxies, which is consistent with the picture that extended emission of this kind represents a different stage in the galaxy formation process. Faint, extended Lya emitters such as these may be lower-mass analogues of the brighter Lya blobs. Our observations provide further, circumstantial evidence that galaxy mergers may promote the production and / or escape of ionizing radiation, and that the halos of interacting galaxies may be significant sources for ionizing photons during the epoch of reionization.

preprint2012arXiv

Neutrino masses and cosmological parameters from a Euclid-like survey: Markov Chain Monte Carlo forecasts including theoretical errors

We present forecasts for the accuracy of determining the parameters of a minimal cosmological model and the total neutrino mass based on combined mock data for a future Euclid-like galaxy survey and Planck. We consider two different galaxy surveys: a spectroscopic redshift survey and a cosmic shear survey. We make use of the Monte Carlo Markov Chains (MCMC) technique and assume two sets of theoretical errors. The first error is meant to account for uncertainties in the modelling of the effect of neutrinos on the non-linear galaxy power spectrum and we assume this error to be fully correlated in Fourier space. The second error is meant to parametrize the overall residual uncertainties in modelling the non-linear galaxy power spectrum at small scales, and is conservatively assumed to be uncorrelated and to increase with the ratio of a given scale to the scale of non-linearity. It hence increases with wavenumber and decreases with redshift. With these two assumptions for the errors and assuming further conservatively that the uncorrelated error rises above 2% at k = 0.4 h/Mpc and z = 0.5, we find that a future Euclid-like cosmic shear/galaxy survey achieves a 1-sigma error on Mnu close to 32 meV/25 meV, sufficient for detecting the total neutrino mass with good significance. If the residual uncorrelated errors indeed rises rapidly towards smaller scales in the non-linear regime as we have assumed here then the data on non-linear scales does not increase the sensitivity to the total neutrino mass. Assuming instead a ten times smaller theoretical error with the same scale dependence, the error on the total neutrino mass decreases moderately from sigma(Mnu) = 18 meV to 14 meV when mildly non-linear scales with 0.1 h/Mpc < k < 0.6 h/Mpc are included in the analysis of the galaxy survey data.

preprint2012arXiv

On the rapid demise of Lyman-alpha emitters at z>7 due to the increasing incidence of optically thick absorption systems

A variety of independent observational studies have now reported a significant decline in the fraction of Lyman-break galaxies which exhibit Ly-a emission over the redshift interval z=6-7. In combination with the strong damping wing extending redward of Ly-a in the spectrum of the bright z=7.085 quasar ULAS 1120+0641, this has strengthened suggestions that the hydrogen in the intergalactic medium (IGM) is still substantially neutral at z~7. Current theoretical models imply HI fractions as large as 40-90 per cent may be required to explain these data assuming there is no intrinsic evolution in the Ly-a emitter population. We propose that such large neutral fractions are not necessary. Based on a hydrodynamical simulation which reproduces the absorption spectra of high-redshift (z~6-7) quasars, we demonstrate that the opacity of the intervening IGM redward of rest-frame Ly-a can rise rapidly in average regions of the Universe simply because of the increasing incidence of absorption systems which are optically thick to Lyman continuum photons as the tail-end of reionisation is approached. Our simulations suggest these data do not require a large change in the IGM neutral fraction by several tens of per cent from z=6-7, but may instead be indicative of the rapid decrease in the typical mean free path for ionising photons expected during the final stages of reionisation.

preprint2011arXiv

Filamentary Infall of Cold Gas and Escape of Lyman Alpha and Hydrogen Ionizing Radiation from an Interacting High-Redshift Galaxy

We present observations of a peculiar Lyman alpha-emitting galaxy at redshift 3.344, discovered in a deep, blind spectroscopic survey for faint Lyman alpha emitters with the Magellan II telescope in the Hubble Ultra Deep Field (HUDF). The galaxy exhibits complex Lyman alpha emission, including an extended, asymmetric component that is partially suppressed by damped Lyman alpha absorption, and two spatially elongated, narrow emission features. Archival HST ACS imaging shows evidence for tidal disruption of the stellar component. This V=27 galaxy appears to give us unprecedented insights into two fundamental stages in the formation of structure at high redshift: the inflow of gas into ordinary galaxies, and the escape of ionizing radiation into the intergalactic medium. Neutral hydrogen, falling in partly in form of a narrow filament, appears to emit fluorescent Lyman alpha photons induced by the stellar ionizing flux escaping from the disturbed galaxy. The in-falling material may represent primary cold accretion or an interaction-triggered inflow. The rate of ionizing photons required by the observed Lyman alpha emission is consistent with the rate of photons produced by the observed stellar population, with roughly 50 percent of ionizing photons escaping from the immediate galaxy and encountering the in-falling gas. The observational properties of the galaxy lend support to a picture where galaxy interactions facilitate the escape of both Lyman alpha and ionizing radiation. We argue that galaxies like the present object may be common at high redshift. This galaxy may therefore be a late example of an interacting population of dwarf galaxies contributing significantly to the reionization of the universe.

preprint2011arXiv

Galactic Winds and Extended Lya Emission from the Host Galaxies of High Column Density QSO Absorption Systems

We present 3D resonant radiative transfer simulations of the spatial and spectral diffusion of the Lya radiation from a central source in the host galaxies of high column density absorption systems at z ~ 3. The radiative transfer simulations are based on a suite of cosmological galaxy formation simulations that reproduce a wide range of observed properties of damped Lya absorption systems. The Lya emission is predicted to be spatially extended up to several arcsec, and the spectral width of the Lya emission is broadened to several hundred (in some case more than thousand) km/s. The distribution and the dynamical state of the gas in the simulated galaxies is complex, the latter with significant contributions from rotation and both in- and out-flows. The emerging Lya radiation extends to gas with column densities of N_HI ~ 10^{18} cm^{-2} and its spectral shape varies strongly with viewing angle. The strong dependence on the central \hi column density and the HI velocity field suggests that the Lya emission will also vary strongly with time on timescales of a few dynamical times of the central region. Such variations with time should be especially pronounced at times where the host galaxy undergoes a major merger and/or starburst. Depending on the pre-dominance of in- or out-flow along a given sightline and the central column density, the spectra show prominent blue peaks, red peaks or double-peaked profiles. Both spatial distribution and spectral shape are very sensitive to details of the galactic wind implementation. Stronger galactic winds result in more spatially extended Lya emission and - somewhat counterintuitively - a narrower spectral distribution.

preprint2011arXiv

How neutral is the intergalactic medium surrounding the redshift z=7.085 quasar ULAS J1120+0641?

The quasar ULAS J1120+0641 at redshift z=7.085 has a highly ionised near zone which is smaller than those around quasars of similar luminosity at z~6. The spectrum also exhibits evidence for a damping wing extending redward of the systemic Lya redshift. We use radiative transfer simulations in a cosmological context to investigate the implications for the ionisation state of the inhomogeneous IGM surrounding this quasar. Our simulations show that the transmission profile is consistent with an IGM in the vicinity of the quasar with a volume averaged HI fraction of f_HI>0.1 and that ULAS J1120+0641 has been bright for 10^6--10^7 yr. The observed spectrum is also consistent with smaller IGM neutral fractions, f_HI ~ 10^-3--10-4, if a damped Lya system in an otherwise highly ionised IGM lies within 5 proper Mpc of the quasar. This is, however, predicted to occur in only ~5 per cent of our simulated sight-lines for a bright phase of 10^6--10^7 yr. Unless ULAS J1120+0641 grows during a previous optically obscured phase, the low age inferred for the quasar adds to the theoretical challenge of forming a 2x10^9 M_sol black hole at this high redshift.

preprint2011arXiv

Improved measurements of the intergalactic medium temperature around quasars: possible evidence for the initial stages of He-II reionisation at z~6

We present measurements of the intergalactic medium (IGM) temperature within ~5 proper Mpc of seven luminous quasars at z~6. The constraints are obtained from the Doppler widths of Lyman-alpha absorption lines in the quasar near-zones and build upon our previous measurement for the z=6.02 quasar SDSS J0818+1722. The expanded data set, combined with an improved treatment of systematic uncertainties, yields an average temperature at the mean density of log (T_0/K) = 4.21 \pm 0.03 (\pm^0.06_0.07) at 68 (95) per cent confidence for a flat prior distribution over 3.2 < log (T_0/K) < 4.8. In comparison, temperatures measured from the general IGM at z~5 are ~0.3 dex cooler, implying an additional source of heating around these quasars which is not yet present in the general IGM at slightly lower redshift. This heating is most likely due to the recent reionisation of He-II in vicinity of these quasars, which have hard, non-thermal ionising spectra. The elevated temperatures may therefore represent evidence for the earliest stages of He-II reionisation in the most biased regions of the high-redshift Universe. The temperature as a function of distance from the quasars is consistent with being constant, log(T_0/K)~4.2, with no evidence for a line-of-sight thermal proximity effect. However, the limited extent of the quasar near-zones prevents the detection of He-III regions larger than ~5 proper Mpc. Under the assumption the quasars have reionised the He-II in their vicinity, we infer the data are consistent with an average optically bright phase of duration in excess of 10^6.5 yr. These measurements represent the highest-redshift IGM temperature constraints to date, and thus provide a valuable data set for confronting models of H-I reionisation.

preprint2010arXiv

A first direct measurement of the intergalactic medium temperature around a quasar at z=6

The thermal state of the intergalactic medium (IGM) provides an indirect probe of both the HI and HeII reionisation epochs. Current constraints on the IGM temperature from the Lya forest are restricted to the redshift range 2<z<4.5, limiting the ability to probe the thermal memory of HI reionisation toward higher redshift. In this work, we present the first direct measurement of the IGM temperature around a z=6 quasar by analysing the Doppler widths of Lya absorption lines in the proximity zone of SDSS J0818+1722. We use a high resolution (R= 40000) Keck/HIRES spectrum in combination with detailed numerical modelling to obtain the temperature at mean density, T_0=23600\pm^5000_6900K (\pm^9200_9300K) at 68 (95) per cent confidence assuming a prior probability 13500K<T_0<38500 K following HI and HeII reionisation. This enables us to place an upper limit on the redshift of HI reionisation, z_H, within 33 comoving Mpc of SDSS J0818+1722. If the quasar reionises the HeII in its vicinity, then in the limit of instantaneous reionisation we infer z_H<9.0 (11.0) at 68 (95) per cent confidence assuming photoheating is the dominant heat source and that HI reionisation is driven by ionising sources with soft spectra, typical of population II stars. If the HI and HeII in the IGM around SDSS J0818+1722 are instead reionised simultaneously by a population of massive metal-free stars, characterised by very hard ionising spectra, we obtain a tighter upper limit of z_H<8.4 (9.4). Initiating reionisation at higher redshifts produces temperatures which are too low with respect to our constraint unless the HI ionising sources or the quasar itself have spectra significantly harder than typically assumed.

preprint2010arXiv

Detection of Extended He II Reionization in the Temperature Evolution of the Intergalactic Medium

We present new measurements of the temperature of the intergalactic medium (IGM) derived from the Lyman-alpha forest over 2.0 < z < 4.8. The small-scale structure in the forest of 61 high-resolution QSO spectra is quantified using a new statistic, the curvature, and the conversion to temperature calibrated using a suite of hydrodynamic simulations. At each redshift we focus on obtaining the temperature at an optimal overdensity probed by the Lyman-alpha forest, T(Delta), where the temperature is nearly a one-to-one function of the curvature regardless of the slope of the temperature-density relation. The median 2-sigma statistical uncertainty in these measurements is 8 per cent, though there may be comparable systematic errors due to the unknown amount of Jeans smoothing in the IGM. We use our T(Delta) results to infer the temperature at the mean density, T0. Even for a maximally steep temperature-density relation, T0 must increase from ~8000 K at z ~ 4.4 to >~12000 K at z ~ 2.8. This increase is not consistent with the monotonic decline in T0 expected in the absence of He II reionization. We therefore interpret the observed rise in temperature as evidence of He II reionization beginning at z >~ 4.4. The evolution of T0 is consistent with an end to He II reionization at z ~ 3, as suggested by opacity measurements of the He II Lyman-alpha forest, although the redshift at which T0 peaks will depend somewhat on the evolution of the temperature-density relation. These new temperature measurements suggest that the heat input due to the reionization of He II dominates the thermal balance of the IGM over an extended period with Delta_z >~ 1.

preprint2010arXiv

Faint Resonantly Scattered Lyman Alpha Emission from the Absorption Troughs of Damped Lyman Alpha Systems at z ~ 3

We demonstrate that the Lyman alpha emission in the absorption troughs of a large sample of stacked damped Lyman alpha absorption systems (DLAS) presented by Rahmani et al (2010) is consistent with the spectral profiles and luminosities of a recently detected population of faint Lyman alpha emitters at z ~ 3. This result supports the suggestion that the faint emitters are to be identified with the host galaxies of DLAS at these redshifts.

preprint2010arXiv

Measurements of the UV background at 4.6 < z < 6.4 using the quasar proximity effect

We present measurements of the ionising ultraviolet background (UVB) at z ~ 5-6 using the quasar proximity effect. The fifteen quasars in our sample cover the range 4.6 < z_q < 6.4, enabling the first proximity effect measurements of the UVB at z > 5. The metagalactic hydrogen ionisation rate, Gamma_bkg, was determined by modelling the combined ionisation field from the quasar and the UVB in the proximity zone on a pixel-by-pixel basis. The optical depths in the spectra were corrected for the expected effect of the quasar until the mean flux in the proximity region equalled that in the average Ly-alpha forest, and from this we make a measurement of Gamma_bkg. A number of systematic effects were tested using synthetic spectra. Noise in the flux was found to be the largest source of bias at z ~ 5, while uncertainties in the mean transmitted Ly-alpha flux are responsible for the largest bias at z ~ 6. The impacts of large-scale overdensities and Lyman limit systems on Gamma_bkg were also investigated, but found to be small at z > 5. We find a decline in Gamma_bkg with redshift, from log(Gamma_bkg) = -12.15 $\pm$ 0.16 at z ~ 5 to log(Gamma_bkg) = -12.84 $\pm$ 0.18 at z ~ 6 (1 sigma errors). Compared to UVB measurements at lower redshifts, our measurements suggest a drop of a factor of five in the HI photoionisation rate between z ~ 4 and z ~ 6. The decline of Gamma_bkg appears to be gradual, and we find no evidence for a sudden change in the UVB at any redshift that would indicate a rapid change in the attenuation length of ionising photons. Combined with recent measurements of the evolution of the mean free path of ionising photons, our results imply decline in the emissivity of ionising photons by roughly a factor of two from z ~ 5 to 6, albeit with significant uncertainty due to the measurement errors in both Gamma_bkg and the mean free path.

preprint2010arXiv

The effect of neutrinos on the matter distribution as probed by the Intergalactic Medium

We present a suite of full hydrodynamical cosmological simulations that quantitatively address the impact of neutrinos on the (mildly non-linear) spatial distribution of matter and in particular on the neutral hydrogen distribution in the Intergalactic Medium (IGM), which is responsible for the intervening Lyman-alpha absorption in quasar spectra. The free-streaming of neutrinos results in a (non-linear) scale-dependent suppression of power spectrum of the total matter distribution at scales probed by Lyman-alpha forest data which is larger than the linear theory prediction by about 25% and strongly redshift dependent. By extracting a set of realistic mock quasar spectra, we quantify the effect of neutrinos on the flux probability distribution function and flux power spectrum. The differences in the matter power spectra translate into a ~2.5% (5%) difference in the flux power spectrum for neutrino masses with Sigma m_ν = 0.3 eV (0.6 eV). This rather small effect is difficult to detect from present Lyman-alpha forest data and nearly perfectly degenerate with the overall amplitude of the matter power spectrum as characterised by sigma_8. If the results of the numerical simulations are normalized to have the same sigma_8 in the initial conditions, then neutrinos produce a smaller suppression in the flux power of about 3% (5%) for Sigma m_ν = 0.6$ eV (1.2 eV) when compared to a simulation without neutrinos. We present constraints on neutrino masses using the Sloan Digital Sky Survey flux power spectrum alone and find an upper limit of Sigma m_ν < 0.9$ eV (2 sigma C.L.), comparable to constraints obtained from the cosmic microwave background data or other large scale structure probes.

preprint2009arXiv

Faint extended Lyalpha emission due to star formation at the centre of high-column density QSO absorption systems

We use detailed Lyalpha radiative transfer calculations to further test the claim of Rauch et al. (2008) that they have detected spatially extended faint Lyalpha emission from the elusive host population of Damped Lyalpha Absorption systems (DLAs) in their recent ultra-deep spectroscopic survey. We investigate the spatial and spectral distribution of Lyalpha emission due to star-formation at the centre of DLAs, and its dependence on the spatial and velocity structure of the gas. Our model simultaneously reproduces the observed properties of DLAs and the faint Lyalpha emitters, including the velocity width and column density distribution of DLAs and the large spatial extent of the emission of the faint emitters. Our modelling confirms previous suggestions that DLAs are predominately hosted by Dark Matter (DM) halos in the mass range 10^{9.5}-10^{12} M_sun, and are thus of significantly lower mass than those inferred for L_* Lyman Break Galaxies (LBGs). Our modelling suggests that DM halos hosting DLAs retain up to 20% of the cosmic baryon fraction in the form of neutral hydrogen, and that star formation at the centre of the halos is responsible for the faint Lyalpha emission. The scattering of a significant fraction of the Lyalpha emission to the observed radii, which can be as large as 50 kpc or more, requires the amplitude of the bulk motions of the gas at the centre of the halos to be moderate. The observed space density and size distribution of the emitters together with the incidence rate of DLAs suggests that the Lyalpha emission due to star formation has a duty cycle of ~ 25%.

preprint2009arXiv

Growing the first bright quasars in cosmological simulations of structure formation

We employ cosmological hydrodynamical simulations to study the growth of massive black holes (BHs) at high redshifts subject to BH merger recoils from gravitational wave emission. We select the most massive dark matter halo at z=6 from the Millennium simulation, and resimulate its formation at much higher resolution including gas physics and a model for BH seeding, growth and feedback. Assuming that the initial BH seeds are relatively massive, of the order of 10^5 Msun, and that seeding occurs around z~15 in dark matter haloes of mass 10^9-10^10 Msun, we find that it is possible to build up supermassive BHs (SMBHs) by z=6 that assemble most of their mass during extended Eddington-limited accretion periods. The properties of the simulated SMBHs are consistent with observations of z=6 quasars in terms of the estimated BH masses and bolometric luminosities, the amount of star formation occurring within the host halo, and the presence of highly enriched gas in the innermost regions of the host galaxy. After a peak in the BH accretion rate at z=6, the most massive BH has become sufficiently massive for the growth to enter into a much slower phase of feedback-regulated accretion. We explore the full range of expected recoils and radiative efficiencies, and also consider models with spinning BHs. In the most `pessimistic' case where BH spins are initially high, we find that the growth of the SMBHs can be potentially hampered if they grow mostly in isolation and experience only a small number of mergers. Whereas BH kicks can expel a substantial fraction of low mass BHs, they do not significantly affect the build up of the SMBHs. On the contrary, a large number of BH mergers has beneficial consequences for the growth of the SMBHs by considerably reducing their spin. [Abridged]

preprint2008arXiv

A joint model for the emission and absorption properties of damped Lyman alpha absorption systems

The recently discovered population of ultra-faint extended line emitters can account for the majority of the incidence rate of Damped Lyman Alpha systems (DLAs) at z ~ 3 if the line emission is interpreted as Ly alpha. We show here that a model similar to that proposed by Haehnelt, Steinmetz, & Rauch (2000), which explains the incidence rate and kinematics of DLAs in the context of $Λ$CDM models for structure formation, also reproduces the size distribution of the new population of faint Ly alpha emitters for plausible parameters. This lends further support to identification of the emitters with the hitherto elusive population of DLA host galaxies. The observed incidence rate of DLAs together with the observed space density and size distribution of the emitters suggest a duty cycle of ~ 0.2 - 0.4 for the Ly alpha emission from DLA host galaxies. We further show that Ly alpha cooling is expected to contribute little to the Ly alpha emission for the majority of emitters. This leaves centrally concentrated star formation at a rate of a few tenths M_sun/yr, surrounded by extended Ly alpha halos with radii up to 30-50 kpc, as the most plausible explanation for the origin of the emission. Both the luminosity function of Ly alpha emission and the velocity width distribution of low ionization absorption require that galaxies inside Dark Matter (DM) halos with virial velocities < 50 - 70 km/s contribute little to the incidence rate of DLAs at z ~ 3, suggesting that energy and momentum input due to star formation efficiently removes gas from these halos. Galaxies with DM halos with virial velocities of 100 - 150 km/s appear to account for the majority of DLA host galaxies. DLA host galaxies at z ~ 3 should thus become the building blocks of typical present-day galaxies.

preprint2008arXiv

Pathways to massive black holes and compact star clusters in pre-galactic dark matter haloes with virial temperatures > 10000K

Large dynamic range numerical simulations of atomic cooling driven collapse of gas in pre-galactic DM haloes with T_vir ~ 10000 K show that the gas loses 90% and more of its angular momentum before rotational support sets in. In a fraction of these haloes where the metallicity is low and UV radiation suppresses H_2 cooling, conditions are thus very favourable for the rapid build-up of massive black holes. Depending on the progression of metal enrichment, the continued suppression of H_2 cooling by external and internal UV radiation and the ability to trap the entropy produced by the release of gravitational energy, the gas at the centre of the halo is expected to form a supermassive star, a stellar-mass black hole accreting at super-Eddington accretion rates or a compact star-cluster undergoing collisional run-away of massive stars at its centre. In all three cases a massive black hole of initially modest mass finds itself at the center of a rapid inflow of gas with inflow rates of ~ 1 M_solar\yr. The massive black hole will thus grow quickly to a mass of 10^5 to 10^6 M_solar until further inflow is halted either by consumption of gas by star formation or by the increasing energy and momentum feedback from the growing massive black hole. Conditions for the formation of massive seed black holes in this way are most favourable in haloes with T_vir ~ 15000 K and V_vir ~ 20 km\s with less massive haloes not allowing collapse of gas by atomic cooling and more massive haloes being more prone to fragmentation. This should imprint a characteristic mass on the mass spectrum of an early population of massive black hole seeds in pre-galactic haloes which will later grow into the observed population of supermassive black holes in galactic bulges.

preprint2005arXiv

Cosmological and astrophysical parameters from the SDSS flux power spectrum and hydrodynamical simulations of the Lyman-alpha forest

(abridged) The flux power spectrum of the Lyman-alpha forest in quasar (QSO) absorption spectra is sensitive to a wide range of cosmological and astrophysical parameters and instrumental effects. Modelling the flux power spectrum in this large parameter space to an accuracy comparable to the statistical uncertainty of large samples of QSO spectra is very challenging. We use here a coarse grid of hydrodynamical simulations run with GADGET-2 to obtain a ``best guess'' model around which we calculate a finer grid of flux power spectra using a Taylor expansion of the flux power spectrum to first order. We find that the SDSS flux power spectrum alone is able to constrain a wide range of parameters including the amplitude of the matter power spectrum sigma_8, the matter density Omega_m, the spectral index of primordial density fluctuations n, the effective optical depth tau_eff and its evolution. The thermal history of the Intergalactic Medium (IGM) is, however, poorly constrained and the SDSS data favour either an unplausibly large temperature or an unplausibly steep temperature-density relation. By enforcing a thermal history of the IGM consistent with that inferred from high-resolution QSO spectra, we find the following values for the best fitting model (assuming a flat Universe with a cosmological constant and zero neutrino mass): Omega_ m=0.28 \pm 0.03, n=0.95\pm0.04, σ_8=0.91\pm0.07 (1σerror bars).We argue that the major uncertainties in this measurement are still systematic rather than statistical.

preprint1997arXiv

High-redshift galaxies, their active nuclei and central black holes

We demonstrate that the luminosity function of the recently detected population of star-forming galaxies and the QSO luminosity function at z=3 can be matched with the mass function of dark matter haloes predicted by hierarchical cosmogonies for lifetimes of optically bright QSOs in the range 10^6 to 10^8 yr. We suggest that the mass of supermassive black holes may be limited by the back-reaction of the emitted energy on the accretion flow in a self-gravitating disc. This would imply a relation of black hole to halo mass of the form M_bh ~(v_halo)^5 and a typical duration of the optically bright QSO phase of a few times 10^7 yr. The high integrated mass density of black holes inferred from recent black hole mass estimates in nearby galaxies may indicate that the overall efficiency of supermassive black holes for producing blue light is smaller than previously assumed. We discuss three possible accretion modes with low optical emission efficiency: (i) accretion at far above the Eddington rate, (ii) accretion obscured by dust, and (iii) accretion below the critical rate leading to an advection dominated accretion flow lasting for a Hubble time. We further argue that accretion with low optical efficiency might be closely related to the origin of the hard X-ray background and that the ionizing background might be progressively dominated by stars rather than QSOs at higher redshift.

preprint1996arXiv

Non-equilibrium effects on line-of-sight size estimates of QSO absorption systems

Estimates of the linear extent of heavy-element absorption systems along the line-of-sight to a QSO often assume that the cloud is photoionized and that the temperature takes the equilibrium value where photo-heating balances line cooling. We show that rather small deviations from this photoionization equilibrium temperature caused by additional heating processes will lead to an overestimate of the neutral hydrogen fraction and thus to an underestimate of the thickness of the absorber by about two orders of magnitude. Such temperature deviations are indicated both by observations and numerical simulations. This interpretation reconciles the discrepancy between the rather small extent of heavy-element-absorption systems parallel to the line-of-sight obtained from a standard photoionization analysis and the much larger transverse sizes estimates inferred from the observation of common absorption in the spectra of close quasar pairs.

preprint1996arXiv

QSO Metal Absorption Systems at High Redshift and the Signature of Hierarchical Galaxy Formation

In hierarchical cosmogonies at redshift three the matter content of a typical present-day galaxy is dispersed over several individual clumps embedded in sheet-like structures, often aligned along filaments. We have used hydrodynamical simulations to investigate the spatial distribution and absorption properties of metal enriched gas in such regions of ongoing galaxy formation. Strong CII and SiIV absorption occurs for lines-of-sight passing the densest regions close to the center of the protogalactic clumps, while CIV is a good tracer of the prominent filamentary structures and OVI becomes the strongest absorption feature for lines-of-sight through low-density regions far away from fully-collapsed objects. Observed column density ratios can be well reproduced if a mean metallicity [Z/H] = -2.5, relative abundances as found in metal-poor stars, a UV background with intensity J_{-22} = 3 at the Lyman limit, and either a power law spectrum (alpha = -1.5) or the spectral shape proposed by Haardt & Madau (1996) are assumed. The observed scatter in [C/H] is about a magnitude larger than that in the simulations suggesting inhomogeneous metal distribution. The detailed small and large scale kinematics observed in CIV and HI are also well reproduced.

preprint1995arXiv

Omega_baryon and the Geometry of Intermediate Redshift Lyman alpha Absorption Systems

Estimates of Omega_baryon from primordial nucleosynthesis together with standard assumptions about the ionization state of low column density Lyman alpha forest clouds can be used to determine an upper limit for the cloud thickness along the line of sight. This upper limit provides significant constraints on the axis ratio of the absorbers and on the overall fraction of the baryonic matter contained in these objects if the recently measured large values for their transversal size are representative. The absorbers have to be considerably flattened structures and may possibly contain a substantial fraction of the baryonic matter at high redshift.

preprint1995arXiv

Using the kinematic Sunyaev-Zeldovich effect to determine the peculiar velocities of clusters of galaxies

We have investigated the possibility of inferring peculiar velocities for clusters of galaxies from the Doppler shift of scattered cosmic microwave background (CMB) photons. We find that if the core radius of the gas distribution or the beam size of the instrument is larger than 3-7 arcminutes, then the maximum attainable signal-to-noise ratio is determined by confusion with primary fluctuations. For smaller angular scales, ``cosmic confusion'' is less important and instrumental noise and/or foreground emission will be the limiting factor. For a cluster with the optical depth of the Coma cluster and for an optimal filtering technique, typical one-sigma errors span the wide range from 400 to 1600 km/s, depending on the cosmological model, the resolution of the instrument and the core radius of the cluster. The results have important implications for the design of future high-resolution surveys of the CMB. Individual peculiar velocities will be measurable only for a few fast moving clusters at intermediate redshift unless cosmic fluctuations are smaller than most standard cosmological scenarios predict. However, a reliable measurement of bulk velocities of ensembles of X-ray bright clusters will be possible on very large scales (100-500 Mpc/h).

preprint1994arXiv

Momentum transfer by an internal source of ionizing radiation - an important feedback process during galaxy formation ?!

The role of momentum transfer (``radiation pressure'') due to an internal source of ionizing radiation for the formation of baryonic structures is investigated. Fully-ionized self-gravitating gaseous objects can be radiation-pressure supported on a characteristic length scale $D_{\rm rp} \sim 100 \pc - 3 \kpc$. On smaller scales momentum transfer due to ionizing radiation will be the dominant force and for spherical collapsing objects of mass $\la 10^{10}\Msol$ a complete bounce is possible. A population of massive stars and/or accretion onto a central compact object are natural sources of ionizing radiation in newly-forming baryonic structures. Radiation pressure is therefore likely to play an important role for the dynamical and thermal evolution of the intergalactic medium and the mass-to-light ratio of small galaxies. The effect will be especially important for hierarchical cosmogonies where galactic structures build up by merging of smaller objects. Radiation pressure due to ionizing radiation might furthermore be responsible for substructures of size $D_{\rm rp}$ in the stellar component of large galaxies.