Using the latest upper limits on the 21-cm power spectrum at $z\approx9.1$ from the Low Frequency Array (LOFAR), we explore regions of parameter space which are inconsistent with the data. We use 21CMMC, a Monte Carlo Markov Chain sampler of 21cmFAST which directly forward models the 3D cosmic 21-cm signal in a fully Bayesian framework. We use the astrophysical parameterisation from 21cmFAST, which includes mass-dependent star formation rates and ionising escape fractions as well as soft-band X-ray luminosities to place limits on the properties of the high-$z$ galaxies. Further, we connect the disfavoured regions of parameter space with existing observational constraints on the Epoch of Reionisation such as ultra-violet (UV) luminosity functions, background UV photoionisation rate, intergalactic medium (IGM) neutral fraction and the electron scattering optical depth. We find that all models exceeding the 21-cm signal limits set by LOFAR at $z\approx9.1$ are excluded at $\gtrsim2\sigma$ by other probes. Finally, we place limits on the IGM spin temperature from LOFAR, disfavouring at 95 per cent confidence spin temperatures below $\sim2.6$ K across an IGM neutral fraction range of $0.15 \lesssim \bar{x}_{H{\scriptscriptstyle I}} \lesssim 0.6$. Note, these limits are only obtained from 141 hrs of data in a single redshift bin. With tighter upper limits, across multiple redshift bins expected in the near future from LOFAR, more viable models will be ruled out. Our approach demonstrates the potential of forward modelling tools such as 21CMMC in combining 21-cm observations with other high-$z$ probes to constrain the astrophysics of galaxies.

The ARCADE2 and LWA1 experiments have claimed an excess over the Cosmic Microwave Background (CMB) at low radio frequencies. If the cosmological high-redshift contribution to this radio background is between 0.1% and 22% of the CMB at $1.42\,$GHz, it could explain the tentative EDGES Low-Band detection of the anomalously deep absorption in the 21-cm signal of neutral hydrogen (Fialkov & Barkana 2019). We use the upper limit on the 21-cm signal from the Epoch of Reionization ($z=9.1$) based on $141\,$hours of observations with LOFAR (Mertens et al. 2020) to evaluate the contribution of the high redshift Universe to the detected radio background. Marginalizing over astrophysical properties of star-forming halos, we find (at 68% C.L.) that the cosmological radio background can be at most 0.7% of the CMB at $1.42\,$GHz. This limit rules out strong contribution of the high-redshift Universe to the ARCADE2 and LWA1 measurements. Even though LOFAR places limit on the extra radio background, excess of $0.1-0.7$% over the CMB (at $1.42\,$GHz) is still allowed and could explain the EDGES Low-Band detection. We also constrain the thermal and ionization state of the gas at $z = 9.1$, and put limits on the properties of the first star-forming objects. We find that, in agreement with the limits from EDGES High-Band data, LOFAR data disfavour scenarios with inefficient X-ray sources and cases where the Universe was ionized by massive halos only.

A new upper limit on the 21-cm signal power spectrum at a redshift of $z \approx 9.1$ is presented, based on 141 hours of data obtained with the Low-Frequency Array (LOFAR). The analysis includes significant improvements in spectrally-smooth gain-calibration, Gaussian Process Regression (GPR) foreground mitigation and optimally-weighted power spectrum inference. Previously seen `excess power' due to spectral structure in the gain solutions has markedly reduced but some excess power still remains with a spectral correlation distinct from thermal noise. This excess has a spectral coherence scale of $0.25 - 0.45$\,MHz and is partially correlated between nights, especially in the foreground wedge region. The correlation is stronger between nights covering similar local sidereal times. A best 2-$\sigma$ upper limit of $\Delta^2_{21} < (73)^2\,\mathrm{mK^2}$ at $k = 0.075\,\mathrm{h\,cMpc^{-1}}$ is found, an improvement by a factor $\approx 8$ in power compared to the previously reported upper limit. The remaining excess power could be due to residual foreground emission from sources or diffuse emission far away from the phase centre, polarization leakage, chromatic calibration errors, ionosphere, or low-level radio-frequency interference. We discuss future improvements to the signal processing chain that can further reduce or even eliminate these causes of excess power.

We derive constraints on the thermal and ionization states of the intergalactic medium (IGM) at redshift $\approx$ 9.1 using new upper limits on the 21-cm power spectrum measured by the LOFAR radio-telescope and a prior on the ionized fraction at that redshift estimated from recent cosmic microwave background (CMB) observations. We have used results from the reionization simulation code GRIZZLY and a Bayesian inference framework to constrain the parameters which describe the physical state of the IGM. We find that, if the gas heating remains negligible, an IGM with ionized fraction $\gtrsim 0.13$ and a distribution of the ionized regions with a characteristic size $\gtrsim 8 ~h^{-1}$ comoving megaparsec (Mpc) and a full width at the half maximum (FWHM) $\gtrsim 16 ~h^{-1}$ Mpc is ruled out. For an IGM with a uniform spin temperature $T_{\rm S} \gtrsim 3$ K, no constraints on the ionized component can be computed. If the large-scale fluctuations of the signal are driven by spin temperature fluctuations, an IGM with a volume fraction $\lesssim 0.34$ of heated regions with a temperature larger than CMB, average gas temperature 7-160 K and a distribution of the heated regions with characteristic size 3.5-70 $h^{-1}$ Mpc and FWHM of $\lesssim 110$ $h^{-1}$ Mpc is ruled out. These constraints are within the 95 per cent credible intervals. With more stringent future upper limits from LOFAR at multiple redshifts, the constraints will become tighter and will exclude an increasingly large region of the parameter space.

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.

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.

The resonance scattering of Ly$\alpha$ 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$\alpha$ 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$\alpha$ background is ruled out. Instead the results favour a weak Ly$\alpha$ 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$\alpha$ heating is taken into account.

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.

We present the prospects of extracting information about the Epoch of Reionization by identifying the remaining neutral regions, referred to as islands, in tomographic observations of the redshifted 21-cm signal. Using simulated data sets we show that at late times the 21-cm power spectrum is fairly insensitive to the details of the reionization process but that the properties of the neutral islands can distinguish between different reionization scenarios. We compare the properties of these islands with those of ionized bubbles. At equivalent volume filling fractions, neutral islands tend to be fewer in number but larger compared to the ionized bubbles. In addition, the evolution of the size distribution of neutral islands is found to be slower than that of the ionized bubbles and also their percolation behaviour differs substantially. Even though the neutral islands are relatively rare, they will be easier to identify in observations with the low-frequency component of the Square Kilometre Array (SKA-Low) due to their larger size and the lower noise levels at lower redshifts. The size distribution of neutral islands at the late stages of reionization is found to depend on the source properties, such as the ionizing efficiency of the sources and their minimum mass. We find the longest line of sight through a neutral region to be more than 100 comoving Mpc until very late stages (90-95 per cent reionized), which may have relevance for the long absorption trough at $z=5.6-5.8$ in the spectrum of quasar ULAS J0148+0600.

Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged]

Observations of the redshifted 21-cm hyperfine line of neutral hydrogen from early phases of the Universe such as Cosmic Dawn and the Epoch of Reionization promise to open a new window onto the early formation of stars and galaxies. We present the first upper limits on the power spectrum of redshifted 21-cm brightness temperature fluctuations in the redshift range $z = 19.8 - 25.2$ ($54-68$ MHz frequency range) using 14 hours of data obtained with the LOFAR-Low Band Antenna (LBA) array. We also demonstrate the application of a multiple pointing calibration technique to calibrate the LOFAR-LBA dual-pointing observations centred on the North Celestial Pole and the radio galaxy 3C220.3. We observe an unexplained excess of $\sim 30-50\%$ in Stokes $I$ noise compared to Stokes $V$ for the two observed fields, which decorrelates on $\gtrsim 12$ seconds and might have a physical origin. We show that enforcing smoothness of gain errors along frequency direction during calibration reduces the additional variance in Stokes $I$ compared Stokes $V$ introduced by the calibration on sub-band level. After subtraction of smooth foregrounds, we achieve a $2\sigma$ upper limit on the 21-cm power spectrum of $\Delta_{21}^2 < (14561\,\text{mK})^2$ at $k\sim 0.038\,h\,\text{cMpc}^{-1}$ and $\Delta_{21}^2 < (14886\,\text{mK})^2$ at $k\sim 0.038 \,h\,\text{cMpc}^{-1}$ for the 3C220 and NCP fields respectively and both upper limits are consistent with each other. The upper limits for the two fields are still dominated by systematics on most $k$ modes.

The upcoming radio interferometer Square Kilometre Array (SKA) is expected to directly detect the redshifted 21-cm signal from the neutral hydrogen present during the Cosmic Dawn. Temperature fluctuations from X-ray heating of the neutral intergalactic medium can dominate the fluctuations in the 21-cm signal from this time. This heating depends on the abundance, clustering, and properties of the X-ray sources present, which remain highly uncertain. We present a suite of three new large-volume, 349\,Mpc a side, fully numerical radiative transfer simulations including QSO-like sources, extending the work previously presented in Ross et al. (2017). The results show that our QSOs have a modest contribution to the heating budget, yet significantly impact the 21-cm signal. Initially, the power spectrum is boosted on large scales by heating from the biased QSO-like sources, before decreasing on all scales. Fluctuations from images of the 21-cm signal with resolutions corresponding to SKA1-Low at the appropriate redshifts are well above the expected noise for deep integrations, indicating that imaging could be feasible for all the X-ray source models considered. The most notable contribution of the QSOs is a dramatic increase in non-Gaussianity of the signal, as measured by the skewness and kurtosis of the 21-cm probability distribution functions. However, in the case of late Lyman-$\alpha$ saturation, this non-Gaussianity could be dramatically decreased particularly when heating occurs earlier. We conclude that increased non-Gaussianity is a promising signature of rare X-ray sources at this time, provided that Lyman-$\alpha$ saturation occurs before heating dominates the 21-cm signal.

Potential small-scale discrepancies in the picture of galaxy formation painted by the $\Lambda$CDM paradigm have led to considerations of modified dark matter models. One such dark matter model that has recently attracted much attention is fuzzy dark matter (FDM). In FDM models, the dark matter is envisaged to be an ultra-light scalar field with a particle mass $m_{\rm FDM} \sim 10^{-22} $ eV. This yields astronomically large de Broglie wavelengths which can suppress small-scale structure formation and give rise to the observed kpc-sized density cores in dwarf galaxies. We investigate the evolution of the 21-cm signal during Cosmic Dawn and the Epoch of Reionization (EoR) in $\Lambda$FDM cosmologies using analytical models. The delay in source formation and the absence of small halos in $\Lambda$FDM significantly postpone the Ly$\alpha$ coupling, heating, as well as the reionization of the neutral hydrogen of the intergalactic medium. As a result, the absorption feature in the evolution of the global 21-cm signal has a significantly smaller full width at half maximum ($\Delta z \lesssim 3$), than $\Lambda$CDM ($\Delta z \simeq 6$). This alone rules out $m_{\rm FDM} < 6 \times 10^{-22}$ eV as a result of the $2\sigma$ lower limit $\Delta z \gtrsim 4$ from EDGES High-Band. As a result, $\Lambda$FDM is not a viable solution to the potential small-scale problems facing $\Lambda$CDM. Finally, we show that any detection of the 21-cm signal at redshifts $z > 14$ by interferometers such as the SKA can also exclude $\Lambda$FDM models.

This white paper highlights the crucial and urgent synergies required between WFIRST, Subaru Hyper Suprime-Cam or other >25m-class telescopes galaxy observations and SKA 21cm measurements to constrain the nature of reionization (ionization history and topology) and its sources.

The 21-cm signal from the epoch of reionization is non-Gaussian. Current radio telescopes are focused on detecting the 21-cm power spectrum, but in the future the Square Kilometre Array is anticipated to provide a first measurement of the bispectrum. Previous studies have shown that the position-dependent power spectrum is a simple and efficient way to probe the squeezed-limit bispectrum. In this approach, the survey is divided into subvolumes and the correlation between the local power spectrum and the corresponding mean density of the subvolume is computed. This correlation is equivalent to an integral of the bispectrum in the squeezed limit, but is much simpler to implement than the usual bispectrum estimators. It also has a clear physical interpretation: it describes how the small-scale power spectrum of tracers such as galaxies and the 21-cm signal respond to a large-scale environment. Reionization naturally couples large and small scales as ionizing radiation produced by galactic sources can travel up to tens of Megaparsecs through the intergalactic medium during this process. Here we apply the position-dependent power spectrum approach to fluctuations in the 21-cm background from reionization. We show that this statistic has a distinctive evolution in time that can be understood with a simple analytic model. We also show that the statistic can easily distinguish between simple "inside-out" and "outside-in" models of reionization. The position-dependent power spectrum is thus a promising method to validate the reionization signal and to extract higher-order information on this process.

We present analysis of the normalised 21-cm bispectrum from fully-numerical simulations of intergalactic-medium heating by stellar sources and high-mass X-ray binaries (HMXB) during the cosmic dawn. Lyman-$\alpha$ coupling is assumed to be saturated, we therefore probe the nature of non-Gaussianities produced by X-ray heating processes. We find the evolution of the normalised bispectrum to be very different from that of the power spectrum. It exhibits a turnover whose peak moves from large to small scales with decreasing redshift, and corresponds to the typical separation of emission regions. This characteristic scale reduces as more and more regions move into emission with time. Ultimately, small-scale fluctuations within heated regions come to dominate the normalised bispectrum, which at the end of the simulation is almost entirely driven by fluctuations in the density field. To establish how generic the qualitative evolution of the normalised bispectrum we see in the stellar + HMXB simulation is, we examine several other simulations - two fully-numerical simulations that include QSO sources, and two with contrasting source properties produced with the semi-numerical simulation 21cmFAST. We find the qualitative evolution of the normalised bispectrum during X-ray heating to be generic, unless the sources of X-rays are, as with QSOs, less numerous and so exhibit more distinct isolated heated profiles. Assuming mitigation of foreground and instrumental effects are ultimately effective, we find that we should be sensitive to the normalised bispectrum during the epoch of heating, so long as the spin temperature has not saturated by $z \approx 19$.

The ability of the future low frequency component of the Square Kilometre Array radio telescope (SKA-Low) to produce tomographic images of the redshifted 21-cm signal will enable direct studies of the evolution of the sizes and shapes of ionized regions during the Epoch of Reionization. However, a reliable identification of ionized regions in noisy interferometric data is not trivial. Here, we introduce an image processing method known as superpixels for this purpose. We compare this method with two other previously proposed ones, one relying on a chosen threshold and the other employing automatic threshold determination using the K-Means algorithm. We use a correlation test and compare power spectra and bubble size distributions to show that the superpixels method provides a better identification of ionized regions, especially in the case of noisy data. We also describe some possible additional applications of the superpixel method, namely the derivation of the ionization history and constraints on the source properties in specific regions.

The epoch of reionization (EoR) 21-cm signal is expected to be highly non-Gaussian in nature and this non-Gaussianity is also expected to evolve with the progressing state of reionization. Therefore the signal will be correlated between different Fourier modes ($k$). The power spectrum will not be able capture this correlation in the signal. We use a higher-order estimator -- the bispectrum -- to quantify this evolving non-Gaussianity. We study the bispectrum using an ensemble of simulated 21-cm signal and with a large variety of $k$ triangles. We observe two competing sources driving the non-Gaussianity in the signal: fluctuations in the neutral fraction ($x_{\rm HI}$) field and fluctuations in the matter density field. We find that the non-Gaussian contribution from these two sources vary, depending on the stage of reionization and on which $k$ modes are being studied. We show that the sign of the bispectrum works as a unique marker to identify which among these two components is driving the non-Gaussianity. We propose that the sign change in the bispectrum, when plotted as a function of triangle configuration $\cos{\theta}$ and at a certain stage of the EoR can be used as a confirmative test for the detection of the 21-cm signal. We also propose a new consolidated way to visualize the signal evolution (with evolving $\overline{x}_{\rm HI}$ or redshift), through the trajectories of the signal in a power spectrum and equilateral bispectrum i.e. $P(k)-B(k, k, k)$ space.

The upcoming radio interferometer Square Kilometre Array is expected to directly detect the redshifted 21-cm signal from the Cosmic Dawn for the first time. In this era temperature fluctuations from X-ray heating of the neutral intergalactic medium can impact this signal dramatically. Previously, in (Ross et al, 2017), we presented the first large-volume, 244 $h^{-1}$Mpc = 349 Mpc a side, fully numerical radiative transfer simulations of X-ray heating. This work is a follow-up where we now also consider QSO-like sources in addition to high mass X-ray binaries. Images of the two cases are clearly distinguishable at SKA1-LOW resolution and have RMS fluctuations above the expected noise. The inclusion of QSOs leads to a dramatic increase in non-Gaussianity of the signal, as measured by the skewness and kurtosis of the 21-cm signal. We conclude that this increased non-Gaussianity is a promising signature of early QSOs.

The upcoming SKA1-Low radio interferometer will be sensitive enough to produce tomographic imaging data of the redshifted 21-cm signal from the Epoch of Reionization. Due to the non-Gaussian distribution of the signal, a power spectrum analysis alone will not provide a complete description of its properties. Here, we consider an additional metric which could be derived from tomographic imaging data, namely the bubble size distribution of ionized regions. We study three methods that have previously been used to characterize bubble size distributions in simulation data for the hydrogen ionization fraction - the spherical-average, mean-free-path and friends-of-friends methods - and apply them to simulated 21-cm data cubes. Our simulated data cubes have the (sensitivity-dictated) resolution expected for the SKA1-Low reionization experiment and we study the impact of both the light-cone and redshift space distortion effects. To identify ionized regions in the 21-cm data we introduce a new, self-adjusting thresholding approach based on the K-Means algorithm. We find that the fraction of ionized cells identified in this way consistently falls below the mean volume-averaged ionized fraction. From a comparison of the three bubble size methods, we conclude that all three methods are useful, but that the mean-free-path method performs best in terms of tracking the progress of reionization and separating different reionization scenarios. The light-cone effect is found to affect data spanning more than about 10~MHz in frequency ($\Delta z\sim0.5$). We find that redshift space distortions only marginally affect the bubble size distributions.

Three-dimensional radiative transfer simulations of the epoch of reionization can produce realistic results, but are computationally expensive. On the other hand, simulations relying on one-dimensional radiative transfer solutions are faster but limited in accuracy due to their more approximate nature. Here, we compare the performance of the reionization simulation codes grizzly and C2-ray which use 1D and 3D radiative transfer schemes respectively. The comparison is performed using the same cosmological density fields, halo catalogues and source properties. We find that the ionization maps, as well as the 21-cm signal maps from these two simulations are very similar even for complex scenarios which include thermal feedback on low mass halos. The comparison between the schemes in terms of the statistical quantities such as the power spectrum of the brightness temperature fluctuation agree with each other within 10% error throughout the entire reionization history. grizzly seems to perform slightly better than the semi-numerical approaches considered in Majumdar et al. (2014) which are based on the excursion set principle. We argue that grizzly can be efficiently used for exploring parameter space, establishing observations strategies and estimating parameters from 21-cm observations.

We introduce a novel technique, called "granulometry", to characterize and recover the mean size and the size distribution of HII regions from 21-cm tomography. The technique is easy to implement, but places the previously not very well defined concept of morphology on a firm mathematical foundation. The size distribution of the cold spots in 21-cm tomography can be used as a direct tracer of the underlying probability distribution of HII region sizes. We explore the capability of the method using large-scale reionization simulations and mock observational data cubes while considering capabilities of SKA1-low and a future extension to SKA2. We show that the technique allows the recovery of the HII region size distribution with a moderate signal-to-noise ratio from wide-field imaging ($\rm SNR\lesssim3$), for which the statistical uncertainty is sample variance dominated. We address the observational requirements on the angular resolution, the field-of-view, and the thermal noise limit for a successful measurement. To achieve a full scientific return from 21-cm tomography and to exploit a synergy with 21-cm power spectra, we suggest an observing strategy using wide-field imaging (several tens of square degrees) by an interferometric mosaicking/multi-beam observation with additional intermediate baselines (~2-4 km).

Upcoming observations of the 21-cm signal from the Epoch of Reionization will soon provide the first direct detection of this era. This signal is influenced by many astrophysical effects, including long range X-ray heating of the intergalactic gas. During the preceding Cosmic Dawn era the impact of this heating on the 21-cm signal is particularly prominent, especially before spin temperature saturation. We present the largest-volume (349\,Mpc comoving=244~$h^{-1}$Mpc) full numerical radiative transfer simulations to date of this epoch that include the effects of helium and multi-frequency heating, both with and without X-ray sources. We show that X-ray sources contribute significantly to early heating of the neutral intergalactic medium and, hence, to the corresponding 21-cm signal. The inclusion of hard, energetic radiation yields an earlier, extended transition from absorption to emission compared to the stellar-only case. The presence of X-ray sources decreases the absolute value of the mean 21-cm differential brightness temperature. These hard sources also significantly increase the 21-cm fluctuations compared the common assumption of temperature saturation. The 21-cm differential brightness temperature power spectrum is initially boosted on large scales, before decreasing on all scales. Compared to the case of the cold, unheated intergalactic medium, the signal has lower rms fluctuations and increased non-Gaussianity, as measured by the skewness and kurtosis of the 21-cm probability distribution functions. Images of the 21-cm signal with resolution around 11~arcmin still show fluctuations well above the expected noise for deep integrations with the SKA1-Low, indicating that direct imaging of the X-ray heating epoch could be feasible.

We present the first limits on the Epoch of Reionization (EoR) 21-cm HI power spectra, in the redshift range $z=7.9-10.6$, using the Low-Frequency Array (LOFAR) High-Band Antenna (HBA). In total 13\,h of data were used from observations centred on the North Celestial Pole (NCP). After subtraction of the sky model and the noise bias, we detect a non-zero $\Delta^2_{\rm I} = (56 \pm 13 {\rm mK})^2$ (1-$\sigma$) excess variance and a best 2-$\sigma$ upper limit of $\Delta^2_{\rm 21} < (79.6 {\rm mK})^2$ at $k=0.053$$h$cMpc$^{-1}$ in the range $z=$9.6-10.6. The excess variance decreases when optimizing the smoothness of the direction- and frequency-dependent gain calibration, and with increasing the completeness of the sky model. It is likely caused by (i) residual side-lobe noise on calibration baselines, (ii) leverage due to non-linear effects, (iii) noise and ionosphere-induced gain errors, or a combination thereof. Further analyses of the excess variance will be discussed in forthcoming publications.

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 $\Delta 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.

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.

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.

A proposed method for dealing with foreground emission in upcoming 21-cm observations from the epoch of reionization is to limit observations to an uncontaminated window in Fourier space. Foreground emission can be avoided in this way, since it is limited to a wedge-shaped region in $k_{\parallel}, k_{\perp}$ space. However, the power spectrum is anisotropic owing to redshift-space distortions from peculiar velocities. Consequently, the 21-cm power spectrum measured in the foreground avoidance window---which samples only a limited range of angles close to the line-of-sight direction---differs from the full spherically-averaged power spectrum which requires an average over \emph{all} angles. In this paper, we calculate the magnitude of this "wedge bias" for the first time. We find that the bias is strongest at high redshifts, where measurements using foreground avoidance will over-estimate the power spectrum by around 100 per cent, possibly obscuring the distinctive rise and fall signature that is anticipated for the spherically-averaged 21-cm power spectrum. In the later stages of reionization, the bias becomes negative, and smaller in magnitude ($\lesssim 20$ per cent). The effect shows only a weak dependence on spatial scale and reionization topology.

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.

This study aims to characterize linear polarization structures in LOFAR observations of the interstellar medium (ISM) in the 3C196 field, one of the primary fields of the LOFAR-Epoch of Reionization key science project. We have used the high band antennas (HBA) of LOFAR to image this region and Rotation Measure (RM) synthesis to unravel the distribution of polarized structures in Faraday depth. The brightness temperature of the detected Galactic emission is $5-15~{\rm K}$ in polarized intensity and covers the range from -3 to +8 ${\rm rad~m^{-2}}$ in Faraday depth. The most interesting morphological feature is a strikingly straight filament at a Faraday depth of $+0.5~{\rm rad~m^{-2}}$ running from north to south, right through the centre of the field and parallel to the Galactic plane. There is also an interesting system of linear depolarization canals conspicuous in an image showing the peaks of Faraday spectra. We used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. For the first time, we see some common morphology in the RM cubes made at 150 and 350~{\rm MHz}. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies and previous LOFAR observations. Based on our results, we determined physical parameters of the ISM and proposed a simple model that may explain the observed distribution of the intervening magneto-ionic medium. The mean line-of-sight magnetic field component, $B_\parallel$, is determined to be $0.3\pm0.1~{\rm \mu G}$ and its spatial variation across the 3C196 field is $0.1~{\rm \mu G}$. The filamentary structure is probably an ionized filament in the ISM, located somewhere within the Local Bubble. This filamentary structure shows an excess in thermal electron density ($n_e B_\parallel>6.2~{\rm cm^{-3}\mu G}$) compared to its surroundings.

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.

We investigate the feasibility of detecting 21cm absorption features in the afterglow spectra of high redshift long Gamma Ray Bursts (GRBs). This is done employing simulations of cosmic reionization, together with the instrumental characteristics of the LOw Frequency ARray (LOFAR). We find that absorption features could be marginally (with a S/N larger than a few) detected by LOFAR at z>7 if the GRB originated from PopIII stars, while the detection would be easier if the noise were reduced by one order of magnitude, i.e. similar to what is expected for the first phase of the Square Kilometer Array (SKA1-low). On the other hand, more standard GRBs are too dim to be detected even with ten times the sensitivity of SKA1-low, and only in the most optimistic case can a S/N larger than a few be reached at z>9.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

Detection of the 21-cm signal coming from the epoch of reionization (EoR) is challenging especially because, even after removing the foregrounds, the residual Stokes $I$ maps contain leakage from polarized emission that can mimic the signal. Here, we discuss the instrumental polarization of LOFAR and present realistic simulations of the leakages between Stokes parameters. From the LOFAR observations of polarized emission in the 3C196 field, we have quantified the level of polarization leakage caused by the nominal model beam of LOFAR, and compared it with the EoR signal using power spectrum analysis. We found that at 134--166 MHz, within the central 4$^\circ$ of the field the $(Q,U)\rightarrow I$ leakage power is lower than the EoR signal at $k<0.3$ Mpc$^{-1}$. The leakage was found to be localized around a Faraday depth of 0, and the rms of the leakage as a fraction of the rms of the polarized emission was shown to vary between 0.2-0.3\%, both of which could be utilized in the removal of leakage. Moreover, we could define an `EoR window' in terms of the polarization leakage in the cylindrical power spectrum above the PSF-induced wedge and below $k_\parallel\sim 0.5$ Mpc$^{-1}$, and the window extended up to $k_\parallel\sim 1$ Mpc$^{-1}$ at all $k_\perp$ when 70\% of the leakage had been removed. These LOFAR results show that even a modest polarimetric calibration over a field of view of $\lesssim 4^\circ$ in the future arrays like SKA will ensure that the polarization leakage remains well below the expected EoR signal at the scales of 0.02-1 Mpc$^{-1}$.

We provide an overview of 21cm tomography of the Cosmic Dawn and Epoch of Reionization as possible with SKA-Low. We show why tomography is essential for studying CD/EoR and present the scales which can be imaged at different frequencies for the different phases of SKA- Low. Next we discuss the different ways in which tomographic data can be analyzed. We end with an overview of science questions which can only be answered by tomography, ranging from the characterization of individual objects to understanding the global processes shaping the Universe during the CD/EoR

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.

The Square Kilometre Array (SKA) will offer an unprecedented view onto the early Universe, using interferometric observations of the redshifted 21cm line. The 21cm line probes the thermal and ionization state of the cosmic gas, which is governed by the birth and evolution of the first structures in our Universe. Here we show how the evolution of the 21cm signal will allow us to study when the first generations of galaxies appeared, what were their properties, and what was the structure of the intergalactic medium. We highlight qualitative trends which will offer robust insights into the early Universe.

We investigate the scale dependence of fluctuations inside a realistic model of an evolving turbulent HII region and to what extent these may be studied observationally. We find that the multiple scales of energy injection from champagne flows and the photoionization of clumps and filaments leads to a flatter spectrum of fluctuations than would be expected from top-down turbulence driven at the largest scales. The traditional structure function approach to the observational study of velocity fluctuations is shown to be incapable of reliably determining the velocity power spectrum of our simulation. We find that a more promising approach is the Velocity Channel Analysis technique of Lazarian & Pogosyan (2000), which, despite being intrinsically limited by thermal broadening, can successfully recover the logarithmic slope of the velocity power spectrum to a precision of +-0.1 from high resolution optical emission line spectroscopy.

We study the prospects for constraining the ionized fraction of the intergalactic medium (IGM) at $z>6$ with the next generation of large Ly$\alpha$ emitter surveys. We make predictions for the upcoming Subaru Hyper Suprime-Cam (HSC) Ly$\alpha$ survey and a hypothetical spectroscopic survey performed with the James Webb Space Telescope (JWST). Considering various scenarios where the observed evolution of the Ly$\alpha$ luminosity function of Ly$\alpha$ 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$\alpha$ 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.

We present radio observations of the Moon between $35$ and $80$ MHz to demonstrate a novel technique of interferometrically measuring large-scale diffuse emission extending far beyond the primary beam (global signal) for the first time. In particular, we show that (i) the Moon appears as a negative-flux source at frequencies $35<\nu<80$ MHz since it is `colder' than the diffuse Galactic background it occults, (ii) using the (negative) flux of the lunar disc, we can reconstruct the spectrum of the diffuse Galactic emission with the lunar thermal emission as a reference, and (iii) that reflected RFI (radio-frequency interference) is concentrated at the center of the lunar disc due to specular nature of reflection, and can be independently measured. Our RFI measurements show that (i) Moon-based Cosmic Dawn experiments must design for an Earth-isolation of better than $80$ dB to achieve an RFI temperature $<1$ mK, (ii) Moon-reflected RFI contributes to a dipole temperature less than $20$ mK for Earth-based Cosmic Dawn experiments, (iii) man-made satellite-reflected RFI temperature exceeds $20$ mK if the aggregate cross section of visible satellites exceeds $80$ m$^2$ at $800$ km height, or $5$ m$^2$ at $400$ km height. Currently, our diffuse background spectrum is limited by sidelobe confusion on short baselines (10-15% level). Further refinement of our technique may yield constraints on the redshifted global $21$-cm signal from Cosmic Dawn ($40>z>12$) and the Epoch of Reionization ($12>z>5$).

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.

This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%.

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 ($\Delta z$) of reionization can be constrained assuming a parametrization. We use an EoR simulation of $z_r = 7.68$ and $\Delta 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 $\Delta 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.

A novel method to characterize the topology of the early-universe intergalactic medium during the epoch of cosmic reionization is presented. The 21-cm radiation background from high redshift is analyzed through the calculation of the 2-dimensional (2D) genus. The radiative transfer of hydrogen-ionizing photons and ionization-rate equations are calculated in a suite of numerical simulations under various input parameters. The 2D genus is calculated from the mock 21-cm images of the high-redshift Universe. We construct the 2D genus curve by varying the threshold differential brightness temperature, and compare this to the 2D genus curve of the underlying density field. We find that (1) the 2D genus curve reflects the evolutionary track of cosmic reionization and (2) the 2D genus curve can discriminate between certain reionization scenarios and thus indirectly probe the properties of radiation-sources. Choosing the right beam shape of a radio antenna is crucial for this analysis. To this end, the Square Kilometre Array (SKA) is found to be a suitable apparatus for this analysis in terms of sensitivity, even though some deterioration of the data for this purpose is unavoidable under the planned size of the antenna core.

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.

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.

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 $\Delta 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)

Faint undetected sources of radio-frequency interference (RFI) might become visible in long radio observations when they are consistently present over time. Thereby, they might obstruct the detection of the weak astronomical signals of interest. This issue is especially important for Epoch of Reionisation (EoR) projects that try to detect the faint redshifted HI signals from the time of the earliest structures in the Universe. We explore the RFI situation at 30-163 MHz by studying brightness histograms of visibility data observed with LOFAR, similar to radio-source-count analyses that are used in cosmology. An empirical RFI distribution model is derived that allows the simulation of RFI in radio observations. The brightness histograms show an RFI distribution that follows a power-law distribution with an estimated exponent around -1.5. With several assumptions, this can be explained with a uniform distribution of terrestrial radio sources whose radiation follows existing propagation models. Extrapolation of the power law implies that the current LOFAR EoR observations should be severely RFI limited if the strength of RFI sources remains strong after time integration. This is in contrast with actual observations, which almost reach the thermal noise and are thought not to be limited by RFI. Therefore, we conclude that it is unlikely that there are undetected RFI sources that will become visible in long observations. Consequently, there is no indication that RFI will prevent an EoR detection with LOFAR.

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.

Using a combination of N-body simulations, semi-analytic models and radiative transfer calculations, we have estimated the theoretical cross power spectrum between galaxies and the 21cm emission from neutral hydrogen during the epoch of reionization. In accordance with previous studies, we find that the 21cm emission is initially correlated with halos on large scales (> 30 Mpc), anti-correlated on intermediate (~ 5 Mpc), and uncorrelated on small (< 3 Mpc) scales. This picture quickly changes as reionization proceeds and the two fields become anti-correlated on large scales. The normalization of the cross power spectrum can be used to set constraints on the average neutral fraction in the intergalactic medium and its shape can be a tool to study the topology of reionization. When we apply a drop-out technique to select galaxies and add to the 21cm signal the noise expected from the LOFAR telescope, we find that while the normalization of the cross power spectrum remains a useful tool for probing reionization, its shape becomes too noisy to be informative. On the other hand, for a Lyalpha Emitter (LAE) survey both the normalization and the shape of the cross power spectrum are suitable probes of reionization. A closer look at a specific planned LAE observing program using Subaru Hyper-Suprime Cam reveals concerns about the strength of the 21cm signal at the planned redshifts. If the ionized fraction at z ~ 7 is lower that the one estimated here, then using the cross power spectrum may be a useful exercise given that at higher redshifts and neutral fractions it is able to distinguish between two toy models with different topologies.

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.

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).

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.

The aim of the LOFAR Epoch of Reionization (EoR) project is to detect the spectral fluctuations of the redshifted HI 21cm signal. This signal is weaker by several orders of magnitude than the astrophysical foreground signals and hence, in order to achieve this, very long integrations, accurate calibration for stations and ionosphere and reliable foreground removal are essential. One of the prospective observing windows for the LOFAR EoR project will be centered at the North Celestial Pole (NCP). We present results from observations of the NCP window using the LOFAR highband antenna (HBA) array in the frequency range 115 MHz to 163 MHz. The data were obtained in April 2011 during the commissioning phase of LOFAR. We used baselines up to about 30 km. With about 3 nights, of 6 hours each, effective integration we have achieved a noise level of about 100 microJy/PSF in the NCP window. Close to the NCP, the noise level increases to about 180 microJy/PSF, mainly due to additional contamination from unsubtracted nearby sources. We estimate that in our best night, we have reached a noise level only a factor of 1.4 above the thermal limit set by the noise from our Galaxy and the receivers. Our continuum images are several times deeper than have been achieved previously using the WSRT and GMRT arrays. We derive an analytical explanation for the excess noise that we believe to be mainly due to sources at large angular separation from the NCP.

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.

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.

We discuss the feasibility of the detection of the 21cm forest in the diffuse IGM with the radio telescope LOFAR. The optical depth to the 21cm line has been derived using simulations of reionization which include detailed radiative transfer of ionizing photons. We find that the spectra from reionization models with similar total comoving hydrogen ionizing emissivity but different frequency distribution look remarkably similar. Thus, unless the reionization histories are very different from each other (e.g. a predominance of UV vs. x-ray heating) we do not expect to distinguish them by means of observations of the 21cm forest. Because the presence of a strong x-ray background would make the detection of 21cm line absorption impossible, the lack of absorption could be used as a probe of the presence/intensity of the x-ray background and the thermal history of the universe. Along a random line of sight LOFAR could detect a global suppression of the spectrum from z>12, when the IGM is still mostly neutral and cold, in contrast with the more well-defined, albeit broad, absorption features visible at lower redshift. Sharp, strong absorption features associated with rare, high density pockets of gas could be detected also at z~7 along preferential lines of sight.

Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills.

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, \sigma = 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.

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.

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.

The first generation of redshifted 21 cm detection experiments, carried out with arrays like LOFAR, MWA and GMRT, will have a very low signal-to-noise ratio per resolution element (\sim 0.2). In addition, whereas the variance of the cosmological signal decreases on scales larger than the typical size of ionization bubbles, the variance of the formidable galactic foregrounds increases, making it hard to disentangle the two on such large scales. The poor sensitivity on small scales on the one hand, and the foregrounds effect on large scales on the other hand, make direct imaging of the Epoch of Reionization of the Universe very difficult, and detection of the signal therefore is expected to be statistical.Despite these hurdles, in this paper we argue that for many reionization scenarios low resolution images could be obtained from the expected data. This is because at the later stages of the process one still finds very large pockets of neutral regions in the IGM, reflecting the clustering of the large-scale structure, which stays strong up to scales of \sim 120 comoving Mpc/h (\sim 1 degree). The coherence of the emission on those scales allows us to reach sufficient S/N (\sim 3) so as to obtain reionization 21 cm images. Such images will be extremely valuable for answering many cosmological questions but above all they will be a very powerful tool to test our control of the systematics in the data. The existence of this typical scale (\sim 120 comoving Mpc/h) also argues for designing future EoR experiments, e.g., with SKA, with a field of view of at least 4 degree.

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)

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.

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]

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.

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)

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.

We present the results of radiation-magnetohydrodynamic simulations of the expansion of HII regions and surrounding photodissociation regions in turbulent, magnetised, molecular clouds on scales of up to 4 parsecs, including the effects of ionising and non-ionising ultraviolet radiation and x rays from young star clusters. We find that HII region expansion reduces the disordered component of the B field, imposing a large-scale order on the field around its border. The field in the neutral gas lies along the ionisation front, while the field in the ionised gas tends to be perpendicular to this. The highest pressure compressed neutral/molecular gas is driven towards approximate equipartition between thermal/magnetic/turbulent energy densities, whereas lower pressure neutral/molecular gas divides into quiescent, magnetically dominated regions, and, on the other hand, turbulent, demagnetised regions. The ionised gas shows approximate thermal/turbulent equipartition, but with magnetic energy densities 1 to 3 orders of magnitude lower. A high velocity dispersion (approx 8 km/s) is maintained in the ionised gas throughout our simulations, despite the mean expansion velocity being significantly lower. The B field does not significantly brake the HII region expansion on the length and timescales accessible to our simulations, but it does tend to suppress the small-scale fragmentation and radiation-driven implosion of neutral/molecular gas that forms globules and pillars at the edge of the HII region. However, the relative luminosity of ionising and non-ionising radiation has a much larger influence than the presence or absence of the B field. When the radiation field is relatively soft (as in the case of a lower mass cluster, with earliest spectral type of B0.5), then fragmentation is less vigorous and a thick, relatively smooth PDR forms. Movies available at http://youtube.com/user/divBequals0

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.

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.

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.

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.

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.

An obstacle to the detection of redshifted 21cm emission from the epoch of reionization (EoR) is the presence of foregrounds which exceed the cosmological signal in intensity by orders of magnitude. We argue that in principle it would be better to fit the foregrounds non-parametrically - allowing the data to determine their shape - rather than selecting some functional form in advance and then fitting its parameters. Non-parametric fits often suffer from other problems, however. We discuss these before suggesting a non-parametric method, Wp smoothing, which seems to avoid some of them. After outlining the principles of Wp smoothing we describe an algorithm used to implement it. We then apply Wp smoothing to a synthetic data cube for the LOFAR EoR experiment. The performance of Wp smoothing, measured by the extent to which it is able to recover the variance of the cosmological signal and to which it avoids leakage of power from the foregrounds, is compared to that of a parametric fit, and to another non-parametric method (smoothing splines). We find that Wp smoothing is superior to smoothing splines for our application, and is competitive with parametric methods even though in the latter case we may choose the functional form of the fit with advance knowledge of the simulated foregrounds. Finally, we discuss how the quality of the fit is affected by the frequency resolution and range, by the characteristics of the cosmological signal and by edge effects.

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.

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.

Detecting redshifted 21cm emission from neutral hydrogen in the early Universe promises to give direct constraints on the epoch of reionization (EoR). It will, though, be very challenging to extract the cosmological signal (CS) from foregrounds and noise which are orders of magnitude larger. Fortunately, the signal has some characteristics which differentiate it from the foregrounds and noise, and we suggest that using the correct statistics may tease out signatures of reionization. We generate mock datacubes simulating the output of the Low Frequency Array (LOFAR) EoR experiment. These cubes combine realistic models for Galactic and extragalactic foregrounds and the noise with three different simulations of the CS. We fit out the foregrounds, which are smooth in the frequency direction, to produce residual images in each frequency band. We denoise these images and study the skewness of the one-point distribution in the images as a function of frequency. We find that, under sufficiently optimistic assumptions, we can recover the main features of the redshift evolution of the skewness in the 21cm signal. We argue that some of these features - such as a dip at the onset of reionization, followed by a rise towards its later stages - may be generic, and give us a promising route to a statistical detection of 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.

Observations of high-redshift Ly-alpha sources are a major tool for studying the high-redshift Universe. We discuss the effect of the reionizing intergalactic medium on the observability of Ly-alpha sources based on large simulations of early structure formation with radiative transfer. This takes into account self-consistently the reionization history, density, velocity and ionization structures and nonlinear source clustering. We find that all fields are highly anisotropic and as a consequence there are very large variations in opacity among the different lines-of-sight. The velocity effects, from both infall and source peculiar velocity are most important for the luminous sources, affecting the line profile and depressing the bright end of the luminosity function. The line profiles are generally asymmetric and the line centers of the luminous sources are always absorbed due to the high density of the local IGM. For both luminous and average sources the damping wing effects are of similar magnitude and remain significant until fairly late. The ionizing flux in the ionized patch surrounding a high density peak is generally strongly dominated, particularly at late times, by the cluster of faint sources, rather than the central massive galaxy. The IGM absorption does not change appreciably the correlation function of sources at high redshift. Our derived luminosity function assuming constant mass-to-light ratio provides an excellent match to the shape of the observed luminosity function at z=6.6 with faint-end slope of alpha=-1.5. The resulting mass-to-light ratio implies that the majority of sources responsible for reionization are too faint to be observed by the current surveys. (abridged)

We present an efficient method to generate large simulations of the Epoch of Reionization (EoR) without the need for a full 3-dimensional radiative transfer code. Large dark-matter-only simulations are post-processed to produce maps of the redshifted 21cm emission from neutral hydrogen. Dark matter haloes are embedded with sources of radiation whose properties are either based on semi-analytical prescriptions or derived from hydrodynamical simulations. These sources could either be stars or power-law sources with varying spectral indices. Assuming spherical symmetry, ionized bubbles are created around these sources, whose radial ionized fraction and temperature profiles are derived from a catalogue of 1-D radiative transfer experiments. In case of overlap of these spheres, photons are conserved by redistributing them around the connected ionized regions corresponding to the spheres. The efficiency with which these maps are created allows us to span the large parameter space typically encountered in reionization simulations. We compare our results with other, more accurate, 3-D radiative transfer simulations and find excellent agreement for the redshifts and the spatial scales of interest to upcoming 21cm experiments. We generate a contiguous observational cube spanning redshift 6 to 12 and use these simulations to study the differences in the reionization histories between stars and quasars. Finally, the signal is convolved with the LOFAR beam response and its effects are analyzed and quantified. Statistics performed on this mock data set shed light on possible observational strategies for LOFAR.

The first, self-consistent calculations of the cosmological H_2 dissociating UV background produced during the epoch of reionization (EOR) by the sources of reionization are presented. 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 solar mass, 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 background in the Lyman-Werner band of H_2 between 11 and 13.6 eV.

The redshifted 21-cm line of distant neutral H atoms provides a probe of the cosmic ``dark ages'' and the epoch of reionization (``EOR'') which ended them. The radio continuum produced by this redshifted line can be seen in absorption or emission against the CMB at meterwaves, yielding information about the thermal and ionization history of the universe and the primordial density perturbation spectrum that led to galaxy and large-scale structure formation. Observing this 21-cm background is a great challenge. A new generation of low-frequency radio arrays is currently under development to search for this background. Accurate theoretical predictions of the spectrum and anisotropy of this background, necessary to guide and interpret future observations, are also quite challenging. It is necessary to model the inhomogeneous reionization of the intergalactic medium and determine the spin temperature of the 21-cm transition and its variations in time and space as it decouples from the temperature of the CMB. Here, we focus on just a few of the predictions for the 21-cm background from the EOR, based on our newest, large-scale simulations of patchy reionization. These simulations are the first with enough N-body particles (from 5 to 29 billion) and radiative transfer rays to resolve the formation of and trace the ionizing radiation from each of the millions of dwarf galaxies believed responsible for reionization, down to 10^8 M_solar, in a cubic volume large enough (90 and 163 comoving Mpc on a side) to make meaningful statistical predictions of the fluctuating 21-cm background. (abridged)

The Cosmic Dark Ages and the Epoch of Reionization constitute a crucial missing link in our understanding of the evolution of the intergalactic medium and the formation and evolution of galaxies. Due to the complex nature of this global process it is best studied through large-scale numerical simulations. This presents considerable computational challenges. The dominant contributors of ionizing radiation were dwarf galaxies. These tiny galaxies must be resolved in very large cosmological volumes in order to derive their clustering properties and the corresponding observational signatures correctly, which makes this one of the most challenging problems of numerical cosmology. We have recently performed the largest and most detailed simulations of the formation of early cosmological large-scale structures and their radiative feedback leading to cosmic reionization. This was achieved by running extremely large (up to 29 billion-particle) N-body simulations of the formation of the Cosmic Web, with enough particles and sufficient force resolution to resolve all the galactic halos with total masses larger than 10^8 Solar masses in computational volumes of up to (163 Mpc)^3. These results were then post-processed by propagating the ionizing radiation from all sources by using fast and accurate ray-tracing radiative transfer method. Both of our codes are parallelized using a combination of MPI and OpenMP and to this date have been run efficiently on up to 2048 cores (N-body) and up to 10000 cores (radiative transfer) on the newly-deployed Sun Constellation Linux Cluster at the Texas Advanced Computing Center. In this paper we describe our codes, parallelization strategies, scaling and some preliminary scientific results. (abridged)

We present a numerical study of rapid, so called type III migration for Jupiter-sized planets embedded in a protoplanetary disc. We limit ourselves to the case of outward migration, and study in detail its evolution and physics, concentrating on the structure of the co-rotation and circumplanetary regions, and processes for stopping migration. We also consider the dependence of the migration behaviour on several key parameters. We perform this study using global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. We find that the outward directed type III migration can be started if the initial conditions support $Z > 1$, that corresponds to initial value $M_\rmn{\Delta} \ga 1.5$. Unlike the inward directed migration, in the outward migration the migration rate increases due to the growing of the volume of the co-orbital region. We find the migration to be strongly dependent on the rate of the mass accumulation in the circumplanetary disc, leading to two possible regimes of migration, fast and slow. The structure of the co-orbital region and the stopping mechanism differ between these two regimes.

Observations of extrasolar planets reveal that planets can be found in close binary systems, where the semi-major axis of the binary orbit is less than 20 AU. The existence of these planets challenges planet formation theory, because the strong gravitational perturbations due to the companion increase encounter velocities between planetesimals and make it difficult for them to grow through accreting collisions. We study planetesimal encounter velocities in binary systems, where the planetesimals are embedded in a circumprimary gas disc that is allowed to evolve under influence of the gravitational perturbations of the companion star. We find that the encounter velocities between planetesimals of different size strongly depend on the gas disc eccentricity. In all cases studied, inclusion of the full gas dynamics increases the encounter velocity compared to the case of a static, circular gas disc. Full numerical parameter exploration is still impossible, but we derive analytical formulae to estimate encounter velocities between bodies of different sizes given the gas disc eccentricity. The gas dynamical evolution of a protoplanetary disc in a binary system tends to make planetesimal accretion even more difficult than in a static, axisymmetric gas disc.

Direct detection of the Dark Ages and the Epoch of Reionization (EOR) is among the main scientific objectives of all current and future low-frequency radio facilities. In this paper we summarize and discuss recent results, based on state-of-the-art numerical simulations, regarding the fundamental EOR properties and its observability with current and future radio arrays, like the Giant Metrewave Radio Telescope (GMRT), the Low Frequency Array (LOFAR), the 21-CM Array (21CMA), the Murchison Widefield Array (MWA) and the Square Kilometre Array (SKA). Results show that the optimal observational frequencies for statistical detection are 140-160 MHz. The signals are strongly non-Gaussian at late times. The correlation widths between 21-cm maps at neighbouring frequencies are short, of order 300-800 kHz, which should help with the cleaning of the strong foregrounds. Direct comparison of the resolutions and expected sensitivities of GMRT and MWA indicate that their optimal sensitivity ranges are similar, at scales k~0.2-0.4 h/Mpc, however, all else being equal the former should require shorter integration times due to its significantly larger collecting area.

We study the interaction of a low-mass planet with a protoplanetary disk with a realistic treatment of the energy balance by doing radiation-hydrodynamical simulations. We look at accretion and migration rates and compare them to isothermal studies. We used a three-dimensional version of the hydrodynamical method RODEO, together with radiative transport in the flux-limited diffusion approach. The accretion rate, as well as the torque on the planet, depend critically on the ability of the disk to cool efficiently. For densities appropriate to 5 AU in the solar nebula, the accretion rate drops by more than an order of magnitude compared to isothermal models, while at the same time the torque on the planet is positive, indicating outward migration. It is necessary to lower the density by a factor of 2 to recover inward migration and more than 2 orders of magnitude to recover the usual Type I migration. The torque appears to be proportional to the radial entropy gradient in the unperturbed disk. These findings are critical for the survival of protoplanets, and they should ultimately find their way into population synthesis models.

Recently, we have presented the first large-scale radiative transfer simulations of reionization. Here we present new simulations which extend the source halo mass range downward to 10^8M_solar, to capture the full range of halo masses thought to be primarily responsible for reionization by their star formation following atomic hydrogen radiative cooling and gravitational collapse. Haloes below about 10^9M_solar, however, are subject to Jeans-mass filtering in the ionized regions, which suppresses their baryonic content and their ability to release ionizing radiation. By including these smaller-mass haloes but accounting for their suppression, too, we find that reionization is ``self-regulating,'' as follows. As the mean ionized fraction rises, so does the fraction of the volume within which suppression occurs. Hence, the degree of suppression is related to the mean ionized fraction. Since low-mass haloes with high emissivity achieve a given mean ionized fraction earlier than do those with low efficiency, Jeans-mass filtering compensates for the difference in the emissivity of the suppressible haloes in these two cases. As a result, in the presence of lower-mass source haloes, reionization begins earlier, but the later stages of reionization and the time of overlap are dictated by the efficiency of the higher-mass haloes, independent of the efficiency of the suppressible, lower-mass haloes. Reionization histories consistent with current observational constraints are shown to be achievable with standard stellar sources in haloes above 10^8M_solar. Neither minihalos nor exotic sources are required, and the phenomenon of ``double reionization'' previously suggested does not occur. (abridged)

A number of large current experiments aim to detect the signatures of the Cosmic Reionization at redshifts z>6. Their success depends crucially on understanding the character of the reionization process and its observable consequences and designing the best strategies to use. We use large-scale simulations of cosmic reionization to evaluate the reionization signatures at redshifted 21-cm and small-scale CMB anisotropies in the best current model for the background universe, with fundamental cosmological parameters given by WMAP 3-year results (WMAP3). We find that the optimal frequency range for observing the ``global step'' of the 21-cm emission is 120-150 MHz, while statistical studies should aim at 140-160 MHz, observable by GMRT. Some strongly-nongaussian brightness features should be detectable at frequencies up to ~190 MHz. In terms of sensitivity-signal trade-off relatively low resolutions, corresponding to beams of at least a few arcminutes, are preferable. The CMB anisotropy signal from the kinetic Sunyaev-Zel'dovich effect from reionized patches peaks at tens of muK at arcminute scales and has an rms of ~1 muK, and should be observable by the Atacama Cosmology Telescope and the South Pole Telescope. We discuss the various observational issues and the uncertainties involved, mostly related to the poorly-known reionization parameters and, to a lesser extend, to the uncertainties in the background cosmology.

Recently the numerical simulations of the process of reionization of the universe at z>6 have made a qualitative leap forward, reaching sufficient sizes and dynamic range to determine the characteristic scales of this process. This allowed making the first realistic predictions for a variety of observational signatures. We discuss recent results from large-scale radiative transfer and structure formation simulations on the observability of high-redshift Ly-alpha sources. We also briefly discuss the dependence of the characteristic scales and topology of the ionized and neutral patches on the reionization parameters.

We investigate the fast (type III) migration regime of high-mass protoplanets orbiting in protoplanetary disks. This type of migration is dominated by corotational torques. We study the details of flow structure in the planet's vicinity, the dependence of migration rate on the adopted disc model, and the numerical convergence of models (independence of certain numerical parameters such as gravitational softening). We use two-dimensional hydrodynamical simulations with adaptive mesh refinement,based on the FLASH code with improved time-stepping scheme. We perform global disk simulations with sufficient resolution close to the planet, which is allowed to freely move throughout the grid. We employ a new type of equation of state in which the gas temperature depends on both the distance to the star and planet, and a simplified correction for self-gravity of the circumplanetary gas. We find that the migration rate in the type III migration regime depends strongly on the gas dynamics inside the Hill sphere (Roche lobe of the planet) which, in turn, is sensitive to the aspect ratio of the circumplanetary disc. Furthermore, corrections due to the gas self-gravity are necessary to reduce numerical artifacts that act against rapid planet migration. Reliable numerical studies of Type III migration thus require consideration of both the thermal andthe self-gravity corrections, as well as a sufficient spatial resolution and the calculation of disk-planet attraction both inside and outside the Hill sphere. With this proviso, we find Type III migration to be a robust mode of migration, astrophysically promising because of a speed much faster than in the previously studied modes of migration.

We present a numerical study of rapid, so called type III migration for Jupitersized planets embedded in a protoplanetary disc. We limit ourselves to the case of inward migration, and study in detail its evolution and physics, concentrating on the structure of the corotation and circumplanetary regions, and processes for stopping migration. We also consider the dependence of the migration behaviour on several key parameters. We perform this study using the results of global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. The initial conditions are chosen to satisfy the condition for rapid inward migration. We find that type III migration can be divided into two regimes, fast and slow. The structure of the coorbital region, mass accumulation rate, and migration behaviour differ between these two regimes. All our simulations show a transition from the fast to the slow regime, ending type III migration well before reaching the star. The stopping radius is found to be larger for more massive planets and less massive discs. A sharp density drop is also found to be an efficient stopping mechanism. In the fast migration limit the migration rate and induced eccentricity are lower for less massive discs, but almost do not depend on planet mass. Eccentricity is damped on the migration time scale.

In recent years there has been considerable progress in our understanding of the nature and properties of the reionization process. In particular, the numerical simulations of this epoch have made a qualitative leap forward, reaching sufficiently large scales to derive the characteristic scales of the reionization process and thus allowing for realistic observational predictions. Our group has recently performed the first such large-scale radiative transfer simulations of reionization, run on top of state-of-the-art simulations of early structure formation. This allowed us to make the first realistic observational predictions about the Epoch of Reionization based on detailed radiative transfer and structure formation simulations. We discuss the basic features of reionization derived from our simulations and some recent results on the observational implications for the high-redshift Ly-alpha sources.

The reionization of the Universe has profound effects on the way galaxies form and on their observed properties at later times. Of particular importance is the relative timing of the reionization history of a region and its halo assembly history, which can affect the nature of the first stars formed in that region, the properties and radial distribution of its stellar halo, globular cluster population and its satellite galaxies. We distinguish two basic cases for the reionization of a halo - internal reionization, whereby the stars forming in situ reionize their host galaxy, and external reionization, whereby the progenitor of a galaxy is reionized by external radiation before its own stars are able to form in sufficient numbers. We use a set of large-scale radiative transfer and structure formation simulations, based on cosmologies derived from both WMAP 1-year and WMAP 3-year data, to evaluate the mean reionization redshifts and the probability of internal/external reionization for Local Group-like systems, galaxies in the field and central cD galaxies in clusters. We find that these probabilities are strongly dependent on the underlying cosmology and the efficiency of photon production, but also on the halo mass. There is a rapid transition between predominantly external and predominantly internal reionization at a mass scale of 1.0e12 Msun (corresponding roughly to L*galaxies), with haloes less massive than this being reionized preferentially from distant sources. We provide a fit for the reionization redshift as a function of halo mass, which could be helpful to parameterize reionization in semi-analytical models of galaxy formation on cosmological scales. We find no statistical correlation between the reionization history of field galaxies and their environment.

The inhomogeneous ionization state of the universe when the first sources of ionizing radiation appeared should lead to anisotropies in the polarization of the cosmic microwave background. We use cosmological simulations of the process by which the first sources ionized the intergalactic medium to study the induced polarization anisotropies. We find that the polarization anisotropies have rms of order ~0.01 \mu K, and local peak values of ~0.1 \mu K, smaller than those due to gravitational lensing on small scales. The polarization direction is highly coherent over degree scales. This directional coherence is not expected from either primary anisotropy or gravitational lensing effects, making the largest signals due to inhomogeneous ionization relatively easy to isolate, should experiments achieve the necessary very low noise levels.

Large fluctuations in the electron column density can occur during the reionization process. We investigate the possibility of deriving the electron density fluctuations through detailed mapping of the redshifted 21-cm emission from the neutral medium during reionization. We find that the electron-scattering optical depth and 21-cm differential brightness temperature are strongly anti-correlated, allowing optical depth estimates based entirely on redshifted 21-cm measurements. This should help isolate the CMB polarization fluctuations that are due to reionization, allowing both cleaning of the patchy reionization polarization signal as a contaminating source of confusion to other signals and a measurement of the primordial quadrupole that would be measured at various locations in the universe at the epoch of reionization. This latter application in principle allows mapping of the primordial density field at z~1100 over a large fraction of the Hubble volume.

Aims: We investigate the effect of including a proper energy balance on the interaction of a low-mass planet with a protoplanetary disk. Methods: We use a three-dimensional version of the RODEO method to perform hydrodynamical simulations including the energy equation. Radiation is included in the flux-limited diffusion approach. Results: The sign of the torque depends sensitively on the ability of the disk to radiate away the energy generated in the immediate surroundings of the planet. In the case of high opacity, corresponding to the dense inner regions of protoplanetary disks, migration is directed \emph{outward}, instead of the usual inward migration that was found in locally isothermal disks. For low values of the opacity we recover inward migration, and we show that torques originating in the coorbital region are responsible for the change in migration direction.

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.

We present detailed predictions for the redshifted 21cm signal from the epoch of reionization. These predictions are obtained from radiative transfer calculations on the results of large scale (100/h Mpc), high dynamic range, cosmological simulations. We consider several scenarios for the reionization history, of both early and extended reionization. From the simulations we construct and analyze a range of observational characteristics, from the global signal, via detailed images and spectra, to statistical representations of rms fluctuations, angular power spectra, and probability distribution functions to characterize the non-gaussianity of the 21cm signal. (abbreviated abstract)

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.

We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analyzed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5% for the Eulerian simulations while the SPH results predict the correct shape of the gap although have less resolution in the low density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10%. The spread is larger in the tidal torques acting on the planet which agree within a factor 2 at the end of the simulation. In the Neptune case the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.

We present the first large-scale radiative transfer simulations of cosmic reionization, in a simulation volume of (100/h Mpc)^3, while at the same time capturing the dwarf galaxies which are primarily responsible for reionization. We achieve this by combining the results from extremely large, cosmological, N-body simulations with a new, fast and efficient code for 3D radiative transfer, C^2-Ray. The resulting electron-scattering optical depth is in good agreement with the first-year WMAP polarization data. We show that reionization clearly proceeded in an inside-out fashion, with the high-density regions being ionized earlier, on average, than the voids. Ionization histories of smaller-size (5 to 10 comoving Mpc) subregions exibit a large scatter about the mean and do not describe the global reionization history well. The minimum reliable volume size for such predictions is ~30 Mpc. We derive the power-spectra of the neutral, ionized and total gas density fields and show that there is a significant boost of the density fluctuations in both the neutral and the ionized components relative to the total at arcminute and larger scales. We find two populations of HII regions according to their size, numerous, mid-sized (~10 Mpc) regions and a few, rare, very large regions tens of Mpc in size. We derive the statistical distributions of the ionized fraction and ionized gas density at various scales and for the first time show that both distributions are clearly non-Gaussian. (abridged)

We present numerical radiation-hydrodynamic simulations of the evolution of HII regions formed in an inhomogeneous medium resulting from turbulence simulations. We find that the filamentary structure of the underlying density distribution produces a highly irregular shape for the ionized region, in which the ionization front escapes to large distances in some directions within 80,000 years. In other directions, on the other hand, neutral gas in the form of dense globules persists within 1 parsec of the central star for the full duration of our simulation (400,000 years). Divergent photoablation flows from these globules maintain a root-mean-squared velocity in the ionized gas that is close to the ionized sound speed. Simulated images in optical emission lines show morphologies that are in strikingly detailed agreement with those observed in real HII regions.

Radiative transfer simulations are now at the forefront of numerical astrophysics. They are becoming crucial for an increasing number of astrophysical and cosmological problems; at the same time their computational cost has come to the reach of currently available computational power. Further progress is retarded by the considerable number of different algorithms (including various flavours of ray-tracing and moment schemes) developed, which makes the selection of the most suitable technique for a given problem a non-trivial task. Assessing the validity ranges, accuracy and performances of these schemes is the main aim of this paper, for which we have compared 11 independent RT codes on 5 test problems: (0) basic physics, (1) isothermal H II region expansion and (2) H II region expansion with evolving temperature, (3) I-front trapping and shadowing by a dense clump, (4) multiple sources in a cosmological density field. The outputs of these tests have been compared and differences analyzed. The agreement between the various codes is satisfactory although not perfect. The main source of discrepancy appears to reside in the multi-frequency treatment approach, resulting in different thicknesses of the ionized-neutral transition regions and different temperature structure. The present results and tests represent the most complete benchmark available for the development of new codes and improvement of existing ones. To this aim all test inputs and outputs are made publicly available in digital form.

We study the dynamics of gas and dust in a protoplanetary disk in the presence of embedded planets. We investigate the conditions for dust-gap formation in terms of particle size and planetary mass. We also monitor the amount of dust that is accreted by the planet relative to the amount of gas, which is an important parameter in determining the enrichment of solids in giant planets compared to the solid content of the central star. We use a new two-fluid hydrodynamics code to solve the flow equations for both gas and dust. For the gas, we use a Godunov-type scheme with an approximate Riemann solver (the Roe solver). The dust is treated as a pressureless fluid by essentially the same numerical method as is used for the gas. We find that it only takes a planet of 0.05 Jupiter masses to open up a gap in a disk with a significant population of mm-sized particles. Dust particles larger than 150 micron participate in gap formation. We also find that the formation of the gap severely slows down dust accretion compared to that in the gas. Therefore, it is not possible to enrich a newly formed giant planet in solids, if these solids are contained in particles with sizes from 150 micron to approximately 10 cm.

In this paper we describe a new method for studying the hydrodynamical problem of a planet embedded in a gaseous disk. We use a finite volume method with an approximate Riemann solver (the Roe solver), together with a special way to integrate the source terms. This new source term integration scheme sheds new light on the Coriolis instability, and we show that our method does not suffer from this instability. The first results on flow structure and gap formation are presented, as well as accretion and migration rates. For Mpl < 0.1 M_J and Mpl > 1.0 M_J (M_J = Jupiter's mass) the accretion rates do not depend sensitively on numerical parameters, and we find that within the disk's lifetime a planet can grow to 3-4 M_J. In between these two limits numerics play a major role, leading to differences of more than 50 % for different numerical parameters. Migration rates are not affected by numerics at all as long as the mass inside the Roche lobe is not considered. We can reproduce the Type I and Type II migration for low-mass and high-mass planets, respectively, and the fastest moving planet of 0.1 M_J has a migration time of only 2.0 10^4 yr.

We present a new numerical method for calculating the transfer of ionizing radiation, called C^2-Ray=Conservative, Causal Ray-tracing method. The method is explicitly photon-conserving, so the depletion of ionizing photons by bound-free opacity is guaranteed to equal the photoionizations these photons caused. As a result, grid cells can be large and very optically-thick without loss of accuracy. The method also uses an analytical relaxation solution for the ionization rate equations for each time step which can accommodate time steps which greatly exceed the characteristic ionization and ionization front crossing times. Together, these features make it possible to integrate the equation of transfer along a ray with many fewer cells and time steps than previous methods. For multi-dimensional calculations, the code utilizes short-characteristics ray tracing. C^2-Ray is well-suited for coupling radiative transfer to gas and N-body dynamics methods, on both fixed and adaptive grids, without imposing additional limitations on the time step and grid spacing. We present several tests of the code involving propagation of ionization fronts in one and three dimensions, in both homogeneous and inhomogeneous density fields. We compare to analytical solutions for the ionization front position and velocity, some of which we derive here for the first time.

The propagation of cosmological ionization fronts (I-fronts) during reionization is strongly influenced by small-scale structure. Here we summarize our recent attemps to understand the effect of this small-scale structure. We present high resolution cosmological N-body simulations at high-z (z>6) which resolve a wide range of halo mass, from minihalos to clusters of large, rare halos. We also study how minihalos affect I-fronts, through simulations of minihalo photoevaporation. Furthermore, we modify the I-front propagation equations to account for evolving small-scale structure, and incorporate these results into a semi-analytical reionization model. When intergalactic medium clumping and minihalo clustering around sources are included, small-scale structure affects reionization by slowing it down and extending it in time. This helps to explain observations by the Wilkinson Microwave Anisotropy Probe, which imply an early and extended reionization epoch. We also study how source clustering affects the evolution and size of H II regions, finding, in agreement with simulations, that H II regions usually expand, rarely shrinking. Hence, "relic H II regions'' are an exception, rather than the rule. When the suppression of small-mass sources in already-ionized regions by Jeans-mass filtering is accounted for, H II regions are smaller, delaying overlap. We also present a new numerical method for radiative transfer which is fast, efficient, and easily coupled to hydrodynamics and N-body codes, along with sample tests and applications.

We have developed a three-dimensional radiative transfer method designed specifically for use with parallel adaptive mesh refinement hydrodynamics codes. This new algorithm, which we call hybrid characteristics, introduces a novel form of ray tracing that can neither be classified as long, nor as short characteristics, but which applies the underlying principles, i.e. efficient execution through interpolation and parallelizability, of both. Primary applications of the hybrid characteristics method are radiation hydrodynamics problems that take into account the effects of photoionization and heating due to point sources of radiation. The method is implemented in the hydrodynamics package FLASH. The ionization, heating, and cooling processes are modelled using the DORIC ionization package. Upon comparison with the long characteristics method, we find that our method calculates the column density with a similarly high accuracy and produces sharp and well defined shadows. We show the quality of the new algorithm in an application to the photoevaporation of multiple over-dense clumps. We present several test problems demonstrating the feasibility of our method for performing high resolution three-dimensional radiation hydrodynamics calculations that span a large range of scales. Initial performance tests show that the ray tracing part of our method takes less time to execute than other parts of the calculation (e.g. hydrodynamics and adaptive mesh refinement), and that a high degree of efficiency is obtained in parallel execution. Although the hybrid characteristics method is developed for problems involving photoionization due to point sources, the algorithm can be easily adapted to the case of more general radiation fields.

We investigate the interaction of gas and dust in a protoplanetary disk in the presence of a massive planet using a new two-fluid hydrodynamics code. In view of future observations of planet-forming disks we focus on the condition for gap formation in the dust fluid. While only planets more massive than 1 Jupiter mass (MJ) open up a gap in the gas disk, we find that a planet of 0.1 MJ already creates a gap in the dust disk. This makes it easier to find lower-mass planets orbiting in their protoplanetary disk if there is a significant population of mm-sized particles.

We present a search for rings or arcs in the haloes of planetary nebulae (PNe). We discovered such structures in eight PNe, tripling the sample of PNe with known rings. This shows that, contrary to what was believed to date, the occurrence of mass loss fluctuations with timescales of ~100 to ~1000 yrs at the end of the asymptotic giant branch phase (AGB) is common. We estimate a lower limit of the occurrence rate of rings in PN haloes to be ~35%. Using these new detections and the cases previously known, we discuss the statistical properties of ring systems in PNe haloes. We estimate that the mass modulation producing the rings takes place during the last 10 or 20 kyrs of AGB evolution. In PNe, the spacing between rings ranges from <0.01 pc to 0.06 pc, significantly larger than those seen in proto-PNe. This, together with the finding of a possible positive correlation of spacing with the post-AGB age of the nebulae, suggests that the spacing of the rings increases with time. These properties, as well as the modest surface brightness amplitudes of rings, are consistent with the predictions of the dust-driven wind instability model explored by Meijerink et al. (2003), but do not immediately exclude other proposed models.

The distances to individual wind-driven bubbles such as Planetary Nebulae (PNe) can be determined using expansion parallaxes: the angular expansion velocity in the sky is compared to the radial velocity of gas measured spectroscopically. Since the one is a pattern velocity, and the other a matter velocity, these are not necessarily the same. Using the jump conditions for both shocks and ionization fronts, I show that for typical PNe the pattern velocity is 20 to 30% larger than the material velocity, and the derived distances are therefore typically 20 to 30% too low. I present some corrected distances and suggest approaches to be used when deriving distances using expansion parallaxes

At the end of their lives low mass stars such as our Sun lose most of their mass. The resulting planetary nebulae show a wide variety of shapes, from spherical to highly bipolar. According to the generalized interacting stellar winds model, these shapes are due to an interaction between a very fast tenuous outflow, and a denser environment left over from an earlier slow phase of mass loss. Previous analytical and numerical work shows that this mechanism can explain cylindrically symmetric nebulae very well. However, many circumstellar nebulae have a multipolar or point-symmetric shape. With two-dimensional calculations, Icke showed that these seemingly enigmatic forms can be easily reproduced by a two-wind model in which the confining disk is warped, as is expected to occur in irradiated disks. Here, we present the extension to fully three-dimensional adaptive mesh refinement simulations of such an interaction.

Previous analytical and numerical work shows that the generalized interacting stellar winds model can explain the observed bipolar shapes of planetary nebulae very well. However, many circumstellar nebulae have a multipolar or point-symmetric shape. With two-dimensional calculations, Icke showed that these seemingly enigmatic forms can be easily reproduced by a two-wind model in which the confining disk is warped, as is expected to occur in irradiated disks. In this contribution we present the extension to fully three-dimensional adaptive mesh refinement simulations of such an interaction.

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.

This letter presents a numerical study of the evolution of an emission line cloud of initial density 10 cm$^{-3}$, temperature $10^4$ K, and size 200 pc, being overtaken by a strong shock wave. Whereas previous simple models proposed that such a cloud would either be completely destroyed, or simply shrink in size, our results show a different and more complex behaviour: due to rapid cooling, the cloud breaks up into many small and dense fragments, which can survive for a long time. We show that such rapid cooling behaviour is in fact expected for a wide range of cloud and shock properties. This process applies to the evolution of emission line clouds being overtaken by the cocoon of a radio jet. The resulting small clouds would be Jeans unstable, and form stars. Our results thus give theoretical credibility to the process of jet induced star formation, one of the explanations for the alignment of the optical/UV and radio axis observed in high redshift radio galaxies.

We study the dynamics of stellar wind bubbles around hydrogen-deficient stars using numerical simulations with time- and ion dependent cooling. We consider two types of hydrogen-deficient stars, massive WR stars, producing Ring Nebulae, and low mass [WR] stars, producing Planetary Nebulae. We show that for the Planetary Nebulae, the different cooling properties of the hydrogen-deficient wind lead to a later transition from momentum- to energy-driven flow, which could explain the observed turbulence of these nebulae. We find that Ring Nebulae should all be energy-driven, and show how comparing the bubble's momentum and kinetic energy to the input wind momentum and kinetic energy, can give misleading information about the dynamics of the bubble.

We investigate the electron temperature of the inner halo and nebular core regions of NGC 6543, using archival Hubble Space Telescope (HST) Wide Field Planetary Camera 2 (WFPC2) images taken through narrow band [O III] filters. Balick et al. (2001) showed that the inner halo consists of a number of spherical shells. We find the temperature of this inner halo to be much higher (~15000 K) than that of the bright core nebula (~8500 K). Photo-ionization models indicate that hardening of the UV radiation from the central star cannot be the main source of the higher temperature in the halo region. Using a radiation hydrodynamic simulation, we show that mass loss and velocity variations in the AGB wind can explain the observed shells, as well as the higher electron temperature.

We report the discovery of a double ring structure in the waist of the nebula surrounding eta Carinae. The rings are detected in the mid-IR dust continuum at wavelengths of 7.9, 11.9, 12.9 and 20 micronmeter. The dust in the rings has a temperature of about 300 K. The orientation of the rings is inclined with respect to the axis of the homunculus by either 37 or 58 degrees. The central star is not in the projected centre of the structure defined by the two rings. This geometry is reminiscent of that seen in SN1987A and some planetary nebulae. We discuss several possible origins for this remarkable geometry and its orientation.

We explain the Ha emission of the cometary knots in the Helix Nebula (NGC 7293) with an analytical model that describes the emission of the head of the globules as a photoevaporated flow produced by the incident ionizing radiation of the central star.We compare these models with the Ha emission obtained from the HST archival images of the Helix Nebula. From a comparison of the Ha emission with the predictions of the analytical model we obtain a rate of ionizing photons from the central star of about 5e45 s^-1, which is consistent with estimates based on the total Hb flux of the nebula. We also model the tails of the cometary knots as a photoevaporated wind from a neutral shadow region produced by the diffuse ionizing photon field of the nebula. A comparison with the HST images allows us to obtain a direct determination of the value of the diffuse ionizing flux. We compare the ratio of diffuse to direct stellar flux as a function of radius inside an HII region with those obtained from the observational data through the analytical tail and head wind model. The agreement of this model with the values determined from the observations of the knots is excellent.

This review discusses the physics of the formation of planetary nebulae around low mass WR stars, or [WR] stars. It especially focuses on the differences which can be expected due to the different character of the fast winds from these [WR] stars. Their fast winds are more massive and are highly H deficient and metal enriched compared to the winds of normal central stars of planetary nebulae. This is expected to lead to faster expansion velocities for the nebulae and a longer momentum-driven phase in the evolution of the wind-driven bubble, leading to more turbulent nebulae. The observational evidence also shows that the process which produces the [WR] stars is unlikely to influence the onset of aspherical mass loss, something which can be used as a test for models for aspherical mass loss from AGB and post-AGB stars. Finally it is shown that the nebular characteristics rule out a very late He shell flash as the origin of most [WR] stars.

In this paper we describe the results of the first optical kinematic study of the planetary nebula BD+30 3639. This system has a central star of the Wolf-Rayet type and is believed to be fairly young. Emission line spectra were obtained at high spectral and spatial resolution using the Utrecht echelle spectrometer at the William Herschel Telescope. These spectra indicate that the main ionized shell of BD+30 3639 appears to be evolving in a more complex way than previously thought.

The first sources of ionizing radiation to condense out of the dark and neutral IGM sent ionization fronts sweeping outward through their surroundings, overtaking other primordial gas-clouds and photoevaporating them. Results are presented of the first gas dynamical simulations of this process, including radiative transfer, along with some observational diagnostics.

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.

Energy released by a small fraction of the baryons in the universe, which condensed out of while the IGM was cold, dark, and neutral,reheated and reionized it, exposing gas clouds within it to the glare of ionizing radiation. The first gas dynamical simulations of the photoevaporation of an intergalactic cloud by a quasar, including radiative transfer, are presented, along with a few observational diagnostics.

We study the evolution of dense neutral clumps located in the outer parts of planetary nebulae. These clumps will be photo-ionized by the ionizing radiation from the central star and change their structure in the process. The main effect of the ionization process is the setting up of a photo-evaporation flow and a shock running through the clump. Once this shock has moved through the entire clump it starts to accelerate because of the `rocket effect'. This continues until the entire clump has been photo-ionized. We present an analytic model for the shock and accelerating phases and also the results of numerical simulations which include detailed microphysics. We find a good match between the analytic description and the numerical results and use the numerical results to produce some of the clump's observational characteristics at different phases of its evolution. We compare the results with the properties of the fast moving low ionization knots (ansae or FLIERs) seen in a number of planetary nebulae. We find that the models match many of the kinematic and emission properties of FLIERs.

We present the results of numerical hydrodynamic models for the collimation of outflows from young stellar objects. We show that the presence of a toroidal environment can lead to efficient formation of jets and bipolar outflows from initially uncollimated central winds. The interaction between the wind and the environment leads to two types of collimation, one which is dominated by radiative cooling effects, and one which works when cooling is less efficient. We describe the two types of jets as they appear in the simulations and we suggest a description for the long term evolution of these structures in more realistic time-dependent wind sources.

Using a radiation-hydrodynamics code I follow the formation of planetary nebulae around stars of different mass. Because a more massive central star evolves much faster than a lower mass one, it is to be expected that this will affect the formation of the PN. For the stars I use the evolutionary tracks for remnants with masses of 0.605 M0 and 0.836 M0, taken from Bloecker (1995). The AGB wind is assumed to be concentrated in a thin disk, which in models without evolving stars leads to the formation of a bipolar nebula. I find that in the case of the 0.836 M0 remnant the nebula indeed acquires a bipolar shape, whereas for the 0.605 M0 remnant the shape is more elliptical. The reason for this is the time it takes to ionize the AGB material; if this happens sufficiently slowly the density distribution in the AGB wind will be smoothed out, leading to more elliptical shapes. If it happens quickly, the original density distribution (in this case a thin disk) is hardly affected. This result suggests that lower mass central stars will less easily produce bipolar nebulae, which is supported by observations.

We have chosen a reduced set of 18 ionization rate equations (for ions of H, C, N, O, S and Ne), which allow us to obtain a moderately accurate estimate of the non-equilibrium radiative cooling function. We evaluate the accuracy of this approach by comparing our cooling function with previous calculations, computed with a more extended set of ions, for the case of gas that cools from a high temperature at constant density. We also compute steady, plane shock models, which we find to compare well with models calculated with much more detailed microphysics. Using our reduced set of rate equations, we present a simulation of a radiative bow shock formed by a supersonic, plane stream impinging on a rigid sphere. The parameters for the calculation are chosen as to approximately represent typical values found for Herbig-Haro objects, and to give a cooling distance to bow shock radius ratio of 1/10. This simulation is done with an adaptive grid code, which allows a reasonably high resolution (with greater than or equal to 25 points) of the standoff distance between the bow shock and the rigid obstacle. Contrary to the standard expectation, we find that the bow shock standoff distance in the stagnation region is considerably shorter than the cooling distance behind a plane-parallel shock. Also, the centrifugal pressure is found to play an important role in determining the structure of the recombin- ation region. This appears to partially invalidate previous bow shock models based on a ``quasi-1D'' approach, at least for the particular parameters chosen for the present simulation. Finally, we present tabulations of the cooling rates that have been used (for the different species), in order to facilitate the inclusion of this treatment in other gasdynamic codes.

We explore the physics of time-dependent hydrodynamic collimation of jets from Young Stellar Objects (YSOs). Using parameters appropriate to YSOs we have carried out high resolution hydrodynamic simulations modeling the interaction of a central wind with an environment characterized by a moderate opening angle toroidal density distribution. The results show that the the wind/environment interaction produces strongly collimated supersonic jets. The jet is composed of shocked wind gas. Using analytical models of wind blown bubble evolution we show that the scenario studied here should be applicable to YSOs and can, in principle, initiate collimation on the correct scales (R ~ 100 AU). The simulations reveal a number of time-dependent non-linear features not anticipated in previous analytical studies including: a prolate wind shock; a chimney of cold swept-up ambient material dragged into the bubble cavity; a plug of dense material between the jet and bow shocks. We find that the collimation of the jet occurs through both de Laval nozzles and focusing of the wind via the prolate wind shock. Using an analytical model for shock focusing we demonstrate that a prolate wind shock can, by itself, produce highly collimated supersonic jets.

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.

We present a numerical method to solve the equations of general relativistic hydrodynamics in a given external gravitational field. The method is based on a generalization of Roe's approximate Riemann solver for the non relativistic Euler equations in Cartesian coordinates. The new method is applied to a set of standard test problems for general relativistic hydrodynamics, and is shown to perform well in comparison to existing numerical schemes. In contrast to existing explicit methods the present method can cope with strong relativistic shocks. By-products are: the characteristic form of the general relativistic Euler equations, a numerical method for special relativity that can deal with strong discontinuities, a numerical scheme for the integration of the Euler equations in an arbitrary coordinate system, possibly under the influence of (external) gravity, and a novel method to incorporate source terms in numerical schemes.

This paper looks into various aspects brought to light by numerical work on the generalized interacting winds model for planetary nebulae. First, a detailed comparison between radiative and non-radiative models is made, showing that one's naive expectations of the effects of radiative heating and cooling are not always true. Secondly, we consider the evolution of the slow wind after it has gotten ionized. It is found that the initial aspherical density distribution is smoothed out after ionization, thus requiring the aspherical nebula to be formed within the first few thousand years of PN evolution. Thirdly, the nature of the flow in the hot bubble is investigated. Both one and two-dimensional numerical models show strong signs of instabilities or turbulent flow in the hot bubble. Although observationally hard to prove or disprove this turbulent structure is critically examined. It is found that although the turbulence is not a numerical artifact, the full three-dimensional picture will most definitely differ from what is found in two dimensions. The implications for the interpretation of the models are discussed. Finally, the issue of the soft X-ray emission from PNe is considered. It is found that soft X-rays originate mainly from the thin interface between the hot bubble and the actual nebula.

The status of numerical hydrodynamical models for Planetary Nebulae is reviewed. Since all of the numerical work is based on the interacting winds model, we start with a description of this model and give an overview of the early analytical and numerical models. Subsequently we address the numerical models which include radiation effects, first of all the ones which neglect any effects of stellar evolution. These `constant environment' models are shown to closely match typical observed nebulae, both in images and kinematic data. This shows that the basic generalized interacting winds model gives a good description of the situation in aspherical PNe. Next we discuss models that do include the effects of stellar and fast wind evolution. This introduces several new effects, the most important of which are the formation of a surrounding attached envelope, and the modification of the expansion of the nebula, which helps in creating aspherical PNe very early on in their evolution. The ionization of the slow wind also leads to a gradual smoothing out of its aspherical character, working against aspherical PNe forming in later stages. Finally we discuss some applications of the model to nebular problems.

Using a radiation-gasdynamics code the evolution of spherical planetary nebulae is followed, while taking into account the evolution of central star and the fast wind. These models show the importance of ionization fronts for the structure of planetary nebulae, especially for the so called multiple shell nebulae (MSPNe). It is shown that the outer shell is formed by the ionization front while the inner shell is swept-up by the fast wind. These models explain the emission profiles of the outer shells as well as their various kinematic properties. Because they are shaped by the ionization front these outer shells only give indirect information on the AGB mass loss history. The models indicate that typical MSPN structures point to mass loss variations during the AGB phase. The ionization also leads to a stalling of the expansion of the nebula, leading to nebulae with expansion ages lower than their evolutionary age. Values for ionized mass and Zanstra temperatures are derived from the models.

We present the results of two-dimensional radiation-gasdynamic simulations of aspherical Planetary Nebulae (PNe) evolution. These simulations were constructed using the Generalized Interacting Stellar Winds (GISW) scenario of Balick (1987) where a fast, tenuous wind from the central star expands into a toroidal, slow, dense wind. We demonstrate that the GISW model can produce a wide range of aspherical flow patterns. We have constructed self-consistent synthetic observations of the models from forbidden line emissivities used in the energy loss term. We present integrated intensity and long-slit spectrum, (Position-Velocity) maps of the models projected at different angles on the sky. These synthetic observations are compared with real intensity and Position-Velocity maps of PNe. We find that there is a very good match between the synthetic and real observations in terms of morphologies, kinematics, and physical conditions. From the results of these simulations we conclude that the GISW scenario can account for most, if not all, PNe morphologies, thus confirming Balick's (1987) conjecture.