Search

The 1D power spectrum P1D of the Ly α forest provides important information about cosmological and astrophysical parameters, including constraints on warm dark matter models, the sum of the masses of the three neutrino species, and the thermal state of the intergalactic medium. We present the first measurement of P1D with the quadratic maximum likelihood estimator (QMLE) from the Dark Energy Spectroscopic Instrument (DESI) survey early data sample. This early sample of 54 600 quasars is already comparable in size to the largest previous studies, and we conduct a thorough investigation of numerous instrumental and analysis systematic errors to evaluate their impact on DESI data with QMLE. We demonstrate the excellent performance of the spectroscopic pipeline noise estimation and the impressive accuracy of the spectrograph resolution matrix with 2D image simulations of raw DESI images that we processed with the DESI spectroscopic pipeline. We also study metal line contamination and noise calibration systematics with quasar spectra on the red side of the Ly α emission line. In a companion paper, we present a similar analysis based on the Fast Fourier Transform estimate of the power spectrum. We conclude with a comparison of these two approaches and discuss the key sources of systematic error that we need to address with the upcoming DESI Year 1 analysis.

We assess the impact of quasi-stellar objects (QSOs) on the high-redshift (⁠|$z\gt 4$|⁠) intergalactic medium using Monte Carlo realizations of QSO populations and the He  |$\rm \scriptstyle III$| regions they generate, applied to the Sherwood–Relics simulations, allowing for uncertainties in the QSO luminosity function, its evolution, and QSO spectra and ages. While QSO luminosity functions based on optical–infrared selection are unable to reproduce the broadening H  |$\rm \scriptstyle I$| Ly  |$\alpha$| optical depth distributions at |$z\gt 5$|⁠ , much broader distributions are found for the higher numbers of QSOs based on X-ray selection, suggesting a large QSO contribution to the ultraviolet background at |$z\gt 5$| may offer an alternative to late reionization models to account for the broad H  |$\rm \scriptstyle I$| Ly  |$\alpha$| optical depth distributions. Realizations using QSOs based on the higher QSO counts also much better recover the measured pixel flux autocorrelation function at |$z\gt 5$|⁠. The He  |$\rm \scriptstyle III$| regions from QSO sources according to both types of luminosity function suppress the pixel flux power spectrum on small scales, |$k\gt 0.02\, \mathrm{s\, km^{-1}}$|⁠ , while enhancing it on larger, both by amounts of up to tens of per cent at |$z\gt 4$|⁠ , with the magnitude increasing with the intergalactic He  |$\rm \scriptstyle III$| filling factor and the boost in temperature within the He  |$\rm \scriptstyle III$| regions. [ABSTRACT FROM AUTHOR]

Being the most prominent H  i line, Ly |$\alpha$| permeates the cosmic web in emission. Despite its potential as a cosmological probe, its detection on large scales remains elusive. We present a new methodology to perform Ly |$\alpha$| intensity mapping with broad-band optical images, by cross-correlating them with Ly |$\alpha$| forest data using a custom one-parameter estimator. We also develop an analytical large-scale Ly |$\alpha$| emission model with two parameters (average luminosity |$\langle L_{\rm Ly\alpha } \rangle$| and bias |$b_{\rm e}$|⁠) that respects observational constraints from quasar (QSO) luminosity functions. We compute a forecast for Dark Energy Camera Legacy Survey (DECaLS)/Beijing–Arizona Sky Survey (BASS) g -band images cross-correlated with Dark Energy Spectroscopic Instrument (DESI) Ly |$\alpha$| forest data, setting guidelines for reducing images into Ly |$\alpha$| intensity maps. Given the transversal scales of our cross-correlation (26.4 arcmin, |$\sim$| 33 cMpc h−1), our study effectively integrates Ly |$\alpha$| emission over all the cosmic volume inside the DESI footprint at |$2.2 \lt z \lt 3.4$| (the g -band Ly |$\alpha$| redshift range). Over the parameter space (⁠|$\langle L_{\rm Ly\alpha } \rangle$|⁠ , |$b_{\rm e}$|⁠) sampled by our forecast, we find a 3 |$\sigma$| of large-scale structure in Ly |$\alpha$| likely, with a probability of detection of 23.95 per cent for DESI–DECaLS/BASS, and 54.93 per cent for a hypothetical DESI phase-II with twice as much Ly |$\alpha$| QSOs. Without a detection, we derive upper bounds on |$\langle L_{\rm Ly\alpha } \rangle$| competitive with optimistic literature estimates [ |$2.3 \pm 1 \times 10^{\rm 41}$| erg s−1 (cMpc |$^3$|⁠)−1 for DESI, and |$\sim$| 35  per cent lower for its hypothetical phase-II]. Extrapolation to the DESI-Rubin overlap shows that a detection of large-scale structure with Ly |$\alpha$| intensity mapping using next-generation imaging surveys is certain. Such detection would allow constraining |$\langle L_{\rm Ly\alpha } \rangle$|⁠ , and explore the constraining power of Ly |$\alpha$| intensity mapping as a cosmological probe. [ABSTRACT FROM AUTHOR]

Strong low-ionization transitions like the Mg  ii |$\lambda 2796,2804$|  Å doublet are believed to produce a detectable 'metal-line forest', if metals pollute the neutral intergalatic medium (IGM). We measure the autocorrelation of the Mg  ii forest transmission using 10 ground-based |$z \ge 6.80$| quasar spectra probing the redshift range |$5.96 \lt z_{\rm Mg\, {\small II}} \lt 7.42$| (⁠|$z_{\rm Mg\, {\small II}, median} = 6.47$|⁠). The correlation function exhibits strong small-scale clustering and a pronounced peak at the doublet velocity (⁠|$\Delta v = 768~{\rm km\, s^{-1}}$|⁠) arising from discrete absorbers in the circumgalactic medium of galaxies. After these strong absorbers are identified and masked the signal is consistent with noise. Our measurements are compared to a suite of models generated by combining a large hydrodynamical simulation with a seminumerical reionization topology, assuming a simple uniform enrichment model. We obtain a 95 per cent credibility upper limit of |$[{\rm Mg/H}] \lt -3.73$| at |$z_{\rm Mg\, {\small II},median} = 6.47$|⁠ , assuming uninformative priors on [Mg/H] and the IGM neutral fraction |$x_{\rm {H\, {\small I}}}$|⁠. Splitting the data into low- z (⁠|$5.96 \lt z_{\rm Mg\, {\small II}} \lt 6.47$|⁠ ; |$z_{\rm Mg\, {\small II},median} = 6.235$|⁠) and high- z (⁠|$6.47 \lt z_{\rm Mg\, {\small II}} \lt 7.42$|⁠ ; |$z_{\rm Mg\, {\small II},median} = 6.72$|⁠) subsamples again yields null detections and 95 per cent upper limits of |$[{\rm Mg/H}] \lt -3.75$| and |$\lt -3.45$|⁠ , respectively. These first measurements set the stage for making the Mg ii forest an emerging tool to precisely constrain the Universe reionization and enrichment history. [ABSTRACT FROM AUTHOR]

Extended Lyα emission is routinely found around single quasars across cosmic time. However, few studies have investigated how such emission changes in fields with physically associated quasar pairs, which should reside in dense environments and are predicted to be linked through intergalactic filaments. We present VLT/MUSE snapshot observations (45 minutes/source) to unveil extended Lyα emission on scales of the circumgalactic medium (CGM) around the largest sample of physically associated quasar pairs to date, encompassing eight pairs (14 observed quasars) at z ∼ 3 with an i-band magnitude between 18 and 22.75, corresponding to absolute magnitudes Mi(z = 2) between −29.6 and −24.9. The pairs are either at close (∼50–100 kpc, five pairs) or wide (∼450–500 kpc, three pairs) angular separation and have velocity differences of Δv ≤ 2000 km s−1. We detected extended emission around 12 of the 14 targeted quasars and investigated the luminosity, size, kinematics, and morphology of these Lyα nebulae. On average, they span about 90 kpc and are 2.8 × 1043 erg s−1 bright. Irrespective of the quasars' projected distance, the nebulae often (∼45%) extend toward the other quasar in the pair, leading to asymmetric emission whose flux-weighted centroid is at an offset position from any quasar location. We show that large nebulae are preferentially aligned with the large-scale structure, as traced by the direction between the two quasars, and conclude that the cool gas (104 K) in the CGM traces well the direction of cosmic web filaments. Additionally, the radial profile of the Lyα surface brightness around quasar pairs can be described by a power law with a shallower slope (∼−1.6) with respect to single quasars (∼−2), indicative of increased CGM densities out to large radii and/or an enhanced contribution from the intergalactic medium (IGM) due to the dense environments expected around quasar pairs. The sample presented in this study contains excellent targets for ultra-deep observations to directly study filamentary IGM structures in emission. This work demonstrates that a large snapshot survey of quasar pairs will pave the way to direct statistical study of the IGM. [ABSTRACT FROM AUTHOR]

The extended Lyα haloes (LAHs) have been found to be prevalent around high-redshift star-forming galaxies. However, the origin of the LAHs is still a subject of debate. The spatially resolved analysis of Lyα profiles provides an important diagnostic. We analyse the average spatial extent and spectral variation of the circumgalactic LAHs by stacking a sample of 155 Lyα emitters (LAEs) at redshifts of 3 < z < 4 in the MUSE Extremely Deep Field. Our analysis reveals that, with respect to the Lyα line of the target LAE, the peak of the Lyα line at large distances becomes increasingly more blueshifted up to a projected distance of 60 kpc (≈3× virial radius), with a velocity offset of ≈250 km/s. This trend is evident in both the mean and median stacks, suggesting that it is a general property of our LAE sample, which typically has a Lyα luminosity of ≈1041.1 erg s−1. However, due to the absence of systemic redshift data, it remains unclear whether the Lyα line peak at large projected distances is less redshifted compared to the inner regions or truly blueshifted with respect to the systemic velocity. We explore various scenarios to explain the large-scale kinematics of the Lyα line. [ABSTRACT FROM AUTHOR]

We present a deep H  i survey at L band conducted with the Five-hundred-metre Aperture Spherical radio Telescope (FAST) over the COSMOS field. This survey is strategically designed to overlap with the MIGHTEE COSMOS field, aiming to combine the sensitivity of the FAST and high-resolution of the MeerKAT. We observed the field with FAST for approximately 11 hours covering |${\sim} 2$| square degrees, and reduced the raw data to H  i spectral cubes over the frequency range 1310–1420 MHz. The FAST-H  i data reach a median 3 |$\sigma$| column density of |$N_{\rm {H\,{{\small I}}}} \sim 2\times 10^{17}$|  cm |$^{-2}$| over a |${\sim} 5\, {\rm km}\, {\rm s}^{-1}$| channel width, allowing for studies of the distribution of H  i gas in various environments, such as in galaxies, the Circum-Galactic Medium (CGM) and Intergalactic Medium (IGM). We visually searched the spectral cubes for H  i sources, and found a total of 80 H  i detections, of which 56 have been cross-matched with the MIGHTEE-H  i catalogue. With the cross-matched sources, we compare their H  i masses and find that the total H  i mass fraction in the IGM and CGM surrounding the galaxy pairs is statistically higher than the H  i fraction surrounding the isolated galaxies by a difference of |${\sim} 13\pm 4$| per cent, indicating that the CGM and IGM associated with interacting systems are richer in neutral hydrogen compared to those around isolated galaxies in the local Universe. We also describe several FAST-MeerKAT synergy projects, highlighting the full potential of exploiting both single-dish and interferometric observations to study the distribution and evolution of the diffuse H  i gas. [ABSTRACT FROM AUTHOR]

We explore the potential of using the low-redshift Lyman |$\alpha$| (Ly  |$\alpha$|⁠) forest surrounding luminous red galaxies (LRGs) as a tool to constrain active galactic nuclei (AGNs) feedback models. Our analysis is based on snapshots from the Illustris and IllustrisTNG simulations at a redshift of |$z=0.1$|⁠. These simulations offer an ideal platform for studying the influence of AGN feedback on the gas surrounding galaxies, as they share the same initial conditions and underlying code but incorporate different feedback prescriptions. Both simulations show significant impacts of feedback on the temperature and density of the gas around massive haloes. Following our previous work, we adjusted the UV background in both simulations to align with the observed number density of Ly  |$\alpha$| lines (⁠|$\rm dN/dz$|⁠) in the intergalactic medium and study the Ly  |$\alpha$| forest around massive haloes hosting LRGs, at impact parameters (⁠|$r_{\perp }$|⁠) ranging from 0.1 to 100 pMpc. Our findings reveal that |$\rm dN/dz$|⁠ , as a function of |$r_{\perp }$|⁠ , is approximately 1.5 to 2 times higher in IllustrisTNG compared to Illustris up to |$r_{\perp }$| of |$\sim 10$| pMpc. To further assess whether existing data can effectively discern these differences, we search for archival data containing spectra of background quasars probing foreground LRGs. Through a feasibility analysis based on these data, we demonstrate that |${\rm dN/dz} (r_{\perp })$| measurements can distinguish between feedback models of IllustrisTNG and Illustris with a precision exceeding 16 |$\sigma$|⁠. This underscores the potential of |${\rm dN/dz} (r_{\perp })$| measurements around LRGs as a valuable benchmark observation for discriminating between different feedback models. [ABSTRACT FROM AUTHOR]

We present initial results from the Polarization Sky Survey of the Universe's Magnetism (POSSUM), analysing 22 817 Faraday rotation measures (RMs) with median uncertainties of 1.2 rad m |$^{-2}$| across 1520 deg2 to study magnetized gas associated with 55 nearby galaxy groups (⁠|$z\lesssim 0.025$|⁠) with halo masses between |$10^{12.5}$| and |$10^{14.0}$| M |$_\odot$|⁠. We identify two distinct gas phases: the intragroup medium (IGrM) within 0–2 splashback radii and the warm-hot intergalactic medium (WHIM) extending from 2 to 7 splashback radii. These phases enhance the standard deviation of residual (i.e. Galactic foreground RM-subtracted) RMs by |$6.9\pm 1.8$| rad m |$^{-2}$| and |$4.2 \pm 1.2$| rad m |$^{-2}$|⁠ , respectively. Estimated magnetic field strengths are several μ G within the IGrM and 0.1–1 μ G in the WHIM. We estimate the plasma |$\beta$| in both phases, and show that magnetic pressure might be more dynamically important than in the ICM of more massive clusters or sparse cosmic web filaments. Our findings indicate that 'missing baryons' in the WHIM likely extend beyond the gravitational radii of group-mass haloes to Mpc scales, consistent with large-scale, outflow-driven 'magnetized bubbles' seen in cosmological simulations. We demonstrate that RM grids are an effective method for detecting magnetized thermal gas at galaxy group interfaces and within the cosmic web. This approach complements X-ray and Sunyaev-Zel'dovich effect methods, and when combined with fast radio burst dispersion measures, data from the full POSSUM survey – comprising approximately a million RMs – will allow direct magnetic field measurements to further our understanding of baryon circulation in these environments and the magnetized universe. [ABSTRACT FROM AUTHOR]

The X-ray column density ( N HX ) of gamma-ray bursts (GRBs) can probe the local environment of their progenitors over a wide redshift range. Previous work has suggested an increasing trend as a function of redshift. The relevance of this current analysis relies on investigating the selection bias method, such as the effect of the X-ray spectrum in high-redshift GRBs, which complicates the measurement of small N HX ; this has yet to be fully evaluated or discussed elsewhere. In this work, we evaluated these effects through simulations to define appropriate observational limits in the N HX versus redshift plane. We then applied a one-sided nonparametric method developed by Efron and Petrosian. Within the framework of this method, we investigated the redshift dependence of N HX and the local distribution function. Our results show that the evolution of N HX with redshift firmly exists with a significance of more than four sigma and follows a power law of (1 + z) 1.39 (+ 0.22 , − 0.27) . Based on these analyses and previous studies, the GRB progenitor mass varies but is more massive in the high-redshift environment and has a higher gas column density. This suggests that part of the luminosity evolution of GRBs, which has been widely reported, may be due to the evolution of the progenitor's mass. Using the same method, we demonstrate that optically dark GRBs show a consistent evolution: (1 + z) 1.15 (+ 0.67 , − 0.83) . By applying the Kolmogorov–Smirnov (KS) test, it is shown that optically dark GRBs have statistically identical flux and photon index distributions compared to normal GRBs, but the N HX is systematically larger. This result suggests that the darkness of some GRB populations is not due to an intrinsic mechanism, but rather because a higher density surrounds them. [ABSTRACT FROM AUTHOR]

We present a new suite of late-end reionization simulations performed with aton-he , a revised version of the GPU-based radiative transfer code aton that includes helium. The simulations are able to reproduce the Ly  |$\alpha$| flux distribution of the E-XQR-30 sample of QSO absorption spectra at |$5 \lesssim z \lesssim 6.2$|⁠ , and show that a large variety of reionization models are consistent with these data. We explore a range of variations in source models and in the early-stage evolution of reionization. Our fiducial reionization history has a mid-point of reionization at |$z = 6.5$|⁠ , but we also explore an 'Early' reionization history with a mid-point at |$z = 7.5$| and an 'Extremely Early' reionization history with a mid-point at |$z = 9.5$|⁠. Haloes massive enough to host observed Ly  |$\alpha$| emitters are highly biased. The fraction of such haloes embedded in ionized bubbles that are large enough to allow high Ly  |$\alpha$| transmission becomes close to unity much before the volume filling factor of ionized regions. For our fiducial reionization history this happens at |$z = 8$|⁠ , probably too late to be consistent with the detection by JWST of abundant Ly  |$\alpha$| emission out to |$z = 11$|⁠. A reionization history in our 'Early' model or perhaps even our 'Extremely Early' model may be required, suggesting a Thomson scattering optical depth in tension with that reported by Planck, but consistent with recent suggestions of a significantly higher value. [ABSTRACT FROM AUTHOR]

The escaping ionising efficiency from galaxies, fescξion, is a crucial ingredient for understanding their contribution to hydrogen reionisation, but both of its components, fesc and ξion, are extremely difficult to measure. We measured the average escaping ionising efficiency ⟨fescξion⟩ of galaxies at z = 5 implied by the mean level of ionisation in the intergalactic medium via the Lyman-α forest. We used the fact that Nion = ρUVfescξion, N ˙ ion = ρ UV f esc ξ ion the product of the ionising output and the UV density ρUV, can be calculated from the known average strength of the UV background and the mean free path of ionising photons. These quantities, as well as ρUV, are robustly measured at z ≤ 6. We calculated the missing factor of ⟨fescξion⟩ at z = 5 during a convenient epoch after hydrogen reionisation had been completed and the intergalactic medium had reached ionisation equilibrium but before bright quasars began to dominate the ionising photon production. Intuitively, our constraint corresponds to the required escaping ionising production from galaxies in order to avoid over- or under-ionising the Lyman-α forest. We obtained a measurement of log⟨fescξion⟩/erg Hz-1 = 24.28+0.21-0.20 − 1 = 24. 28 − 0.20 + 0.21 at z = 5 when integrating the ρUV down to a limiting magnitude Mlim = −11. Our measurement of the escaping ionising efficiency of galaxies is in rough agreement with both observations of early galaxies and with most models. [ABSTRACT FROM AUTHOR]

Context. The recent discovery of so-called mega radio halos as a new class of diffuse, steep-spectrum radio sources in clusters of galaxies has raised questions about the origin and the evolution of cluster-wide radio emission. Aims. We investigate whether the formation mechanisms of radio halos and mega radio halos differ, or whether they can be produced by different modalities of the same (re-)acceleration mechanism. Here we present results of a cosmological simulation of a disturbed galaxy cluster, with the aim to study the origin of mega radio halos. Methods. We analysed the evolution of cosmic-ray electrons (CRe), subject to gains and losses using a Fokker-Planck solver. In particular, we included the effects of adiabatic stochastic acceleration (ASA) which is caused by the stochastic interaction of cosmic rays with diffusing magnetic field lines in super-Alfvenic turbulence. Moreover, we included shock acceleration and the seeding of CRe by galaxies. Results. Our simulations generate cluster-scale radio sources during mergers, with properties that are in agreement with those observed for real radio halos. Furthermore, we find evidence of additional emission on larger scales. This emission resembles the radial distribution and the spectrum of a mega radio halo, but only when viewed close to the merger axis. Conclusions. In our simulation, the mechanism responsible for the formation of diffuse radio emission, both in the form of classical and mega radio halos, is cosmic-ray re-acceleration by turbulence. This turbulence is more solenoidal and more subsonic in the classical radio halo region, than in the mega radio halo region. [ABSTRACT FROM AUTHOR]

Context. The amplification of astrophysical magnetic fields takes place via dynamo instability in turbulent environments. Vorticity is usually present in any dynamo, but its role is not yet fully understood. Aims. This work is an extension of previous research on the effect of an irrotational subsonic forcing on a magnetized medium in the presence of rotation or a differential velocity profile. We aim to explore a wider parameter space in terms of Reynolds numbers, the magnetic Prandtl number, the forcing scale, and the cooling timescale in a Newtonian cooling. We studied the effect of imposing that either the acceleration or the velocity forcing function be curl-free and evaluated the terms responsible for the evolution vorticity. Methods. We used direct numerical simulations to solve the fully compressible, resistive magnetohydrodynamic equations with the Pencil Code. We studied both isothermal and non-isothermal regimes and addressed the relative importance of different vorticity source terms. Results. We report no small-scale dynamo for the models that do not include shear. We find a hydro instability, followed by a magnetic one, when a shearing velocity profile is applied. The vorticity production is found to be numerical in the purely irrotational acceleration case. Non-isothermality, rotation, shear, and density-dependent forcing, when included, contribute to increasing the vorticity. Conclusions. As in our previous study, we find that turbulence driven by subsonic expansion waves can amplify the vorticity and magnetic field only in the presence of a background shearing profile. The presence of a cooling function makes the instability occur on a shorter timescale. We estimate critical Reynolds and magnetic Reynolds numbers of 40 and 20, respectively. [ABSTRACT FROM AUTHOR]

We present a search for galaxies in the local Universe with extremely low oxygen abundance, that is, more than 25 times lower than solar, which corresponds to 12 + log(O/H) < 7.3. To determine the oxygen abundance, we apply the direct Te method for objects where the [O III]λ4363 line is detected. We identified 21 extremely metal-poor galaxies in the early data release of the Dark Energy Spectroscopic Instrument (DESI EDR), for some of which we also derived N/O, Ne/O, Ar/O, and S/O ratios. We find that many DESI galaxies with extremely low oxygen abundance exhibit a higher N/O ratio in comparison to the reference low-metallicity sample collected from the literature. We suggest that the elevation in N/O ratio may be explained by a contamination with metal-rich gas caused by gas inflow or a merger event. Moreover, contrary to some recent studies, we find that Ar/O and S/O ratios are enhanced as well, while the Ne/O ratio does not show such elevation. One of the galaxies, J0713+5608, has a remarkably low oxygen abundance of 6.978 ± 0.095 dex. This measurement aligns with the lowest known oxygen abundances in galaxies to date. Given the relatively high uncertainty, this galaxy may have the lowest oxygen abundance ever found. Additionally, J0713+5608 exhibited an enhanced N/O ratio compared to the typical N/O ratio observed in metal-poor galaxies within the local Universe. [ABSTRACT FROM AUTHOR]

Aims. Our aim is to present a robust parameter estimation with simulated Lyα forest spectra from Sherwood-Relics simulations suite by using an information-maximizing neural network (IMNN) to extract maximal information from Lyα 1D-transmitted flux in the Fourier space. Methods. We performed 1D estimations using IMNN for intergalactic medium (IGM) thermal parameters T0 and γ at z = 2 − 4, and cosmological parameters σ8 and ns at z = 3 − 4. We compared our results with estimates from the power spectrum using the posterior distribution from a Markov chain Monte Carlo (MCMC). We then checked the robustness of IMNN estimates against deviation in spectral noise levels, continuum uncertainties, and instrumental smoothing effects. Using mock Lyα forest sightlines from the publicly available CAMELS project, we also checked the robustness of the trained IMNN on a different simulation. As a proof of concept, we demonstrated a 2D-parameter estimation for T0 and H I photoionization rates, ΓHI. Results. We obtain improved estimates of T0 and γ using IMNN over the standard MCMC approach. These estimates are also more robust against signal-to-noise deviations at z = 2 and 3. At z = 4, the sensitivity to noise deviations is on par with MCMC estimates. The IMNN also provides T0 and γ estimates that are robust against continuum uncertainties by extracting small-scale continuum-independent information from the Fourier domain. In the cases of σ8 and ns, the IMNN performs on par with MCMC but still offers a significant speed boost in estimating parameters from a new dataset. The improved estimates with IMNN are seen for high instrumental resolution (FWHM = 6 km s−1). At medium or low resolutions, the IMNN performs similarly to MCMC, suggesting an improved extraction of small-scale information with IMNN. We also find that IMNN estimates are robust against the choice of simulation. By performing a 2D-parameter estimation for T0 and ΓHI, we also demonstrate how to take forward this approach observationally in the future. [ABSTRACT FROM AUTHOR]

Context. Cosmic filaments are observationally hard to detect. However, hydrodynamical cosmological simulations are ideal laboratories where the evolution of the cosmic web can be studied, and they allow for easier insight into the nature of the filaments. Aims. We investigate how the intrinsic properties of filaments are evolving in areas extracted from a larger cosmological simulation. We aim to identify significant trends in the properties of the warm-hot intergalactic medium (WHIM) and suggest possible explanations. Methods. To study the filaments and their contents, we selected a subset of regions from the Dianoga simulation. We analysed these regions that were simulated with different baryon physics, namely with and without AGN feedback. We constructed the cosmic web using the subspace constrained mean shift (SCMS) algorithm and the sequential chain algorithm for resolving filaments (SCARF). We examined the basic physical properties of filaments (length, shape, mass, radius) and analysed different gas phases (hot, WHIM, and colder gas components) within those structures. The evolution of the global filament properties and the properties of the gas phases were studied in the redshift range 0 < z < 1.48. Results. Within our simulations, the detected filaments have, on average, lengths below 9 Mpc. The filaments' shape correlates with their length, as the longer they are, the more likely they are curved. We find that the scaling relation between mass M and length L of the filaments is well described by the power law M ∞ L1.7. The radial density profile widens with redshift, meaning that the radius of the filaments becomes larger over time. The fraction of gas mass in the WHIM phase does not depend on the model and rises towards lower redshifts. However, the included baryon physics has a strong impact on the metallicity of gas in filaments, indicating that the AGN feedback impacts the metal content already at redshifts of z ~ 2. [ABSTRACT FROM AUTHOR]

Context. Fast radio bursts (FRBs) are poised to become important cosmological tools in the near future, as the number of observed FRBs is increasing rapidly with multiple surveys underway. A large sample of FRBs will soon have available dispersion measures (DMs) and rotation measures (RMs), which can be used to study the cosmic baryon density and the intergalactic magnetic field. However, the observed DM and RM of FRBs consists of multiple contributions that must be quantified to estimate the DM and RM of the intergalactic medium (IGM). Aims. In this paper, we estimate one such contribution to DM and RM, namely, of FRB host galaxies. We show how this contribution changes with redshift, galaxy type, and the stellar mass of the galaxies. We also investigate its dependence on galaxy inclination and on an FRB's offset from the center of the galaxy. Methods. Using the TNG50 simulation of the IllustrisTNG project, we selected 16 500 galaxies at redshifts of 0≤ ɀ ≤2, with stellar masses in the range of 9 ≤ log(M*/M⊙) ≤ 12. In each galaxy, we calculated the DM and RM contributions of 1000 sightlines; from these, we constructed the DM and RM probability density functions (PDFs). Results. We find that the rest frame DM distributions of all galaxies at a given redshift can be fitted by a log normal function and its median and width increase as a function of redshift. The rest-frame RM distribution is symmetric, with a median RMhost,rf=0 rad m–2 and it can be fitted by a combination of a Lorentzian and two Gaussian functions. The redshift evolution of the distribution width can be fitted by a curved power law. The parameters of these functions change for different subsets of galaxies with different stellar mass, inclination, and FRB offset. These changes are due to an increasing ne with redshift, SFR, and stellar mass. We do find a more ordered B field at lower ɀ compared to higher ɀ, as suggested by the presence of more galaxies with B field reversals and B fields dominated by random B field at higher ɀ. Conclusions. We estimated the FRB host DM and RM contributions, which can be used in the future to isolate the IGM contribution from the observed DM and RM of FRBs. We predict that to constrain a σRM,IGM of 2 rad m–2 to the 95% confidence level, we would need to observe 95 000 FRBs at ɀ = 0.5, but only 9 500 FRBs at ɀ = 2. [ABSTRACT FROM AUTHOR]

Enlightening our understanding of the first galaxies responsible for driving reionization requires detecting the 21-cm signal from neutral hydrogen. Interpreting the wealth of information embedded in this signal requires Bayesian inference. Parameter inference from the 21-cm signal is primarily restricted to the spherically averaged power spectrum (1D PS) owing to its relatively straightforward derivation of an analytic likelihood function enabling traditional Monte Carlo Markov Chain approaches. However, in recent years, simulation-based inference (SBI) has become feasible which removes the necessity of having an analytic likelihood, enabling more complex summary statistics of the 21-cm signal to be used for Bayesian inference. In this work, we use SBI, specifically marginal neural ratio estimation to learn the likelihood-to-evidence ratio with swyft , to explore parameter inference using the cylindrically averaged 2D PS. Since the 21-cm signal is anisotropic, the 2D PS should yield more constraining information compared to the 1D PS which isotropically averages the signal. For this, we consider a mock 1000 h observation of the 21-cm signal using the Square Kilometre Array and compare the performance of the 2D PS relative to the 1D PS. Additionally, we explore two separate foreground mitigation strategies, perfect foreground removal and wedge avoidance. We find the 2D PS outperforms the 1D PS by improving the marginalized uncertainties on individual astrophysical parameters by up to |$\sim 30\!-\!40$|  per cent irrespective of the foreground mitigation strategy. Primarily, these improvements stem from how the 2D PS distinguishes between the transverse, |$k_{\perp }$|⁠ , and redshift-dependent, |$k_{\parallel }$|⁠ , information which enables greater sensitivity to the complex reionization morphology. [ABSTRACT FROM AUTHOR]

Detecting the 21-cm signal at |$z\gtrsim 6$| will reveal insights into the properties of the first galaxies responsible for driving reionization. To extract this information, we perform parameter inference with three-dimensional simulations of the 21-cm signal embedded within a Bayesian inference pipeline. Presently, when performing inference, we must choose which sources of uncertainty to sample and which to hold fixed. Since the astrophysics of galaxies is much more uncertain than that of the underlying halo-mass function (HMF), we typically parametrize and model the former while fixing the latter. However, doing so may bias our inference of the galaxy properties. In this work, we explore the consequences of assuming an incorrect HMF and quantify the relative biases on our inferred astrophysical model parameters when considering the wrong HMF. We then relax this assumption by constructing a generalized five parameter HMF model and simultaneously recover it with our underlying astrophysical model. For this, we use 21cmfast and perform simulation-based inference using marginal neural ratio estimation to learn the likelihood-to-evidence ratio with Swyft. Using a mock 1000-h observation of the 21-cm power spectrum from the forthcoming Square Kilometre Array, conservatively assuming foreground wedge avoidance, we find that assuming the incorrect HMF can bias the recovered astrophysical parameters by up to |$\sim 3{\!-\!}4\sigma$| even when including independent information from observed luminosity functions. Using our generalized HMF model, although we recover our astrophysical parameters with a factor of |$\sim 2{\!-\!}4$| larger marginalized uncertainties, the constraints are unbiased, agnostic to the underlying HMF and therefore more conservative. [ABSTRACT FROM AUTHOR]

Multiple observations now suggest that the hydrogen reionization may have ended well below redshift six. While there has previously been no conclusive proof of extended neutral islands in the |$z \lt 6$| intergalactic medium, it is possible that such islands give rise to the giant Ly  |$\alpha$| absorption troughs seen in the spectra of high-redshift quasars. Here, we present evidence that the deepest and longest known Ly  |$\alpha$| trough at |$z \,\lt\, 6$|⁠ , towards ULAS J0148 + 0600 (J0148), is associated with damping wing absorption. The evidence comes from a window of strong Ly  |$\alpha$| transmission at the edge of the J0148 proximity zone. We show that the relatively smooth profile of this transmission window is highly unlikely to arise from resonant absorption alone, but is consistent with the presence of a damping wing. We further argue that the damping wing is unlikely to arise from a compact source due to the lack of associated metal lines, and is more likely to arise from an extended neutral island associated with the giant Ly  |$\alpha$| trough. We investigate the physical conditions that may give rise to the strong transmission window, and speculate that it may signal an usually deep void, nearby ionizing sources, and/or the recent passage of an ionization front. [ABSTRACT FROM AUTHOR]

Recent quasar absorption line observations suggest that reionization may end as late as |$z \approx 5.3$|⁠. As a means to search for large neutral hydrogen islands at |$z\ \lt\ 6$|⁠ , we revisit long dark gaps in the Ly  |$\beta$| forest in Very Large Telescope/X-Shooter and Keck/Echellette Spectrograph and Imager quasar spectra. We stack the Ly  |$\alpha$| forest corresponding to both edges of these Ly  |$\beta$| dark gaps and identify a damping wing-like extended absorption profile. The average redshift of the stacked forest is |$z=5.8$|⁠. By comparing these observations with reionization simulations, we infer that such a damping wing-like feature can be naturally explained if these gaps are at least partially created by neutral islands. Conversely, simulated dark gaps lacking neutral hydrogen struggle to replicate the observed damping wing features. Furthermore, this damping wing-like profile implies that the volume-averaged neutral hydrogen fraction must be |$\langle x_{\rm H\,{\small {I}}} \rangle \ge 6.1 \pm 3.9~{{\ \rm per\ cent}}$| at |$z = 5.8$|⁠. Our results offer robust evidence that reionization extends below |$z=6$|⁠. [ABSTRACT FROM AUTHOR]

Measurements of the mean free path of Lyman-continuum photons in the intergalactic medium during the epoch of reionization can help constrain the nature of the sources as well as sinks of hydrogen-ionizing radiation. A recent approach to this measurement has been to utilize composite spectra of multiple quasars at |$z\sim 6$|⁠ , and infer the mean free path after correcting the spectra for the presence of quasar proximity zones. This has revealed not only a steep drop in the mean free path from |$z=5$| to |$z=6$|⁠ , but also potentially a mild tension with reionization simulations. We critically examine such direct measurements of the mean free path for biases due to quasar environment, incomplete reionization, and quasar proximity zones. Using cosmological radiative transfer simulations of reionization combined with one-dimensional radiative transfer calculations of quasar proximity zones, we find that the bias in the mean free path due to overdensities around quasars is minimal at |$z\sim 6$|⁠. Patchiness of reionization at this redshift also does not affect the measurements significantly. Fitting our model to the data results in a mean free path of |$\lambda _{\mathrm{mfp}}=1.49^{+0.47}_{-0.52}$|  pMpc at |$z=6$|⁠ , which is consistent with the recent measurements in the literature, indicating robustness with respect to the modelling of quasar proximity zones. We also compare various ways in which the mean free path has been defined in simulations before the end of reionization. Overall, our finding is that recent measurements of the mean free path appear to be robust relative to several sources of potential bias. [ABSTRACT FROM AUTHOR]

ABSTRACT: At z ≲ 1, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas (T ∼ 104 K) in the intergalactic medium (IGM) into warm–hot phase (WHIM) with T > 105 K, resulting in a significant deviation from the previously tight power-law IGM temperature–density relationship, $T=T_0 (\rho / {\bar{\rho }})^{\gamma -1}$. This study explores the impact of the WHIM on measurements of the low-z IGM thermal state, [T0, γ], based on the b–$N_{{\rm H\,{\small I}}} $ distribution of the Ly α forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [T0, γ] can still be reliably measured from the b–$N_{{\rm H\,{\small I}}} $ distribution at z = 0.1, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at z = 0.1. The results suggest that the underlying [T0, γ] of both simulations can be recovered with biases as low as |Δlog (T0/K)| ≲ 0.05 dex, |Δγ| ≲ 0.1, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38 per cent, IllustrisTNG 10 per cent, and Nyx 4 per cent), we conclude that the b–$N_{{\rm H\,{\small I}}} $ distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Ly α absorbers, and discover that although their T and Δ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.

ABSTRACT: The Ly α forest (LAF) at z > 5 probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-z, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting for a late end to reionization, small structures down to the sub-kpc level may survive to later times than conventionally thought due to the reduced time for pressure smoothing to impact the gas, further tightening simulation resolution requirements. We perform a suite of simulations using the Eulerian cosmological hydrodynamics code Nyx, spanning domain sizes of 1.25 − 10 h−1 Mpc  and 5 − 80 h−1 kpc  cells, and explore the interaction of these variables with the timing of reionization on the properties of the matter distribution and the simulated LAF at z = 5.5. In observable Ly α power, convergence within 10 per cent is achieved for k < 0.1 s km–1, but larger k shows deviation of up to 20 per cent. While a later reionization retains more small structure in the density field, because of the greater thermal broadening there is little difference in the convergence of LAF power between early (z = 9) and later (z = 6) reionizations. We conclude that at z ∼ 5.5, resolutions of 10 kpc are necessary for convergence of LAF power at k < 0.1 s km–1, while higher-k modes require higher resolution, and that the timing of reionization does not significantly impact convergence given realistic photoheating.

The intergalactic medium in the vicinity of galaxy protoclusters are interesting testbeds to study complex baryonic effects such as gravitational shocks and feedback. Here, we utilize hydrodynamical simulations from the SIMBA and The Three Hundred suites to study the mechanisms influencing large-scale Lyman- |$\alpha$| transmission in |$2\lt z\lt 2.5$| protoclusters. We focus on the matter overdensity-Lyman- |$\alpha$| transmission relation |$(\delta _m-\delta _F)$| on Megaparsec-scales in these protoclusters, which is hypothesized to be sensitive to the feedback implementations. The lower density regions represented by the SIMBA-100 cosmological volume trace the power-law |$\delta _m-\delta _F$| relationship often known as the fluctuating Gunn–Peterson approximation. This trend is continued into higher density regions covered by simulations that implement stellar feedback only. Simulations with active galactic nucleus (AGN) thermal and AGN jet feedback, however, exhibit progressively more Lyman- |$\alpha$| transmission at fixed matter overdensity. Compared with the seven protoclusters observed in the COSMOS field, only two display the excess absorption expected from protoclusters. The others exhibit deviations: four show some increased transparency suggested by AGN X-ray thermal feedback models while the highly transparent COSTCO-I protocluster appears to reflect intense jet feedback. Discrepancies with the stellar-feedback-only model suggests processes at play beyond gravitational heating and/or stellar feedback as the cause of the protocluster transparencies. Some form of AGN feedback is likely at play in the observed protoclusters, and possibly long-ranged AGN jets in the case of COSTCO-I. While more detailed and resolved simulations are required to move forward, our findings open new avenues for probing AGN feedback at Cosmic Noon. [ABSTRACT FROM AUTHOR]

Recent evidence supporting reionization ending at redshift |$z \lesssim 6$| includes the rapid redshift evolution of the mean free path, |$\lambda _{\rm mfp}$|⁠ , for Lyman-limit photons through the intergalactic medium (IGM) at |$5\lt z\lt 6$|⁠. Here, we investigate |$\lambda _{\rm mfp}$| predicted by the Sherwood–Relics suite of hybrid radiation hydrodynamical simulations. Simulations with comoving volumes of |$40^{3}\, h^{-3}\rm \, cMpc^{3}$| (⁠|$160^{3}\, h^{-3}\rm \, cMpc^{3}$|⁠), calibrated to match the observed Ly |$\alpha$| forest transmission with a late end to reionization at |$z\lt 6$|⁠ , are consistent with recent |$\lambda _{\rm mfp}$| measurements at |$z\lt 5.9$|⁠ , and are |$1.2\sigma$| (⁠|$1.8\sigma$|⁠) above the highest redshift |$\lambda _{\rm mfp}$| measurement at |$z=5.93$|⁠. The majority of the Lyman-limit opacity at the end of reionization is attributable to highly ionized Ly |$\alpha$| forest absorbers with neutral hydrogen column densities |$N_{\rm HI} \leq 10^{16}$| – |$10^{17}\rm \, cm^{-2}$|⁠. Resolving these systems is critical for capturing the redshift evolution of |$\lambda _{\rm mfp}$| during the final stages of reionization. After reionization completes, overdense gas will reduce |$\lambda _{\rm mfp}$| by up to 20 per cent around haloes with masses |$M_{\rm h}\sim 10^{9}$| – |$10^{11}h^{-1}\, {\rm M}_{\odot }$|⁠ , but during reionization ionized bubbles will instead boost |$\lambda _{\rm mfp}$| around haloes by up to an order of magnitude when the IGM is as much as 90 per cent neutral by volume. This effect will play an important role in the visibility of Ly |$\alpha$| emitting galaxies at |$z\gtrsim 10$| discovered with JWST. [ABSTRACT FROM AUTHOR]

We present new measurements of the mean transmitted flux in the hydrogen |$\rm {Ly\,\alpha }$| and a relative transmitted flux measurement in |$\rm {Ly\,\beta }$| using 27 008 quasar spectra from the Fourteenth Data Release (DR14) of the Extended Baryon Oscillation Spectroscopic Survey. Individual spectra are first combined into 16 composites with mean redshifts in the range of |$2.8\lt z\lt 4.9$|⁠. We then apply Markov Chain Monte Carlo (MCMC) inference to produce a piecewise fit of the effective |$\tau _{\rm {Ly\,\alpha }}$| (corrected for metal lines and optically thick absorption) assuming a spline point distribution. We also perform a relative |$\Delta \tau _{\rm {Ly\,\beta }}$| measurement with the same data set, finding |$\Delta \tau _{\rm {Ly\,\beta }}\lt 0.35$| at |$z\lt 4.8$|⁠. The 6–8  per cent precision measurements in the rest frame 1075–1150 Å at |${\it z} \ \lt \rm {4.0}$| and 10–12 per cent precision measurements in the same region at |${\it z} \ \gt \rm {4.0}$| on |$\tau _{\rm {Ly\,\alpha }}$|⁠ , and our determinations of |$\Delta \tau _{\rm {Ly\,\beta }}$|⁠ , are dominated by systematic errors, likely arising from bias and uncertainties in estimates of the quasar continuum. Our |$\tau _{\rm {Ly\,\alpha }}$| values show a smooth increase by a factor of 5 over the redshift range |$z=2.4-4.4$|⁠. [ABSTRACT FROM AUTHOR]

Spectra of the highest redshift galaxies taken with JWST are now allowing us to see into the heart of the reionization epoch. Many of these observed galaxies exhibit strong damping wing absorption redward of their Lyman- |$\alpha$| emission. These observations have been used to measure the redshift evolution of the neutral fraction of the intergalactic medium and sizes of ionized bubbles. However, these estimates have been made using a simple analytic model for the intergalactic damping wing. We explore the recent observations with models of inhomogeneous reionization from the Sherwood–Relics simulation suite. We carry out a comparison between the damping wings calculated from the simulations and from the analytic model. We find that although the agreement is good on the red side of the Lyman- |$\alpha$| emission, there is a discrepancy on the blue side due to residual neutral hydrogen present in the simulations, which saturates the intergalactic absorption. For this reason, we find that it is difficult to reproduce the claimed observations of large bubble sizes at |$z \sim 7$|⁠ , which are driven by a detection of transmitted flux blueward of the Lyman- |$\alpha$| emission. We suggest instead that the observations can be explained by a model with smaller ionized bubbles and larger intrinsic Lyman- |$\alpha$| emission from the host galaxy. [ABSTRACT FROM AUTHOR]

In order to theoretically understand dust properties in the circumgalactic medium (CGM), we construct a dust evolution model that incorporates the evolution of grain size distribution. We treat the galaxy and the CGM as separate one-zone objects, and consider the mass exchange between them. We take into account dust production and interstellar dust processing for the galaxy based on our previous models, and newly incorporate sputtering in the hot phase and shattering in the cool phase for the CGM. We find that shattering increases the dust destruction (sputtering) efficiency in the CGM. The functional shape of the grain size distribution in the CGM evolves following that in the galaxy, but it is sensitive to the balance between sputtering and shattering in the CGM. For an observational test, we discuss the wavelength dependence of the reddening in the CGM traced by background quasar colors, arguing that, in order to explain the observed reddening level, a rapid inflow from the CGM to the galaxy is favored because of quick dust/metal enrichment. Small grain production by shattering in the CGM also helps to explain the rise of dust extinction toward short wavelengths. [ABSTRACT FROM AUTHOR]

Recently, models of the quasar luminosity function (QLF) rooted on large observational compilations have been produced that, unlike their predecessors, feature a smooth evolution with time. This bypasses the need to assume an ionizing emissivity evolution when simulating helium reionization with observations-based QLF, thus yielding more robust constraints. We combine one such QLF with a cosmological hydrodynamical simulation and 3D multifrequency radiative transfer. The simulated reionization history is consistently delayed in comparison to most other models in the literature. The predicted intergalactic medium temperature is larger than the observed one at |$z \lesssim 3$|⁠. Through forward modelling of the He  ii Lyman  |$\alpha$| forest, we show that our model produces an extended helium reionization and successfully matches the bulk of the observed effective optical depth distribution, although it overionizes the Universe at |$z\lesssim 2.8$| as the effect of small-scale Lyman Limit Systems not being resolved. We thoroughly characterize transmission regions and dark gaps in He  ii Lyman  |$\alpha$| forest sightlines. We quantify their sensitivity to the helium reionization, opening a new avenue for further observational studies of this epoch. Finally, we explore the implications for helium reionization of the large number of active galactic nuclei revealed at |$z\gtrsim 5$| by JWST. We find that such modifications do not affect any observable at |$z\le 4$|⁠ , except in our most extreme model, indicating that the observed abundance of high- z AGNs does not bear consequences for helium reionization. [ABSTRACT FROM AUTHOR]

Different star formation models at Cosmic Dawn produce detectable signatures in the observables of upcoming 21-cm experiments. In this work, we consider the physical scenario of feedback-free starbursts (FFB), according to which the star formation efficiency (SFE) is enhanced in sufficiently massive haloes at early enough times, thus explaining the indication from the JWST for an excess of bright galaxies at |$z \ge 10$|⁠. We model the contribution of FFBs to popII SFE and compute the impact these have on the 21-cm global signal and power spectrum. We show that FFBs affect the evolution of the brightness temperature and the 21-cm power spectrum, but they only have a limited effect on the neutral hydrogen fraction. We investigate how the observables are affected by changes in the underlying star formation model and by contribution from popIII stars. Finally, we forecast the capability of next-generation Hydrogen Epoch of Reionization Array (HERA) to detect the existence of FFB galaxies via power spectrum measurements. Our results show the possibility of a significant detection, provided that popII stars are the main drivers of lowering the spin temperature. Efficient popIII star formation will make the detection more challenging. [ABSTRACT FROM AUTHOR]

We study of the properties of a new class of circumgalactic medium absorbers identified in the Ly α forest: 'Strong, Blended Lyman-α' (or SBLA) absorption systems. We study SBLAs at 2.4 < z < 3.1 in SDSS-IV/eBOSS spectra by their strong extended Ly α absorption complexes covering 138 |$\, \, {\rm km}\, {\rm s}^{-1}$| with an integrated |$\log (N_{\rm H\, {\small I}}/\mathrm{cm}^{-2}) =16.04$| |$\substack{+0.05 \\ -0.06}$| and Doppler parameter b = 18.1 |$\substack{+0.7 \\ -0.4}$| |$\, \, {\rm km}\, {\rm s}^{-1}$|⁠. Clustering with the Ly α forest provides a large-scale structure bias of b = 2.34 ± 0.06 and halo mass estimate of |$M_h \approx 10^{12}\, h^{-1}\, {\rm M_{\odot }}$| for our SBLA sample. We measure the ensemble mean column densities of 22 metal features in the SBLA composite spectrum and find that no single-population multiphase model for them is viable. We therefore explore the underlying SBLA population by forward modelling the SBLA absorption distribution. Based on covariance measurements and favoured populations we find that ≈25 per cent of our SBLAs have stronger metals. Using silicon only we find that our strong metal SBLAs trace gas with a log (n H/cm−3) > −2.40 for T = 103.5 K and show gas clumping on <210 parsec scales. We fit multiphase models to this strong subpopulation and find a low ionization phase with n H = 1 cm−3, T = 103.5 K, and [X/H] = 0.8, an intermediate ionization phase with log (n H/cm−3) = −3.05, T = 103.5 K and [X/H] = −0.8, and a poorly constrained higher ionization phase. We find that the low ionization phase favours cold, dense super-solar metallicity gas with a clumping scale of just 0.009 parsecs. [ABSTRACT FROM AUTHOR]

Context. The redshifted 21 cm signal from the Epoch of Reionization (EoR) directly probes the ionization and thermal states of the intergalactic medium during that period. In particular, the distribution of the ionized regions around the radiating sources during EoR introduces scale-dependent features in the spherically averaged EoR 21 cm signal power spectrum. Aims. The goal is to study these scale-dependent features at different stages of reionization using numerical simulations and to build a source model-independent framework to probe the properties of the intergalactic medium using EoR 21 cm signal power spectrum measurements. Methods. Under the assumption of high spin temperature, we modeled the redshift evolution of the ratio of the EoR 21 cm brightness temperature power spectrum to the corresponding density power spectrum using an ansatz consisting of a set of redshift and scale-independent parameters. This set of eight parameters probes the redshift evolution of the average ionization fraction and the quantities related to the morphology of the ionized regions. Results. We tested this ansatz on different reionization scenarios generated using different simulation algorithms and found that it is able to recover the redshift evolution of the average neutral fraction within an absolute deviation ≲0.1. Conclusions. Our framework allows us to interpret 21 cm signal power spectra in terms of parameters related to the state of the IGM. This source model-independent framework is able to efficiently constrain reionization scenarios using multi-redshift power spectrum measurements with ongoing and future radio telescopes such as LOFAR, MWA, HERA, and SKA. This will add independent information regarding the EoR IGM properties. [ABSTRACT FROM AUTHOR]

Absorption line spectroscopy using background quasars can provide strong constraints on galactic outflows. In this paper we investigate possible scaling relations between outflow properties, namely outflow velocity Vout, mass ejection rate Ṁout, and mass loading factor η, and the host galaxy properties, such as star formation rate (SFR), SFR surface density, redshift, and stellar mass, using galactic outflows probed by background quasars from MEGAFLOW and other surveys. We find that Vout (η) is (anti-)correlated with SFR and SFR surface density. We extend the formalism of momentum-driven outflows from a previous study to show that it applies not only to "down-the-barrel" studies, but also to winds probed by background quasars, suggesting a possible universal wind formalism. Under this formalism, we find a clear distinction between strong and weak outflows where strong outflows seem to have tighter correlations with galaxy properties (SFR or galaxy stellar mass) than weak outflows. [ABSTRACT FROM AUTHOR]

The optical depth parameterization is typically used to study the 21-cm signals associated with the properties of the neutral hydrogen (H  i) gas and the ionization morphology during the Epoch of Reionization (EoR), without solving the radiative transfer equation. To assess the uncertainties resulting from this simplification, we conduct explicit radiative transfer calculations using the cosmological 21-cm line radiative transfer (C21LRT) code and examine the imprints of ionization structures on the 21-cm spectrum. We consider a globally averaged reionization history and implement fully ionized cavities (H  ii bubbles) of diameters d ranging from 0.01 to 10 Mpc at epochs within the emission and the absorption regimes of the 21-cm global signal. The single-ray C21LRT calculations show that the shape of the imprinted spectral features are primarily determined by d and the 21-cm line profile, which is parametrized by the turbulent velocity of the H  i gas. It reveals the spectral features tied to the transition from ionized to neutral regions that calculations based on the optical depth parametrization were unable to capture. We also present analytical approximations of the calculated spectral features of the H  ii bubbles. The multiple-ray calculations show that the apparent shape of an H  ii bubble (of d  = 5 Mpc at z  = 8), because of the finite speed of light, differs depending on whether the bubble's ionization front is stationary or expanding. Our study shows the necessity of properly accounting for the effects of line-continuum interaction, line broadening, and cosmological expansion to correctly predict the EoR 21-cm signals. [ABSTRACT FROM AUTHOR]

We present the first observational measurements of the Lyman-α (Ly α) forest flux autocorrelation functions in ten redshift bins from 5.1 ≤ z ≤ 6.0. We use a sample of 35 quasar sightlines at z > 5.7 from the extended XQR-30 data set; these data have signal-to-noise ratios of >20 per spectral pixel. We carefully account for systematic errors in continuum reconstruction, instrumentation, and contamination by damped Ly α systems. With these measurements, we introduce software tools to generate autocorrelation function measurements from any simulation. Our measurements of the smallest bin of the autocorrelation function increase with redshift when normalizing by the mean flux, 〈 F 〉. This increase may come from decreasing 〈 F 〉 or increasing mean free path of hydrogen-ionizing photons, λmfp. Recent work has shown that the autocorrelation function from simulations at z > 5 is sensitive to λmfp, a quantity that contains vital information on the ending of reionization. For an initial comparison, we show our autocorrelation measurements with simulation models for recently measured λmfp values and find good agreements. Further work in modelling and understanding the covariance matrices of the data is necessary to get robust measurements of λmfp from this data. [ABSTRACT FROM AUTHOR]

We present the serendipitous discovery of (1) a large double radio relic associated with the galaxy cluster PSZ2 G277.93 + 12.34 and (2) a new odd radio circle, ORC J1027–4422, both found in the same deep MeerKAT 1.3 GHz wide-band radio continuum image. The angular separation of the two arc-shaped cluster relics is ∼16 arcmin or ∼2.6 Mpc for a cluster redshift of z ≈ 0.158. The thin southern relic, which shows several ridges/shocks including one possibly moving inwards, has a linear extent of ∼1.64 Mpc. In contrast, the northern relic is about twice as wide, twice as bright, but only has a largest linear size of ∼0.66 Mpc. Complementary SRG/eROSITA X-ray images reveal extended emission from hot intracluster gas between the two relics and around the narrow-angle tail (NAT) radio galaxy PMN J1033–4335 (z ≈ 0.153) located just east of the northern relic. The radio morphologies of the NAT galaxy and the northern relic, which are also detected with the Australian Square Kilometer Array Pathfinder (ASKAP) at 888 MHz, suggest both are moving in the same outward direction. The discovery of ORC J1027–4422 in a different part of the same MeerKAT image makes it the fourth known single ORC. It has a diameter of ∼90 arcsec corresponding to 400 kpc at a tentative redshift of z ≈ 0.3 and remains undetected in X-ray emission. Supported by simulations, we discuss similarities between outward moving galaxy and cluster merger shocks as the formation mechanisms for ORCs and radio relics, respectively. [ABSTRACT FROM AUTHOR]

We present a theoretical framework for linking quasar properties, such as quasar age, to the surrounding Ly α emission intensity. In particular, we focus on a method for mapping the large-scale structure of Ly α emission intensity with galaxy spectra from wide-field spectroscopic surveys, e.g. the Subaru Prime Focus Spectrograph (PFS) or the Dark Energy Spectroscopic Instrument (DESI), and consider the quasar-induced Ly α emission from the intergalactic medium (IGM). To do this, we construct a theoretical model based on two physical processes: resonant scattering of quasar Ly α photons and fluorescence due to quasar ionizing photons, finding that the fluorescence contribution due to optically thick gas clouds is dominant. Taking into account the light-cone effect and assuming a typical quasar spectrum, we calculate the fluorescence contribution to the spectrum stacked within each bin of the separation angle from the quasar as a function of quasar age. Furthermore, we compute the quasar-Ly α emission cross-correlation and its SNR for the planned PFS survey. The predicted signal can only account for |$\sim 10~{{\ \rm per\ cent}}$| of the measurements indicated from the BOSS and eBOSS surveys in the outer region of |$\gtrsim 10\ \rm {cMpc}\ \rm {h}^{-1}$|⁠ , and the predicted SNR is not sufficient to detect the quasar-induced contribution. However, we found that our model, combined with the contribution of star-forming galaxies, is not in conflict with these measurements. We discuss other possible contributions to the Ly α emission excess around quasars, the efficiency of using spectroscopic fibres, and the redshift dependence of our model. [ABSTRACT FROM AUTHOR]

The 21-cm hyperfine line of neutral hydrogen is a useful tool to probe the conditions of the Universe during the Dark Ages, Cosmic Dawn, and the Epoch of Reionization. In most of the current calculations, the 21-cm line signals at given frequencies are computed, using an integrated line-of-sight line opacity, with the correction for cosmological expansion. These calculations have not fully captured the line and continuum interactions in the radiative transfer, in response to evolution of the radiation field and the variations of thermal and dynamic properties of the line-of-sight medium. We construct a covariant formulation for the radiative transfer of the 21-cm line and derive the cosmological 21-cm line radiative transfer (C21LRT) equation. The formulation properly accounts for local emission and absorption processes and the interaction between the line and continuum when the radiation propagates across the expanding Universe to the present observer. Our C21LRT calculations show that methods simply summing the line optical depth could lead to error of 5  per cent in the 21-cm signals for redshift z ∼ 12–35 and of |$\gt 10~{{\ \rm per\ cent}}$| for redshift z ≲ 8. Proper covariant radiative transfer is therefore necessary for producing correct theoretical templates for extracting information of the structural evolution of the Universe through the Epoch of Reionization from the 21-cm tomographic data. [ABSTRACT FROM AUTHOR]

Confronting measurements of the Lyman-α forest with cosmological hydrodynamical simulations has produced stringent constraints on models of particle dark matter and the thermal and ionization state of the intergalactic medium. We investigate the robustness of such models of the Lyman-α forest, focusing on the effect of particle initial conditions on the Lyman-α forest statistics in cosmological SPH simulations. We study multiple particle initialization algorithms in simulations that are designed to be identical in other respects. In agreement with the literature, we find that the correct linear theory evolution is obtained when a glass-like configuration is used for initial unperturbed gas particle positions alongside a regular grid configuration for dark matter particles and the use of non-identical initial density perturbations for gas and dark matter. However, we report that this introduces a large scale-dependent distortion in the 1D Lyman-α transmission power spectrum at small scales (k > 0.05 s km−1). The effect is close to 50  per cent at k ∼ 0.1 s km−1, and persists at higher resolution. This can severely bias inferences in parameters such as the dark matter particle mass. By considering multiple initial conditions codes and their variations, we also study the impact of a variety of other assumptions and algorithmic choices, such as adaptive softening, background radiation density, particle staggering, and perturbation theory accuracy, on the matter power spectrum, the Lyman-α flux power spectrum, and the Lyman-α flux PDF. This work reveals possible pathways towards more accurate theoretical models of the Lyman-α forest to match the quality of upcoming measurements. [ABSTRACT FROM AUTHOR]

Recent measurements of the ionizing photon mean free path (MFP) based on composite quasar spectra may point to reionization ending at z < 6. These measurements are challenging because they rely on assumptions about the proximity zones of the quasars. For example, some quasars might have been close to neutral patches where reionization was still ongoing ('neutral islands'), and it is unclear how they would affect the measurements. We address this question with mock MFP measurements from radiative transfer simulations. We find that, even in the presence of neutral islands, our mock MFP measurements agree to within |$30~{{\ \rm per\ cent}}$| with the true spatially averaged MFP in our simulations, which includes opacity from both the ionized medium and the islands. The inferred MFP is sensitive at the |$\lt ~50~{{\ \rm per\ cent}}$| level to assumptions about quasar environments and lifetimes for realistic models. We demonstrate that future analyses with improved data may require explicitly modelling the effects of neutral islands on the composite spectra, and we outline a method for doing this. Lastly, we quantify the effects of neutral islands on Lyman-series transmission, which has been modelled with optically thin simulations in previous MFP analyses. Neutral islands can suppress transmission at λrest < 912 Å significantly, up to a factor of 2 for z qso = 6 in a plausible reionization scenario, owing to absorption by many closely spaced lines as quasar light redshifts into resonance. However, the suppression is almost entirely degenerate with the spectrum normalization and thus does not significantly bias the inferred MFP. [ABSTRACT FROM AUTHOR]

Ly α damping wings in the spectra of bright objects at high redshift are a useful probe of the ionization state of the intergalactic medium during the reionization epoch. It has recently been noted that, despite the inhomogeneous nature of reionization, these damping wings have a characteristic shape which is a strong function of the volume-weighted average neutral hydrogen fraction of the intergalactic medium. We present here a closer examination of this finding using a simulation of patchy reionization from the Sherwood–Relics simulation suite. We show that the characteristic shape and scatter of the damping wings are determined by the average neutral hydrogen density along the line of sight, weighted by its contribution to the optical depth producing the damping wing. We find that there is a redshift dependence in the characteristic shape due to the expansion of the Universe. Finally, we show that it is possible to differentiate between the shapes of damping wings in galaxies and young (or faint) quasars at different points in the reionization history at large velocity offsets from the point where the transmission first reaches zero. [ABSTRACT FROM AUTHOR]

The abundance, temperature, and clustering of metals in the intergalactic medium are important parameters for understanding their cosmic evolution and quantifying their impact on cosmological analysis with the Ly α forest. The properties of these systems are typically measured from individual quasar spectra redward of the quasar's Ly α emission line, yet that approach may provide biased results due to selection effects. We present an alternative approach to measure these properties in an unbiased manner with the two-point statistics commonly employed to quantify large-scale structure. Our model treats the observed flux of a large sample of quasar spectra as a continuous field and describes the one-dimensional, two-point statistics of this field with three parameters per ion: the abundance (column density distribution), temperature (Doppler parameter), and clustering (cloud-cloud correlation function). We demonstrate this approach on multiple ions (e.g. C IV, Si IV, and Mg II) with early data from the Dark Energy Spectroscopic Instrument (DESI) and high-resolution spectra from the literature. Our initial results show some evidence that the C IV abundance is higher than previous measurements and evidence for abundance evolution over time. The first full year of DESI observations will have over an order of magnitude more quasar spectra than this study. In a future paper, we will use those data to measure the growth of clustering and its impact on the Ly α forest, as well as test other DESI analysis infrastructure such as the pipeline noise estimates and the resolution matrix.

We explore the effect of ionizing ultraviolet background (UVB) on the redshift space clustering of low- z (z ≤ 0.5) O  vi absorbers using the Sherwood simulations incorporating 'WIND' (i.e. outflows driven by stellar feedback)-only and 'WIND + AGN (active galactic nucleus)' feedback. These simulations show positive clustering signals up to a scale of 3 Mpc. We find that the effect of feedback is restricted to small scales (i.e. ≤2 Mpc or ≈200 km s−1 at z ∼ 0.3) and 'WIND'-only simulations produce a stronger clustering signal compared to simulations incorporating 'WIND + AGN' feedback. How the clustering signal is affected by the assumed UVB depends on the feedback processes assumed. For the simulations considered here, the effect of the UVB is confined to even smaller scales (i.e. <1 Mpc or ≈100 km s−1 at z ∼ 0.3). These scales are also affected by exclusion caused by line blending. Therefore, our study suggests that clustering at intermediate scales (i.e. 1–2 Mpc for simulations considered here) together with the observed column density distribution can be used to constrain the effect of feedback in simulations. [ABSTRACT FROM AUTHOR]

The attenuation of Lyα photons by neutral hydrogen in the intergalactic medium (IGM) at z ≳ 5 continues to be a powerful probe for studying the epoch of reionization. Given a framework to estimate the intrinsic (true) Lyα emission of high- z sources, one can infer the ionization state of the IGM during reionization. In this work, we use the enlarged XQR-30 sample of 42 high-resolution and high signal-to-noise quasar spectra between |$5.8\lesssim \, z\lesssim \, 6.6$| obtained with VLT/X-shooter to place constraints on the IGM neutral fraction. This is achieved using our existing Bayesian QSO reconstruction framework which accounts for uncertainties such as the: (i) posterior distribution of predicted intrinsic Lyα emission profiles (obtained via covariance matrix reconstruction of the Lyα and N  v emission lines from unattenuated high-ionization emission line profiles; C  iv , Si  iv   + O  iv ], and C  iii ]) and (ii) distribution of ionized regions within the IGM using synthetic damping wing profiles drawn from a 1.63 Gpc3 reionization simulation. Following careful quality control, we used 23 of the 42 available QSOs to obtain constraints/limits on the IGM neutral fraction during the tail-end of reionization. Our median and 68th percentile constraints on the IGM neutral fraction are: |$0.20\substack{+0.14 -0.12}$| and |$0.29\substack{+0.14 -0.13}$| at z  = 6.15 and 6.35. Further, we also report 68th percentile upper limits of |$\bar{x}_{\mathrm{H\, {\small I}}{}} \lt 0.21$|⁠ , 0.20, 0.21, and 0.18 at z  = 5.8, 5.95, 6.05, and 6.55. These results imply reionization is still ongoing at |$5.8\lesssim \, z\lesssim \, 6.55$|⁠ , consistent with previous results from XQR-30 (dark fraction and Lyα forest) along with other observational probes considered in the literature. [ABSTRACT FROM AUTHOR]

The early stages of the epoch of reionization, probed by the 21 cm line, are sensitive to the detailed properties and formation histories of the first galaxies. We use 21cmfast and a simple, self-consistent galaxy model to examine the redshift evolution of the large-scale cross-power spectrum between the 21 cm field and line-emitting galaxies. A key transition in redshift occurs when the 21 cm field shifts from being positively correlated with the galaxy distribution to being negatively correlated. Importantly, this transition redshift is insensitive to the properties of the galaxy tracers but depends sensitively on the thermal and ionization histories traced through the 21 cm field. Specifically, we show that the transition occurs when both ionization fluctuations dominate over 21 cm spin temperature fluctuations and when the average spin temperature exceeds the temperature of the cosmic microwave background. We illustrate this with three different 21 cm models which have largely the same neutral fraction evolution but different heating histories. We find that the transition redshift has a scale dependence, and that this can help disentangle the relative importance of heating and ionization fluctuations. The best prospects for constraining the transition redshift occur in scenarios with late X-ray heating, where the transition occurs at redshifts as low as z ∼ 6–8. In our models, this requires high-redshift galaxy surveys with sensitivities of |$\sim 10^{-18}~\rm erg\,s^{-1}\,cm^{-2}$| for optical lines and |$\sim 10^{-19}~\rm erg\,s^{-1}\,cm^{-2}$| for far-infrared lines. Future measurements of the transition redshift can help discriminate between 21 cm models and will benefit from reduced systematics. [ABSTRACT FROM AUTHOR]

Observed accretion rates onto the Milky Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is a rotating inflow in the volume-filling hot phase (⁠|$\sim 10^6\, {\rm K}$|⁠) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution valid in the slow-rotation limit, we show that a hot inflow spins up as it approaches the galaxy, while remaining hot, subsonic, and quasi-spherical. Within the radius of angular momentum support (⁠|$\sim 15\, {\rm kpc}$| for the Milky Way) the hot flow flattens into a disc geometry and then cools from |$\sim 10^6$| to |$\sim 10^4\, {\rm K}$| at the disc–halo interface. Cooling affects all hot gas, rather than just a subset of individual gas clouds, implying that accretion via hot inflows does not rely on local thermal instability in contrast with 'precipitation' models for galaxy accretion. Prior to cooling and accretion the inflow completes ≈ t cool/ t ff radians of rotation, where t cool/ t ff is the cooling time to free-fall time ratio in hot gas immediately outside the galaxy. The ratio t cool/ t ff may thus govern the development of turbulence and enhancement of magnetic fields in gas accreting onto low-redshift spirals. We show that if rotating hot inflows are common in Milky-Way-size disc galaxies, as predicted, then signatures of the expected hot gas rotation profile should be observable with X-ray telescopes and fast radio burst surveys. [ABSTRACT FROM AUTHOR]

The Epoch of Reionization (EoR) began when galaxies grew in abundance and luminosity, so their escaping Lyman continuum (LyC) radiation started ionizing the surrounding neutral intergalactic medium (IGM). Despite significant recent progress, the nature and role of cosmic reionizers are still unclear: in order to define them, it would be necessary to directly measure their LyC escape fraction (fesc). However, this is impossible during the EoR due to the opacity of the IGM. Consequently, many efforts at low and intermediate redshift have been made to determine measurable indirect indicators in high-redshift galaxies so that their fesc can be predicted. This work presents the analysis of the indirect indicators of 62 spectroscopically confirmed star-forming galaxies at 6 ≤ z ≤ 9 from the Cosmic Evolution Early Release Science (CEERS) survey, combined with 12 sources with public data from other JWST-ERS campaigns. From the NIRCam and NIRSpec observations, we measured their physical and spectroscopic properties. We discovered that on average 6 < z < 9 star-forming galaxies are compact in the rest-frame UV (re ∼ 0.4 kpc), are blue sources (UV-β slope ∼ − 2.17), and have a predicted fesc of about 0.13. A comparison of our results to models and predictions as well as an estimation of the ionizing budget suggests that low-mass galaxies with UV magnitudes fainter than M1500 = −18 that we currently do not characterize with JWST observations probably played a key role in the process of reionization. [ABSTRACT FROM AUTHOR]