A new model framework for circumgalactic Lyα\alpha radiative transfer constrained by galaxy-Lyα\alpha forest clustering

We present a new perturbative approach to "constrained Ly$\alpha$ radiative transfer'" (RT) through the circum- and inter-galactic medium (CGM and IGM). We constrain the HI content and kinematics of the CGM and IGM in a physically motivated model, using the galaxy-Ly$\alpha$ forest clustering data from spectroscopic galaxy surveys in QSO fields at $z\sim2-3$. This enables us to quantify the impact of the CGM/IGM on Ly$\alpha$ emission in an observationally constrained, realistic cosmological environment. Our model predicts that the CGM and IGM at these redshifts transmit $\approx80~\%$ of Ly$\alpha$ photons after having escaped from galaxies. This implies that while the inter-stellar medium primarily regulates Ly$\alpha$ escape, the CGM has a non-negligible impact on the observed Ly$\alpha$ line properties and the inferred Ly$\alpha$ escape fraction, even at $z\sim 2-3$. Ly$\alpha$ scattering in the CGM and IGM further introduces an environmental dependence in the (apparent) Ly$\alpha$ escape fraction, and the observed population of Ly$\alpha$ emitting galaxies: the CGM/IGM more strongly suppresses direct Ly$\alpha$ emission from galaxies in overdense regions in the Universe, and redistributes this emission into brighter Ly$\alpha$ haloes. The resulting mean surface brightness profile of the Ly$\alpha$ haloes is generally found to be a power-law $\propto r^{-2.4}$. Although our model still contains arbitrariness, our results demonstrate how (integral field) spectroscopic surveys of galaxies in QSO fields constrain circumgalactic Ly$\alpha$ RT, and we discuss the potential of these models for studying CGM physics and cosmology.Comment: 20 pages, 14 figures, the version accepted in MNRA

The circumgalactic medium in Lyman-alpha: a new constraint on galactic outflow models

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<30$kpc). Our results suggest that galactic outflows affect Ly$\alpha$ absorption and emission around galaxies mostly at impact parameters $b<50$ kpc, while cold accretion flows dominate at larger distances. We discuss the implications of this result, and underline the potential constraining power of CGM observations - in emission and absorption - on galactic outflow models.Comment: 14 pages, 12 figure

The Epoch of IGM heating by early sources of X-rays

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.Comment: 18 pages, 9 figures. Accepted for publication in MNRA

Opening Reionization: Quantitative Morphology of the Epoch of Reionization and Its Connection to the Cosmic Density Field

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$.Comment: 18 pages, 15 figure, as accepted for publication by MNRA

Probing the high-z IGM with the hyperfine transition of 3^3He+^+

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 MNRA

Observing the Redshifted 21 cm Signal around a Bright QSO at z ∼ 10

We use hydrodynamics and radiative transfer simulations to study the 21 cm signal around a bright QSO at z ~ 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 ionizing and heating the surrounding gas. As expected, a longer QSO lifetime, t QSO, results in a 21 cm signal in emission located at increasingly larger angular radii, θ, and covering a wider range of θ. Similar features can be obtained with a higher galactic emissivity efficiency, f UV, such that determining the origin of a large ionized bubble (i.e., QSO versus stars) is not straightforward. Such degeneracy could be reduced by taking advantage of the finite light travel time effect, which is expected to affect an H ii 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 QSO could be detected by Square Kilometre Array1-low instrument with a high signal-to-noise ratio (S/N). As a reference, for t QSO = 10 Myr, a S/N ~ 8 is expected assuming that no pre-heating of the intergalactic medium has taken place due to high-z energetic sources, while it can reach values above 10 in cases 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 UV and t QSO

Cosmic Galaxy-IGM HI Relation at z∼2−3{\it{z}}\sim 2-3 Probed in the COSMOS/UltraVISTA 1.61.6 deg2^2 Field

We present spatial correlations of galaxies and IGM HI in the COSMOS/UltraVISTA 1.62 deg$^2$ field. Our data consist of 13,415 photo-$z$ galaxies at $z\sim2-3$ with $K_s<23.4$ and the Ly$\alpha$ forest absorptions in the background quasar spectra selected from SDSS data with no signature of damped Ly$\alpha$ system contamination. We estimate a galaxy overdensity $\delta_{gal}$ in an impact parameter of 2.5 pMpc, and calculate the Ly$\alpha$ forest fluctuations $\delta_{\langle F\rangle}$ whose negative values correspond to the strong Ly$\alpha$ forest absorptions. We identify weak evidence of an anti-correlation between $\delta_{gal}$ and $\delta_{\langle F\rangle}$ with a Spearman's rank correlation coefficient of $-0.39$ suggesting that the galaxy overdensities and the Ly$\alpha$ forest absorptions positively correlate in space at the $\sim90\%$ confidence level. This positive correlation indicates that high-$z$ galaxies exist around an excess of HI gas in the Ly$\alpha$ forest. We find four cosmic volumes, dubbed $A_{obs}$-$D_{obs}$, that have extremely large (small) values of $\delta_{gal} \simeq0.8$ ($-1$) and $\delta_{\langle F\rangle}$ $\simeq0.1$ ($-0.4$), three out of which, $B_{obs}$-$D_{obs}$, significantly depart from the correlation, and weaken the correlation signal. We perform cosmological hydrodynamical simulations, and compare with our observational results. Our simulations reproduce the correlation, agreeing with the observational results. Moreover, our simulations have model counterparts of $A_{obs}$-$D_{obs}$, and suggest that the observations pinpoint, by chance, a galaxy overdensity like a proto-cluster, gas filaments lying on the sightline, a large void, and orthogonal low-density filaments. Our simulations indicate that the significant departures of $B_{obs}$-$D_{obs}$ are produced by the filamentary large-scale structures and the observation sightline effects.Comment: 14 pages, 12 figures. Accepted for publication in Ap