Researcher profile

Stuart Wyithe

Stuart Wyithe contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

43 published item(s)

preprint2015arXiv

Connecting Faint End Slopes of the Lyman-$α$ emitter and Lyman-break Galaxy Luminosity Functions

We predict Lyman-$α$ (Ly$α$) luminosity functions (LFs) of Ly$α$-selected galaxies (Ly$α$ emitters, or LAEs) at $z=3-6$ using the phenomenological model of Dijkstra & Wyithe (2012). This model combines observed UV-LFs of Lyman-break galaxies (LBGs, or drop out galaxies), with constraints on their distribution of Ly$α$ line strengths as a function of UV-luminosity and redshift. Our analysis shows that while Ly$α$ LFs of LAEs are generally not Schechter functions, these provide a good description over the luminosity range of $\log_{10}( L_α/{\rm erg}\,{\rm s}^{-1})=41-44$. Motivated by this result, we predict Schechter function parameters at $z=3-6$. Our analysis further shows that (i) the faint end slope of the Ly$α$ LF is steeper than that of the UV-LF of Lyman-break galaxies, (with a median $α_{Lyα} < -2.0$ at $z\gtrsim 4$), and (ii) a turn-over in the Ly$α$ LF of LAEs at Ly$α$ luminosities $10^{40}$ erg s$^{-1}<L_α\lesssim 10^{41}$ erg s$^{-1}$ may signal a flattening of UV-LF of Lyman-break galaxies at $-12>M_{\rm UV}>-14$. We discuss the implications of these results - which can be tested directly with upcoming surveys - for the Epoch of Reionization.

preprint2015arXiv

Correcting the z~8 Galaxy Luminosity Function for Gravitational Lensing Magnification Bias

We present a Bayesian framework to account for the magnification bias from both strong and weak gravitational lensing in estimates of high-redshift galaxy luminosity functions. We illustrate our method by estimating the $z\sim8$ UV luminosity function using a sample of 97 Y-band dropouts (Lyman break galaxies) found in the Brightest of Reionizing Galaxies (BoRG) survey and from the literature. We find the luminosity function is well described by a Schechter function with characteristic magnitude of $M^\star = -19.85^{+0.30}_{-0.35}$, faint-end slope of $α= -1.72^{+0.30}_{-0.29}$, and number density of $\log_{10} Ψ^\star [\textrm{Mpc}^{-3}] = -3.00^{+0.23}_{-0.31}$. These parameters are consistent within the uncertainties with those inferred from the same sample without accounting for the magnification bias, demonstrating that the effect is small for current surveys at $z\sim8$, and cannot account for the apparent overdensity of bright galaxies compared to a Schechter function found recently by Bowler et al. (2014a,b) and Finkelstein et al. (2014). We estimate that the probability of finding a strongly lensed $z\sim8$ source in our sample is in the range $\sim 3-15 \%$ depending on limiting magnitude. We identify one strongly-lensed candidate and three cases of intermediate lensing in BoRG (estimated magnification $μ>1.4$) in addition to the previously known candidate group-scale strong lens. Using a range of theoretical luminosity functions we conclude that magnification bias will dominate wide field surveys -- such as those planned for the Euclid and WFIRST missions -- especially at $z>10$. Magnification bias will need to be accounted for in order to derive accurate estimates of high-redshift luminosity functions in these surveys and to distinguish between galaxy formation models.

preprint2015arXiv

Imaging HII Regions from Galaxies and Quasars During Reionisation with SKA

The ionisation structure of the Intergalactic Medium (IGM) during reionisation is sensitive to the unknown galaxy formation physics that prevailed at that time. This structure introduces non-Gaussian statistics into the redshifted 21 cm fluctuation amplitudes that can only be studied through tomographic imaging, which will clearly discriminate between different galaxy formation scenarios. Imaging the ionisation structure and cosmological HII regions during reionisation is therefore a key goal for the SKA. For example, the SKA1-LOW baseline design with a 1 km diameter core will resolve HII regions expected from galaxy formation models which include strong feedback on low-mass galaxy formation. Imaging the smaller HII regions that result from galaxy formation in the absence of SNe feedback will also be possible for SKA1-LOW in the later stages of reionisation, but may require the greater sensitivity of SKA early in the reionisation era. In addition to having baselines long enough to resolve the HII regions, the field of view for SKA1-LOW reionisation experiments should be at least several degrees in order to image the largest HI structures towards the end of reionisation. The baseline design with 35 meter diameter stations has a field of view within a single primary pointing which is sufficient for this purpose.

preprint2015arXiv

The stellar mass function and star formation rate-stellar mass relation of galaxies at z ~ 4 - 7

We investigate the evolution of the star formation rate-stellar mass relation (SFR-M*) and Galaxy Stellar Mass Function (GSMF) of z ~ 4-7 galaxies, using cosmological simulations run with the smoothed particle hydrodynamics code P-GADGET3(XXL). We explore the effects of different feedback prescriptions (supernova driven galactic winds and AGN feedback), initial stellar mass functions and metal cooling. We show that our fiducial model, with strong energy-driven winds and early AGN feedback, is able to reproduce the observed stellar mass function obtained from Lyman-break selected samples of star forming galaxies at redshift 6 < z < 7. At z ~ 4, observed estimates of the GSMF vary according to how the sample was selected. Our simulations are more consistent with recent results from K-selected samples, which provide a better proxy of stellar masses and are more complete at the high mass end of the distribution. We find that in some cases simulated and observed SFR-M* relations are in tension, and this can lead to numerical predictions for the GSMF in excess of the GSMF observed. By combining the simulated SFR(M*) relationship with the observed star formation rate function at a given redshift, we argue that this disagreement may be the result of the uncertainty in the SFR-M* (Luv-M*) conversion. Our simulations predict a population of faint galaxies not seen by current observations.

preprint2014arXiv

Evolution in the Escape Fraction of Ionizing Photons and the Decline in Strong Lya Emission from z>6 Galaxies

The rapid decline in the number of strong Lyman Alpha (Lya) emitting galaxies at z > 6 provides evidence for neutral hydrogen in the IGM, but is difficult to explain with plausible models for reionization. We demonstrate that the observed reduction in Lya flux from galaxies at z > 6 can be explained by evolution in the escape fraction of ionizing photons, f_esc. We find that the median observed drop in the fraction of galaxies showing strong Lya emission, as well as the observed evolution of the Lya luminosity function both follow from a small increase in f_esc of Delta f_esc ~ 0.1 from f_esc ~ 0.6 at z ~ 6. This high escape fraction may be at odds with current constraints on the ionising photon escape fraction, which favor smaller values of f_esc < 20%. However, models that invoke a redshift evolution of f_ esc that is consistent with these constraints can suppress the z~7 Lya flux to the observed level, if they also include a small evolution in global neutral fraction of Delta x_HI ~ 0.2. Thus, an evolving escape fraction of ionising photons can be a plausible part of the explanation for evolution in the Lya emission of high redshift galaxies. More generally, our analysis also shows that the drop in the Lya fraction is quantitatively consistent with the observed evolution in the Lya luminosity functions of Lya Emitters.

preprint2014arXiv

Simulated star formation rate functions at z~4-7, and the role of feedback in high-z galaxies

We study the role of feedback from supernovae and black holes in the evolution of the star formation rate function (SFRF) of z~4-7 galaxies. We use a new set of cosmological hydrodynamic simulations, ANGUS (AustraliaN GADGET-3 early Universe Simulations), run with a modified and improved version of the parallel TreePM-smoothed particle hydrodynamics code GADGET-3 called P-GADGET3(XXL), that includes a self-consistent implementation of stellar evolution and metal enrichment. In our simulations both Supernova (SN) driven galactic winds and Active Galactic Nuclei (AGN) act simultaneously in a complex interplay. The SFRF is insensitive to feedback prescription at z>5, meaning that it cannot be used to discriminate between feedback models during reionisation. However, the SFRF is sensitive to the details of feedback prescription at lower redshift. By exploring different SN driven wind velocities and regimes for the AGN feedback, we find that the key factor for reproducing the observed SFRFs is a combination of &#34;strong&#34; SN winds and early AGN feedback in low mass galaxies. Conversely, we show that the choice of initial mass function and inclusion of metal cooling have less impact on the evolution of the SFRF. When variable winds are considered, we find that a non-aggressive wind scaling is needed to reproduce the SFRFs at z>4. Otherwise, the amount of objects with low SFRs is greatly suppressed and at the same time winds are not effective enough in the most massive systems.

preprint2013arXiv

A Predicted New Population of UV-faint Galaxies at z>4

We show that a bursty model of high redshift star formation explains several puzzling observations of the high redshift galaxy population. We begin by pointing out that the observed specific star formation rate requires a duty-cycle of ~10%, which is much lower than found in many hydro-dynamical simulations. This value follows directly from the fact that the observed star formation rate in galaxies integrated over a Hubble time exceeds the observed stellar mass by an order of magnitude. We use the large observed specific star formation rate to calibrate the efficiency of feedback in a model for the high redshift star formation rate which includes merger driven star formation regulated by SNe feedback. This model reproduces the star formation rate density function and the stellar mass function of galaxies at 4<z<7. A prediction of the model is that the specific star formation rate does not evolve very rapidly with either mass or redshift, in agreement with observation. This is in contrast to results from hydrodynamical simulations where the star formation closely follows the accretion rate, and so increases strongly towards high redshift. The bursty star formation model naturally explains the observation that there is not enough stellar mass at z~2-4 to account for all of the star-formation observed, without invoking properties like an evolving initial mass function of stars. The finding of a duty cycle that is ~10% implies that there should be ten times the number of known galaxies at fixed stellar mass that have not yet been detected through standard UV selection at high redshift. We therefore predict the existence of a large undetected population of UV-faint galaxies that accounts for most of the stellar mass density at z~4-8.

preprint2012arXiv

A suppressed contribution of low mass galaxies to reionization due to supernova feedback

Motivated by recent observations of the star formation rate density function out to z~7, we describe a simple model for the star formation rate density function at high redshift based on the extended Press-Schechter formalism. This model postulates a starburst following each major merger, lasting for a time t_SF and converting at most f_star of galactic gas into stars. We include a simple physical prescription for supernovae feedback that suppresses star formation in low mass galaxies. Constraining t_SF and f_star to describe the observed star formation rate density at high redshifts, we find that individual starbursts were terminated after a time t_SF~10^7 years. This is comparable to the main-sequence lifetimes of supernova progenitors, indicating that high redshift starbursts are quenched once supernovae feedback had time to develop. High redshift galaxies convert ~10% of their mass into stars for galaxies with star formation rates above ~1 solar mass per year, but a smaller fraction for lower luminosity galaxies. Our best fit model successfully predicts the observed relation between star formation rate and stellar mass at z>~4, while our deduced relation between stellar mass and halo mass is also consistent with data on the dwarf satellites of the Milky Way. We find that supernovae feedback lowers the efficiency of star formation in the lowest mass galaxies and makes their contribution to reionization small. As a result, photo-ionization feedback on low mass galaxy formation does not significantly affect the reionization history. Using a semi-analytic model for the reionization history, we infer that approximately half of the ionizing photons needed to complete reionization have already been observed in star-forming galaxies.

preprint2012arXiv

Fast Holographic Deconvolution: a new technique for precision radio interferometry

We introduce the Fast Holographic Deconvolution method for analyzing interferometric radio data. Our new method is an extension of A-projection/software-holography/forward modeling analysis techniques and shares their precision deconvolution and widefield polarimetry, while being significantly faster than current implementations that use full direction-dependent antenna gains. Using data from the MWA 32 antenna prototype, we demonstrate the effectiveness and precision of our new algorithm. Fast Holographic Deconvolution may be particularly important for upcoming 21 cm cosmology observations of the Epoch of Reionization and Dark Energy where foreground subtraction is intimately related to the precision of the data reduction.

preprint2011arXiv

An Empirical Study of the Relationship between Lyα and UV selected Galaxies: Do Theorists and Observers `Select&#39; the Same Objects?

Lyman Alpha Emitters (LAEs) are galaxies that have been selected on the basis of a strong Lyα emission line in their spectra. Observational campaigns have dramatically increased the sample of known LAEs, which now extends out to z=7. These discoveries have motivated numerous theoretical studies on the subject, which usually define LAEs in their models based on sharp Lyα luminosity and equivalent width (EW) cuts. While broadly representative, this procedure does not mimic the selection from observational programs in detail, which instead use cuts in various colour-spaces. We investigate what implications this disjoint may have for studies that aim to model LAEs. We construct an empirical model for the number density of star forming galaxies as a function of their UV and Lyα luminosity, utilising measured constraints on the luminosity functions (LFs) of drop-out galaxies, and their luminosity dependent probability distribution function of Lyα EW. In particular, we investigate whether the LAE LFs can be reproduced by defining LAEs using a (z-dependent) Lyα luminosity and EW threshold. While we are able to reproduce the observed distribution of Lyα EW among LAEs out to restframe EW 200 A, we find that our formalism over-predicts both the UV and Lyα LFs of LAEs by a factor of 2-3, and is inconsistent with observations at the ~95% level. This tension is partially resolved if we assume the Lyα EW-distribution of drop-out galaxies to be truncated at restframe EW>150 A. However the overprediction indicates that modeling LAEs with simple REW and luminosity cuts does not accurately mimic observed selection criteria, and can lead to uncertainties in the predicted number density of LAEs. On the other hand, the predicted z-evolution is not affected. We apply our formalism to drop-out galaxies at z>6, and predict the LFs of LAEs at z=7-9.

preprint2011arXiv

Non-Gravitational Contributions to the Clustering of Ly-alpha Selected Galaxies: Implications for Cosmological Surveys

We show that the dependence of Ly-alpha absorption on environment leads to significant non-gravitational features in the redshift space power-spectrum of Ly-alpha selected galaxies. We derive a physically motivated fitting formula that can be included in clustering analyses, and use this to discuss the predicted features in the Ly-alpha galaxy power-spectrum based on detailed models in which Ly-alpha absorption is influenced by gas infall and/or by strong galactic outflows. We show that power-spectrum measurements could be used to study the astrophysics of the galaxy-IGM connection, and to measure the properties of outflows from star-forming galaxies. Applying the modified redshift space power-spectrum to a Ly-alpha survey with parameters corresponding to the planned Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), we find that the dependence of observed Ly-alpha flux on velocity gradient and ionising background may compromise the ability of Ly-alpha selected galaxy redshift surveys to constrain cosmology using information from the full power-spectrum. This is because the effects of fluctuating ionizing background and velocity gradients effect the shape of the observed power-spectrum in ways that are similar to the shape of the primordial power-spectrum and redshift space distortions respectively. We use the Alcock-Paczynski test to show that without prior knowledge of the details of Ly-alpha absorption in the IGM, the precision of line-of-sight and transverse distance measurements for HETDEX will be ~1.3-1.7%, decreased by a factor of ~1.5-2 relative to the best case precision of ~0.8% available in a traditional galaxy redshift survey. We specify the precision with which modelling of Ly-alpha radiative transfer must be understood in order for HETDEX to achieve distance measurements that are better than 1%.

preprint2011arXiv

Photon Trapping Enables Super-Eddington Growth of Black-Hole Seeds in Galaxies at High Redshift

We identify a physical mechanism that would have resulted in rapid, obscured growth of seed super-massive black-holes in galaxies at z>6. Specifically, we find that the density at the centre of typical high redshift galaxies was at a level where the Bondi accretion rate implies a diffusion speed of photons that was slower than the gravitational infall velocity, resulting in photons being trapped within the accretion flow and advected into the black-hole. We show that there is a range of black-hole masses (M_bh ~ 10^3-10^5 solar masses) where the accretion flow traps radiation, corresponding to black-holes that were massive enough to generate a photon trapping accretion flow, but small enough that their Bondi radii did not exceed the isothermal scale height of self-gravitating gas. Under these conditions we find that the accretion reaches levels far in excess of the Eddington rate. A prediction of this scenario is that X-ray number counts of active galactic nuclei at z>6 would exhibit a cutoff at the low luminosities corresponding to black-hole masses below ~10^5 solar masses. At low redshifts we find photon trapping to be unimportant because it could only occur in rare low spin halos, and would require black-hole masses in excess of expectations from the observed black-hole - halo mass relation. The super-Eddington growth of ~10^5 solar mass seed black-holes at high redshift may have provided a natural acceleration towards the growth of super-massive black-holes at z~6-7, less than a billion years after the Big Bang.

preprint2011arXiv

The Detectability of Lyman Alpha Emission from Galaxies during the Epoch of Reionization

We study the visibility of the Lyman Alpha (Lya) emission line during the Epoch of Reionization (EoR). Combining galactic outflow models with large-scale semi-numeric simulations of reionization, we quantify the probability distribution function (PDF) of the fraction of Lya photons transmitted through the intergalactic medium (IGM), T_IGM. Our study focusses on galaxies populating dark matter halos with masses of M_halo=1e10 M_sun at z=8.6, which is inspired by the recent reported discovery of a galaxy at z=8.6 with strong Lya line emission. For reasonable model assumptions, we find that winds cause T_IGM>10% [50%], for the majority of galaxies, even when the Universe is ~80% [60%] neutral by volume. Thus, the observed strong Lya emission from the reported z=8.6 galaxy is consistent with a highly neutral IGM. We also investigate the implications of the recent tentative evidence for a observed decrease in the `LAE fraction&#39; among drop-out galaxies between z=6 and z=7. If confirmed, we show that a rapid evolution in x_HI will be required to explain this observation via the effects of reionization.

preprint2011arXiv

The Shocking Truth: The small contribution to hydrogen reionization from gravitational infall

It is commonly thought that stars are responsible for reionizing the Universe. However, deep searches for star-forming galaxies during the epoch of reionization have not yet found sufficient galaxies to provide the necessary ionising flux. Fast accretion shocks associated with gravitational infall of baryons during the formation of galaxies have recently been proposed as an alternative method of generating the required ionising photons. In this Letter we assess the contribution to hydrogen reionization from shocked gas associated with gravitational infall. We find that shocks can ionize no more than a few percent of the cosmic hydrogen by z~6. However, the small fraction of ionizing radiation produced by fast accretion shocks would be significantly more biased than that associated with stars, leading to a modification of the luminosity weighted source clustering by ~10%. This modification of the bias may be measurable with future precision experiments utilising the redshifted 21cm line to study the distribution of hydrogen during the reionization era.

preprint2010arXiv

A gravitational lensing explanation for the excess of strong Mg-II absorbers in GRB afterglow spectra

GRB afterglows offer a probe of the intergalactic medium out to high redshift which complements observations along more abundant quasar lines-of-sight. Although both quasars and GRB afterglows should provide a-priori random sight-lines through the intervening IGM, it has been observed that strong Mg-II absorbers are twice as likely to be found along sight-lines toward GRBs. Several proposals to reconcile this discrepancy have been put forward, but none has been found sufficient to explain the magnitude of the effect. In this paper we estimate the effect of gravitational lensing by galaxies and their surrounding mass distributions on the statistics of Mg-II absorption. We find that the multi-band magnification bias could be very strong in the spectroscopic GRB afterglow population and that gravitational lensing can explain the discrepancy in density of absorbers, for plausibly steep luminosity functions. The model makes the prediction that approximately 20%-60% of the spectroscopic afterglow sample (i.e. ~ 5-15 of 26 sources) would have been multiply imaged, and hence result in repeating bursts. We show that despite this large lensing fraction it is likely that none would yet have been identified by chance owing to the finite sky coverage of GRB searches. We predict that continued optical monitoring of the bright GRB afterglow locations in the months and years following the initial decay would lead to identification of lensed GRB afterglows. A confirmation of the lensing hypothesis would allow us to constrain the GRB luminosity function down to otherwise inaccessibly faint levels, with potential consequences for GRB models.

preprint2010arXiv

Extrapolating the Evolution of Galaxy Sizes to the Epoch of Reionization

We use data on the high-redshift evolution of the size distribution and luminosity function of galaxies to constrain the relationship between their star formation efficiency and starburst lifetime. Based on the derived scaling relations, we predict the angular sizes and average surface brightnesses of faint galaxies that will be discovered with JWST. We find that JWST will be able to resolve galaxies at the magnitude limit m<31 out to a redshift of z~14. The next generation of large ground-based telescopes will resolve all galaxies discovered with JWST, provided they are sufficiently clumpy to enable detection above the bright thermal sky. We combine our constraints with simple models for self regulation of star formation, and show that feedback from supernovae at redshifts z>3 is likely mediated through momentum transfer, with the starburst timescale set by the lifetime of the massive stars rather than the dynamical time in the host galactic disk.

preprint2010arXiv

Near-zone sizes and the rest frame extreme ultra-violet spectral index of the highest redshift quasars

The discovery of quasars with redshifts higher than six has prompted a great deal of discussion in the literature regarding the role of quasars, both as sources of reionization, and as probes of the ionization state of the IGM. However the extreme ultra-violet (EUV) spectral index cannot be measured directly for high redshift quasars owing to absorption at frequencies above the Lyman limit, and as a result, studies of the impact of quasars on the intergalactic medium during reionization must assume a spectral energy distribution in the extreme ultra-violet based on observations at low redshift, z<1. In this paper we use regions of high Ly-alpha transmission (near-zones) around the highest redshift quasars to measure the quasar EUV spectral index at z~6. We jointly fit the available observations for variation of near-zone size with both redshift and luminosity, and propose that the observed relation provides evidence for an EUV spectral index that varies with absolute magnitude in the high redshift quasar sample, becoming softer at higher luminosity. Using a large suite of detailed numerical simulations we find that the typical value of spectral index for a luminous quasar at z~6 is constrained to be alpha=1.3+/-0.4 for a specific luminosity of the form L\proptoν^{-alpha}. We find the scatter in spectral index among individual quasars to be in the range ~0.75-1.25. These values are in agreement with direct observations at low redshift, and indicate that there has been no significant evolution in the EUV spectral index of quasars over 90% of cosmic time.

preprint2010arXiv

Seeing Through the Trough: Outflows and the Detectability of Lyman Alpha Emission from the First Galaxies

The next generation of telescopes aim to directly observe the first generation of galaxies that initiated the reionization process in our Universe. The Lyman Alpha (Lya) emission line is robustly predicted to be the most prominent intrinsic spectral feature of these galaxies, making it an ideal target to search for and study high redshift galaxies. Unfortunately the large Gunn-Peterson optical depth of the surrounding neutral intergalactic medium (IGM) is thought to render this line extremely difficult to detect prior to reionization. In this paper we demonstrate that the radiative transfer effects in the interstellar medium (ISM), which cause Lya flux to emerge from galaxies at frequencies where the Gunn-Peterson optical depth is reduced, can substantially enhance the prospects for detection of the Lya line at high redshift. In particular, scattering off outflows of interstellar HI gas can modify the Lya spectral line shape such that >5% of the emitted Lya radiation is transmitted directly to the observer, even through a fully neutral IGM. It may therefore be possible to directly observe `strong&#39; Lya emission lines (EW > 50 Angstrom rest frame) from the highest redshift galaxies that reside in the smallest HII `bubbles&#39; early in the reionization era with JWST. In addition, we show that outflows can boost the fraction of Lya radiation that is transmitted through the IGM during the latter stages of reionization, and even post-reionization. Coupled with the fact that the first generation of galaxies are thought to have very large intrinsic equivalent Lya equivalent widths, our results suggest that the search for galaxies in their redshifted Lya emission line can be competitive with the drop-out technique out to the highest redshifts that can be probed in the JWST era.

preprint2009arXiv

Determining the escape fraction of ionizing photons during reionization with the GRB derived star-formation rate

The fraction of ionizing photons that escape their host galaxies and so are able to ionize hydrogen in the inter-galactic medium (IGM) is a critical parameter in analyses of the reionization era and early galaxy formation. Studies of the reionization history normally suffer from a degeneracy between the unknown values for the efficiency with which high redshift galaxies turn mass into stars and the escape fraction of ionizing photons. Recent gamma-ray burst (GRB) measurements of the star formation rate density during reionization provide the first opportunity to break this degeneracy. We confront a semi-analytic model for reionization with the GRB-derived star formation rate, as well as observations of the Ly-alpha forest and the CMB. Assuming that UV photons produced in star-forming galaxies dominate the reionization process, we show that the escape fraction of ionizing photons from high redshift galaxies is ~5% [log(f_esc)=-1.35+/-0.15 (68%)] for our fiducial model. This value is reasonably stable against uncertainties in the modeling, including the implementation of radiative feedback, the possibility of an evolving escape fraction, and the unknown shape of the IMF, which in sum contribute ~0.2 dex of additional systematic uncertainty on the value of escape fraction.

preprint2009arXiv

Modification of the 21-cm power spectrum by X-rays during the epoch of reionisation

We incorporate a contribution to reionization from X-rays within analytic and semi-numerical simulations of the 21-cm signal arising from neutral hydrogen during the epoch of reionization. We explore the impact that X-ray ionizations have on the power spectrum (PS) of 21-cm fluctuations by varying both the average X-ray MFP and the fractional contribution of X-rays to reionization. In general, prior to the epoch when the intergalactic medium is dominated by ionized regions (H {\sevensize II} regions), X-ray-induced ionization enhances fluctuations on spatial scales smaller than the X-ray MFP, provided that X-ray heating does not strongly supress galaxy formation. Conversely, at later times when \H2 regions dominate, small-scale fluctuations in the 21-cm signal are suppressed by X-ray ionization. Our modelling also shows that the modification of the 21-cm signal due to the presence of X-rays is sensitive to the relative scales of the X-ray MFP, and the characteristic size of \H2 regions. We therefore find that X-rays imprint an epoch and scale-dependent signature on the 21-cm PS, whose prominence depends on fractional X-ray contribution. The degree of X-ray heating of the IGM also determines the extent to which these features can be discerned. We show that the MWA will have sufficient sensitivity to detect this modification of the PS, so long as the X-ray photon MFP falls within the range of scales over which the array is most sensitive ($\sim0.1$ Mpc$^{-1}$). In cases in which this MFP takes a much smaller value, an array with larger collecting area would be required.

preprint2009arXiv

On merger bias and the clustering of quasars

We use the large catalogues of haloes available for the Millennium Simulation to test whether recently merged haloes exhibit stronger large-scale clustering than other haloes of the same mass. This effect could help to understand the very strong clustering of quasars at high redshift. However, we find no statistically significant excess bias for recently merged haloes over the redshift range 2 < z < 5, with the most massive haloes showing an excess of at most ~5%. We also consider galaxies extracted from a semianalytic model built on the Millennium Simulation. At fixed stellar mass, we find an excess bias of ~ 20-30% for recently merged objects, decreasing with increasing stellar mass. The fact that recently-merged galaxies are found in systematically more massive haloes than other galaxies of the same stellar mass accounts for about half of this signal, and perhaps more for high-mass galaxies. The weak merger bias of massive systems suggests that objects of merger-driven nature, such as quasars, do not cluster significantly differently than other objects of the same characteristic mass. We discuss the implications of these results for the interpretation of clustering data with respect to quasar duty cycles, visibility times, and evolution in the black hole-host mass relation.

preprint2009arXiv

The 21cm Power Spectrum After Reionization

We discuss the 21cm power spectrum (PS) following the completion of reionization. In contrast to the reionization era, this PS is proportional to the PS of mass density fluctuations, with only a small modulation due to fluctuations in the ionization field on scales larger than the mean-free-path of ionizing photons. We derive the form of this modulation, and demonstrate that its effect on the 21cm PS will be smaller than 1% for physically plausible models of damped Ly-alpha systems. In contrast to the 21cm PS observed prior to reionization, in which HII regions dominate the ionization structure, the simplicity of the 21cm PS after reionization will enhance its utility as a cosmological probe by removing the need to separate the PS into physical and astrophysical components. As a demonstration, we consider the Alcock-Paczynski test and show that the next generation of low-frequency arrays could measure the angular distortion of the PS at the percent level for z~3-5.

preprint2009arXiv

The Halo Occupation Distribution of HI from 21cm Intensity Mapping at Moderate Redshifts

The spatial clustering properties of HI galaxies can be studied using the formalism of the halo occupation distribution (HOD). The resulting parameter constraints describe properties like gas richness verses environment. Unfortunately, clustering studies based on individual HI galaxies will be restricted to the local Universe for the foreseeable future, even with the deepest HI surveys. Here we discuss how clustering studies of the HI HOD could be extended to moderate redshift, through observations of fluctuations in the combined 21cm intensity of unresolved galaxies. In particular we make an analytic estimate for the clustering of HI in the HOD. Our joint goals are to estimate i) the amplitude of the signal, and ii) the sensitivity of telescopes like the Australian SKA Pathfinder to HOD parameters. We find that the power spectrum of redshifted 21cm intensity could be used to study the distribution of HI within dark matter halos at z>0.5 where individual galaxies cannot be detected. In addition to the HOD of HI, the amplitude of the 21cm power spectrum would also yield estimates of the cosmic HI content at epochs between the local Universe, and redshifts probed by damped Ly-alpha absorbers.

preprint2008arXiv

A Method to Measure the Mass of Damped Ly-alpha Absorber Host Galaxies Using Fluctuations in 21cm Emission

Observations of damped Ly-alpha absorbers (DLA) indicate that the fraction of hydrogen in its neutral form (HI) is significant by mass at all redshifts. This gas represents the reservoir of material that is available for star formation at late times. As a result, observational identification of the systems in which this neutral hydrogen resides is an important missing ingredient in models of galaxy formation. Precise identification of DLA host mass via traditional clustering studies is not practical owing to the small numbers of known systems being spread across sparsely distributed sight lines. However following the completion of reionization, 21cm surface brightness fluctuations will be dominated by neutral hydrogen in DLAs. Observations of these fluctuations will measure the combined clustering signal from all DLAs within a large volume. We show that measurement of the spherically averaged power-spectrum of 21cm intensity fluctuations due to DLAs could be used to measure the galaxy bias for DLA host galaxies when combined with an independent measurement of the cosmological HI mass density from quasar absorption studies. Utilising this technique, the low frequency arrays now under construction could measure the characteristic DLA host mass with a statistical precision as low as 0.3 dex at z~4. In addition, high signal-to-noise observations of the peculiar-motion induced anisotropy of the power-spectrum would facilitate measurement of both the DLA host mass and the cosmic HI density directly from 21cm fluctuations. By exploiting this anisotropy, a second generation of low frequency arrays with an order of magnitude increase in collecting area could measure the values of cosmic HI density and DLA host mass, with uncertainties of a few percent and a few tens of percent respectively.

preprint2008arXiv

Evidence for Merger-Driven Activity in the Clustering of High Redshift Quasars

Recently, a very large clustering length has been measured for quasars at a redshift of z~4. In combination with the observed quasar luminosity function we assess the implications of this clustering for the relationship between quasar luminosity and dark matter halo mass. Our analysis allows for non-linearity and finite scatter in the relation between quasar luminosity and halo mass, as well as a luminosity dependent quasar lifetime. The additional novel ingredient in our modelling is the allowance for an excess in the observed bias over the underlying halo bias owing to the merger driven nature of quasar activity. We find that the observations of clustering and luminosity function can be explained only if both of the following conditions hold: (i) The luminosity to halo mass ratio increases with halo mass; (ii) The observed clustering amplitude is in excess of that expected solely from halo bias. The latter result is statistically significant at the 99% level. Taken together, the observations provide compelling evidence for merger driven quasar activity, with a black-hole growth that is limited by feedback. In difference from previous analyses, we show that there could be scatter in the luminosity halo mass relation of up to 1 dex, and that quasar clustering can not be used to estimate the quasar lifetime.

preprint2008arXiv

The escape fraction of ionizing photons from high redshift galaxies

The fraction of ionizing photons which escape their host galaxy and so are able to ionize hydrogen in the inter-galactic medium (IGM) is a critical parameter in studies of the reionization era and early galaxy formation. In this paper we combine observations of Lyman-alpha absorption towards high redshift quasars with the measured UV luminosity function of high redshift galaxies to constrain the escape fraction (f_esc) of ionizing photons from galaxies at z ~ 5.5-6. The observed Lyman-alpha transmission constrains the escape fraction to lie in the range f_esc ~ 10-25 % (at z ~ 5.5-6). Excluding halos with M< 10^10 M_sun (as might be expected if galaxy formation is suppressed due to the reionization of the IGM) implies a larger escape fraction of f_esc ~ 20-45 %. Using the numerical results to calibrate an analytic relation between the escape fraction and minimum galaxy halo mass we also extrapolate our results to a mass (M~10^8 M_sun) corresponding to the hydrogen cooling threshold. In this case we find f_esc ~ 5-10 %, consistent with observed estimates at lower redshift. We find that the escape fraction of high redshift galaxies must be greater than 5 % irrespepctive of galaxy mass. Based on these results we use a semi-analytic description to model the reionization history of the IGM, assuming ionizing sources with escape fractions suggested by our numerical simulations. We find that the IBG observed at z ~ 5.5-6 implies a sufficient number of ionizing photons to have reionized the Universe by z ~ 6. However, if the minimum mass for star-formation were greater than 10^9 M_sun, the IBG would be over-produced at redshifts less than z ~ 5. In summary, our results support a scenario in which the IGM was reionized by low mass galaxies.

preprint2008arXiv

The Impact of HI in Galaxies on 21-cm Intensity Fluctuations During the Reionisation Epoch

We investigate the impact of neutral hydrogen (HI) in galaxies on the statistics of 21-cm fluctuations using analytic and semi-numerical modelling. Following the reionisation of hydrogen the HI content of the Universe is dominated by damped absorption systems (DLAs), with a cosmic density in HI that is observed to be constant at a level equal to ~2% of the cosmic baryon density from z~1 to z~5. We show that extrapolation of this constant fraction into the reionisation epoch results in a reduction of 10-20% in the amplitude of 21-cm fluctuations over a range of spatial scales. The assumption of a different percentage during the reionisation era results in a proportional change in the 21-cm fluctuation amplitude. We find that consideration of HI in galaxies/DLAs reduces the prominence of the HII region induced shoulder in the 21-cm power spectrum (PS), and hence modifies the scale dependence of 21-cm fluctuations. We also estimate the 21cm-galaxy cross PS, and show that the cross PS changes sign on scales corresponding to the HII regions. From consideration of the sensitivity for forthcoming low-frequency arrays we find that the effects of HI in galaxies/DLAs on the statistics of 21-cm fluctuations will be significant with respect to the precision of a PS or cross PS measurement. In addition, since overdense regions are reionised first we demonstrate that the cross-correlation between galaxies and 21-cm emission changes sign at the end of the reionisation era, providing an alternative avenue to pinpoint the end of reionisation. The sum of our analysis indicates that the HI content of the galaxies that reionise the universe will need to be considered in detailed modelling of the 21-cm intensity PS in order to correctly interpret measurements from forthcoming low-frequency arrays.

preprint2007arXiv

Baryonic Acoustic Oscillations in 21cm Emission: A Probe of Dark Energy out to High Redshifts

Low-frequency observatories are currently being constructed with the goal of detecting redshifted 21cm emission from the epoch of reionization. These observatories will also be able to detect intensity fluctuations in the cumulative 21cm emission after reionization, from hydrogen in unresolved damped Ly-alpha absorbers (such as gas rich galaxies) down to a redshift z~3.5. The inferred power spectrum of 21cm fluctuations at all redshifts will show acoustic oscillations, whose co-moving scale can be used as a standard ruler to infer the evolution of the equation of state for the dark energy. We find that the first generation of low-frequency experiments (such as MWA or LOFAR) will be able to constrain the acoustic scale to within a few percent in a redshift window just prior to the end of the reionization era, provided that foregrounds can be removed over frequency band-passes of >8MHz. This sensitivity to the acoustic scale is comparable to the best current measurements from galaxy redshift surveys, but at much higher redshifts. Future extensions of the first generation experiments (involving an order of magnitude increase in the antennae number of the MWA) could reach sensitivities below one percent in several redshift windows and could be used to study the dark energy in the unexplored redshift regime of 3.5<z<12. Moreover, new experiments with antennae designed to operate at higher frequencies would allow precision measurements (<1%) of the acoustic peak to be made at more moderate redshifts (1.5<z<3.5), where they would be competitive with ambitious spectroscopic galaxy surveys covering more than 1000 square degrees. Together with other data sets, observations of 21cm fluctuations will allow full coverage of the acoustic scale from the present time out to z~12.

preprint2007arXiv

Biased Reionisation and Non-Gaussianity in Redshifted 21cm Intensity Maps of the Reionisation Epoch

Spatial dependence in the statistics of redshifted 21cm fluctuations promises to provide the most powerful probe of the reionisation epoch. In this paper we consider the second and third moments of the redshifted 21cm intensity distribution using a simple model that accounts for galaxy bias during the reionisation process. We demonstrate that skewness in redshifted 21cm maps should be substantial throughout the reionisation epoch and on all angular scales, owing to the effects of galaxy bias which leads to early reionisation in over-dense regions of the IGM. The variance (or power-spectrum) of 21cm fluctuations will exhibit a minimum in redshift part way through the reionisation process, when the global ionisation fraction is around 50%. This minimum is generic, and is due to the transition from 21cm intensity being dominated by over-dense too under-dense regions as reionisation progresses. We show that the details of the reionisation history, including the presence of radiative feedback are encoded in the evolution of the auto-correlation and skewness functions with redshift and mean IGM neutral fraction. We discuss the detection of skewness by first generation instruments, and conclude that the Mileura Widefield Array - Low Frequency Demonstrator will have sufficient sensitivity to detect skewness on a range of angular scales at redshifts near the end of reionisation, while a subsequent instrument of 10 times the collecting area could map out the evolution of skewness in detail. The observation of a minimum in variance during the reionisation history, and the detection of skewness would both provide important confirmation of the cosmological origin of redshifted 21cm intensity fluctuations.

preprint2007arXiv

Fluctuations in 21cm Emission After Reionization

The fluctuations in the emission of redshifted 21cm photons from neutral inter-galactic hydrogen will provide an unprecedented probe of the reionization era. Conventional wisdom assumes that this 21cm signal disappears as soon as reionization is complete, when little atomic hydrogen is left through most of the volume of the inter-galactic medium (IGM). However observations of damped Ly-alpha absorbers indicate that the fraction of hydrogen in its neutral form is significant by mass at all redshifts. Here we use a physically-motivated model to show that residual neutral gas, confined to dense regions in the IGM with a high recombination rate, will generate a significant post-reionization 21cm signal. We show that the power-spectrum of fluctuations in this signal will be detectable by the first generation of low-frequency observatories at a signal-to-noise that is comparable to that achievable in observations of the reionization era. The statistics of 21cm fluctuations will therefore probe not only the pre-reionization IGM, but rather the entire process of HII region overlap, as well as the appearance of the diffuse ionized IGM.

preprint2007arXiv

Reionization Bias in High Redshift Quasar Near-Zones

Absorption spectra of high redshift quasars exhibit an increasingly thick Ly-alpha forest towards z~6. However, the interpretation of these spectra is complicated by the fact that the Ly-alpha optical depth is already large for neutral hydrogen fractions in excess of 10^-4, and also because quasars are expected to reside in dense regions of the IGM. We present a model for the evolution of the ionization state of the IGM which is applicable to the dense, biased regions around high-redshift quasars as well as more typical regions in the IGM, and combine this with numerical radiative transfer simulations. Our model is able to simultaneously reproduce the observed Ly-alpha forest opacity at 4<z<6, the ionizing photon mean-free-path at z~4 and the rapid evolution of highly ionized near-zone sizes around high-redshift quasars at 5.8<z<6.4. We find that within 5 physical Mpc of a high redshift quasar, the evolution of the ionization state of the IGM precedes that in more typical regions by around 0.3 redshift units. More importantly, when combined with the rapid increase in the ionizing photon mean-free-path expected shortly after overlap, this offset results in an ionizing background near the quasar which exceeds the value in the rest of the IGM by a factor of ~2-3. We further find that in the post-overlap phase of reionization the size of the observed quasar near-zones is not directly sensitive to the neutral hydrogen fraction of the IGM. Instead, these sizes probe the level of the background ionization rate and the temperature of the surrounding IGM. The observed rapid evolution of the quasar near-zone sizes at 5.8<z<6.4 can thus be explained by the rapid evolution of the ionizing background, which in our model is caused by the completion of overlap at the end of reionization by 6<z<7.

preprint2006arXiv

Constraining the Quasar Contribution to the Reionisation of Cosmic Hydrogen

Absorption spectra of high redshift quasars suggest that the reionisation of cosmic hydrogen was complete near z~6. The dominant sources of ionising photons responsible for this reionisation are generally thought to be stars and quasars. In this paper we make a quantitative estimate of the relative contributions made by these sources. Our approach is to compute the evolution of the post overlap ionising background radiation by combining semi-analytic descriptions of reionisation in a clumpy medium with a model for the quasar luminosity function. Our overall model has two free parameters, the star formation efficiency and the minimum quasar luminosity. By adjusting these parameters, we constrain the relative contributions made by stars and quasars through comparison with reported observations (Fan et al. 2005). We find that the relative quasar contribution (at z=5.7) to the ionising background was between 1.4% and 14.5%. The range of uncertainty is dominated by the unknown minimum quasar luminosity.

preprint2006arXiv

Dwarf Galaxy Formation Was Suppressed By Cosmic Reionization

A large number of faint galaxies, born less than a billion years after the big bang, have recently been discovered. The fluctuations in the distribution of these galaxies contributed to a scatter in the ionization fraction of cosmic hydrogen on scales of tens of Mpc, as observed along the lines of sight to the earliest known quasars. Theoretical simulations predict that the formation of dwarf galaxies should have been suppressed after cosmic hydrogen was reionized, leading to a drop in the cosmic star formation rate. Here we present evidence for this suppression. We show that the post-reionization galaxies which produced most of the ionizing radiation at a redshift z~5.5, must have had a mass in excess of ~10^{10.6+/-0.4} solar masses or else the aforementioned scatter would have been smaller than observed. This limiting mass is two orders of magnitude larger than the galaxy mass that is thought to have dominated the reionization of cosmic hydrogen (~10^8 solar masses). We predict that future surveys with space-based infrared telescopes will detect a population of smaller galaxies that reionized the Universe at an earlier time, prior to the epoch of dwarf galaxy suppression.

preprint2006arXiv

Properties of High Redshift Quasars-II: What does the quasar luminosity function tell us about super-massive black-hole evolution?

In the local universe, the masses of Super-Massive Black-Holes (SMBH) appear to correlate with the physical properties of their hosts, including the mass of the dark-matter halo. Using these clues as a starting point many studies have produced models that can explain phenomena like the quasar luminosity function. The shortcoming of this approach is that working models are not unique, and as a result it is not always clear what input physics is being constrained. Here we take a different approach. We identify critical parameters that describe the evolution of SMBHs at high redshift, and constrain their parameter space based on observations of high redshift quasars from the Sloan Digital Sky Survey. We find that the luminosity function taken in isolation is somewhat limited in its ability to constrain SMBH evolution due to some strong degeneracies. This explains the presence in the literature of a range of equally successful models based on different physical hypotheses. Including the constraint of the local SMBH to halo mass ratio breaks some of the degeneracies, and our results suggest halo masses at z~4.8 of 10^{12.5+/-0.3}M_solar (with 90% confidence), with a SMBH to halo mass ratio that decreases with time (>99%). We also find a quasar luminosity to halo mass ratio that increases with halo mass (>99%). These features need to be incorporated in all successful models of SMBH evolution. On the other hand current observations do not permit any conclusions regarding the evolution of quasar lifetime, or the SMBH occupation fraction in dark matter halos.

preprint2006arXiv

Smooth Boundaries to Cosmological HII Regions from Galaxy Clustering

The HII regions around quasars and galaxies at redshifts beyond the epoch of reionisation will provide prime targets for upcoming 21cm campaigns using a new generation of low-frequency radio observatories. Here we show that the boundaries of these HII regions will not be sharp. Rather, the clustering of sources near massive galaxies results in a neutral fraction that rises gradually towards large radii from an interior value near zero. A neutral fraction corresponding to the global background value is typically reached at a distance of 2-5 times the radius of the HII region around the central massive galaxy.

preprint2006arXiv

The Extended Starformation History of the First Generation of Stars, and the Reionization of Cosmic Hydrogen

Population-III (Pop-III) starformation (SF) is thought to be quenched when the metallicity of the star-forming gas reaches a critical level. At high z, when the general intergalactic medium (IGM) was enriched with metals, the fraction of primordial gas already collapsed in minihalos was significantly larger than the fraction of primordial gas that had already been involved in Pop-III SF. We argue that this minihalo gas remained largely in a metal-free state, until these minihalos merged into large systems and formed stars. As a result, the era of Pop-III SF was significantly prolonged, leading to an integrated Pop-III SF an order of magnitude larger than expected for an abrupt transition redshift. The contribution of Pop-III SF to the reionization of hydrogen could have been significant until z~10 and may have extended to z~6. Our modeling allows for gradual enrichment of the IGM, feedback from photo-ionization and screening of reionization by minihalos. Nevertheless, extended Pop-III SF can result in complex, multi-peaked reionization histories. The contribution of Pop-III stars to reionization will be tested by the three-year WMAP results: (1) if Pop-III stars do not contribute to reionization, tau_es<0.05-0.06 and a rapid reionization at z~6 is expected; (2) if the product of star formation efficiency and escape fraction for Pop-III stars is significantly larger than for Pop-II stars, then a maximum tau_es=0.21 is achievable; (3) in a scenario where the product of star formation efficiency and escape fraction for Pop-III stars is comparable to that for Pop-II stars, tau_es=0.09-0.12 would be observed, with reionization histories characterized by an extended ionization plateau from z=7-12. This result holds regardless of the redshift where the IGM becomes enriched with metals.

preprint2005arXiv

Constraints on the Process that Regulates the Growth of Supermassive Black Holes Based on the Intrinsic Scatter in the M_bh-sigma Relation

We show that the observed scatter in the relations between the mass of supermassive black holes (SMBHs), M_bh, and the velocity dispersion sigma or mass M_sph of their host spheroid, place interesting constraints on the process that regulates SMBH growth in galaxies. When combined with the observed properties of early-type SDSS galaxies, the observed intrinsic scatters imply that SMBH growth is regulated by the spheroid velocity dispersion rather than its mass. The M_bh-M_sph relation is therefore a by-product of a more fundamental M_bh-sigma relation. We construct a theoretical model for the scatter among baryon modified dark matter halo profiles, out of which we generate a population of spheroid hosts and show that these naturally lead to a relation between effective radius and velocity dispersion of the form R_sph ~ sigma^1.5 with a scatter of ~0.2dex, in agreement with the corresponding projection of the fundamental plane for early type galaxies in SDSS. At the redshift of formation, our model predicts the minimum scatter that SMBHs can have at fixed velocity dispersion or spheroid mass under different formation scenarios. We also estimate the additional scatter that is introduced into these relations through collisionless mergers of purely stellar spheroids at z<1. We find that the observed scatter in the M_bh-sigma and M_bh-M_sph relations preclude the properties of dark matter halos from being the governing factor in SMBH growth. Finally, we show that SMBH growth governed by the properties of the host spheroid can lead to the observed values of scatter in the M_bh-sigma and M_bh-M_sph relations, only if the SMBH growth is limited by momentum or energy feedback over the dynamical time of the host spheroid.

preprint2004arXiv

A Size of ~10 Mpc for the Ionized Bubbles at the End of Cosmic Reionization

The first galaxies to appear in the universe at redshifts z>20 created ionized bubbles in the intergalactic medium of neutral hydrogen left over from the Big-Bang. It is thought that the ionized bubbles grew with time, surrounded clusters of dwarf galaxies and eventually overlapped quickly throughout the universe over a narrow redshift interval near z~6. This event signaled the end of the reionization epoch when the universe was a billion years old. Measuring the hitherto unknown size distribution of the bubbles at their final overlap phase is a focus of forthcoming observational programs aimed at highly redshifted 21cm emission from atomic hydrogen. Here we show that the combined constraints of cosmic variance and causality imply an observed bubble size at the end of the overlap epoch of ~10 physical Mpc, and a scatter in the observed redshift of overlap along different lines-of-sight of ~0.15. This scatter is consistent with observational constraints from recent spectroscopic data on the farthest known quasars. Our novel result implies that future radio experiments should be tuned to a characteristic angular scale of ~0.5 degrees and have a minimum frequency band-width of ~8 MHz for an optimal detection of 21cm flux fluctuations near the end of reionization.

preprint2004arXiv

Cosmic Hydrogen Was Significantly Neutral a Billion Years After the Big Bang

The ionization fraction of cosmic hydrogen, left over from the big bang, provides crucial fossil evidence for when the first stars and quasar black holes formed in the infant universe. Spectra of the two most distant quasars known show nearly complete absorption of photons with wavelengths shorter than the Ly-alpha transition of neutral hydrogen, indicating that hydrogen in the intergalactic medium (IGM) had not been completely ionized at a redshift z~6.3, about a billion years after the big bang. Here we show that the radii of influence of ionizing radiation from these quasars imply that the surrounding IGM had a neutral hydrogen fraction of tens of percent prior to the quasar activity, much higher than previous lower limits of ~0.1%. When combined with the recent inference of a large cumulative optical depth to electron scattering after cosmological recombination from the WMAP data, our result suggests the existence of a second peak in the mean ionization history, potentially due to an early formation episode of the first stars.

preprint2003arXiv

Detection of Gravitational Waves from the Coalescence of Population-III Remnants with Advanced LIGO

The comoving mass density of massive black hole (MBH) remnants from pre-galactic star formation could have been similar in magnitude to the mass-density of supermassive black holes (SMBHs) in the present-day universe. We show that the fraction of MBHs that coalesce during the assembly of SMBHs can be extracted from the rate of ring-down gravitational waves that are detectable by Advanced LIGO. Based on the SMBH formation history inferred from the evolution of the quasar luminosity function, we show that an observed event rate of 1 per year will constrain the SMBH mass fraction that was contributed by MBHs coalescence down to a level of ~10^-6 for 20 solar mass MBH remnants (or ~10^-4 for 260 solar mass remnants).

preprint2003arXiv

Was the Universe Reionized by Massive Population-III Stars?

The WMAP satellite has measured a large optical depth to electron scattering after cosmological recombination of 0.17+-0.04, implying significant reionization of the primordial gas only ~200 million years after the big bang. However, the most recent overlap of intergalactic HII regions must have occured at z<9 based on the Lyman-alpha forest constraint on the thermal history of the intergalactic medium. Here we argue that a first generation of metal-free stars with a heavy (rather than Salpeter) mass function is therefore required to account for much of the inferred optical depth. This conclusion holds if feedback regulates star formation in early dwarf galaxies as observed in present-day dwarfs.

preprint2002arXiv

Gravitational lenses magnify up to one third of the most distant quasars

Exceptionally bright quasars with redshifts up to z=6.28 have recently been discovered. Quasars are thought to be powered by the accretion of gas onto supermassive black holes at the centres of galaxies. Their maximum (Eddington) luminosity is proportional to the mass of the black hole, and so these bright quasars are inferred to have black holes with masses of more than a few billion solar masses. The existence of such massive black holes poses a challenge to models for the formation of structures in the early Universe, as it requires that the black holes would grow so massive in less than a billion years after the Big Bang. Here we show that up to a third of known quasars with z~6 will have their observed flux magnified by a factor of 10 or more through gravitational lensing by galaxies along the line of sight. The inferred abundance of quasar host galaxies, as well as the luminosity density provided by the quasars, are therefore substantially overestimated.

preprint2002arXiv

Measuring the Size of Quasar Broad-Line Clouds Through Time Delay Light-Curve Anomalies of Gravitational Lenses

Intensive monitoring campaigns have recently attempted to measure the time delays between multiple images of gravitational lenses. Some of the resulting light-curves show puzzling low-level, rapid variability which is unique to individual images, superimposed on top of (and concurrent with) longer time-scale intrinsic quasar variations which repeat in all images. We demonstrate that both the amplitude and variability time-scale of the rapid light-curve anomalies, as well as the correlation observed between intrinsic and microlensed variability, are naturally explained by stellar microlensing of a smooth accretion disk which is occulted by optically-thick broad-line clouds. The rapid time-scale is caused by the high velocities of the clouds (~5x10^3 km/s), and the low amplitude results from the large number of clouds covering the magnified or demagnified parts of the disk. The observed amplitudes of variations in specific lenses implies that the number of broad-line clouds that cover ~10% of the quasar sky is ~10^5 per 4 pi steradian. This is comparable to the expected number of broad line clouds in models where the clouds originate from bloated stars.