Source author record

Garrelt Mellema

Garrelt Mellema 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

48works
7topics
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

48 published item(s)

preprint2022arXiv

Astrophysical information from the Rayleigh-Jeans Tail of the CMB

One of the explanations for the recent EDGES-LOW band 21-cm measurements of a strong absorption signal around 80~MHz is the presence of an excess radio background to the Cosmic Microwave Background (CMB). Such excess can be produced by the decay of unstable particles into small mass dark photons which have a non-zero mixing angle with electromagnetism. We use the EDGES-LOW band measurements to derive joint constraints on the properties of the early galaxies and the parameters of such a particle physics model for the excess radio background. A Bayesian analysis shows that a high star formation efficiency and X-ray emission of $4-7 \times 10^{48} ~\rm erg$ per solar mass in stars are required along with a suppression of star formation in halos with virial temperatures $\lesssim 2\times 10^4$ K. The same analysis also suggests a 68 percent credible intervals for the mass of the decaying dark matter particles, it's lifetime, dark photon mass and the mixing angle of the dark and ordinary photon oscillation of $[10^{-3.5}, 10^{-2.4}]$ eV, $[10^{1.1}, 10^{2.7}]\times 13.8 ~\rm Gyr$, $[10^{-12.2}, 10^{-10}]$ eV and $[10^{-7}, 10^{-5.6}]$ respectively. This implies an excess radio background which is $\approx 5.7$ times stronger than the CMB around 80~MHz. This value is a factor $\sim 3$ higher than the previous predictions which used a simplified model for the 21-cm signal.

preprint2022arXiv

The large-scale 21-cm power spectrum from reionization

Radio interferometers, such as the Low-Frequency Array and the future Square Kilometre Array, are attempting to measure the spherically averaged 21-cm power spectrum from the Epoch of Reionization. Understanding of the dominant physical processes which influence the power spectrum at each length-scale is therefore crucial for interpreting any future detection. We study a decomposition of the 21-cm power spectrum and quantify the evolution of its constituent terms for a set of numerical and semi-numerical simulations of a volume of $(714~\mathrm{Mpc})^3$, focusing on large scales with $k\lesssim 0.3$~Mpc$^{-1}$. We find that after $\sim 10$ per cent of the Universe has been ionized, the 21-cm power spectrum follows the power spectrum of neutral hydrogen fluctuations, which itself beyond a certain scale follows the matter power spectrum. Hence the signal has a two-regime form where the large-scale signal is a biased version of the cosmological density field, and the small-scale power spectrum is determined by the astrophysics of reionization. We construct a bias parameter to investigate the relation between the large-scale 21-cm signal and the cosmological density field. We find that the transition scale between the scale-independent and scale-dependent bias regimes is directly related to the value of the mean free path of ionizing photons ($λ_{\mathrm{MFP}}$), and is characterised by the empirical formula $k_{\mathrm{trans}} \approx 2/λ_{\mathrm{MFP}}$. Furthermore, we show that the numerical implementation of the mean free path effect has a significant impact on the shape of this transition. Most notably, the transition is more gradual if the mean free path effect is implemented as an absorption process rather than as a barrier.

preprint2022arXiv

The multi-frequency angular power spectrum in parameter studies of the cosmic 21-cm signal

The light-cone effect breaks the periodicity and statistical homogeneity (ergodicity) along the line-of-sight direction of cosmological emission/absorption line surveys. The spherically averaged power spectrum (SAPS), which by definition assumes ergodicity and periodicity in all directions, can only quantify some of the second-order statistical information in the 3D light-cone signals and therefore gives a biased estimate of the true statistics. The multi-frequency angular power spectrum (MAPS), by extracting more information from the data, does not rely on these assumptions. It is therefore better aligned with the properties of the signal. We have compared the performance of the MAPS and SAPS metrics for parameter estimation for a mock 3D light-cone observation of the 21-cm signal from the Epoch of Reionization. Our investigation is based on a simplified 3-parameter 21cmFAST model. We find that the MAPS produces parameter constraints which are a factor of $\sim 2$ more stringent than when the SAPS is used. The significance of this result does not change much even in the presence of instrumental noise expected for 128 hours of SKA-Low observations. Our results therefore suggest that a parameter estimation framework based on the MAPS metric would yield superior results over one using the SAPS metric.

preprint2021arXiv

Redshifted 21-cm bispectrum II: Impact of the spin temperature fluctuations and redshift space distortions on the signal from the Cosmic Dawn

We present a study of the 21-cm signal bispectrum (which quantifies the non-Gaussianity in the signal) from the Cosmic Dawn (CD). For our analysis, we have simulated the 21-cm signal using radiative transfer code GRIZZLY, while considering two types of sources (mini-QSOs and HMXBs) for Ly$α$ coupling and the X-ray heating of the IGM. Using this simulated signal, we have, for the first time, estimated the CD 21-cm bispectra for all unique $k$-triangles and for a range of $k$ modes. We observe that the redshift evolution of the bispectra magnitude and sign follow a generic trend for both source models. However, the redshifts at which the bispectra magnitude reach their maximum and minimum values and show their sign reversal depends on the source model. When the Ly$α$ coupling and the X-ray heating of the IGM occur simultaneously, we observe two consecutive sign reversals in the bispectra for small $k$-triangles (irrespective of the source models). One arising at the beginning of the IGM heating and the other at the end of Ly$α$ coupling saturation. This feature can be used in principle to constrain the CD history and/or to identify the specific CD scenarios. We also quantify the impact of the spin temperature ($T_{\rm S}$) fluctuations on the bispectra. We find that $T_{\rm S}$ fluctuations have maximum impact on the bispectra magnitude for small $k$-triangles and at the stage when Ly$α$ coupling reaches saturation. Furthermore, we are also the first to quantify the impact of redshift space distortions (RSD), on the CD bispectra. We find that the impact of RSD on the CD 21-cm bispectra is significant ($> 20\%$) and the level depends on the stages of the CD and the $k$-triangles for which the bispectra are being estimated.

preprint2020arXiv

Bubble mapping with the Square Kilometer Array -- I. Detecting galaxies with Euclid, JWST, WFIRST and ELT within ionized bubbles in the intergalactic medium at z>6

The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts.

preprint2020arXiv

Predictions for the 21cm-galaxy cross-power spectrum observable with SKA and future galaxy surveys

In this paper we use radiative transfer + N-body simulations to explore the feasibility of measurements of cross-correlations between the 21cm field observed by the Square Kilometer Array (SKA) and high-z Lyman Alpha Emitters (LAEs) detected in galaxy surveys with the Subaru Hyper Supreme Cam (HSC), Subaru Prime Focus Spectrograph (PFS) and Wide Field Infrared Survey Telescope (WFIRST). 21cm-LAE cross-correlations are in fact a powerful probe of the epoch of reionization as they are expected to provide precious information on the progress of reionization and the typical scale of ionized regions at different redshifts. The next generation observations with SKA will have a noise level much lower than those with its precursor radio facilities, introducing a significant improvement in the measurement of the cross-correlations. We find that an SKA-HSC/PFS observation will allow to investigate scales below ~10 Mpc/h and ~60 Mpc/h at z=7.3 and 6.6, respectively. WFIRST will allow to access also higher redshifts, as it is expected to observe spectroscopically ~900 LAEs per square degree and unit redshift in the range 7.5<z<8.5. Because of the reduction of the shot noise compared to HSC and PFS, observations with WFIRST will result in more precise cross-correlations and increased observable scales.

preprint2019arXiv

Impact of Ly$α$ heating on the global 21-cm signal from the Cosmic Dawn

The resonance scattering of Ly$α$ photons with neutral hydrogen atoms in the intergalactic medium not only couples the spin temperature to the kinetic temperature but also leads to a heating of the gas. We investigate the impact of this heating on the average brightness temperature of the 21-cm signal from the Cosmic Dawn in the context of the claimed detection by the EDGES low-band experiment. We model the evolution of the global signal taking into account the Ly$α$ coupling and heating and a cooling which can be stronger than the Hubble cooling. Using the claimed detection of a strong absorption signal at $z\approx 17$ as a constraint, we find that a strong Ly$α$ background is ruled out. Instead the results favour a weak Ly$α$ background combined with an excess cooling mechanism which is substantially stronger than previously considered. We also show that the cooling mechanism driven by the interaction between millicharged baryons and dark matter particles no longer provides a viable explanation for the EDGES result when Ly$α$ heating is taken into account.

preprint2019arXiv

The impact of inhomogeneous subgrid clumping on cosmic reionization

Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, ${\left<n^2\right>/\left<n\right>^2}$, which corrects for unresolved fluctuations in gas density ${n}$ on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution $N$-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially-varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale $N$-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally-averaged, time-evolving clumping factor from the same high-resolution $N$-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20% error, suggesting it may be a useful approximation.

preprint2016arXiv

The Effects of Lyman-Limit Systems on the Evolution and Observability of the Epoch of Reionization

We present the first large-scale, full radiative transfer simulations of the reionization of the intergalactic medium in the presence of Lyman-limit systems (LLSs). To illustrate the impact of LLS opacity, possibly missed by previous simulations, we add either a uniform or spatially-varying hydrogen bound-free opacity. This opacity, implemented as the mean free path (mfp) of the ionizing photons, extrapolates the observed, post-reionization redshift dependence into the epoch of reionization. In qualitative agreement with previous studies, we find that at late times the presence of LLSs slows down the ionization fronts, and alters the size distribution of H II regions. We quantitatively characterize the size distribution and morphological evolution of H II regions and examine the effects of the LLSs on the redshifted 21-cm signal from the patchy reionization. The presence of LLSs extends the ionization history by $Δz \sim 0.8$. The LLS absorbers significantly impede the late-time growth of the H II regions. The position dependent LLS distribution slows reionization further and additionally limits the late growth of the ionized regions. However, there is no "freeze out" of the H II regions and the largest regions grow to the size of the simulation volume. The 21-cm power spectra show that at large scales the power drops by a factor of 2 for 50% and 75% ionization stages (at $k = 0.1$ $\text{h} \, \text{Mpc}^{-1} $) reflecting the limiting effect of the LLSs on the growth of ionized patches. The statistical observables such as the RMS of the brightness temperature fluctuations and the peak amplitudes of the 21-cm power spectra at large-scales ($k = 0.05 - 0.1$ $\text{h} \, \text{Mpc}^{-1} $) are diminished by the presence of LLS.

preprint2015arXiv

Effects of the sources of reionization on 21-cm redshift-space distortions

The observed 21-cm signal from the epoch of reionization will be distorted along the line-of-sight by the peculiar velocities of matter particles. These redshift-space distortions will affect the contrast in the signal and will also make it anisotropic. This anisotropy contains information about the cross-correlation between the matter density field and the neutral hydrogen field, and could thus potentially be used to extract information about the sources of reionization. In this paper, we study a collection of simulated reionization scenarios assuming different models for the sources of reionization. We show that the 21-cm anisotropy is best measured by the quadrupole moment of the power spectrum. We find that, unless the properties of the reionization sources are extreme in some way, the quadrupole moment evolves very predictably as a function of global neutral fraction. This predictability implies that redshift-space distortions are not a very sensitive tool for distinguishing between reionization sources. However, the quadrupole moment can be used as a model-independent probe for constraining the reionization history. We show that such measurements can be done to some extent by first-generation instruments such as LOFAR, while the SKA should be able to measure the reionization history using the quadrupole moment of the power spectrum to great accuracy.

preprint2015arXiv

Efficient Photo-heating Algorithms in Time-dependent Photo-ionization Simulations

We present an extension to the time-dependent photo-ionization code C$^2$-Ray to calculate photo-heating in an efficient and accurate way. In C$^2$-Ray, the thermal calculation demands relatively small time-steps for accurate results. We describe two novel methods to reduce the computational cost associated with small time-steps, namely, an adaptive time-step algorithm and an asynchronous evolution approach. The adaptive time-step algorithm determines an optimal time-step for the next computational step. It uses a fast ray-tracing scheme to quickly locate the relevant cells for this determination and only use these cells for the calculation of the time-step. Asynchronous evolution allows different cells to evolve with different time-steps. The asynchronized clocks of the cells are synchronized at the times where outputs are produced. By only evolving cells which may require short time-steps with these short time-steps instead of imposing them to the whole grid, the computational cost of the calculation can be substantially reduced. We show that our methods work well for several cosmologically relevant test problems and validate our results by comparing to the results of another time-dependent photo-ionization code.

preprint2015arXiv

Epoch of Reionization modelling and simulations for SKA

In this chapter we provide an overview of the current status of the simulations and modelling of the Cosmic Dawn and Epoch of Reionization. We discuss the modelling requirements as dictated by the characteristic scales of the problem and the SKA instrumental properties and the planned survey parameters. Current simulations include most of the relevant physical processes. They can follow the full nonlinear dynamics and are now reaching the required scale and dynamic range, although small-scale physics still needs to be included at sub-grid level. However, despite a significant progress in developing novel numerical methods for efficient utilization of current hardware they remain quite computationally expensive. In response, a number of alternative approaches, particularly semi-analytical/semi-numerical methods, have been developed. While necessarily more approximate, if appropriately constructed and calibrated on simulations they could be used to quickly explore the vast parameter space available. Further work is still required on including some physical processes in both simulations and semi-analytical modelling. This hybrid approach of fast, approximate modelling calibrated on numerical simulations can then be used to construct large libraries of reionization models for reliable interpretation of the observational data.

preprint2015arXiv

Predictions for the 21cm-galaxy cross-power spectrum observable with LOFAR and Subaru

The 21cm-galaxy cross-power spectrum is expected to be one of the promising probes of the Epoch of Reionization (EoR), as it could offer information about the progress of reionization and the typical scale of ionized regions at different redshifts. With upcoming observations of 21cm emission from the EoR with the Low Frequency Array (LOFAR), and of high redshift Lyalpha emitters (LAEs) with Subaru's Hyper Suprime Cam (HSC), we investigate the observability of such cross-power spectrum with these two instruments, which are both planning to observe the ELAIS-N1 field at z=6.6. In this paper we use N-body + radiative transfer (both for continuum and Lyalpha photons) simulations at redshift 6.68, 7.06 and 7.3 to compute the 3D theoretical 21cm-galaxy cross-power spectrum, as well as to predict the 2D 21cm-galaxy cross-power spectrum expected to be observed by LOFAR and HSC. Once noise and projection effects are accounted for, our predictions of the 21cm-galaxy cross-power spectrum show clear anti-correlation on scales larger than ~ 60 h$^{-1}$ Mpc (corresponding to k ~ 0.1 h Mpc$^{-1}$), with levels of significance p=0.04 at z=6.6 and p=0.048 at z=7.3. On smaller scales, instead, the signal is completely contaminated.

preprint2015arXiv

The Large-Scale Observational Signatures of Low-Mass Galaxies During Reionization

Observations of the epoch of reionization give us clues about the nature and evolution of the sources of ionizing photons, or early stars and galaxies. We present a new suite of structure formation and radiative transfer simulations from the PRACE4LOFAR project designed to investigate whether the mechanism of radiative feedback, or the suppression of star formation in ionized regions from UV radiation, can be inferred from these observations. Our source halo mass extends down to $10^8 M_\odot$, with sources in the mass range $10^8$ to $10^9 M_\odot$ expected to be particularly susceptible to feedback from ionizing radiation, and we vary the aggressiveness and nature of this suppression. Not only do we have four distinct source models, we also include two box sizes (67 Mpc and 349 Mpc), each with two grid resolutions. This suite of simulations allows us to investigate the robustness of our results. All of our simulations are broadly consistent with the observed electron-scattering optical depth of the cosmic microwave background and the neutral fraction and photoionization rate of hydrogen at $z\sim6$. In particular, we investigate the redshifted 21-cm emission in anticipation of upcoming radio interferometer observations. We find that the overall shape of the 21-cm signal and various statistics are robust to the exact nature of source suppression, the box size, and the resolution. There are some promising model discriminators in the non-Gaussianity and small-scale power spectrum of the 21-cm signal.

preprint2014arXiv

Constraining the epoch of reionization with the variance statistic: simulations of the LOFAR case

Several experiments are underway to detect the cosmic redshifted 21-cm signal from neutral hydrogen from the Epoch of Reionization (EoR). Due to their very low signal-to-noise ratio, these observations aim for a statistical detection of the signal by measuring its power spectrum. We investigate the extraction of the variance of the signal as a first step towards detecting and constraining the global history of the EoR. Signal variance is the integral of the signal's power spectrum, and it is expected to be measured with a high significance. We demonstrate this through results from a simulation and parameter estimation pipeline developed for the Low Frequency Array (LOFAR)-EoR experiment. We show that LOFAR should be able to detect the EoR in 600 hours of integration using the variance statistic. Additionally, the redshift ($z_r$) and duration ($Δz$) of reionization can be constrained assuming a parametrization. We use an EoR simulation of $z_r = 7.68$ and $Δz = 0.43$ to test the pipeline. We are able to detect the simulated signal with a significance of 4 standard deviations and extract the EoR parameters as $z_r = 7.72^{+0.37}_{-0.18}$ and $Δz = 0.53^{+0.12}_{-0.23}$ in 600 hours, assuming that systematic errors can be adequately controlled. We further show that the significance of detection and constraints on EoR parameters can be improved by measuring the cross-variance of the signal by cross-correlating consecutive redshift bins.

preprint2014arXiv

Light cone effect on the reionization 21-cm signal II: Evolution, anisotropies and observational implications

Measurements of the HI 21-cm power spectra from the reionization epoch will be influenced by the evolution of the signal along the line-of-sight direction of any observed volume. We use numerical as well as semi-numerical simulations of reionization in a cubic volume of 607 Mpc across to study this so-called light cone effect on the HI 21-cm power spectrum. We find that the light cone effect has the largest impact at two different stages of reionization: one when reionization is $\sim 20\%$ and other when it is $\sim 80\%$ completed. We find a factor of $\sim 4$ amplification of the power spectrum at the largest scale available in our simulations. We do not find any significant anisotropy in the 21-cm power spectrum due to the light cone effect. We argue that for the power spectrum to become anisotropic, the light cone effect would have to make the ionized bubbles significantly elongated or compressed along the line-of-sight, which would require extreme reionization scenarios. We also calculate the two-point correlation functions parallel and perpendicular to the line-of-sight and find them to differ. Finally, we calculate an optimum frequency bandwidth below which the light cone effect can be neglected when extracting power spectra from observations. We find that if one is willing to accept a $10 \%$ error due to the light cone effect, the optimum frequency bandwidth for $k= 0.056 \, \rm{Mpc}^{-1}$ is $\sim 7.5$ MHz. For $k = 0.15$ and $0.41 \, \rm{Mpc}^{-1}$ the optimum bandwidth is $\sim 11$ and $\sim 16$ MHz respectively.

preprint2014arXiv

On the use of semi-numerical simulations in predicting the 21-cm signal from the epoch of reionization

We present a detailed comparison of three different simulations of the epoch of reionization (EoR). The radiative transfer simulation (${\rm C}^2$-RAY) among them is our benchmark. Radiative transfer codes can produce realistic results, but are computationally expensive. We compare it with two semi-numerical techniques: one using the same halos as ${\rm C}^2$-RAY as its sources (Sem-Num), and one using a conditional Press-Schechter scheme (CPS+GS). These are vastly more computationally efficient than ${\rm C}^2$-RAY, but use more simplistic physical assumptions. We evaluate these simulations in terms of their ability to reproduce the history and morphology of reionization. We find that both Sem-Num and CPS+GS can produce an ionization history and morphology that is very close to ${\rm C}^2$-RAY, with Sem-Num performing slightly better compared to CPS+GS. We also study different redshift space observables of the 21-cm signal from EoR: the variance, power spectrum and its various angular multipole moments. We find that both semi-numerical models perform reasonably well in predicting these observables at length scales relevant for present and future experiments. However, Sem-Num performs slightly better than CPS+GS in producing the reionization history, which is necessary for interpreting the future observations.

preprint2014arXiv

Stars and Reionization: The Cross-Correlation of the 21cm Line and the Near Infrared Background

With improving telescopes, it may now be possible to observe the Epoch of Reionization in multiple ways. We examine two of these observables - the excess light in the near-infrared background that may be due to high redshift stars and ionized HII bubbles, and the 21 cm emission from neutral hydrogen. Because these two forms of emission should result from different, mutually exclusive regions, an anticorrelation should exist between them. We discuss the strengths of using cross-correlations between these observations to learn more about high redshift star formation and reionization history. In particular, we create simulated maps of emission from both the near-infrared background and 21 cm emission. We find that these observations are anticorrelated, with the strongest anticorrelation originating from times when the universe is half ionized. This result is robust and does not depend on the properties of the stars themselves. Rather, it depends on the ionization history. Cross-correlations can provide redshift information, which the near-infrared background cannot provide alone. In addition, cross-correlations can help separate foreground emission from the true high redshift component, making it possible to say with greater certainty that we are indeed witnessing the Epoch of Reionization.

preprint2014arXiv

Studying reionization with the next generation of Ly-alpha emitter surveys

We study the prospects for constraining the ionized fraction of the intergalactic medium (IGM) at $z>6$ with the next generation of large Ly$α$ emitter surveys. We make predictions for the upcoming Subaru Hyper Suprime-Cam (HSC) Ly$α$ survey and a hypothetical spectroscopic survey performed with the James Webb Space Telescope (JWST). Considering various scenarios where the observed evolution of the Ly$α$ luminosity function of Ly$α$ emitters at $z>6$ is explained partly by an increasingly neutral IGM and partly by intrinsic galaxy evolution, we show how clustering measurements will be able to distinguish between these scenarios. We find that the HSC survey should be able to detect the additional clustering induced by a neutral IGM if the global IGM neutral fraction is greater than $\sim$20 per cent at $z=6.5$. If measurements of the Ly$α$ equivalent widths (EWs) are also available, neutral fractions as small as 10 per cent may be detectable by looking for correlation between the EW and the local number density of objects. In this case, if it should turn out that the IGM is significantly neutral at $z=6.5$ and the intrinsic EW distribution is relatively narrow, the observed EWs can also be used to construct a map of the locations and approximate sizes of the largest ionized regions. For the JWST survey, the results appear a bit less optimistic. Since such surveys probe a large range of redshifts, the effects of the IGM will be mixed up with any intrinsic galaxy evolution that is present, making it difficult to disentangle the effects. However, we show that a survey with the JWST will have a possibility of observing a large group of galaxies at $z\sim7$, which would be a strong indication of a partially neutral IGM.

preprint2013arXiv

Probing reionization with LOFAR using 21-cm redshift space distortions

One of the most promising ways to study the epoch of reionization (EoR) is through radio observations of the redshifted 21-cm line emission from neutral hydrogen. These observations are complicated by the fact that the mapping of redshifts to line-of-sight positions is distorted by the peculiar velocities of the gas. Such distortions can be a source of error if they are not properly understood, but they also encode information about cosmology and astrophysics. We study the effects of redshift space distortions on the power spectrum of 21-cm radiation from the EoR using large scale $N$-body and radiative transfer simulations. We quantify the anisotropy introduced in the 21-cm power spectrum by redshift space distortions and show how it evolves as reionization progresses and how it relates to the underlying physics. We go on to study the effects of redshift space distortions on LOFAR observations, taking instrument noise and foreground subtraction into account. We find that LOFAR should be able to directly observe the power spectrum anisotropy due to redshift space distortions at spatial scales around $k \sim 0.1$ Mpc$^{-1}$ after $\gtrsim$ 1000 hours of integration time. At larger scales, sample errors become a limiting factor, while at smaller scales detector noise and foregrounds make the extraction of the signal problematic. Finally, we show how the astrophysical information contained in the evolution of the anisotropy of the 21-cm power spectrum can be extracted from LOFAR observations, and how it can be used to distinguish between different reionization scenarios.

preprint2013arXiv

Reionization and the Cosmic Dawn with the Square Kilometre Array

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.

preprint2013arXiv

Simulating cosmic reionization: How large a volume is large enough?

We present the largest-volume (425 Mpc/h=607 Mpc on a side) full radiative transfer simulation of cosmic reionization to date. We show that there is significant additional power in density fluctuations at very large scales. We systematically investigate the effects this additional power has on the progress, duration and features of reionization, as well as on selected reionization observables. We find that comoving simulation volume of ~100 Mpc/h per side is sufficient for deriving a convergent mean reionization history, but that the reionization patchiness is significantly underestimated. We use jackknife splitting to quantify the convergence of reionization properties with simulation volume for both mean-density and variable-density sub-regions. We find that sub-volumes of ~100 Mpc/h per side or larger yield convergent reionization histories, except for the earliest times, but smaller volumes of ~50 Mpc/h or less are not well converged at any redshift. Reionization history milestones show significant scatter between the sub-volumes, of Delta z=0.6-1 for ~50 Mpc/h volumes, decreasing to Delta z=0.3-0.5 for ~100 Mpc/h volumes, and $Δz$~0.1 for ~200 Mpc/h volumes. If we only consider mean-density sub-regions the scatter decreases, but remains at Delta z~0.1-0.2 for the different size sub-volumes. Consequently, many potential reionization observables like 21-cm rms, 21-cm PDF skewness and kurtosis all show good convergence for volumes of ~200 Mpc/h, but retain considerable scatter for smaller volumes. In contrast, the three-dimensional 21-cm power spectra at large scales (k<0.25 h/Mpc) do not fully converge for any sub-volume size. These additional large-scale fluctuations significantly enhance the 21-cm fluctuations, which should improve the prospects of detection considerably, given the lower foregrounds and greater interferometer sensitivity at higher frequencies. (abridged)

preprint2013arXiv

The Kinetic Sunyaev-Zel'dovich effect as a probe of the physics of cosmic reionization: the effect of self-regulated reionization

We calculate the angular power spectrum of the Cosmic Microwave Background (CMB) temperature fluctuations induced by the kinetic Sunyaev-Zel'dovich (kSZ) effect from the epoch of reionization (EOR). We use detailed N-body+radiative transfer simulations to follow inhomogeneous reionization of the intergalactic medium (IGM). For the first time we take into account the "self-regulation" of reionization: star formation in low-mass dwarf galaxies (10^8 M_\sun \lesssim M \lesssim 10^9 M_\sun) or minihalos (10^5 M_\sun \lesssim M \lesssim 10^8 M_\sun) is suppressed if these halos form in the regions that were already ionized or Lyman-Werner dissociated. Some previous work suggested that the amplitude of the kSZ power spectrum from the EOR can be described by a two-parameter family: the epoch of half ionization and the duration of reionization. However, we argue that this picture applies only to simple forms of the reionization history which are roughly symmetric about the half-ionization epoch. In self-regulated reionization, the universe begins to be ionized early, maintains a low level of ionization for an extended period, and then finishes reionization as soon as high-mass atomically-cooling halos dominate. While inclusion of self-regulation affects the amplitude of the kSZ power spectrum only modestly (\sim 10 %), it can change the duration of reionization by a factor of more than two. We conclude that the simple two-parameter family does not capture the effect of a physical, yet complex, reionization history caused by self-regulation. When added to the post-reionization kSZ contribution, our prediction for the total kSZ power spectrum is below the current upper bound from the South Pole Telescope. Therefore, the current upper bound on the kSZ effect from the EOR is consistent with our understanding of the physics of reionization.

preprint2013arXiv

Will Nonlinear Peculiar Velocity and Inhomogeneous Reionization Spoil 21cm Cosmology from the Epoch of Reionization?

The 21cm background from the epoch of reionization is a promising cosmological probe: line-of-sight velocity fluctuations distort redshift, so brightness fluctuations in Fourier space depend upon angle, which linear theory shows can separate cosmological from astrophysical information. Nonlinear fluctuations in ionization, density and velocity change this, however. The validity and accuracy of the separation scheme are tested here for the first time, by detailed reionization simulations. The scheme works reasonably well early in reionization (< 40% ionized), but not late (> 80% ionized).

preprint2012arXiv

Detecting the Rise and Fall of the First Stars by Their Impact on Cosmic Reionization

The intergalactic medium was reionized before redshift z~6, most likely by starlight which escaped from early galaxies. The very first stars formed when hydrogen molecules (H2) cooled gas inside the smallest galaxies, minihalos of mass between 10^5 and 10^8 solar masses. Although the very first stars began forming inside these minihalos before redshift z~40, their contribution has, to date, been ignored in large-scale simulations of this cosmic reionization. Here we report results from the first reionization simulations to include these first stars and the radiative feedback that limited their formation, in a volume large enough to follow the crucial spatial variations that influenced the process and its observability. We show that, while minihalo stars stopped far short of fully ionizing the universe, reionization began much earlier with minihalo sources than without, and was greatly extended, which boosts the intergalactic electron-scattering optical depth and the large-angle polarization fluctuations of the cosmic microwave background significantly. Although within current WMAP uncertainties, this boost should be readily detectable by Planck. If reionization ended as late as z_ov<~7, as suggested by other observations, Planck will thereby see the signature of the first stars at high redshift, currently undetectable by other probes.

preprint2012arXiv

Light cone effect on the reionization 21-cm power spectrum

Observations of redshifted 21-cm radiation from neutral hydrogen during the epoch of reionization (EoR) are considered to constitute the most promising tool to probe that epoch. One of the major goals of the first generation of low frequency radio telescopes is to measure the 3D 21-cm power spectrum. However, the 21-cm signal could evolve substantially along the line of sight (LOS) direction of an observed 3D volume, since the received signal from different planes transverse to the LOS originated from different look-back times and could therefore be statistically different. Using numerical simulations we investigate this so-called light cone effect on the spherically averaged 3D 21-cm power spectrum. For this version of the power spectrum, we find that the effect mostly `averages out' and observe a smaller change in the power spectrum compared to the amount of evolution in the mean 21-cm signal and its rms variations along the LOS direction. Nevertheless, changes up to 50% at large scales are possible. In general the power is enhanced/suppressed at large/small scales when the effect is included. The cross-over mode below/above which the power is enhanced/suppressed moves toward larger scales as reionization proceeds. When considering the 3D power spectrum we find it to be anisotropic at the late stages of reionization and on large scales. The effect is dominated by the evolution of the ionized fraction of hydrogen during reionization and including peculiar velocities hardly changes these conclusions. We present simple analytical models which explain qualitatively all the features we see in the simulations.

preprint2012arXiv

On the Use of Ly-alpha Emitters as Probes of Reionization

We use numerical simulations to study the effects of the patchiness of a partly reionized intergalactic medium (IGM) on the observability of Ly-alpha emitters (LAEs) at high redshifts (z ~ 6). We present a new model that divides the Ly-alpha radiative transfer into a (circum-)galactic and an extragalactic (IGM) part, and investigate how the choice of intrinsic line model affects the IGM transmission results. We use our model to study the impact of neutral hydrogen on statistical observables such as the Ly-alpha restframe equivalent width (REW) distribution, the LAE luminosity function and the two-point correlation function. We find that if the observed changes in LAE luminosity functions and equivalent width distributions between z ~ 6 and z ~ 7 are to be explained by an increased IGM neutral fraction alone, we require an extremely late and rapid reionization scenario, where the Universe was ~ 40 % ionized at z = 7, ~ 50 % ionized at z = 6.5 and ~ 100 % ionized at z = 6. This is in conflict with other observations, suggesting that intrinsic LAE evolution at z > 6 cannot be completely neglected. We show how the two-point correlation function can provide more robust constraints once future observations obtain larger LAE samples, and provide predictions for the sample sizes needed to tell different reionization scenarios apart.

preprint2012arXiv

Prospects of observing a quasar HII region during the Epoch of Reionization with redshifted 21cm

We present a study of the impact of a bright quasar on the redshifted 21cm signal during the Epoch of Reionization (EoR). Using three different cosmological radiative transfer simulations, we investigate if quasars are capable of substantially changing the size and morphology of the H II regions they are born in. We choose stellar and quasar luminosities in a way that is favourable to seeing such an effect. We find that even the most luminous of our quasar models is not able to increase the size of its native H II region substantially beyond those of large H II regions produced by clustered stellar sources alone. However, the quasar H II region is found to be more spherical. We next investigate the prospects of detecting such H II regions in the redshifted 21cm data from the Low Frequency Array (LOFAR) by means of a matched filter technique. We find that H II regions with radii ~ 25 comoving Mpc or larger should have a sufficiently high detection probability for 1200 hours of integration time. Although the matched filter can in principle distinguish between more and less spherical regions, we find that when including realistic system noise this distinction can no longer be made. The strong foregrounds are found not to pose a problem for the matched filter technique. We also demonstrate that when the quasar position is known, the redshifted 21cm data can still be used to set upper limits on the ionizing photon rate of the quasar. If both the quasar position and its luminosity are known, the redshifted 21 cm data can set new constrains on quasar lifetimes.

preprint2012arXiv

Radiative transfer of energetic photons: X-rays and helium ionization in C2-Ray

We present an extension to the short-characteristic ray-tracing and non-equilibrium photon-ionization code C2Ray. The new version includes the effects of helium and improved multi-frequency heating. The motivation for this work is to be able to deal with harder ionizing spectra, such as for example from quasar-like sources during cosmic reionization. We review the basic algorithmic ingredients of C2-Ray before describing the changes implemented, which include a treatment of the full on the spot (OTS) approximation, secondary ionization, and multi-frequency photo-ionization and heating. We performed a series of tests against equilibrium solutions from CLOUDY as well as comparisons to the hydrogen only solutions by C2-Ray in the extensive code comparison in Iliev et al. (2006). We show that the full, coupled OTS approximation is more accurate than the simplified, uncoupled one. We find that also with helium and a multi-frequency set up, long timesteps (up to ~10% of the recombination time) still give accurate results for the ionization fractions. On the other hand, accurate results for the temperature set strong constrains on the timestep. The details of these constraints depend however on the optical depth of the cells. We use the new version of the code to confirm that the assumption made in many reionization simulations, namely that helium is singly ionized everywhere were hydrogen is, is indeed valid when the sources have stellar-like spectra.

preprint2012arXiv

Redshift Space Distortion of the 21cm Background from the Epoch of Reionization I: Methodology Re-examined

The peculiar velocity of the intergalactic gas responsible for the cosmic 21cm background from the epoch of reionization and beyond introduces an anisotropy in the three-dimensional power spectrum of brightness temperature fluctuations. Measurement of this anisotropy by future 21cm surveys is a promising tool for separating cosmology from 21cm astrophysics. However, previous attempts to model the signal have often neglected peculiar velocity or only approximated it crudely. This paper re-examines the effects of peculiar velocity on the 21cm signal in detail, improving upon past treatment and addressing several issues for the first time. (1) We show that properly accounting for finite optical depth eliminates the unphysical divergence of 21cm brightness temperature in overdense regions of the IGM found by previous work that employed the usual optically-thin approximation. (2) The approximation made previously to circumvent the diverging brightness temperature problem by capping velocity gradient can misestimate the power spectrum on all scales. (3) The observed power spectrum in redshift-space remains finite even in the optically-thin approximation if one properly accounts for the redshift-space distortion. However, results that take full account of finite optical depth show that this approximation is only accurate in the limit of high spin temperature. (4) The linear theory for redshift-space distortion results in ~30% error in the observationally relevant wavenumber range, at the 50% ionized epoch. (5) We describe and test two numerical schemes to calculate the 21cm signal from reionization simulations to incorporate peculiar velocity effects in the optically-thin approximation accurately. One is particle-based, the other grid-based, and while the former is most accurate, we demonstrate that the latter is computationally more efficient and can achieve sufficient accuracy. [Abridged]

preprint2012arXiv

Simulating Cosmic Reionization and the Radiation Backgrounds from the Epoch of Reionization

Large-scale reionization simulations are described which combine the results of cosmological N-body simulations that model the evolving density and velocity fields and identify the galactic halo sources, with ray-tracing radiative transfer calculations which model the nonequilibrium ionization of the intergalactic medium. These simulations have been used to predict some of the signature effects of reionization on cosmic radiation backgrounds, including the CMB, near-IR, and redshifted 21cm backgrounds. We summarize some of our recent progress in this work, and address the question of whether observations of such signature effects can be used to distinguish the relative contributions of galaxies of different masses to reionization.

preprint2012arXiv

The inhomogenous reionization of the inter-galactic medium by metal-poor globular clusters

We present detailed radiative transfer simulations of the reionization history of the Milky Way by metal-poor globular clusters. We identify potential metal-poor globular cluster candidates within the Aquarius simulation using dark matter halo velocity dispersions. We calculate the local ionization fields via a photon-conserving, three dimensional non-equilibrium chemistry code and allow the model to propagate through to the present day. The key feature of the model is that globular cluster formation is suppressed if the local gas is ionized. We find that our spatial treatment of the ionization field leads to drastically different numbers and spatial distributions when compared to models where globular cluster formation is simply truncated at a given redshift. We find that it is possible for metal-poor globular clusters to have formed via the dark matter halo formation channel as our secondary model (delayed formation) combined with truncation at z = 10 produces radial distributions statistically consistent with that of the Milky Way metal-poor globular clusters. If globular clusters do indeed form within high-redshift dark matter halos, if only in-part, their contributions to the reionization of the local (i.e. 2^3 h^-3 Mpc^3 centred on the host galaxy) volume and mass by redshift 10 could be as high as 98% and 90%, respectively. In our photon poorest model, this contribution drops to 60% and 50%. The surviving clusters in all models have a narrow average age range (mean = 13.34 Gyr, σ= 0.04 Gyr) consistent with current ages estimates of the Milky Way metal-poor globular clusters. We also test a simple dynamical destruction model and estimate that ~60% of all metal-poor globular clusters formed at high redshift have since been destroyed via tidal interactions with the host galaxy.

preprint2012arXiv

Topology and Sizes of HII Regions during Cosmic Reionization

We use the results of large-scale simulations of reionization to explore methods for characterizing the topology and sizes of HII regions during reionization. We use four independent methods for characterizing the sizes of ionized regions. Three of them give us a full size distribution: the friends-of-friends (FOF) method, the spherical average method (SPA) and the power spectrum (PS) of the ionized fraction. These latter three methods are complementary: While the FOF method captures the size distribution of the small scale H II regions, which contribute only a small amount to the total ionization fraction, the spherical average method provides a smoothed measure for the average size of the H II regions constituting the main contribution to the ionized fraction, and the power spectrum does the same while retaining more details on the size distribution. Our fourth method for characterizing the sizes of the H II regions is the average size which results if we divide the total volume of the H II regions by their total surface area, (i.e. 3V/A), computed in terms of the ratio of the corresponding Minkowski functionals of the ionized fraction field. To characterize the topology of the ionized regions, we calculate the evolution of the Euler Characteristic. We find that the evolution of the topology during the first half of reionization is consistent with inside-out reionization of a Gaussian density field. We use these techniques to investigate the dependence of size and topology on some basic source properties, such as the halo mass-to-light ratio, susceptibility of haloes to negative feedback from reionization, and the minimum halo mass for sources to form. We find that suppression of ionizing sources within ionized regions slows the growth of H II regions, and also changes their size distribution. Additionally, the topology of simulations including suppression is more complex. (abridged)

preprint2011arXiv

Can 21-cm observations discriminate between high-mass and low-mass galaxies as reionization sources?

The prospect of detecting the first galaxies by observing their impact on the intergalactic medium as they reionized it during the first billion years leads us to ask whether such indirect observations are capable of diagnosing which types of galaxies were most responsible for reionization. We attempt to answer this by considering a set of large-scale radiative transfer simulations of reionization in sufficiently large volumes to make statistically meaningful predictions of observable signatures, while also directly resolving all atomically-cooling halos down to 10^8 M_solar. We focus here on predictions of the 21-cm background, to see if upcoming observations are capable of distinguishing a universe ionized primarily by high-mass halos from one in which both high-mass and low-mass halos are responsible, and to see how these results depend upon the uncertain source efficiencies. We find that 21-cm fluctuation power spectra observed by the first generation EoR/21-cm radio interferometer arrays should be able to distinguish the case of reionization by high-mass halos alone from that by both high- and low-mass halos, together. Some reionization scenarios yield very similar power spectra and rms evolution and thus can only be discriminated by their different mean reionization history and 21-cm PDF distributions. We find that the skewness of the 21-cm PDF distribution smoothed over LOFAR-like window shows a clear feature correlated with the rise of the rms due to patchiness. Measurements of the mean photoionization rates are sensitive to the average density of the regions being studied and therefore could be strongly skewed in certain cases. (abridged)

preprint2011arXiv

Observational constraints on supermassive dark stars

Some of the first stars could be cooler and more massive than standard stellar models would suggest, due to the effects of dark matter annihilation in their cores. It has recently been argued that such objects may attain masses in the 10^4--10^7 solar mass range, and that such supermassive dark stars should be within reach of the upcoming James Webb Space Telescope. Notwithstanding theoretical difficulties with this proposal, we argue here that some of these objects should also be readily detectable with both the Hubble Space Telescope and ground-based 8--10 m class telescopes. Existing survey data already place strong constraints on 10^7 solar mass dark stars at z~10. We show that such objects must be exceedingly rare or short-lived to have avoided detection.

preprint2010arXiv

A Physical Model of Lyman Alpha Emitters

We present a simple physical model for populating dark matter halos with Lyman Alpha Emiiters(LAEs) and predict the physical properties of LAEs at z~3-7. The central tenet of this model is that the Ly-alpha luminosity is proportional to the star formation rate (SFR) which is directly related to the halo mass accretion rate. The only free parameter in our model is then the star-formation efficiency (SFE). An efficiency of 2.5% provides the best-fit to the Ly-alpha luminosity function (LF) at redshift z=3.1, and we use this SFE to construct Ly-alpha LFs at other redshifts. Our model reproduce the Ly-alpha LFs, stellar ages, SFR ~1-10; Msun/yr, stellar masses ~ 10^7-10^8 Msun and the clustering properties of LAEs at z~3-7. We find the spatial correlation lengths ro ~ 3-6 Mpc/h, in agreement with the observations. Finally, we estimate the field-to-field variation ~ 30% for current volume and flux limited surveys, again consistent with observations. Our results suggest that the star formation, and hence Ly-alpha emission in LAEs is powered by the accretion of new material, and that the physical properties of LAEs do not evolve significantly over a wide range of redshifts. Relating the accreted mass, rather than the total mass of halos, to the Ly-alpha luminosity of LAEs naturally gives rise to the duty cycle of LAEs.

preprint2010arXiv

Power spectrum extraction for redshifted 21-cm epoch of reionization experiments: the LOFAR case

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

preprint2010arXiv

Reionization of the Local Group of Galaxies

We present the first detailed structure formation and radiative transfer simulations of the reionization history of our cosmic neighbourhood. To this end, we follow the formation of the Local Group of galaxies and nearby clusters by means of constrained simulations, which use the available observational constraints to construct a representation of those structures which reproduces their actual positions and properties at the present time. We find that the reionization history of the Local Group is strongly dependent on the assumed photon production efficiencies of the ionizing sources, which are still poorly constrained. If sources are relatively efficient, i.e. the process is 'photon-rich', the Local Group is primarily ionized externally by the nearby clusters. Alternatively, if the sources are inefficient, i.e. reionization is 'photon-poor' the Local Group evolves largely isolated and reionizes itself. The mode of reionization, external vs. internal, has important implications for the evolution of our neighbourhood, in terms of e.g. its satellite galaxy populations and primordial stellar populations. This therefore provides an important avenue for understanding the young universe by detailed studies of our nearby structures.

preprint2009arXiv

Cosmological Radiative Transfer Comparison Project II: The Radiation-Hydrodynamic Tests

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

preprint2009arXiv

Cross-correlation Study between the Cosmological 21-cm Signal and the kinetic Sunyaev-Zel'dovich effect

The Universe's Epoch of Reionization (EoR) can be studied using a number of observational probes that provide complementary or corroborating information. Each of these probes suffers from its own systematic and statistical uncertainties. It is therefore useful to consider the mutual information that these data sets contain. In this paper we present a cross-correlation study between the kinetic Sunyaev-Zel'dovich effect (kSZ) -- produced by the scattering of CMB photons off free electrons produced during the reionization process -- and the cosmological 21cm signal -- which reflects the neutral hydrogen content of the Universe, as a function of redshift. The study is carried out using a simulated reionization history in 100Mpc/h scale N-body simulations with radiative transfer. In essence we find that the two probes anti-correlate. The significance of the anti-correlation signal depends on the extent of the reionization process, wherein extended histories result in a much stronger signal compared to instantaneous cases. Unfortunately however, once the primary CMB fluctuations are included into our simulation they serve as a source of large correlated noise that renders the cross-correlation signal insignificant, regardless of the reionization scenario.

preprint2009arXiv

Measuring the History of Cosmic Reionization using the 21-cm PDF from Simulations

The 21-cm PDF (i.e., distribution of pixel brightness temperatures) is expected to be highly non-Gaussian during reionization and to provide important information on the distribution of density and ionization. We measure the 21-cm PDF as a function of redshift in a large simulation of cosmic reionization and propose a simple empirical fit. Guided by the simulated PDF, we then carry out a maximum likelihood analysis of the ability of upcoming experiments to measure the shape of the 21-cm PDF and derive from it the cosmic reionization history. Under the strongest assumptions, we find that upcoming experiments can measure the reionization history in the mid to late stages of reionization to 1-10% accuracy. Under a more flexible approach that allows for four free parameters at each redshift, a similar accuracy requires the lower noise levels of second-generation 21-cm experiments.

preprint2009arXiv

Radiation-magnetohydrodynamic simulations of the photoionisation of magnetised globules

We present the first three-dimensional radiation-magnetohydrodynamic simulations of the photoionisation of a dense, magnetised molecular globule by an external source of ultraviolet radiation. We find that, for the case of a strong ionising field, significant deviations from the non-magnetic evolution are seen when the initial magnetic field threading the globule has an associated magnetic pressure that is greater than one hundred times the gas pressure. In such a strong-field case, the photoevaporating globule will adopt a flattened or "curled up" shape, depending on the initial field orientation, and magnetic confinement of the ionised photoevaporation flow can lead to recombination and subsequent fragmentation during advanced stages of the globule evolution. We find suggestive evidence that such magnetic effects may be important in the formation of bright, bar-like emission features in H II regions. We include simple but realistic fits to heating and cooling rates in the neutral and molecular gas in the vicinity of a high-mass star cluster and show that the frequently used isothermal approximation can lead to an overestimate of the importance of gravitational instability in the radiatively imploded globule. For globules within 2 parsecs of a high-mass star cluster, we find that heating by stellar x rays prevents the molecular gas from cooling below 50 K.

preprint2009arXiv

The Inhomogeneous Background of Hydrogen-Molecule Dissociating Radiation during Cosmic Reionization

The first, self-consistent calculations are presented of the cosmological, H2-dissociating UV background produced during the epoch of reionization by the sources of reionization. Large-scale radiative transfer simulations of reionization trace the impact of all the ionizing starlight on the IGM from all the sources in our simulation volume down to dwarf galaxies of mass ~10^8 Msun, identified by very high-resolution N-body simulations, including the self-regulating effect of IGM photoheating on dwarf galaxy formation. The UV continuum emitted below 13.6 eV by each source is then transferred through the same IGM, attenuated by atomic H Lyman series resonance lines, to predict the evolution of the inhomogeneous radiation background in the Lyman-Werner bands of H2 between 11 and 13.6 eV. On average, the intensity of this Lyman-Werner background is found to rise to the threshold level at which dissociation suppresses H2 cooling and star formation inside minihalos, long before reionization is complete. Spatial variations in the Lyman-Werner background are found which result from the clustering of sources associated with large-scale structure formation, such that intensity fluctuations correlate with matter density fluctuations. As a result, the Lyman-Werner background rises to the threshold level for H2 suppression earlier in the vicinity of the reionization sources and their H II regions.

preprint2006arXiv

The Kinetic Sunyaev-Zel'dovich Effect from Patchy Reionization: the View from the Simulations

We present the first calculation of the kinetic Sunyaev-Zel'dovich (kSZ) effect due to the inhomogeneus reionization of the universe based on detailed large-scale radiative transfer simulations of reionization. The resulting sky power spectra peak at l=2000-8000 with maximum values of l(l+1)C_l/(2pi)_{max}~4-7x10^{-13}. The scale roughly corresponds to the typical ionized bubble sizes observed in our simulations, of ~5-20 Mpc. The kSZ anisotropy signal from reionization dominates the primary CMB signal above l=3000. At large scales the patchy kSZ signal depends only on the source efficiencies. It is higher when sources are more efficient at producing ionizing photons, since such sources produce larger ionized regions, on average, than less efficient sources. The introduction of sub-grid gas clumping in the radiative transfer simulations produce significantly more power at small scales, but has little effect at large scales. The patchy reionization kSZ signal is dominated by the post-reionization signal from fully-ionized gas, but the two contributions are of similar order at scales l~3000-10^4, indicating that the kSZ anisotropies from reionization are an important component of the total kSZ signal at these scales.

preprint2006arXiv

The Kinetic Sunyaev-Zel'dovich Effect from Radiative Transfer Simulations of Patchy Reionization

We present the first calculation of the kinetic Sunyaev-Zel'dovich (kSZ) effect due to the inhomogeneous reionization of the universe based on detailed large-scale radiative transfer simulations of reionization. The resulting sky power spectra peak at l=2000-8000 with maximum values of l^2C_l~1\times10^{-12}. The peak scale is determined by the typical size of the ionized regions and roughly corresponds to the ionized bubble sizes observed in our simulations, ~5-20 Mpc. The kSZ anisotropy signal from reionization dominates the primary CMB signal above l=3000. This predicted kSZ signal at arcminute scales is sufficiently strong to be detectable by upcoming experiments, like the Atacama Cosmology Telescope and South Pole Telescope which are expected to have ~1' resolution and ~muK sensitivity. The extended and patchy nature of the reionization process results in a boost of the peak signal in power by approximately one order of magnitude compared to a uniform reionization scenario, while roughly tripling the signal compared with that based upon the assumption of gradual but spatially uniform reionization. At large scales the patchy kSZ signal depends largely on the ionizing source efficiencies and the large-scale velocity fields: sources which produce photons more efficiently yield correspondingly higher signals. The introduction of sub-grid gas clumping in the radiative transfer simulations produces significantly more power at small scales, and more non-Gaussian features, but has little effect at large scales. The patchy nature of the reionization process roughly doubles the total observed kSZ signal for l~3000-10^4 compared to non-patchy scenarios with the same total electron-scattering optical depth.

preprint2003arXiv

The post-AGB evolution of AGB mass loss variations

We present new numerical hydrodynamical modelling of the evolution of Asymptotic Giant Branch (AGB) mass loss fluctuations during the post-AGB/Planetary Nebula phase. These models show that after ionization, the observable effects of the mass loss fluctuations disappear in a few thousand years, consistent with the fact that only few PNe have been found to be surrounded by `rings'. We derive the observational characteristics of these rings, and compare them to reported observations of the rings around NGC 6543, finding a good match of emission properties and line shapes. We predict small variations in the observable electron temperatures.

preprint1997arXiv

Outflow collimation in young stellar objects

In this paper we explore the effect of radiative losses on purely hydrodynamic jet collimation models applicable to Young Stellar Objects (YSOs). In our models aspherical bubbles form from the interaction of a central YSO wind with an aspherical circum-protostellar density distribution. Building on a previous non-radiative study (Frank & Mellema 1996) we demonstrate that supersonic jets are a natural and robust consequence of aspherical wind-blown bubble evolution. The simulations show that the addition of radiative cooling makes the hydrodynamic collimation mechanisms studied by Frank & Mellema (1996) more effective. An unexpected result is the production of cool, dense jets forming through conical converging flows at the poles of the bubbles. For steady winds the formation of these jets occurs early in the bubble evolution. At later times we find that the dynamical and cooling time scales for the jet material become similar. Since observations of YSO jets show considerable variability in the jet beam we present a simple one-dimensional (1-D) model for the evolution of a variable wind interacting with an accreting environment. These models indicate that, in the presence of a varying protostellar wind, the hydrodynamic collimation processes studied in our simulations can produce cool jets with sizes and time scales consistent with observations.

preprint1995arXiv

Radiation gasdynamics of planetary nebulae -- VI. The evolution of aspherical planetary nebulae

This paper reports the results of the numerical study of the formation of aspherical Planetary Nebulae through the Generalised Interacting Winds Model, taking into account the effects caused by the evolving central star and fast wind. The results show for the first time that aspherical nebulae do form within the required time scale. Consideration of the development of the nebula shows that in the early stages it is the ionization of the aspherical AGB wind that contributes considerably to the shaping of the nebula. Furthermore the passing through of the ionization front may modifies the density distribution in the slow wind, leading to the formation of a surrounding envelope, and sometimes a different morphology for the nebula than was to be expected from the initial conditions. I consider how the different phases of ionization fronts and wind swept bubbles can be observationally distinguished.