Source author record

Steven Furlanetto

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

9works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

9 published item(s)

preprint2010arXiv

21cmFAST: A Fast, Semi-Numerical Simulation of the High-Redshift 21-cm Signal

We introduce a powerful semi-numeric modeling tool, 21cmFAST, designed to efficiently simulate the cosmological 21-cm signal. Our code generates 3D realizations of evolved density, ionization, peculiar velocity, and spin temperature fields, which it then combines to compute the 21-cm brightness temperature. Although the physical processes are treated with approximate methods, we compare our results to a state-of-the-art large-scale hydrodynamic simulation, and find good agreement on scales pertinent to the upcoming observations (>~ 1 Mpc). The power spectra from 21cmFAST agree with those generated from the numerical simulation to within 10s of percent, down to the Nyquist frequency. We show results from a 1 Gpc simulation which tracks the cosmic 21-cm signal down from z=250, highlighting the various interesting epochs. Depending on the desired resolution, 21cmFAST can compute a redshift realization on a single processor in just a few minutes. Our code is fast, efficient, customizable and publicly available, making it a useful tool for 21-cm parameter studies.

preprint2010arXiv

Constraints on Quasar Lifetimes and Beaming from the HeII Lyman-alpha Forest

We show that comparisons of HeII Lyman-alpha forest lines of sight to nearby quasar populations can strongly constrain the lifetimes and emission geometry of quasars. By comparing the HeII and HI Lyman-alpha forests along a particular line of sight, one can trace fluctuations in the hardness of the radiation field (which are driven by fluctuations in the HeII ionization rate). Because this high-energy background is highly variable - thanks to the rarity of the bright quasars that dominate it and the relatively short attenuation lengths of these photons - it is straightforward to associate features in the radiation field with their source quasars. Here we quantify how finite lifetimes and beamed emission geometries affect these expectations. Finite lifetimes induce a time delay that displaces the observed radiation peak relative to the quasar. For beamed emission, geometry dictates that sources invisible to the observer can still create a peak in the radiation field. We show that both these models produce substantial populations of "bare" peaks (without an associated quasar) for reasonable parameter values (lifetimes ~10^6-10^8 yr and beaming angles <90 degrees). A comparison to existing quasar surveys along two HeII Lyman-alpha forest lines of sight rules out isotropic emission and infinite lifetime at high confidence; they can be accommodated either by moderate beaming or lifetimes ~10^7-10^8 yr. We also show that the distribution of radial displacements between peaks and their quasars can unambiguously distinguish these two models, although larger statistical samples are needed.

preprint2010arXiv

The distribution of Lyman-limit absorption systems during and after reionization

Lyman-limit absorption systems can play many important roles during and after cosmological reionization. Unfortunately, due to the prohibitively large dynamic range required, it is impossible to self-consistently include these systems in cosmological simulations. Using fast and versatile semi-numeric simulations, we systematically explore the spatial distribution of absorption systems during and following reionization. We self-calibrate the resulting number of absorbers to the mean free path (mfp) of the ionizing ultraviolet background (UVB), and present results at a given mfp and neutral hydrogen fraction. We use a simple optical depth criterion to identify the locations of absorbers. Our approach is fairly robust to uncertainties such as missing subgrid structure. Unlike at lower redshifts where the UVB is relatively uniform, at higher redshifts the fluctuations in the UVB and the HII morphology of reionization can drive the large-scale distribution of absorption systems. Specifically, we find that absorbers are highly correlated with the density field on small scales, and then become anti-correlated with the UVB on large scales. After reionization, the large-scale power spectrum of the absorbers traces the UVB power spectrum, which can be predicted with a simple analytic extension of the halo-model. During reionization, absorbers tend to preferentially lie inside overdensities (i.e. filaments) of the recently-ionized intergalactic medium (IGM). Absorbers may also dominate the small-scale (k > 1/Mpc) 21-cm power during and after reionization. Conversely, they smooth the contrast on moderate scales. Once the HII regions grow to surpass the mfp, the absorbers add to the large-scale 21-cm power. Our results should prove useful in interpreting future observations of the reionization epoch.

preprint2009arXiv

Large-Scale Fluctuations in the HeII Lyman-alpha Forest and HeII Reionization

We examine large-scale fluctuations in the HeII Lyman-alpha forest transmission during and after HeII reionization. We use a simple Monte Carlo model to distribute quasars throughout a large volume and compute the resulting radiation field along one-dimensional skewers. In agreement with previous studies, we find that the rarity of these sources induces order unity fluctuations in the mean optical depth after reionization, even when averaged over large segments (~10-100 Mpc across). We compare our models to existing data along five HeII Lyman-alpha forest lines of sight spanning z~2-3.2. The large cosmic variance contained in our model plausibly explains many of the observed fluctuations at z<2.5. But our models cannot accommodate the large fluctuations toward high optical depths on ~30 Mpc scales observed at z~2.7-2.9, and the measured optical depths (tau>4) at z>2.9 are difficult to explain with a smoothly-evolving mean radiation field. In order to better understand this data, we construct a toy model of HeII reionization, in which we assume that regions with the smallest radiation fields in a post-reionization Universe (or farthest from strong ionizing sources) are completely dark during reionization. The observed fluctuations fit much more comfortably into this model, and we therefore argue that, according to present data, HeII reionization does not complete until z<2.9.

preprint2009arXiv

Secondary ionization and heating by fast electrons

We examine the fate of fast electrons (with energies E>10 eV) in a thermal gas of primordial composition. To follow their interactions with the background gas, we construct a Monte Carlo model that includes: (1) electron-electron scattering (which transforms the electron kinetic energy into heat), (2) collisional ionization of hydrogen and helium (which produces secondary electrons that themselves scatter through the medium), and (3) collisional excitation (which produces secondary photons, whose fates we also follow approximately). For the last process, we explicitly include all transitions to upper levels n<=4, together with a well-motivated extrapolation to higher levels. In all cases, we use recent calculated cross-sections at E<1 keV and the Bethe approximation to extrapolate to higher energies. We compute the fractions of energy deposited as heat, ionization (tracking HI and the helium species separately), and excitation (tracking HI Lyman-alpha separately) under a broad range of conditions appropriate to the intergalactic medium. The energy deposition fractions depend on both the background ionized fraction and the electron energy but are nearly independent of the background density. We find good agreement with some, but not all, previous calculations at high energies. Electronic tables of our results are available on request.

preprint2009arXiv

The Equation of State of the Intergalactic Medium After Hydrogen Reionization

We use an analytic model to study how inhomogeneous hydrogen reionization affects the temperature distribution of the intergalactic medium (IGM). During this process, the residual energy of each ionizing photon is deposited in the IGM as heat, increasing its temperature to 20,000-30,000 K; subsequent expansion of the Universe then cools the gas. Because reionization most likely proceeds from high to low densities, underdense voids are ionized last, have less time to cool, and are (on average) warmer than mean-density gas immediately after reionization is complete (an "inverted" equation of state). From this initial configuration, the low-density gas cools quickly and eventually returns to a more normal equation of state. The rapidly evolving temperature introduces systematic uncertainties in measurements of the ionizing background at z~6. For example, late reionization implies rapid cooling, so that the ionizing background would have to evolve even more rapidly at z ~5-6 than typically claimed. This degeneracy is difficult to disentangle, because the Lyman-alpha forest probes only a narrow range in densities (over which the gas is nearly isothermal). However, higher Lyman-series transitions probe wider density ranges, sampling different effective temperatures, and offer a new way to measure the IGM equation of state that should work where nearly saturated absorption precludes other methods. This will help to separate evolution in temperature from that in the ionizing background. While more detailed study with hydrodynamic simulations is needed, we show that such measurements could potentially distinguish early and late reionization using only a handful of lines of sight.

preprint2009arXiv

The Inhomogeneous Ionizing Background Following Reionization

We study the spatial fluctuations in the hydrogen ionizing background in the epoch following reionization (z ~ 5--6). The rapid decrease with redshift in the photon mean free path (m.f.p.), combined with the clustering of increasingly rare ionizing sources, can result in a very inhomogenous ionizing background during this epoch. We systematically investigate the probability density functions (PDFs) and power spectra of ionizing flux, by varying several parameters such as the m.f.p., minimum halo mass capable of hosting stars, and halo duty cycle. In order to be versatile, we make use of analytic, semi-numeric and numeric approaches. Our models show that the ionizing background indeed has sizable fluctuations during this epoch sourced by the clustering of sources, with the PDFs being a factor of few wide at half of the maximum likelihood. The distributions also show marked asymmetries, with a high-value tail set by clustering on small scales, and a shorter low-value tail which is set by the mean free path. The power spectrum of the ionizing background is much more sensitive to source properties than the PDF and can be well-understood analytically with a framework similar to the halo model (usually used to describe dark matter clustering). Nevertheless, we find that Lya forest spectra are extremely insensitive to the details of the UVB, despite marked differences in the PDFs and power spectra of our various ionizing backgrounds. Assuming a uniform ionizing background only underestimates the value of the mean ionization rate inferred from the Lya forest by a few percent. Instead, analysis of the Lya forest is dominated by the uncertainties in the density field. Thus, our results justify the common assumption of a uniform ionizing background in Lya forest analysis even during this epoch.

preprint2007arXiv

Reionization and the large-scale 21 cm-cosmic microwave background cross correlation

Of the many probes of reionization, the 21 cm line and the cosmic microwave background (CMB) are among the most effective. We examine how the cross-correlation of the 21 cm brightness and the CMB Doppler fluctuations on large angular scales can be used to study this epoch. We employ a new model of the growth of large scale fluctuations of the ionized fraction as reionization proceeds. We take into account the peculiar velocity field of baryons and show that its effect on the cross correlation can be interpreted as a mixing of Fourier modes. We find that the cross-correlation signal is strongly peaked toward the end of reionization and that the sign of the correlation should be positive because of the inhomogeneity inherent to reionization. The signal peaks at degree scales (l~100) and comes almost entirely from large physical scales (k~0.01 Mpc). Since many of the foregrounds and noise that plague low frequency radio observations will not correlate with CMB measurements, the cross correlation might appear to provide a robust diagnostic of the cosmological origin of the 21 cm radiation around the epoch of reionization. Unfortunately, we show that these signals are actually only weakly correlated and that cosmic variance dominates the error budget of any attempted detection. We conclude that the detection of a cross-correlation peak at degree-size angular scales is unlikely even with ideal experiments.

preprint2004arXiv

Polarization Signals of the 21 cm Background from the Era of Reionization

While emission and absorption lines of the 21 cm spin-flip transition of neutral Hydrogen are intrinsically unpolarized, a magnetic field creates left- and right-handed polarized components through the Zeeman effect. Here we consider the resulting polarization of the redshifted 21 cm background from the intergalactic medium before reionization. The polarization is detectable in regions with a strong gradient in the mean brightness temperature. In principle, this can open a new window on the evolution of intergalactic magnetic fields. One possible approach is an extended integration of an individual target during this era, such as the Mpc-scale HII regions inferred to surround quasars at z ~ 6.5. The differential intensity between the two polarization states can be used as a probe of the magnetic field at the edge of the HII region. We estimate that the SKA could (ignoring systematics) detect B \~ 200 (10) \muG coherent over several kiloparsecs with an observational bandwidth of 100 (2) kHz. Beyond individual sources, the statistical properties of wide-field 21 cm polarization maps, such as the angular power spectrum, can be used to constrain the large-scale magnetic field. In this case, the SKA can detect $B \sim 100 μ$G fields coherent over many megaparsecs. The magnetic field can be measured in any epoch over which the 21 cm background changes rapidly (for example because the ionized fraction or spin temperature change). Although the resulting constraints with SKA are relatively weak compared to theoretical expectations, they nevertheless offer a unique direct probe of magnetic fields in the high-redshift universe.