The correlation of high-redshift galaxies with the thermal Sunyaev-Zel’dovich effect traces reionization

We explore a potential new probe of reionization: the cross-correlation of high-redshift galaxies with maps of the thermal Sunyaev-Zel’dovich (tSZ) effect. We consider two types of high redshift galaxies: Lyman break galaxies (LBGs) and Lyman-α emitters (LAEs). LBGs and LAEs will be detected in large numbers at high redshift (z ≈ 4 – 7) by ongoing and future surveys. We consider a future LBG sample from The Rubin Observatory Legacy Survey of Space and Time (LSST), and a selection of LAEs modelled after the Subaru SILVERRUSH program, but covering a much larger sky fraction. The tSZ effect is sensitive to a line-of-sight integral of the ionized gas pressure, and is measured across large patches of sky using multi-frequency CMB surveys. We consider forecast tSZ maps from CMB Stage 4 and more futuristic observations. Using a suite of hydrodynamical simulations, we show that LBGs and LAEs are correlated with the tSZ signal from reionization. The cross-spectra between LBGs/LAEs with tSZ maps contain information about the reionization history of the Universe, such as the distribution of bubble sizes, and could be used to directly measure the timing of reionization. The amplitude of the signal is small, however, and its detectability is hindered by low-redshift contributions to tSZ maps and by instrumental noise. If the low-redshift contribution to the observed tSZ signal is suppressed by masking of massive halos, a combination of overlapping futuristic CMB and galaxy surveys could probe this signal

The Lyman α forest power spectrum from the XQ-100 legacy survey

We present the Lyman α flux power spectrum measurements of the XQ-100 sample of quasar spectra obtained in the context of the European Southern Observatory Large Programme ‘Quasars and their absorption lines: a legacy survey of the high redshift universe with VLT/XSHOOTER’. Using 100 quasar spectra with medium resolution and signal-to-noise ratio, we measure the power spectrum over a range of redshifts z = 3–4.2 and over a range of scales k = 0.003–0.06 km−1 s. The results agree well with the measurements of the one-dimensional power spectrum found in the literature. The data analysis used in this paper is based on the Fourier transform and has been tested on synthetic data. Systematic and statistical uncertainties of our measurements are estimated, with a total error (statistical and systematic) comparable to the one of the BOSS data in the overlapping range of scales, and smaller by more than 50 per cent for higher redshift bins (z > 3.6) and small scales (k > 0.01 km−1 s). The XQ-100 data set has the unique feature of having signal-to-noise ratios and resolution intermediate between the two data sets that are typically used to perform cosmological studies, i.e. BOSS and high-resolution spectra (e.g. UVES/VLT or HIRES). More importantly, the measured flux power spectra span the high-redshift regime that is usually more constraining for structure formation models

Deep Learning of Dark Energy Spectroscopic Instrument Mock Spectra to Find Damped Ly alpha Systems

We have updated and applied a convolutional neural network (CNN) machine-learning model to discover and characterize damped Lyα systems (DLAs) based on Dark Energy Spectroscopic Instrument (DESI) mock spectra. We have optimized the training process and constructed a CNN model that yields a DLA classification accuracy above 99% for spectra that have signal-to-noise ratios (S/N) above 5 per pixel. The classification accuracy is the rate of correct classifications. This accuracy remains above 97% for lower S/N ≈1 spectra. This CNN model provides estimations for redshift and H i column density with standard deviations of 0.002 and 0.17 dex for spectra with S/N above 3 pixel-1. Also, this DLA finder is able to identify overlapping DLAs and sub-DLAs. Further, the impact of different DLA catalogs on the measurement of baryon acoustic oscillations (BAO) is investigated. The cosmological fitting parameter result for BAO has less than 0.61% difference compared to analysis of the mock results with perfect knowledge of DLAs. This difference is lower than the statistical error for the first year estimated from the mock spectra: above 1.7%. We also compared the performances of the CNN and Gaussian Process (GP) models. Our improved CNN model has moderately 14% higher purity and 7% higher completeness than an older version of the GP code, for S/N > 3. Both codes provide good DLA redshift estimates, but the GP produces a better column density estimate by 24% less standard deviation. A credible DLA catalog for the DESI main survey can be provided by combining these two algorithms

The wide-field, multiplexed, spectroscopic facility WEAVE: Survey design, overview, and simulated implementation

WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator.Comment: 41 pages, 27 figures, accepted for publication by MNRA

The wide-field, multiplexed, spectroscopic facility WEAVE : survey design, overview, and simulated implementation

Funding for the WEAVE facility has been provided by UKRI STFC, the University of Oxford, NOVA, NWO, Instituto de Astrofísica de Canarias (IAC), the Isaac Newton Group partners (STFC, NWO, and Spain, led by the IAC), INAF, CNRS-INSU, the Observatoire de Paris, Région Île-de-France, CONCYT through INAOE, Konkoly Observatory (CSFK), Max-Planck-Institut für Astronomie (MPIA Heidelberg), Lund University, the Leibniz Institute for Astrophysics Potsdam (AIP), the Swedish Research Council, the European Commission, and the University of Pennsylvania.WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366-959 nm at R ∼ 5000, or two shorter ranges at R ∼ 20,000. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for ∼ 3 million stars and detailed abundances for ∼ 1.5 million brighter field and open-cluster stars; (ii) survey ∼ 0.4 million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey  ∼ 400 neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in z 1 million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at z > 2. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator.PostprintPeer reviewe

The wide-field, multiplexed, spectroscopic facility WEAVE: Survey design, overview, and simulated implementation

WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366−959\,nm at R∼5000, or two shorter ranges at R∼20000. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for ∼3 million stars and detailed abundances for ∼1.5 million brighter field and open-cluster stars; (ii) survey ∼0.4 million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey ∼400 neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in z1 million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at z>2. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator

The correlation of high-redshift galaxies with the thermal Sunyaev-Zel’dovich effect traces reionization

We explore a potential new probe of reionization: the cross-correlation of high-redshift galaxies with maps of the thermal Sunyaev-Zel’dovich (tSZ) effect. We consider two types of high redshift galaxies: Lyman break galaxies (LBGs) and Lyman-α emitters (LAEs). LBGs and LAEs will be detected in large numbers at high redshift (z ≈ 4 – 7) by ongoing and future surveys. We consider a future LBG sample from The Rubin Observatory Legacy Survey of Space and Time (LSST), and a selection of LAEs modelled after the Subaru SILVERRUSH program, but covering a much larger sky fraction. The tSZ effect is sensitive to a line-of-sight integral of the ionized gas pressure, and is measured across large patches of sky using multi-frequency CMB surveys. We consider forecast tSZ maps from CMB Stage 4 and more futuristic observations. Using a suite of hydrodynamical simulations, we show that LBGs and LAEs are correlated with the tSZ signal from reionization. The cross-spectra between LBGs/LAEs with tSZ maps contain information about the reionization history of the Universe, such as the distribution of bubble sizes, and could be used to directly measure the timing of reionization. The amplitude of the signal is small, however, and its detectability is hindered by low-redshift contributions to tSZ maps and by instrumental noise. If the low-redshift contribution to the observed tSZ signal is suppressed by masking of massive halos, a combination of overlapping futuristic CMB and galaxy surveys could probe this signal

Small-scale clumping of dark matter and the mean free path of ionizing photons at z=6z=6

Recently, the mean free path of ionizing photons in the $z = 6$ intergalactic medium (IGM) was measured to be very short, presenting a challenge to existing reionization models. At face value, the measurement can be interpreted as evidence that the IGM clumps on scales $M\lesssim 10^8$ M$_\odot$, a key but largely untested prediction of the cold dark matter (CDM) paradigm. Motivated by this possibility, we study the role that the underlying dark matter cosmology plays in setting the $z > 5$ mean free path. We use two classes of models to contrast against the standard CDM prediction: (1) thermal relic warm dark matter (WDM), representing models with suppressed small-scale power; (2) an ultralight axion exhibiting a white noise-like power enhancement. Differences in the mean free path between the WDM and CDM models are subdued by pressure smoothing and the possible contribution of neutral islands to the IGM opacity. For example, comparing late reionization scenarios with a fixed volume-weighted mean neutral fraction of $20\%$ at $z=6$, the mean free path is $19~(45)~\%$ longer in a WDM model with $m_x = 3~(1)$ keV. The enhanced power in the axion-like model produces better agreement with the short mean free path measured at $z = 6$. However, drawing robust conclusions about cosmology is hampered by large uncertainties in the reionization process, extragalactic ionizing background, and thermal history of the Universe. This work highlights some key open questions about the IGM opacity during reionization.Comment: 23 pages, 11 figure

Forecasts for WEAVE-QSO: 3D clustering of critical points with Lyman-alpha tomography

The upcoming WEAVE-QSO survey will target a high density of quasars over a large area, enabling the reconstruction of the 3D density field through Lyman-훼 tomography over unprecedented volumes smoothed on intermediate cosmological scales (≈ 16 Mpc/h). We produce mocks of the Lyman-훼 forest using LyMAS, and reconstruct the 3D density field between sightlines through Wiener filtering in a configuration compatible with the future WEAVE-QSO observations. The fidelity of the reconstruction is assessed by measuring one- and two-point statistics from the distribution of critical points in the cosmic web. In addition, initial Lagrangian statistics are predicted from first principles, and measurements of the connectivity of the cosmic web are performed. The reconstruction captures well the expected features in the auto- and cross-correlations of the critical points. This remains true after a realistic noise is added to the synthetic spectra, even though sparsity of sightlines introduces systematics, especially in the cross-correlations of points with mixed signature. Specifically, the most striking clustering features involving filaments and walls could be measured with up to 4 sigma of significance with a WEAVE-QSO-like survey. Moreover, the connectivity of each peak identified in the reconstructed field is globally consistent with its counterpart in the original field, indicating that the reconstruction preserves the geometry of the density field not only statistically, but also locally. Hence the critical points relative positions within the tomographic reconstruction could be used as standard rulers for dark energy by WEAVE-QSO and similar surveys