BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Steven R. Furlanetto

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

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-ph.IMastro-ph

Catalog footprint

What is connected

32works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

32 published item(s)

preprint2022arXiv

A galaxy-free phenomenological model for the 21-cm power spectrum during reionization

Upper limits from the current generation of interferometers targeting the 21-cm signal from high redshifts have recently begun to rule out physically realistic, though still extreme, models of the Epoch of Reionization (EoR). While inferring the detailed properties of the first galaxies is one of the most important motivations for measuring the high-$z$ 21-cm signal, they can also provide useful constraints on the properties of the intergalactic medium (IGM). Motivated by this, we build a simple, phenomenological model for the 21-cm power spectrum that works directly in terms of IGM properties, which bypasses the computationally expensive 3-D semi-numerical modeling generally employed in inference pipelines and avoids explicit assumptions about galaxy properties. The key simplifying assumptions are that (i) the ionization field is binary, and composed of spherical bubbles with an abundance described well by a parametric bubble size distribution, and (ii) that the spin temperature of the ``bulk'' IGM outside bubbles is uniform. Despite the simplicity of the model, the mean ionized fraction and spin temperature of the IGM recovered from mock 21-cm power spectra generated with \textsc{21cmfast} are generally in good agreement with the true input values. This suggests that it is possible to obtain comparable constraints on the IGM using models with very different assumptions, parameters, and priors. Our approach will thus be complementary to semi-numerical models as upper limits continue to improve in the coming years.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A joint measurement of galaxy luminosity functions and large-scale field densities during the Epoch of Reionization

One of the most exciting advances of the current generation of telescopes has been the detection of galaxies during the epoch of reionization, using deep fields that have pushed these instruments to their limits. It is essential to optimize our analyses of these fields in order to extract as much information as possible from them. In particular, standard methods of measuring the galaxy luminosity function discard information on large-scale dark matter density fluctuations, even though this large-scale structure drives galaxy formation and reionization during the Cosmic Dawn. Measuring these densities would provide a bedrock observable, connecting galaxy surveys to theoretical models of the reionization process and structure formation. Here, we use existing Hubble deep field data to simultaneously fit the universal luminosity function and measure large-scale densities for each Hubble deep field at $z =$ 6--8 by directly incorporating priors on the large-scale density field and galaxy bias. Our fit of the universal luminosity function is consistent with previous methods but differs in the details. For the first time, we measure the underlying densities of the survey fields, including the most over/under-dense Hubble fields. We show that the distribution of densities is consistent with current predictions for cosmic variance. This analysis on just 17 fields is a small sample of what will be possible with the James Webb Space Telescope, which will measure hundreds of fields at comparable (or better) depths and at higher redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Automated Detection of Antenna Malfunctions in Large-N Interferometers: A Case Study with the Hydrogen Epoch of Reionization Array

We present a framework for identifying and flagging malfunctioning antennas in large radio interferometers. We outline two distinct categories of metrics designed to detect outliers along known failure modes of large arrays: cross-correlation metrics, based on all antenna pairs, and auto-correlation metrics, based solely on individual antennas. We define and motivate the statistical framework for all metrics used, and present tailored visualizations that aid us in clearly identifying new and existing systematics. We implement these techniques using data from 105 antennas in the Hydrogen Epoch of Reionization Array (HERA) as a case study. Finally, we provide a detailed algorithm for implementing these metrics as flagging tools on real data sets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improved treatments of the ionizing photon mean free path in semi-numerical simulations of reionization

Efficient and accurate simulations of the reionization epoch are crucial to exploring the vast uncharted parameter space that will soon be constrained by measurements of the 21 cm power spectrum. One of these parameters, $R_{\rm max}$, is meant to characterize the absorption of photons by residual neutral gas inside of ionized regions, but has historically been implemented in a very simplistic fashion acting only as a maximum filtering scale. We leverage the correspondence between excursion set methods and the integrated flux from ionizing sources to define two physically-motivated prescriptions of the mean free path of ionizing photons that smoothly attenuate the contribution from distant sources. Implementation of these methods in semi-numerical reionization codes requires only modest additional computational effort due to the fact that spatial filtering is still performed on scales larger than the characteristic absorption distance. We find that our smoothly-defined mean free path prescriptions more effectively suppress large-scale structures in the ionization field in semi-numerical reionization simulations compared to the standard $R_{\rm max}$ approach, and the magnitude of the mean free path modulates the power spectrum in a much smoother manner. We show that this suppression of large-scale power is significant enough to be relevant for upcoming 21 cm power spectrum observations. Finally, we show that in our model the mean free path plays a larger role in regulating the reionization history than in models using $R_{\rm max}$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Framework for Simultaneously Measuring Field Densities and the High-z Luminosity Function

Cosmic variance from large-scale structure will be a major source of uncertainty for galaxy surveys at z > 6, but that same structure will also provide an opportunity to identify and study dense environments in the early Universe. Using a robust model for galaxy clustering, we directly incorporate large-scale densities into an inference framework that simultaneously measures the high-z (z > 6) UV luminosity function and the average matter density of each distinct volume in a survey. Through this framework, we forecast the performance of several major upcoming James Webb Space Telescope (JWST) galaxy surveys. We find that they can constrain field matter densities down to the theoretical limit imposed by Poisson noise and unambiguously identify over-dense (and under-dense) regions on transverse scales of tens of comoving Mpc. We also predict JWST will measure the luminosity function with a precision at z = 12 comparable to existing Hubble Space Telescope's constraints at z = 8 (and even better for the faint-end slope). We also find that wide-field surveys are especially important in distinguishing luminosity function models.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Bursty star formation during the Cosmic Dawn driven by delayed stellar feedback

In recent years, several analytic models have demonstrated that simple assumptions about halo growth and feedback-regulated star formation can match the (limited) existing observational data on galaxies at z>6. By extending such models, we demonstrate that imposing a time delay on stellar feedback (as inevitably occurs in the case of supernova explosions) induces burstiness in small galaxies. Although supernova progenitors have short lifetimes (~5-30 Myr), the delay exceeds the dynamical time of galaxies at such high redshifts. As a result, star formation proceeds unimpeded by feedback for several cycles and "overshoots" the expectations of feedback-regulated star formation models. We show that such overshoot is expected even in atomic cooling halos, with masses up to ~10^10.5 Msun at z>6. However, these burst cycles damp out quickly in massive galaxies, because large haloes are more resistant to feedback so retain a continuous gas supply. Bursts in small galaxies - largely beyond the reach of existing observations - induce a scatter in the luminosity of these haloes (of ~1 mag) and increase the time-averaged star formation efficiency by up to an order of magnitude. This kind of burstiness can have substantial effects on the earliest phases of star formation and reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Deep Realistic Extragalactic Model (DREaM) Galaxy Catalogs: Predictions for a Roman Ultra-Deep Field

In the next decade, deep galaxy surveys from telescopes such as the James Webb Space Telescope and Roman Space Telescope will provide transformational data sets that will greatly enhance the understanding of galaxy formation during the epoch of reionization (EoR). In this work, we present the Deep Realistic Extragalactic Model (DREaM) for creating synthetic galaxy catalogs. Our model combines dark matter simulations, subhalo abundance matching and empirical models, and includes galaxy positions, morphologies, and spectral energy distributions (SEDs). The resulting synthetic catalog extends to redshifts $z \sim 12 $, and galaxy masses $\log_{10}(M/M_{\odot}) = 5 $ covering an area of $1\, {\rm deg}^2$ on the sky. We use DREaM to explore the science returns of a $1\, {\rm deg}^2$ \emph{Roman} UDF, and to provide a resource for optimizing ultra-deep survey designs. We find that a \emph{Roman} UDF to $\sim 30\,m_{\rm AB}$ will potentially detect more than $10^6$ $M_{\rm UV}<-17$ galaxies, with more than $10^4$ at redshifts $z>7$, offering an unparalleled dataset for constraining galaxy properties during the EoR. Our synthetic catalogs and simulated images are made publicly available to provide the community with a tool to prepare for upcoming data.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

First Results from HERA Phase I: Upper Limits on the Epoch of Reionization 21 cm Power Spectrum

We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground power. This yields a 95% confidence upper limit on the 21 cm power spectrum of $Δ^2_{21} \le (30.76)^2\ {\rm mK}^2$ at $k=0.192\ h\ {\rm Mpc}^{-1}$ at $z=7.9$, and also $Δ^2_{21} \le (95.74)^2\ {\rm mK}^2$ at $k=0.256\ h\ {\rm Mpc}^{-1}$ at $z=10.4$. At $z=7.9$, these limits are the most sensitive to-date by over an order of magnitude. While we find evidence for residual systematics at low line-of-sight Fourier $k_\parallel$ modes, at high $k_\parallel$ modes we find our data to be largely consistent with thermal noise, an indicator that the system could benefit from deeper integrations. The observed systematics could be due to radio frequency interference, cable sub-reflections, or residual instrumental cross-coupling, and warrant further study. This analysis emphasizes algorithms that have minimal inherent signal loss, although we do perform a careful accounting in a companion paper of the small forms of loss or bias associated with the pipeline. Overall, these results are a promising first step in the development of a tuned, instrument-specific analysis pipeline for HERA, particularly as Phase II construction is completed en route to reaching the full sensitivity of the experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Validation of the HERA Phase I Epoch of Reionization 21 cm Power Spectrum Software Pipeline

We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 sigma level) for k > 0.2 h/Mpc for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong `detection' at k ~0.2 h/Mpc -- at the level of ~25 sigma -- with foregrounds dominating on larger scales, and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of > 10^7. Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Detection of Cosmic Structures using the Bispectrum Phase. II. First Results from Application to Cosmic Reionization Using the Hydrogen Epoch of Reionization Array

Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1σ$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $κ_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $κ_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Effects of Population III X-ray and Radio Backgrounds on the Cosmological 21-cm Signal

We investigate the effects of Population III (Pop III) stars and their remnants on the cosmological 21-cm global signal. By combining a semi-analytic model of Pop III star formation with a global 21-cm simulation code, we investigate how X-ray and radio emission from accreting Pop III black holes may affect both the timing and depth of the 21-cm absorption feature that follows the initial onset of star formation during the Cosmic Dawn. We compare our results to the findings of the EDGES experiment, which has reported the first detection of a cosmic 21-cm signal. In general, we find that our fiducial Pop III models, which have peak star formation rate densities of $\sim 10^{-4}$ M$_\odot$ yr$^{-1}$ Mpc$^{-3}$ between $z \sim 10$ and $z \sim 15$, are able to match the timing of the EDGES signal quite well, in contrast to models that ignore Pop III stars. To match the unexpectedly large depth of the EDGES signal without recourse to exotic physics, we vary the parameters of emission from accreting black holes (formed as Pop III remnants) including the intrinsic strength of X-ray and radio emission as well as the local column density of neutral gas. We find that models with strong radio emission and relatively weak X-ray emission can self-consistently match the EDGES signal, though this solution requires fine-tuning. We are only able to produce signals with sharp features similar to the EDGES signal if the Pop~III IMF is peaked narrowly around $140 \, M_\odot$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constraints on the Star Formation Efficiency of Galaxies During the Epoch of Reionization

Reionization is thought to have occurred in the redshift range of $6 < z < 9$, which is now being probed by both deep galaxy surveys and CMB observations. Using halo abundance matching over the redshift range $5<z<8$ and assuming smooth, continuous gas accretion, we develop a model for the star formation efficiency $f_{\star}$ of dark matter halos at $z>6$ that matches the measured galaxy luminosity functions at these redshifts. We find that $f_{\star}$ peaks at $\sim 30\%$ at halo masses $M \sim 10^{11}$--$10^{12}$~M$_\odot$, in qualitative agreement with its behavior at lower redshifts. We then investigate the cosmic star formation histories and the corresponding models of reionization for a range of extrapolations to small halo masses. We use a variety of observations to further constrain the characteristics of the galaxy populations, including the escape fraction of UV photons. Our approach provides an empirically-calibrated, physically-motivated model for the properties of star-forming galaxies sourcing the epoch of reionization. In the case where star formation in low-mass halos is maximally efficient, an average escape fraction $\sim0.1$ can reproduce the optical depth reported by Planck, whereas inefficient star formation in these halos requires either about twice as many UV photons to escape, or an escape fraction that increases towards higher redshifts. Our models also predict how future observations with JWST can improve our understanding of these galaxy populations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Reionization Through the Lens of Percolation Theory

The reionization of intergalactic hydrogen has received intense theoretical scrutiny over the past two decades. Here, we approach the process formally as a percolation process and phase transition. Using semi-numeric simulations, we demonstrate that an infinitely-large ionized region abruptly appears at an ionized fraction of ~0.1 and quickly grows to encompass most of the ionized gas: by an ionized fraction of 0.3, nearly ninety percent of the ionized material is part of this region. Throughout most of reionization, nearly all of the intergalactic medium is divided into just two regions, one ionized and one neutral, and both infinite in extent. We also show that the discrete ionized regions that exist before and near this transition point follow a near-power law distribution in volume, with equal contributions to the total filling factor per logarithmic interval in size up to a sharp cutoff in volume. These qualities are generic to percolation processes, with the detailed behavior a result of long-range correlations in the underlying density field. These insights will be crucial to understanding the distribution of ionized and neutral gas during reionization and provide precise meaning to the intuitive description of reionization as an "overlap" process.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Distribution of Bubble Sizes During Reionization

A key physical quantity during reionization is the size of HII regions. Previous studies found a characteristic bubble size which increases rapidly during reionization, with apparent agreement between simulations and analytic excursion set theory. Using four different methods, we critically examine this claim. In particular, we introduce the use of the watershed algorithm -- widely used for void finding in galaxy surveys -- which we show to be an unbiased method with the lowest dispersion and best performance on Monte-Carlo realizations of a known bubble size PDF. We find that a friends-of-friends algorithm declares most of the ionized volume to be occupied by a network of volume-filling regions connected by narrow tunnels. For methods tuned to detect the volume-filling regions, previous apparent agreement between simulations and theory is spurious, and due to a failure to correctly account for the window function of measurement schemes. The discrepancy is already obvious from visual inspection. Instead, HII regions in simulations are significantly larger (by factors of 10-1000 in volume) than analytic predictions. The size PDF is narrower, and evolves more slowly with time, than predicted. It becomes more sharply peaked as reionization progresses. These effects are likely caused by bubble mergers, which are inadequately modeled by analytic theory. Our results have important consequences for high-redshift 21cm observations, the mean free path of ionizing photons, and the visibility of Ly-alpha emitters, and point to a fundamental failure in our understanding of the characteristic scales of the reionization process.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Global 21-cm Signal in the Context of the High-z Galaxy Luminosity Function

Motivated by recent progress in studies of the high-$z$ Universe, we build a new model for the global 21-cm signal that is explicitly calibrated to measurements of the galaxy luminosity function (LF) and further tuned to match the Thomson scattering optical depth of the cosmic microwave background, $τ_e$. Assuming that the $z \lesssim 8$ galaxy population can be smoothly extrapolated to higher redshifts, the recent decline in best-fit values of $τ_e$ and the inefficient heating induced by X-ray binaries (HMXBs; the presumptive sources of the X-ray background at high-$z$) imply that the entirety of cosmic reionization and reheating occurs at redshifts $z \lesssim 12$. In contrast to past global 21-cm models, whose $z \sim 20$ ($ν\sim 70$ MHz) absorption features and strong $\sim 25$ mK emission features were driven largely by the assumption of efficient early star-formation and X-ray heating, our new fiducial model peaks in absorption at $ν\sim 110$ MHz at a depth of $\sim -160$ mK and has a negligible emission component. As a result, a strong emission signal would provide convincing evidence that HMXBs are not the only drivers of cosmic reheating. Shallow absorption troughs should accompany strong heating scenarios, but could also be caused by a low escape fraction of Lyman-Werner photons. Generating signals with troughs at $ν\lesssim 95$ MHz requires a floor in the star-formation efficiency in halos below $\sim 10^{9} M_{\odot}$, which is equivalent to steepening the faint-end of the galaxy LF. These findings demonstrate that the global 21-cm signal is a powerful complement to current and future galaxy surveys and efforts to better understand the interstellar medium in high-$z$ galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Cosmic Reionization and Early Star-Forming Galaxies: A Joint Analysis of New Constraints from Planck and Hubble Space Telescope

We discuss new constraints on the epoch of cosmic reionization and test the assumption that most of the ionizing photons responsible arose from high redshift star-forming galaxies. Good progress has been made in charting the end of reionization through spectroscopic studies of z~6-8 QSOs, gamma-ray bursts and galaxies expected to host Lyman-alpha emission. However, the most stringent constraints on its duration have come from the integrated optical depth, tau, of Thomson scattering to the cosmic microwave background. Using the latest data on the abundance and luminosity distribution of distant galaxies from Hubble Space Telescope imaging, we simultaneously match the reduced value tau=0.066 +/- 0.012 recently reported by the Planck collaboration and the evolving neutrality of the intergalactic medium with a reionization history within 6 <~ z <~ 10, thereby reducing the requirement for a significant population of very high redshift (z>>10) galaxies. Our analysis strengthens the conclusion that star-forming galaxies dominated the reionization process and has important implications for upcoming 21cm experiments and searches for early galaxies with James Webb Space Telescope.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Large fluctuations in the hydrogen-ionizing background and mean free path following the epoch of reionization

Extremely large opaque troughs in the Lyman-alpha forest have been interpreted as a sign of an extended reionization process below z~6. Such features are impossible to reproduce with simple models of the intergalactic ionizing background that assume a uniform mean free path of ionizing photons. We build a self-consistent model of the ionizing background that includes fluctuations in the mean free path due to the varying strength of the ionizing background and large-scale density field. The dominant effect is the suppression of the ionizing background in large-scale voids due to "self-shielding" by an enhanced number of optically thick absorbers. Our model results in a distribution of 50 Mpc/h Lyman-alpha forest effective optical depths that significantly improves agreement with the observations at z~5.6. Extrapolation to z~5.4 and z~5.8 appears promising, but matching the mean background evolution requires evolution in the absorber population beyond the scope of the present model. We also demonstrate the need for extremely large volumes (>400 Mpc on a side) to accurately determine the incidence of rare large-scale features in the Lyman-alpha forest.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

PAPER-64 Constraints On Reionization II: The Temperature Of The z=8.4 Intergalactic Medium

We present constraints on both the kinetic temperature of the intergalactic medium (IGM) at z=8.4, and on models for heating the IGM at high-redshift with X-ray emission from the first collapsed objects. These constraints are derived using a semi-analytic method to explore the new measurements of the 21 cm power spectrum from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), which were presented in a companion paper, Ali et al. (2015). Twenty-one cm power spectra with amplitudes of hundreds of mK^2 can be generically produced if the kinetic temperature of the IGM is significantly below the temperature of the Cosmic Microwave Background (CMB); as such, the new results from PAPER place lower limits on the IGM temperature at z=8.4. Allowing for the unknown ionization state of the IGM, our measurements find the IGM temperature to be above ~5 K for neutral fractions between 10% and 85%, above ~7 K for neutral fractions between 15% and 80%, or above ~10 K for neutral fractions between 30% and 70%. We also calculate the heating of the IGM that would be provided by the observed high redshift galaxy population, and find that for most models, these galaxies are sufficient to bring the IGM temperature above our lower limits. However, there are significant ranges of parameter space that could produce a signal ruled out by the PAPER measurements; models with a steep drop-off in the star formation rate density at high redshifts or with relatively low values for the X-ray to star formation rate efficiency of high redshift galaxies are generally disfavored. The PAPER measurements are consistent with (but do not constrain) a hydrogen spin temperature above the CMB temperature, a situation which we find to be generally predicted if galaxies fainter than the current detection limits of optical/NIR surveys are included in calculations of X-ray heating.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The 21-cm Line as a Probe of Reionization

One of the most exciting probes of the early phases of structure formation and reionization is the spin-flip line of neutral hydrogen, with a rest wavelength of 21 cm. This chapter introduces the physics of this transition and the astrophysical parameters upon which it depends, including discussions of the radiation fields that permeate the intergalactic medium that fix the brightness of this transition. We describe the critical points in the evolution of the 21-cm background and focus on the sky-averaged brightness and the power spectrum as representative measurements. Finally, we include a discussion of observations and the challenges they face in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Flatness and Sudden Evolution of the Intergalactic Ionizing Background

The ionizing background of cosmic hydrogen is an important probe of the sources and absorbers of ionizing radiation in the post-reionization universe. Previous studies show that the ionization rate should be very sensitive to changes in the source population: as the emissivity rises, absorbers shrink in size, increasing the ionizing mean free path and, hence, the ionizing background. By contrast, observations of the ionizing background find a very flat evolution from z~2-5, before falling precipitously at z~6. We resolve this puzzling discrepancy by pointing out that, at z~2-5, optically thick absorbers are associated with the same collapsed halos that host ionizing sources. Thus, an increasing abundance of galaxies is compensated for by a corresponding increase in the absorber population, which moderates the instability in the ionizing background. However, by z~5-6, gas outside of halos dominates the absorption, the coupling between sources and absorbers is lost, and the ionizing background evolves rapidly. Our halo based model reproduces observations of the ionizing background, its flatness and sudden decline, as well as the redshift evolution of the ionizing mean free path. Our work suggests that, through much of their history, both star formation and photoelectric opacity in the universe track halo growth.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Molecular Cloud Properties and CO Line Emission in z >~ 6 Galaxies

We explore molecular cloud properties and the physics of CO transition lines in z >~ 6 Lyman-break galaxies and predict their CO fluxes using an analytic formalism built from global models of star formation in high-redshift galaxies that minimizes our reliance on local observations. Our model includes a new approach to calculating the molecular gas fraction that takes the total gas density, the star formation rate, and the star formation efficiency in clouds as the principal inputs. This method agrees with chemical equilibrium calculations of the molecular fraction based on local chemistry if galaxies at z >~ 6 have metallicities of order a few percent of solar. Such low metallicities in turn imply that much of the carbon in these systems exists in ionized form rather than as CO. Moreover, we find that the higher-order CO transitions observable at high redshift with ALMA will typically be sub-thermally populated but that the details depend sensitively on the presence of turbulent clumps within molecular clouds. Ultimately, we expect current facilities will only be able to observe the CO signal from reionization epoch galaxies with great difficulty. We estimate that at least ~100 hours of integration time with ALMA will be required to detect the CO(6--5) transition in z = 6 systems with rest-frame UV magnitudes of -20.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

New Constraints on Cosmic Reionization from the 2012 Hubble Ultra Deep Field Campaign

Understanding cosmic reionization requires the identification and characterization of early sources of hydrogen-ionizing photons. The 2012 Hubble Ultra Deep Field (UDF12) campaign has acquired the deepest infrared images with the Wide Field Camera 3 aboard Hubble Space Telescope and, for the first time, systematically explored the galaxy population deep into the era when cosmic microwave background (CMB) data indicates reionization was underway. The UDF12 campaign thus provides the best constraints to date on the abundance, luminosity distribution, and spectral properties of early star-forming galaxies. We synthesize the new UDF12 results with the most recent constraints from CMB observations to infer redshift-dependent ultraviolet (UV) luminosity densities, reionization histories, and electron scattering optical depth evolution consistent with the available data. Under reasonable assumptions about the escape fraction of hydrogen ionizing photons and the intergalactic medium clumping factor, we find that to fully reionize the universe by redshift z~6 the population of star-forming galaxies at redshifts z~7-9 likely must extend in luminosity below the UDF12 limits to absolute UV magnitudes of M_UV\sim -13 or fainter. Moreover, low levels of star formation extending to redshifts z~15-25, as suggested by the normal UV colors of z\simeq7-8 galaxies and the smooth decline in abundance with redshift observed by UDF12 to z\simeq10, are additionally likely required to reproduce the optical depth to electron scattering inferred from CMB observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Semi-numeric simulations of helium reionization and the fluctuating radiation background

Recent He II Lyman-alpha forest observations from 2.0 < z < 3.2 show large fluctuations in the optical depth at z > 2.7. These results point to a fluctuating He-ionizing background, which may be due to the end of helium reionization of this era. We present a fast, semi-numeric procedure to approximate detailed cosmological simulations. We compute the distribution of dark matter halos, ionization state of helium, and density field at z = 3 in broad agreement with recent simulations. Given our speed and flexibility, we investigate a range of ionizing source and active quasar prescriptions. Spanning a large area of parameter space, we find order-of-magnitude fluctuations in the He II ionization rate in the post-reionization regime. During reionization, the fluctuations are even stronger and develop a bimodal distribution, in contrast to semi-analytic models and the hydrogen equivalent. These distributions indicate a low-level ionizing background even at significant He II fractions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Effect of Fluctuations on the Helium-Ionizing Background

Interpretation of He II Lyα absorption spectra after the epoch of He II reionization requires knowledge of the He II ionizing background. While past work has modelled the evolution of the average background, the standard cosmological radiative transfer technique assumes a uniform radiation field despite the discrete nature of the (rare) bright quasars that dominate the background. We implement a cosmological radiative transfer model that includes the most recent constraints on the ionizing spectra and luminosity function of quasars and the distribution of IGM absorbers. We also estimate, for the first time, the effects of fluctuations on the evolving continuum opacity in two ways: by incorporating the complete distribution of ionizing background amplitudes into the standard approach, and by explicitly treating the quasars as discrete -- but isolated -- sources. Our model results in a He II ionization rate that evolves steeply with redshift, increasing by a factor ~2 from z=3.0 to z=2.5. This causes rapid evolution in the mean He II Lyα optical depth -- as recently observed -- without appealing to the reionization of He II. The observed behaviour could instead result from rapid evolution in the mean free path of ionizing photons as the helium in higher H I column density absorbers becomes fully ionized.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evolution of the Sizes of Galaxies over 7<z<12 Revealed by the 2012 Hubble Ultra Deep Field Campaign

We analyze the redshift- and luminosity-dependent sizes of dropout galaxy candidates in the redshift range z~7-12 using deep images from the UDF12 campaign, data which offers two distinct advantages over that used in earlier work. Firstly, we utilize the increased S/N ratio offered by the UDF12 imaging to provide improved size measurements for known galaxies at z=6.5-8 in the HUDF. Specifically, we stack the new deep F140W image with the existing F125W data in order to provide improved measurements of the half-light radii of z-dropouts. Similarly we stack this image with the new deep UDF12 F160W image to obtain new size measurements for a sample of Y-dropouts. Secondly, because the UDF12 data have allowed the construction of the first robust galaxy sample in the HUDF at z>8, we have been able to extend the measurement of average galaxy size out to significantly higher redshifts. Restricting our size measurements to sources which are now detected at >15sigma, we confirm earlier indications that the average half-light radii of z~7-12 galaxies are extremely small, 0.3-0.4 kpc, comparable to the sizes of giant molecular associations in local star-forming galaxies. We also confirm that there is a clear trend of decreasing half-light radius with increasing redshift, and provide the first evidence that this trend continues beyond z~8. Modeling the evolution of the average half-light radius as a power-law (1+z)^s, we obtain a best-fit index of s=-1.28+/-0.13 over the redshift range z~4-12, mid-way between the physically expected evolution for baryons embedded in dark halos of constant mass (s=-1) and constant velocity (s=-1.5). A clear size-luminosity relation, such as that found at lower redshift, is also evident in both our z- and Y-dropout sample. This relation can be interpreted in terms of a constant surface density of star formation over a range in luminosity of 0.05-1.0L*_z=3.(abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Faint AGN in z>~6 Lyman-break Galaxies Powered by Cold Accretion and Rapid Angular Momentum Transport

We develop a radiation pressure-balanced model for the interstellar medium of high-redshift galaxies that describes many facets of galaxy formation at z>~6, including star formation rates and distributions and gas accretion onto central black holes. We first show that the vertical gravitational force in the disk of such a model is dominated by the disk self-gravity supported by the radiation pressure of ionizing starlight on gas. Constraining our model to reproduce the UV luminosity function of Lyman-break galaxies (LBGs), we limit the available parameter-space to wind mass-loading factors 1--4 times the canonical value for momentum-driven winds. We then focus our study by exploring the effects of different angular momentum transport mechanisms in the galactic disk and find that accretion driven by gravitational torques, such as from linear spiral waves or non-linear orbit crossings, can build up black hole masses by z=6 consistent with the canonical M-sigma relation with a duty cycle of unity, while accretion mediated by a local viscosity such as in an alpha-disk results in negligible BH accretion. Both gravitational torque models produce X-ray emission from active galactic nuclei (AGN) in high-redshift LBGs in excess of the estimated contribution from high-mass X-ray binaries. Using a recent analysis of deep Chandra observations by Cowie et al., we can already begin to rule out the most extreme regions of our parameter-space: the inflow velocity of gas through the disk must either be less than one percent of the disk circular velocity or the X-ray luminosity of the AGN must be substantially obscured. Moderately deeper future observations or larger sample sizes will be able to probe the more reasonable range of angular momentum transport models and obscuring geometries.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The UV Luminosity Function of star-forming galaxies via dropout selection at redshifts z ~ 7 and 8 from the 2012 Ultra Deep Field campaign

We present a catalog of high redshift star-forming galaxies selected to lie within the redshift range z ~ 7-8 using the Ultra Deep Field 2012 (UDF12), the deepest near-infrared (near-IR) exposures yet taken with the Hubble Space Telescope. As a result of the increased near-infrared exposure time compared to previous HST imaging in this field, we probe 0.65 (0.25) mag fainter in absolute UV magnitude, at z ~ 7 (8), which increases confidence in a measurement of the faint end slope of the galaxy luminosity function. Through a 0.7 mag deeper limit in the key F105W filter that encompasses or lies just longward of the Lyman break, we also achieve a much-refined color-color selection that balances high redshift completeness and a low expected contamination fraction. We improve the number of drop-out selected UDF sources to 47 at z ~ 7 and 27 at z ~ 8. Incorporating brighter archival and ground-based samples, we measure the z ~ 7 UV luminosity function to an absolute magnitude limit of M_UV = -17 and find a faint end Schechter slope of α= -1.87+/- 0.18. Using a similar color-color selection at z ~ 8 that takes account of our newly-added imaging in the F140W filter, and incorporating archival data from the HIPPIES and BoRG campaigns, we provide a robust estimate of the faint end slope at z ~ 8, α= -1.94 +/- 0.23. We briefly discuss our results in the context of earlier work and that derived using the same UDF12 data but with an independent photometric redshift technique (McLure et al 2012).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fluctuations in the High-Redshift Lyman-Werner and Lyman-alpha Radiation Backgrounds

We use a new method to model fluctuations of the Lyman-Werner (LW) and Lyman-alpha radiation backgrounds at high redshift. At these early epochs the backgrounds are symptoms of a universe newly lit with its first stars. LW photons (11.5-13.6 eV) are of particular interest because they dissociate molecular hydrogen, the primary coolant in the first minihalos. By using a variation of the halo model, we efficiently generate power spectra for any choice of radiation background. We find that the LW power spectrum typically traces the matter power spectrum at large scales but turns over at the scale corresponding to the effective `horizon' of LW photons (~100 comoving Mpc), unless the sources are extremely rare. The series of horizons that characterize the Lyman-alpha flux profile shape the fluctuations of that background in a similar fashion, though those imprints are washed out once one considers fluctuations in the brightness temperature of the 21-cm signal. The Lyman-alpha background strongly affects the redshifted 21-cm signal at just about the time the LW background begins to dissociate molecular hydrogen, so measuring that background's properties will reveal important information about the transition from early Population III stars to more normal stars. Around this time we find that fluctuations in the LW background are weak; the fractional standard deviation is less than ~0.5 on scales > 10 cMpc, only rising to be of order unity on scales < 1 cMpc. This should not lead to substantial spatial fluctuations in molecular hydrogen content, except at the earliest times. Even then, most halos form far from other sources, so the transition from star formation in low-mass to high-mass halos is rather homogeneous across the universe.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Intensity Mapping with Carbon Monoxide Emission Lines and the Redshifted 21 cm Line

We quantify the prospects for using emission lines from rotational transitions of the CO molecule to perform an `intensity mapping' observation at high redshift during the Epoch of Reionization (EoR). The aim of CO intensity mapping is to observe the combined CO emission from many unresolved galaxies, to measure the spatial fluctuations in this emission, and use this as a tracer of large scale structure at very early times in the history of our Universe. This measurement would help determine the properties of molecular clouds -- the sites of star formation -- in the very galaxies that reionize the Universe. We further consider the possibility of cross-correlating CO intensity maps with future observations of the redshifted 21 cm line. The cross spectrum is less sensitive to foreground contamination than the auto power spectra, and can therefore help confirm the high redshift origin of each signal. Furthermore, the cross spectrum measurement would help extract key information about the EoR, especially regarding the size distribution of ionized regions. We discuss uncertainties in predicting the CO signal at high redshift, and discuss strategies for improving these predictions. Under favorable assumptions, and feasible specifications for a CO survey mapping the CO(2-1) and CO(1-0) lines, the power spectrum of CO emission fluctuations and its cross spectrum with future 21 cm measurements from the MWA are detectable at high significance.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulations and Analytic Calculations of Bubble Growth During Hydrogen Reionization

We present results from a large volume simulation of Hydrogen reionization. We combine 3d radiative transfer calculations and an N-body simulation, describing structure formation in the intergalactic medium, to detail the growth of HII regions around high redshift galaxies. Our simulation tracks 1024^3 dark matter particles, in a box of co-moving side length 65.6 Mpc/h. This large volume allows us to accurately characterize the size distribution of HII regions throughout most of the reionization process. At the same time, our simulation resolves many of the small galaxies likely responsible for reionization. It confirms a picture anticipated by analytic models: HII regions grow collectively around highly-clustered sources, and have a well-defined characteristic size, which evolves from a sub-Mpc scale at the beginning of reionization to R > 10 Mpc towards the end. We present a detailed statistical description of our results, and compare them with a numerical scheme based on the analytic model by Furlanetto, Zaldarriaga, and Hernquist. We find that the analytic calculation reproduces the size distribution of HII regions and the 21 cm power spectrum of the radiative transfer simulation remarkably well. The ionization field from the simulation, however, has more small scale structure than the analytic calculation, owing to Poisson scatter in the simulated abundance of galaxies on small scales. We propose and validate a simple scheme to incorporate this scatter into our calculations. Our results suggest that analytic calculations are sufficiently accurate to aid in predicting and interpreting the results of future 21 cm surveys. In particular, our fast numerical scheme is useful for forecasting constraints from future 21 cm surveys, and in constructing mock surveys to test data analysis procedures.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The evolution of the Lyman-alpha forest effective optical depth following HeII reionisation

Three independent observational studies have now detected a narrow (Δz ~ 0.5) dip centred at z=3.2 in the otherwise smooth redshift evolution of the Lya forest effective optical depth. This feature has previously been interpreted as an indirect signature of rapid photo-heating in the IGM during the epoch of HeII reionisation. We examine this interpretation using a semi-analytic model of inhomogeneous HeII reionisation and high resolution hydrodynamical simulations of the Lya forest. We instead find that a rapid (Δz ~ 0.2) boost to the IGM temperature (ΔT ~ 10^4 K) beginning at z=3.4 produces a well understood and generic evolution in the Lya effective optical depth, where a sudden reduction in the opacity is followed by a gradual, monotonic recovery driven largely by adiabatic cooling in the low density IGM. This behaviour is inconsistent with the narrow feature in the observational data. If photo-heating during HeII reionisation is instead extended over several redshift units, as recent theoretical studies suggest, then the Lya opacity will evolve smoothly with redshift. We conclude that the sharp dip observed in the Lya forest effective optical depth is instead most likely due to a narrow peak in the hydrogen photo-ionisation rate around z=3.2, and suggest that it may arise from the modulation of either reprocessed radiation during HeII reionisation, or the opacity of Lyman limit systems.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Cluster Merger Variance and the Luminosity Gap Statistic

The presence of multiple luminous galaxies in clusters can be explained by the finite time over which a galaxy sinks to the center of the cluster and merges with the the central galaxy. The simplest measurable statistic to quantify the dynamical age of a system of galaxies is the luminosity (magnitude) gap, which is the difference in photometric magnitude between the two most luminous galaxies. We present a simple analytical estimate of the luminosity gap distribution in groups and clusters as a function of dark matter halo mass. The luminosity gap is used to define "fossil" groups; we expect the fraction of fossil systems to exhibit a strong and model-independent trend with mass: ~1-3% of massive clusters and ~5-40% of groups should be fossil systems. We also show that, on cluster scales, the observed intrinsic scatter in the central galaxy luminosity-halo mass relation can be ascribed to dispersion in the merger histories of satellites within the cluster. We compare our predictions to the luminosity gap distribution in a sample of 730 clusters in the Sloan Digital Sky Survey C4 Catalog and find good agreement. This suggests that theoretical excursion set merger probabilities and the standard theory of dynamical segregation are valid on cluster scales.

0 citations0 reviewsNo code linkedOpen access