Using a suite of detailed numerical simulations we estimate the level of anisotropy generated by the time evolution along the light cone of the 21cm signal from the epoch of reionization. Our simulations include the physics necessary to model the signal during both the late emission regime and the early absorption regime, namely X-ray and Lyman-band 3D radiative transfer in addition to the usual dynamics and ionizing UV transfer. The signal is analysed using correlation functions perpendicular and parallel to the line of sight (LOS). We reproduce general findings from previous theoretical studies: the overall amplitude of the correlations and the fact that the light cone anisotropy is visible only on large scales (100 cMpc). However, the detailed behaviour is different. At 3 different epochs, the amplitude of the correlations along and perpendicular to the LOS differ from each other, indicating anisotropy. These 3 epochs are associated with 3 events of the global reionization history: the overlap of ionized bubbles, the onset of mild heating by X-rays in regions around the sources, and the onset of efficient Lyman-alpha coupling in regions around the sources. A 20x20 deg^2 survey area may be necessary to mitigate sample variance when we use the directional correlation functions. On a 100 cMpc scale the light cone anisotropy dominates over the anisotropy generated by peculiar velocity gradients computed in the linear regime. By modelling instrumental noise and limited resolution, we find that the anisotropy should be easily detectable by the SKA, assuming perfect foreground removal, the limiting factor being a large enough survey size. In the case of the LOFAR, it is likely that only first anisotropy episode will fall in the observing frequency range and will be detectable only if sample variance is much reduced (i.e. a larger than 20x20 deg^2 survey, which is not presently planned).

It is predicted that sources emitting UV radiation in the Lyman band during the epoch of reionization (EoR) showed a series of discontinuities in their Ly-alpha flux radial profile as a consequence of the thickness of the Lyman line series in the primeval intergalactic medium. Through unsaturated Wouthuysen-Field coupling, these spherical discontinuities are also present in the 21 cm emission of the neutral IGM. In this article, we study the effects these discontinuities have on the differential brightness temperature of the 21 cm signal of neutral hydrogen in a realistic setting including all other sources of fluctuations. We focus on the early phases of the EoR, and we address the question of the detectability by the planned Square Kilometre Array. Such a detection would be of great interest, because these structures could provide an unambiguous diagnostic for the cosmological origin of the signal remaining after the foreground cleaning procedure. Also, they could be used as a new type of standard rulers. We determine the differential brightness temperature of the 21 cm signal in the presence of inhomogeneous Wouthuysen-Field effect using simulations which include (hydro)dynamics and both ionizing and Lyman lines 3D radiative transfer with the code LICORICE. We find that the Lyman horizons are clearly visible on the maps and radial profiles around the first sources of our simulations, but for a limited time interval, typically \Delta z \approx 2 at z \sim 13. Stacking the profiles of the different sources of the simulation at a given redshift results in extending this interval to \Delta z \approx 4. When we take into account the implementation and design planned for the SKA (collecting area, sensitivity, resolution), we find that detection will be challenging. It may be possible with a 10 km diameter for the core, but will be difficult with the currently favored design of a 5 km core.

Primordial molecules were formed during the Dark Ages, i.e. the time between recombination and reionization in the early Universe. The purpose of this article is to analyze the formation of primordial molecules based on heavy elements during the Dark Ages, with elemental abundances taken from different nucleosynthesis models. We present calculations of the full non-linear equation set governing the primordial chemistry. We consider the evolution of 45 chemical species and use an implicit multistep method of variable order of precision with an adaptive stepsize control. We find that the most abundant Dark Ages molecules based on heavy elements are CH and OH. Non-standard nucleosynthesis can lead to higher heavy element abundances while still satisfying the observed primordial light abundances. In that case, we show that the abundances of molecular species based on C, N, O and F can be enhanced by two orders of magnitude compared to the standard case, leading to a CH relative abundance higher than that of HD+ or H2D+.

Some models of Big Bang nucleosynthesis suggest that very high baryon density regions were formed in the early Universe, and generated the production of heavy elements other than lithium such as fluorine F. We present a comprehensive chemistry of fluorine in the post-recombination epoch. Calculation of F, F- and HF abundances, as a function of redshift z, are carried out. The main result is that the chemical conditions in the early Universe can lead to the formation of HF. The final abundance of the diatomic molecule HF is predicted to be close to 3.75 10(-17) when the initial abundance of neutral fluorine F is 10(-15). These results indicate that molecules of fluorine HF were already present during the dark age. This could have implications on the evolution of proto-objects and on the anisotropies of cosmic microwave background radiation. Hydride of fluorine HF may affect enhancement of the emission line intensity from the proto-objects and could produce spectral-spatial fluctuations.