Search

To self-consistently model galactic properties, reionization of the intergalactic medium, and the associated 21-cm signal, we have developed the algorithm polar by integrating the one-dimensional radiative transfer code grizzly with the semi-analytical galaxy formation code L-Galaxies 2020. Our proof-of-concept results are consistent with observations of the star formation rate history, UV luminosity function and the CMB Thomson scattering optical depth. We then investigate how different galaxy formation models affect UV luminosity functions and 21-cm power spectra, and find that while the former are most sensitive to the parameters describing the merger of halos, the latter have a stronger dependence on the supernovae feedback parameters, and both are affected by the escape fraction model., Comment: MNRAS accepted

The presence of charged dust grains is known to have a profound impact on the physical evolution of the multiphase interstellar medium (ISM). Despite its importance, this process is still poorly explored in numerical simulations due to its complex physics and the tight dependence on the environment. Here we introduce a novel implementation of grain charging in the cosmological radiative transfer code CRASH. We first benchmark the code predictions on a series of idealized dusty HII regions created by a single star, in order to assess the impact of grain properties on the resulting spatial distribution of charges. Second, we perform a realistic radiative transfer simulation of a star forming region extracted from a dusty galaxy evolving in the Epoch of Reionization. We find that $\sim 13$ % of the total dust mass gets negatively charged, mainly silicate and graphite grains of radius micron. A complex spatial distribution of grain charges is also found, primarily depending on the exposure to stellar radiation and strongly varying along different lines of sight, as a result of radiative transfer effects. We finally assess the impact of grain properties (both chemical composition and size) on the resulting charge distribution. The new implementation described here will open up a wide range of possible studies investigating the physical evolution of the dusty ISM, nowadays accessible to observations of high and low redshift galaxies., Comment: MNRAS in press

Overdense, metal-rich regions, shielded from stellar radiation might remain neutral throughout reionization and produce metal as well as 21 cm absorption lines. Simultaneous absorption from metals and 21 cm can complement each other as probes of underlying gas properties. We use Aurora, a suite of high resolution radiation-hydrodynamical simulations of galaxy formation, to investigate the occurrence of such "aligned" absorbers. We calculate absorption spectra for 21 cm, OI, CII, SiII and FeII. We find velocity windows with absorption from at least one metal and 21 cm, and classify the aligned absorbers into two categories: 'aligned and cospatial absorbers' and 'aligned but not cospatial absorbers'. While 'aligned and cospatial absorbers' originate from overdense structures and can be used to trace gas properties, 'aligned but not cospatial absorbers' are due to peculiar velocity effects. The incidence of absorbers is redshift dependent, as it is dictated by the interplay between reionization and metal enrichment, and shows a peak at $z \approx 8$ for the aligned and cospatial absorbers. While aligned but not cospatial absorbers disappear towards the end of reionization because of the lack of an ambient 21 cm forest, aligned and cospatial absorbers are associated with overdense pockets of neutral gas which can be found at lower redshift. We produce mock observations for a set of sightlines for the next generation of telescopes like the ELT and SKA1-LOW, finding that given a sufficiently bright background quasar, these telescopes will be able to detect both types of absorbers throughout reionization., Accepted in MNRAS. Comments welcome

The Murchison Widefield Array (MWA) team has derived new upper limits on the spherically averaged power spectrum of the 21-cm signal at six redshifts in the range $z \approx 6.5-8.7$. We use these upper limits and a Bayesian inference framework to derive constraints on the ionization and thermal state of the intergalactic medium (IGM) as well as on the strength of a possible additional radio background. We do not find any constraints on the state of the IGM for $z\gtrsim 7.8$ if no additional radio background is present. In the presence of such a radio background, the 95 per cent credible intervals of the disfavoured models at redshift $\gtrsim 6.5 $ correspond to an IGM with a volume averaged fraction of ionized regions below 0.6 and an average gas temperature $\lesssim 10^3$ K. In these models, the heated regions are characterised by a temperature larger than that of the radio background, and by a distribution with characteristic size $\lesssim 10$ $h^{-1}$ Mpc and a full width at half maximum (FWHM) of $\lesssim 30$ $h^{-1}$ Mpc. Within the same credible interval limits, we exclude an additional radio background of at least $0.008\%$ of the CMB at 1.42 GHz., Comment: 13 Pages, 6 Figures, 5 Tables, Accepted for publication in MNRAS

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., Comment: 12 pages, 4 pages, 1 table. Submitted to MNRAS, comments welcome

We use a set of 3D radiative transfer simulations to study the effect that a large fraction of binary stars in galaxies during the epoch of reionization has on the physical properties of the intergalactic medium (i.e. the gas temperature and the ionization state of hydrogen and helium), on the topology of the ionized bubbles and on the 21 cm power spectra. Consistently to previous literature, we find that the inclusion of binary stars can speed up the reionization process of HI and HeI, while HeII reionization is still dominated by more energetic sources, especially accreting black holes. The earlier ionization attained with binary stars allows for more time for cooling and recombination, so that gas fully ionized by binary stars is typically colder than that ionized by single stars at any given redshift. With the same volume averaged ionization fraction, the inclusion of binary stars results in fewer small ionized bubbles and more large ones, with visible effects also on the large scales of the 21 cm power spectrum., Comment: 14 pages, 11 figures, MNRAS accepted

Full list of authors: Morabito, L. K.; Jackson, N. J.; Mooney, S.; Sweijen, F.; Badole, S.; Kukreti, P.; Venkattu, D.; Groeneveld, C.; Kappes, A.; Bonnassieux, E.; Drabent, A.; Iacobelli, M.; Croston, J. H.; Best, P. N.; Bondi, M.; Callingham, J. R.; Conway, J. E.; Deller, A. T.; Hardcastle, M. J.; McKean, J. P.; Miley, G. K.; Moldon, J.; Röttgering, H. J. A.; Tasse, C.; Shimwell, T. W.; van Weeren, R. J.; Anderson, J. M.; Asgekar, A.; Avruch, I. M.; van Bemmel, I. M.; Bentum, M. J.; Bonafede, A.; Brouw, W. N.; Butcher, H. R.; Ciardi, B.; Corstanje, A.; Coolen, A.; Damstra, S.; de Gasperin, F.; Duscha, S.; Eislöffel, J.; Engels, D.; Falcke, H.; Garrett, M. A.; Griessmeier, J.; Gunst, A. W.; van Haarlem, M. P.; Hoeft, M.; van der Horst, A. J.; Jütte, E.; Kadler, M.; Koopmans, L. V. E.; Krankowski, A.; Mann, G.; Nelles, A.; Oonk, J. B. R.; Orru, E.; Paas, H.; Pandey, V. N.; Pizzo, R. F.; Pandey-Pommier, M.; Reich, W.; Rothkaehl, H.; Ruiter, M.; Schwarz, D. J.; Shulevski, A.; Soida, M.; Tagger, M.; Vocks, C.; Wijers, R. A. M. J.; Wijnholds, S. J.; Wucknitz, O.; Zarka, P.; Zucca, P., The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz. However, it is technically and logistically challenging to process LOFAR data at this resolution. To date only a handful of publications have exploited this capability. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. It is implemented in a pipeline using mostly dedicated LOFAR software tools and the same processing framework as the LOFAR Two-metre Sky Survey (LoTSS). We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LoTSS pointing with an 8 h observation and 13 international stations. We perform in-field delay calibration, solution referencing to other calibrators in the field, self-calibration of these calibrators, and imaging of example directions of interest in the field. We find that for this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5° away, while phase solution transferral works well over ~1°. We also demonstrate a check of the astrometry and flux density scale with the in-field delay calibrator source. Imaging in 17 directions, we find the restoring beam is typically ~0.3′′ ×0.2′′ although this varies slightly over the entire 5 deg2 field of view. We find we can achieve ~80–300 μJy bm−1 image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 μJy bm−1 for the 8 h observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image roughly 900 sources per LoTSS pointing. This equates to ~ 3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution makes this estimate a lower limit. © ESO 2022., This work made use of the Dutch national e-infrastructure with the support of the SURF Cooperative using grant no. EINF-262 LKM is grateful for support from the Medical Research Council (grant MR/T042842/1). S.M. acknowledges support from the Governmentof Ireland Postgraduate Scholarship Programme. E.B. acknowledges support from the ERC-ERG grant DRANOEL, n.714245. A.D. acknowledges support by the BMBF Verbundforschung under the grant 052020. J.H.C. acknowledges support from the UK Science and Technology Facilities Council (ST/R000794/1). P.N.B. is grateful for support from the UK STFC via grant ST/R000972/1. J.R.C. thanks the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO) for support via the Talent Programme Veni grant. M.J.H. acknowledges support from the UK Science and Technology Facilities Council (ST/R000905/1). J.P.M. acknowledges support from the NetherlandsOrganization for Scientific Research (NWO, project number 629.001.023) and the Chinese Academy of Sciences (CAS, project number 114A11KYSB20170054). J.M. acknowledges financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award to the Instituto de Astrofísicade Andalucía (SEV-2017-0709) and from the grant RTI2018-096228-B-C31 (MICIU/FEDER, EU). R.J.v.W. acknowledges support from the ERC Starting Grant ClusterWeb 804208. D.J.S. acknowledges support by the GermanFederal Ministry for Science and Research BMBF-Verbundforschungsprojekt D-LOFAR 2.0 (grant numbers 05A20PB1). LOFAR (van Haarlem et al. 2013) is the Low Frequency Array designed and constructed by ASTRON. It has observing, data processing, and data storage facilities in several countries, that are owned by various parties (each with their own funding sources), and that are collectively operated by the ILT foundation under a joint scientific policy. The ILT resources have benefitted from the following recent major funding sources: CNRS-INSU, Observatoire de Paris and Université d’Orléans, France; BMBF, MIWF-NRW, MPG, Germany; Science Foundation Ireland (SFI), Department of Business, Enterprise and Innovation (DBEI), Ireland; NWO, The Netherlands; The Science and Technology Facilities Council, UK; Ministry of Science and Higher Education, Poland.

Cosmic gas makes up about 90% of baryonic matter in the Universe and H$_2$ is the closest molecule to star formation. In this work we study cold neutral gas and its H$_2$ component at different epochs, exploiting state-of-the-art hydrodynamic simulations that include time-dependent atomic and molecular non-equilibrium chemistry coupled to star formation, feedback effects, different UV backgrounds presented in the recent literature and a number of additional processes - such as gas self-shielding, H$_2$ dust grain catalysis, photoelectric and cosmic-ray heating - occurring during structure formation (ColdSIM). We find neutral-gas mass density parameters $ \Omega_{\rm neutral} \simeq $10$^{-3}$ and increasing from lower to higher redshift, in agreement with available HI data. Resulting H$_2$ fractions can be as high as $\sim $50% at $z\sim $4-8, in line with the latest high-$z$ measurements. Albeit dependent on the adopted UV background, derived $ \Omega_{\rm H_2} $ values agree with observations up to $z\sim$7 and both HI and H$_2$ trends are better reproduced by our non-equilibrium H$_2$-based star formation modelling. The predicted gas depletion timescales decrease towards lower $z$, with H$_2$ depletion times remaining below the Hubble time and comparable to the dynamical time at all considered redshifts. This implies that non-equilibrium molecular cooling is efficient at driving cold-gas collapse in a broad variety of environments and since the very early cosmic epochs. In appendix, we show detailed analyses of individual processes, as well as simple numerical parameterizations and fits to account for them. Our findings suggest that, in addition to HI, non-equilibrium H$_2$ observations are pivotal probes for assessing cold-gas abundances and the role of UV background radiation - Abridged, Comment: 15 pages + appendices; accepted for publication on A&A; comments are welcome

We present simulations of cosmic reionization and reheating from $z=18$ to $z=5$, investigating the role of stars (emitting soft UV-photons), nuclear black holes (BHs, with power-law spectra), X-ray binaries (XRBs, with hard X-ray dominated spectra), and the supernova-associated thermal bremsstrahlung of the diffuse interstellar medium (ISM, with soft X-ray spectra). We post-process the hydrodynamical simulation Massive-Black II (MBII) with multifrequency ionizing radiative transfer. The source properties are directly derived from the physical environment of MBII, and our only real free parameter is the ionizing escape fraction $f_{\rm esc}$. We find that, among the models explored here, the one with an escape fraction that decreases with decreasing redshift yields results most in line with observations, such as of the neutral hydrogen fraction and the Thomson scattering optical depth. Stars are the main driver of hydrogen reionization and consequently of the thermal history of the intergalactic medium (IGM). We obtain $\langle x_{\rm HII} \rangle = 0.99998$ at $z=6$ for all source types, with volume averaged temperatures $\langle T \rangle \sim 20,000~{\rm K}$. BHs are rare and negligible to hydrogen reionization, but conversely they are the only sources which can fully ionize helium, increasing local temperatures by $\sim 10^4~{\rm K}$. The thermal and ionization state of the neutral and lowly ionized hydrogen differs significantly with different source combinations, with ISM and (to a lesser extent) XRBs, playing a significant role and, as a consequence, determining the transition from absorption to emission of the 21 cm signal from neutral hydrogen., Comment: 17 pages, 19 figures, accepted for publication in MNRAS

We introduce a versatile and spatially resolved morphological characterisation of binary fields, rooted in the opening transform of mathematical morphology. We subsequently apply it to the thresholded ionization field in simulations of cosmic reionization and study the morphology of ionized regions. We find that an ionized volume element typically resides in an ionized region with radius $\sim8\,h^{-1}\mathrm{cMpc}$ at the midpoint of reionization ($z\approx7.5$) and follow the bubble size distribution even beyond the overlap phase. We find that percolation of the fully ionized component sets in when 25% of the universe is ionized and that the resulting infinite cluster incorporates all ionized regions above $\sim8\,h^{-1}\mathrm{cMpc}$. We also quantify the clustering of ionized regions of varying radius with respect to matter and on small scales detect the formation of superbubbles in the overlap phase. On large scales we quantify the bias values of the centres of ionized and neutral regions of different sizes and not only show that the largest ones at the high-point of reionization can reach $b\approx 30$, but also that early small ionized regions are positively correlated with matter and large neutral regions and late small ionized regions are heavily anti-biased with respect to matter, down to $b\lesssim-20$., 18 pages, 15 figure, as accepted for publication by MNRAS

We investigate systematic effects in direction dependent gain calibration in the context of the Low-Frequency Array (LOFAR) 21-cm Epoch of Reionization (EoR) experiment. The LOFAR EoR Key Science Project aims to detect the 21-cm signal of neutral hydrogen on interferometric baselines of $50-250 \lambda$. We show that suppression of faint signals can effectively be avoided by calibrating these short baselines using only the longer baselines. However, this approach causes an excess variance on the short baselines due to small gain errors induced by overfitting during calibration. We apply a regularised expectation-maximisation algorithm with consensus optimisation (sagecal-co) to real data with simulated signals to show that overfitting can be largely mitigated by penalising spectrally non-smooth gain solutions during calibration. This reduces the excess power with about a factor 4 in the simulations. Our results agree with earlier theoretical analysis of this bias-variance trade off and support the gain-calibration approach to the LOFAR 21-cm signal data., Comment: This is a pre-copyedited, author-produced PDF of an article accepted for publication in MNRAS following peer review

We present a definitive assessment of the role of Inverse Compton scattering of Cosmic Microwave Background photons (IC/CMB) in the context of radio galaxies. Owing to the steep increase of the CMB radiation energy density, IC/CMB is supposed to become progressively more important with respect to radio synchrotron cooling as the redshift increases. For typical energies at play, this process will up-scatter the CMB photons into the X-ray band, and is thus expected to yield a redshift-dependent, concurrent X-ray brightening and radio dimming of the jet-powered structures. Here we show how a conclusive proof of this effect hinges on high-resolution imaging data in which the extended lobes can be distinguished from the compact hot spots where synchrotron-self-Compton dominates the X-ray emission regardless of redshift. We analyze Chandra and Very Large Array data of 11 radio galaxies between 1.33 radio galaxies below, Comment: 16 pages, 5 figures, 5 tables. Accepted for publication in MNRAS

High-precision pulsar timing requires accurate corrections for dispersive delays of radio waves, parametrized by the dispersion measure (DM), particularly if these delays are variable in time. In a previous paper we studied the Solar-wind (SW) models used in pulsar timing to mitigate the excess of DM annually induced by the SW, and found these to be insufficient for high-precision pulsar timing. Here we analyze additional pulsar datasets to further investigate which aspects of the SW models currently used in pulsar timing can be readily improved, and at what levels of timing precision SW mitigation is possible. Our goals are to verify: a) whether the data are better described by a spherical model of the SW with a time-variable amplitude rather than a time-invariant one as suggested in literature, b) whether a temporal trend of such a model's amplitudes can be detected. We use the pulsar-timing technique on low-frequency pulsar observations to estimate the DM and quantify how this value changes as the Earth moves around the Sun. Specifically, we monitor the DM in weekly to monthly observations of 14 pulsars taken with LOFAR across time spans of up to 6 years. We develop an informed algorithm to separate the interstellar variations in DM from those caused by the SW and demonstrate the functionality of this algorithm with extensive simulations. Assuming a spherically symmetric model for the SW density, we derive the amplitude of this model for each year of observations. We show that a spherical model with time-variable amplitude models the observations better than a spherical model with constant amplitude, but that both approaches leave significant SW induced delays uncorrected in a number of pulsars in the sample. The amplitude of the spherical model is found to be variable in time, as opposed to what has been previously suggested., 14 pages, 9 figures. Accepted for publication in Astronomy and Astrophysics

To investigate the potential abundance and impact of nuclear black holes (BHs) during reionization, we generate a neural network that estimates their masses and accretion rates by training it on 23 properties of galaxies harbouring them at $z=6$ in the cosmological hydrodynamical simulation Massive-Black II. We then populate all galaxies in the simulation from $z=18$ to $z=5$ with BHs from this network. As the network allows to robustly extrapolate to BH masses below those of the BH seeds, we predict a population of faint BHs with a turnover-free luminosity function, while retaining the bright (and observed) BHs, and together they predict a Universe in which intergalactic hydrogen is $15\%$ ionized at $z=6$ for a clumping factor of 5. Faint BHs may play a stronger role in H reionization without violating any observational constraints. This is expected to have an impact also on pre-heating and -ionization, which is relevant to observations of the 21 cm line from neutral H. We also find that BHs grow more efficiently at higher $z$, but mainly follow a redshift-independent galaxy-BH relation. We provide a power law parametrisation of the hydrogen ionizing emissivity of BHs., 8 pages, 7 figures, accepted for publication in MNRAS

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, Comment: 8 pages, 8 figures, accepted for publication in MNRAS

We use hydrodynamics and radiative transfer simulations to study the 21~cm signal around a bright QSO at $z \sim 10$. Due to its powerful UV and X-ray radiation, the QSO quickly increases the extent of the fully ionized bubble produced by the pre-existing stellar type sources, in addition to partially ionize and heat the surrounding gas. As expected, a longer QSO lifetime, $t_{\rm QSO}$, results in a 21~cm signal in emission located at increasingly larger angular radii, $\theta$, and covering a wider range of $\theta$. Similar features can be obtained with a higher galactic emissivity efficiency, $f_{\rm UV}$, so that determining the origin of a large ionized bubble (i.e. QSO vs stars) is not straightforward. Such degeneracy could be reduced by taking advantage of the finite light traveltime effect, which is expected to affect an HII region produced by a QSO differently from one created by stellar type sources. From an observational point of view, we find that the 21 cm signal around a QSO at various $t_{\rm QSO}$ could be detected by SKA1-low with a high signal-noise ratio (S/N). As a reference, for $t_{\rm QSO} = 10\,\rm Myr$, a S/N $\sim 8$ is expected assuming that no pre-heating of the IGM has taken place due to high-$z$ energetic sources, while it can reach value above 10 in case of pre-heating. Observations of the 21~cm signal from the environment of a high-$z$ bright QSO could then be used to set constraints on its lifetime, as well as to reduce the degeneracy between $f_{\rm UV}$ and $t_{\rm QSO}$., Comment: 11 pages, 10 figures, apj accepted

We report on deep Chandra X-ray Telescope imaging observations of 4C 63.20, one of the few known radio galaxies at z>3.5. The X-ray counterpart is resolved into a core plus two off-nuclear sources that (combined) account for close to 30% of the total X-ray flux. Their morphology and orientation are consistent with a diffuse, lobe-like nature, albeit compact hotspots cannot be ruled out. The broadband spectral energy distribution of 4C 63.20 can be reproduced with a jet model where the majority of the radio flux can be ascribed to synchrotron emission from the hotspots, whereas the (non-nuclear) X-ray emission is produced via Inverse Compton (IC) off of Cosmic Microwave Background (CMB) photons within the extended lobes. This scenario is broadly consistent with the expectation from highly magnetized lobes in a hotter CMB, and supports the view that IC/CMB may quench less extreme radio lobes at high redshifts., 10 pages, 8 figures, accepted to MNRAS

This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR "core" reveals two different velocities in the scintillation pattern: a primary velocity of ~30m/s with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110m/s with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller--scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported., Comment: 24 pages, 16 figures. Accepted for open-access publication in the Journal of Space Weather and Space Climate. For associated movie file, see https://www.swsc-journal.org/10.1051/swsc/2020010/olm

The hyperfine transition of $^3$He$^+$ at 3.5cm has been thought as a probe of the high-z IGM since it offers a unique insight into the evolution of the helium component of the gas, as well as potentially give an independent constraint on the 21cm signal from neutral hydrogen. In this paper, we use radiative transfer simulations of reionization driven by sources such as stars, X-ray binaries, accreting black holes and shock heated interstellar medium, and simulations of a high-z quasar to characterize the signal and analyze its prospects of detection. We find that the peak of the signal lies in the range 1-50 $\mu$K for both environments, but while around the quasar it is always in emission, in the case of cosmic reionization a brief period of absorption is expected. As the evolution of HeII is determined by stars, we find that it is not possible to distinguish reionization histories driven by more energetic sources. On the other hand, while a bright QSO produces a signal in 21cm that is very similar to the one from a large collection of galaxies, its signature in 3.5cm is very peculiar and could be a powerful probe to identify the presence of the QSO. We analyze the prospects of the signal's detectability using SKA1-mid as our reference telescope. We find that the noise power spectrum dominates over the power spectrum of the signal, although a modest S/N ratio can be obtained when the wavenumber bin width and the survey volume are sufficiently large., Comment: 10 pages, 13 figures, accepted for publication in MNRAS

2MASS J050051.85$-$093054.9 is the closest known low-mass helium-core white dwarf in a binary system. We used three high-band international LOFAR stations to perform a targeted search for a pulsar companion, reaching sensitivities of ~3 mJy for a 10-ms pulsar at a DM = 1 pc$\cdot$cm$^{-3}$. No pulsed signal was detected, confidently excluding the presence of a detectable radio pulsar in the system., Accepted for publication in the Research Notes of the AAS

Recent results have suggested that active galactic nuclei (AGN) could provide enough photons to reionise the Universe. We assess the viability of this scenario using a semi-numerical framework for modeling reionisation, to which we add a quasar contribution by constructing a Quasar Halo Occupation Distribution (QHOD) based on Giallongo et al. observations. Assuming a constant QHOD, we find that an AGN-only model cannot simultaneously match observations of the optical depth $\tau_e$, neutral fraction, and ionising emissivity. Such a model predicts $\tau_e$ too low by $\sim 2\sigma$ relative to Planck constraints, and reionises the Universe at $z\lesssim 5$. Arbitrarily increasing the AGN emissivity to match these results yields a strong mismatch with the observed ionising emissivity at $z\sim 5$. If we instead assume a redshift-independent AGN luminosity function yielding an emissivity evolution like that assumed in Madau & Haardt model, then we can match $\tau_e$ albeit with late reionisation, however such evolution is inconsistent with observations at $z\sim 4-6$ and poorly motivated physically. These results arise because AGN are more biased towards massive halos than typical reionising galaxies, resulting in stronger clustering and later formation times. AGN-dominated models produce larger ionising bubbles that are reflected in $\sim\times 2$ more 21cm power on all scales. A model with equal parts galaxies and AGN contribution is still (barely) consistent with observations, but could be distinguished using next-generation 21cm experiments HERA and SKA-low. We conclude that, even with recent claims of more faint AGN than previously thought, AGN are highly unlikely to dominate the ionising photon budget for reionisation., Comment: 16 pages, 9 figures, matches the accepted version for publication in MNRAS, 2017

We use numerical N-body hydrodynamical simulations with varying PopIII stellar models to investigate the possibility of detecting first star signatures with observations of high-redshift damped Lyα absorbers (DLAs). The simulations include atomic and molecular cooling, star formation, energy feedback and metal spreading due to the evolution of stars with a range of masses and metallicities. Different initial mass functions (IMFs) and corresponding metal-dependent yields and lifetimes are adopted to model primordial stellar populations. The DLAs in the simulations are selected according to either the local gas temperature (temperature selected) or the host mass (mass selected). We find that 3 per cent (40 per cent) of mass (temperature)-selected high-z (z ≥ 5.5) DLAs retain signatures of pollution from PopIII stars, independent of the first star model. Such DLAs have low halo mass ( < 109.6 M☉), metallicity ( < 10- 3 Z☉) and star formation rate ( < 10- 1.5 M☉ yr- 1). Metal abundance ratios of DLAs imprinted in the spectra of quasi-stellar object can be useful as tools to infer the properties of the polluting stellar generation and to constrain the first star mass ranges. Comparing the abundance ratios derived from our simulations to those observed in DLAs at z ≥ 5, we find that most of these DLAs are consistent within errors with PopII star dominated enrichment and strongly disfavour the pollution pattern of very massive first stars (i.e. 100-500 M☉). However, some of them could still result from the pollution of first stars in the mass range [0.1, 100] M☉. In particular, we find that the abundance ratios from SDSS J1202+3235 are consistent with those expected from PopIII enrichment dominated by massive (but not extreme) first stars.

We study statistics of the 21-cm signal from the epoch of reionization. We propose to use [OIII] line emitting galaxies to cross-correlate with the 21-cm signal from $z = 7 - 10$. To this aim, we employ simulations of reionization obtained post-processing the high-resolution cosmological hydrodynamic simulation Massive Black-II with the three-dimensional (3D) radiative transfer code CRASH to follow the propagation of ionizing photons from a variety of sources. We show that, during the early phases of reionization, the 21-cm signal is positively correlated with the spatial distribution of the [OIII] emitters on large scales $(k < 1h~\rm cMpc^{-1})$. This positive correlation is generated by the temperature-galaxy correlation and it is a few times larger than when we assume that the heating is saturated. As the reionized regions expand, the correlation changes its sign to negative from $z = 10$ to 8. The signals at this epoch can be detected by combining the Square Kilometre Array (SKA) and a wide-field [OIII] emitter survey. We also calculate the cross-power spectrum with a 3D [OIII] intensity field, aiming at exploiting future intensity mapping observations. We conclude that high-redshift [OIII] line emitters can be used to probe the reionization process when the intergalactic medium is largely neutral., 7 pages, 5 figures, published in MNRAS

Galactic outflows are critical to our understanding of galaxy formation and evolution. However the details of the underlying feedback process remain unclear. We compare Ly$\alpha$ observations of the circumgalactic medium (CGM) of Lyman Break Galaxies (LBGs) with mock observations of their simulated CGM. We use cosmological hydrodynamical `zoom-in' simulations of an LBG which contains strong, momentum-driven galactic outflows. Simulation snapshots at $z=2.2$ and $z=2.65$ are used, corresponding to the available observational data. The simulation is post-processed with the radiative transfer code \textsc{crash} to account for the impact of ionising photons on hydrogen gas surrounding the simulated LBG. We generate mock absorption line maps for comparison with data derived from observed close galaxy-galaxy pairs. We perform calculations of Ly$\alpha$ photons scattering through the CGM with our newly developed Monte-Carlo code \textsc{slaf}, and compare to observations of diffuse Ly$\alpha$ halos around LBGs. Our fiducial galactic outflow model comes closer to reproducing currently observed characteristics of the CGM in Ly$\alpha$ than a reference inefficient feedback model used for comparison. Nevertheless, our fiducial model still struggles to reproduce the observed data of the inner CGM (at impact parameter $b, Comment: 14 pages, 12 figures

We present cosmological hydrodynamical simulations including atomic and molecular non-equilibrium chemistry, multi-frequency radiative transfer (0.7-100 eV sampled over 150 frequency bins) and stellar population evolution to investigate the host candidates of the seeds of supermassive black holes coming from direct collapse of gas in primordial haloes (direct-collapse black holes, DCBHs). We consistently address the role played by atomic and molecular cooling, stellar radiation and metal spreading of C, N, O, Ne, Mg, Si, S, Ca, Fe, etc. from primordial sources, as well as their implications for nearby quiescent proto-galaxies under different assumptions for early source emissivity, initial mass function and metal yields. We find that putative DCBH host candidates need powerful primordial stellar generations, since common solar-like stars and hot OB-type stars are neither able to determine the conditions for direct collapse nor capable of building up a dissociating Lyman-Werner background radiation field. Thermal and molecular features of the identified DCBH host candidates in the scenario with very massive primordial stars seem favourable, with illuminating Lyman-Werner intensities featuring values of 1-50 J21. Nevertheless, additional non-linear processes, such as merger events, substructure formation, rotational motions and photo-evaporation, should inhibit pure DCBH formation in 2/3 of the cases. Local turbulence may delay gas direct collapse almost irrespectively from other environmental conditions. The impact of large Lyman-Werner fluxes at distances smaller than 5 kpc is severely limited by metal pollution., references added

We present a new dust extension to the Monte Carlo radiative transfer code crash, which enables it to simulate the propagation of ionizing radiation through mixtures of gas and dust. The new code is applied to study the impact of dust absorption on idealized galactic H II regions and on small scale reionization. We find that H II regions are reduced in size by the presence of dust, while their inner temperature and ionization structure remain largely unaffected. In the small scale reionization simulation, dust hardens ionization fronts and delays the overlap of ionized bubbles. This effect is found to depend only weakly on the assumed abundance of dust in underdense regions., 17 pages, 14 figures. Accepted for publication in MNRAS

Publisher's version (útgefin grein), We use high-resolution hydrodynamical simulations to study the contribution to the X-ray background from high-z energetic sources, such as X-ray binaries, accreting nuclear black holes, and shock heated interstellar medium. Adopting the model discussed in Eide et al. (2018), we find that these X-ray sources during the Epoch of Reionization (EoR) contribute less than a few per cent of the unresolved X-ray background. The same sources contribute to less than ∼2 per cent of the measured angular power spectrum of the fluctuations of the X-ray background. The outputs of radiative transfer simulations modeling the EoR are used to evaluate the cross-correlations of X-ray background with the 21 cm signal from neutral hydrogen. Such correlation could be used to confirm the origin of the 21 cm signal, as well as give information on the properties of the X-ray sources during the EoR. We find that the correlations are positive during the early stages of reionization when most of the hydrogen is neutral, while they become negative when the intergalactic medium gets highly ionized, with the transition from positive to negative depending on both the X-ray model and the scale under consideration. With SKA as the reference instrument for the 21 cm experiment, the predicted S/N for such correlations is 10−15ergcm−2s−1⁠, while the cumulative S/N from l = 1000 to 104 at xHI=0.5 is ∼5 if sources with observed flux >10−17ergcm−2s−1 are detected., The authors are grateful to James Aird, Andrea Merloni, Arne Rau, and Mara Salvato for inputs on planned X-ray surveys, Ruben Salvaterra for useful comments, and an anonimous referee for his or her comments. The tools for bibliographic research are offered by the NASA Astrophysics Data Systems and by the JSTOR archive. QM is supported by the National Basic Research Program (“973” Program) of China (grant No. 2014CB845800), the National Natural Science Foundation of China (grant Nos 11673068 and 11725314), the Youth Innovation Promotion Association (2011231), the Key Research Program of Frontier Sciences (QYZDB-SSW-SYS005), the Strategic Priority Research Program ‘Multiwaveband gravitational wave Universe’ (grant No. XDB23000000) of the Chinese Academy of Sciences. MBE thanks the Astronomy and Astrophysics department at UCSC for their kind hospitality and is a fellow of the U.S.-Norway Fulbright Foundation. KH acknowledges support from the Icelandic Research Fund, grant number 173728-051.

Observations of the 21 cm line from neutral hydrogen indicate that an Epoch of Heating (EoH) might have preceded the later Epoch of Reionization (EoR). Here we study the effects on the ionization state and the thermal history of the Intergalactic Medium (IGM) during the EoH induced by different assumptions on ionizing sources in the high redshift Universe: (i) stars, (ii) X-ray binaries (XRBs), (iii) thermal bremsstrahlung of the hot Interstellar Medium (ISM), and (iv) accreting nuclear black holes (BHs). To this aim, we post-process outputs from the ($100 h^{-1}$ cMpc)$^3$ hydrodynamical simulation MassiveBlack-II with the cosmological 3D radiative transfer code CRASH, which follows the propagation of UV and X-ray photons, computing the thermal and ionization state of hydrogen and helium through the EoH. We find that stars determine the fully ionized morphology of the IGM, while the spectrally hard XRBs pave way for efficient subsequent heating and ionization by the spectrally softer ISM. With the seeding prescription in MassiveBlack-II, BHs do not contribute significantly to either ionization or heating. With only stars, most of the IGM remains in a cold state (with a median $T=11$ K at $z=10$), however, the presence of more energetic sources raises the temperature of regions around the brightest and more clustered sources above that of the CMB, opening the possibility to observing the 21 cm signal in emission., 18 pages, 9 figures. Accepted for publication in MNRAS

In an effort to understand the cause of the apparent depletion in the number density of radio-loud AGNs at $z>3$, this work investigates the viability of the so-called Cosmic Microwave Background (CMB) quenching mechanism of intrinsically jetted, high-z AGNs, whereby Inverse Compton scattering of CMB photons off electrons within the extended lobes results in a substantial dimming of the lobe synchrotron emission at GHz frequencies, while simultaneously boosting their diffuse X-ray signal. We focus on five $z>3.5$ radio galaxies that have sufficiently deep Chandra exposure (> 50 ks) to warrant a meaningful investigation of any extended X-ray emission. For those objects with evidence for statistically significant extended X-ray lobes (4C 41.17 and 4C 03.24), we combine the Chandra measurements with literature data at lower frequencies to assemble the systems' Spectral Energy Distributions (SEDs), and utilize state-of-the-art SED modelling (Ghisellini et al. 2015) -- including emission from the disk, torus, jet, hotspots, and lobes -- to infer their physical parameters. For both radio galaxies, the magnetic energy density in the hotspots is found to exceed the energy density in CMB photons, wheres the opposite is true for the lobes. This implies that any extended synchrotron emission likely originates from the hotspots themselves, rather than the lobes. Conversely, Inverse Compton scattering of CMB photons dominates the extended X-ray emission from the lobes, which are effectively "radio-quenched". As a result, CMB quenching is effective in these systems in spite of the fact that the observed X-ray to radio luminosity ratio does not bear the signature $(1+z)^4$ dependence of the CMB energy density., MNRAS in press; 17 pages, 8 figures, 4 tables

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)., 19 pages, 19 figures, the version accepted in MNRAS

We present a detailed analysis of the ionization and thermal structure of the intergalactic medium (IGM) around a high-redshift QSO using a large suite of cosmological, multi-frequency radiative transfer (RT) simulations, exploring the contribution from galaxies as well as the QSO, and the effect of X-rays and secondary ionization. We show that in high-z QSO environments both the central QSO and the surrounding galaxies concertedly control the reionization morphology of hydrogen and helium and have a non-linear impact on the thermal structure of the IGM. A QSO imprints a distinctive morphology on H II regions if its total ionizing photon budget exceeds that of the surrounding galaxies since the onset of hydrogen reionization; otherwise, the morphology shows little difference from that of H II regions produced only by galaxies. In addition, the spectral shape of the collective radiation field from galaxies and QSOs controls the thickness of the I-fronts. While a UV-obscured QSO can broaden the I-front, the contribution from other UV sources, either galaxies or unobscured QSO, is sufficient to maintain a sharp I-front. X-rays photons from the QSO are responsible for a prominent extended tail of partial ionization ahead of the I-front. QSOs leave a unique imprint on the morphology of He II / He III regions. We suggest that, while the physical state of the IGM is modified by QSOs, the most direct test to understand the role of galaxies and QSOs during reionization is to perform galaxy surveys in a region of sky imaged by 21 cm tomography., 19 pages, 22 figures, submitted

We study the possibility to detect and distinguish signatures of enrichment from PopIII stars in observations of PopII GRBs (GRBIIs) at high redshift by using numerical N-body/hydrodynamical simulations including atomic and molecular cooling, star formation and metal spreading from stellar populations with different initial mass functions (IMFs), yields and lifetimes. PopIII and PopII star formation regimes are followed simultaneously and both a top-heavy and a Salpeter-like IMF for pristine PopIII star formation are adopted. We find that the fraction of GRBIIs hosted in a medium previously enriched by PopIII stars (PopIII-dominated) is model independent. Typical abundance ratios, such as [Si/O] vs [C/O] and [Fe/C] vs [Si/C], can help to disentangle enrichment from massive and intermediate PopIII stars, while low-mass first stars are degenerate with regular PopII generations. The properties of galaxies hosting PopIII-dominated GRBIIs are not very sensitive to the particular assumption on the mass of the first stars., 9 pages, 4 figures

We present the first theoretical study of metals in damped Ly-alpha (DLA) systems at redshift z=7. The features of cold, primordial gas are studied by means of N-body, hydro, chemistry simulations, including atomic and molecular non-equilibrium chemistry, cooling, star formation for population III and population II-I regimes, stellar evolution, cosmic metal spreading according to proper yields (for He, C, O, Si, Fe, Mg, S, etc.) and lifetimes, and feedback effects. Theoretical expectations are then compared to recently available constraints from DLA observations. We find that DLA galaxies at z=7 account for 10 per cent of the whole galaxy population and for most of the metal-poor galaxies at these epochs. About 7 per cent of these DLA galaxies contain purely pristine material and 34 per cent of them consist of very weakly polluted gas, being, therefore, suitable candidate s as population III sites. The remaining 59 per cent are enriched above ~10^{-4} Zsun. Additionally, DLA candidates appear to have: gas masses ~3x10^{20} cm^{-2}. They present no or weak correlations between their gas mass and Z, NHI, or x_mol; a moderate correlation between x_mol and Z, linked to the ongoing molecular-driven star formation and metal pollution processes; a mild anti-correlation between NHI and x_mol, due to H depletion into molecules; and a chemical content that is subject to environmental dependencies., Comment: 9 pages, 4 figures; MNRAS, in Press - "Rex sedet in vertice caveat ruinam nam sub axe legimus Hecubam reginam"

We propose to use high-redshift long $\gamma$-ray bursts (GRBs) as cosmological tools to constrain the amount of primordial non-Gaussianity in the density field. By using numerical, N-body, hydrodynamic, chemistry simulations of different cosmological volumes with various Gaussian and non-Gaussian models, we self-consistently relate the cosmic star formation rate density to the corresponding GRB rate. Assuming that GRBs are fair tracers of cosmic star formation, we find that positive local non-Gaussianities, described in terms of the non-linear parameter, \fnl, might boost significantly the GRB rate at high redshift, $z \gg 6$. Deviations with respect to the Gaussian case account for a few orders of magnitude if \fnl$\sim 1000$, one order of magnitude for \fnl$\sim 100$, and a factor of $\sim 2$ for \fnl$\sim 50$. These differences are found only at large redshift, while at later times the rates tend to converge. Furthermore, a comparison between our predictions and the observed GRB data at $z > 6$ allows to exclude large negative \fnl, consistently with previous works. Future detections of any long GRB at extremely high redshift ($z\sim 15-20$) could favor non-Gaussian scenarios with positive \fnl. More stringent constraints require much larger high-$z$ GRB complete samples, currently not available in literature. By distinguishing the contributions to the GRB rate from the metal-poor population III regime, and the metal-enriched population II-I regime, we conclude that the latter is a more solid tracer of the underlying matter distribution, while the former is strongly dominated by feedback mechanisms from the first, massive, short-lived stars, rather than by possible non-Gaussian fluctuations. This holds quite independently of the assumed population III initial mass function.

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/)., Comment: 38 pages, 11 figures, 5 tables, A&A in press, updated with editorial corrections

The 21 cm-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 21 cm emission from the EoR with the Low Frequency Array (LOFAR), and of high-redshift Ly α emitters 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 Ly α photons) simulations at redshift 6.68, 7.06 and 7.3 to compute the 3D theoretical 21 cm-galaxy cross-power spectrum and cross-correlation function, as well as to predict the 2D 21 cm-galaxy cross-power spectrum and cross- correlation function expected to be observed by LOFAR and HSC. Once noise and projection effects are accounted for, our predictions of the 21 cm-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.003 at z = 6.6 and p = 0.08 at z = 7.3. On smaller scales, instead, the signal is completely contaminated. On the other hand, our 21 cm-galaxy cross-correlation function is strongly contaminated by noise on all scales, since the noise is no longer being separated by its k modes.

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., Comment: 17 pages, 12 figures, 1 table. Accepted for publication in MNRAS. Replaced to match the accepted version

The redshifted 21 cm line of neutral hydrogen is a promising probe of the Epoch of Reionization (EoR). However, its detection requires a thorough understanding and control of the systematic errors. We study two systematic biases observed in the LOFAR EoR residual data after calibration and subtraction of bright discrete foreground sources. The first effect is a suppression in the diffuse foregrounds, which could potentially mean a suppression of the 21 cm signal. The second effect is an excess of noise beyond the thermal noise. The excess noise shows fluctuations on small frequency scales, and hence it can not be easily removed by foreground removal or avoidance methods. Our analysis suggests that sidelobes of residual sources due to the chromatic point spread function and ionospheric scintillation can not be the dominant causes of the excess noise. Rather, both the suppression of diffuse foregrounds and the excess noise can occur due to calibration with an incomplete sky model containing predominantly bright discrete sources. We show that calibrating only on bright sources can cause suppression of other signals and introduce an excess noise in the data. The levels of the suppression and excess noise depend on the relative flux of sources which are not included in the model with respect to the flux of modeled sources. We discuss possible solutions such as using only long baselines to calibrate the interferometric gain solutions as well as simultaneous multi-frequency calibration along with their benefits and shortcomings., Submitted to MNRAS, 14 pages, 11 figures

We attempt to interpret existing data on the evolution of the UV luminosity function and UV colours, $\beta$, of galaxies at $5 \leq z \leq 8$, to improve our understanding of their dust content and ISM properties. To this aim, we post-process the results of a cosmological hydrodynamical simulation with a chemical evolution model, which includes dust formation by supernovae and intermediate mass stars, dust destruction in supernova shocks, and grain growth by accretion of gas-phase elements in dense gas. We find that observations require a steep, Small Magellanic Cloud-like extinction curve and a clumpy dust distribution, where stellar populations younger than 15 Myr are still embedded in their dusty natal clouds. Investigating the scatter in the colour distribution and stellar mass, we find that the observed trends can be explained by the presence of two populations: younger, less massive galaxies where dust enrichment is mainly due to stellar sources, and massive, more chemically evolved ones, where efficient grain growth provides the dominant contribution to the total dust mass. Computing the IR-excess - UV color relation we find that all but the dustiest model galaxies follow a relation shallower than the Meurer et al. (1999) one, usually adopted to correct the observed UV luminosities of high-$z$ galaxies for the effects of dust extinction. As a result, their total star formation rates might have been over-estimated. Our study illustrates the importance to incorporate a proper treatment of dust in simulations of high-$z$ galaxies, and that massive, dusty, UV-faint galaxies might have already appeared at $z \lesssim 7$., Comment: accepted in MNRAS

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 perform a systematic study of how the inhomogeneities in the Inter-Galactic Medium (IGM) affect the observability of Lyman Alpha Emitters (LAEs) around the Epoch of Reionization. We focus on the IGM close to the galaxies as the detailed ionization distribution and velocity fields of this region could significantly influence the scattering of Ly-alpha photons off neutral H atoms as they traverse the IGM after escaping from the galaxy. We simulate the surface brightness (SB) maps and spectra of more than 100 LAEs at z=7.7 as seen by an observer at z=0. To achieve this, we extract the source properties of galaxies and their surrounding IGM from cosmological simulations of box sizes 5-30 Mpc/h and follow the coupled radiative transfer of ionizing and Ly-alpha radiation through the IGM using CRASH-alpha. We find that the simulated SB profiles are extended and their detailed structure is affected by inhomogeneities in the IGM, especially at high neutral fractions. The detectability of LAEs and the fraction of the flux observed depend heavily on the shape of the SB profile and the SB threshold (SB_th) of the observational campaign. Only ultradeep observations (e.g. SB_th ~ 10^-23 ergs/s/cm^2/arcsec^2) would be able to obtain the true underlying mass-luminosity relation and luminosity functions of LAEs. The details of our results depend on whether Ly-alpha photons are significantly shifted in the galaxy to longer wavelengths, the mean ionization fraction in the IGM and the clustering of ionizing sources. These effects can lead to an easier escape of Ly-alpha photons with less scattering in the IGM and a concentrated SB profile similar to the one of a point source. Finally, we show that the SB profiles are steeper at high ionization fraction for the same LAE sample which can potentially be observed from the stacked profile of a large number of LAEs.

The report that follows gives glimpses into the enormous amount of work in cosmology over the past three years, both in theory and in observations. It is impossible to do justice, within the space allotted, to all the exciting work and some important topics like structure formation could not be covered. Nevertheless, the contributors to different sections have done a conscientious job and I thank them all. Some of them had to step in to fill the gaps at short notice.