Revealing the z~2.5 Cosmic Web With 3D Lyman-Alpha Forest Tomography: A Deformation Tensor Approach

Studies of cosmological objects should take into account their positions within the cosmic web of large-scale structure. Unfortunately, the cosmic web has only been extensively mapped at low-redshifts ($z<1$), using galaxy redshifts as tracers of the underlying density field. At $z>1$, the required galaxy densities are inaccessible for the foreseeable future, but 3D reconstructions of Lyman-$\alpha$ forest absorption in closely-separated background QSOs and star-forming galaxies already offer a detailed window into $z\sim2-3$ large-scale structure. We quantify the utility of such maps for studying the cosmic web by using realistic $z=2.5$ Ly$\alpha$ forest simulations matched to observational properties of upcoming surveys. A deformation tensor-based analysis is used to classify voids, sheets, filaments and nodes in the flux, which is compared to those determined from the underlying dark matter field. We find an extremely good correspondence, with $70\%$ of the volume in the flux maps correctly classified relative to the dark matter web, and $99\%$ classified to within 1 eigenvalue. This compares favorably to the performance of galaxy-based classifiers with even the highest galaxy densities at low-redshift. We find that narrow survey geometries can degrade the cosmic web recovery unless the survey is $\gtrsim 60\,h^{-1}\,\mathrm{Mpc}$ or $\gtrsim 1\,\mathrm{deg}$ on the sky. We also examine halo abundances as a function of the cosmic web, and find a clear dependence as a function of flux overdensity, but little explicit dependence on the cosmic web. These methods will provide a new window on cosmological environments of galaxies at this very special time in galaxy formation, "high noon", and on overall properties of cosmological structures at this epoch.Comment: 10 pages, 8 figures. Accepted by Ap

Joint Bayesian Estimation of Quasar Continua and the Lyman-Alpha Forest Flux Probability Distribution Function

We present a new Bayesian algorithm making use of Markov Chain Monte Carlo sampling that allows us to simultaneously estimate the unknown continuum level of each quasar in an ensemble of high-resolution spectra, as well as their common probability distribution function (PDF) for the transmitted Ly$\alpha$ forest flux. This fully automated PDF regulated continuum fitting method models the unknown quasar continuum with a linear Principal Component Analysis (PCA) basis, with the PCA coefficients treated as nuisance parameters. The method allows one to estimate parameters governing the thermal state of the intergalactic medium (IGM), such as the slope of the temperature-density relation $\gamma-1$, while marginalizing out continuum uncertainties in a fully Bayesian way. Using realistic mock quasar spectra created from a simplified semi-numerical model of the IGM, we show that this method recovers the underlying quasar continua to a precision of $\simeq7\%$ and $\simeq10\%$ at $z=3$ and $z=5$, respectively. Given the number of principal component spectra, this is comparable to the underlying accuracy of the PCA model itself. Most importantly, we show that we can achieve a nearly unbiased estimate of the slope $\gamma-1$ of the IGM temperature-density relation with a precision of $\pm8.6\%$ at $z=3$, $\pm6.1\%$ at $z=5$, for an ensemble of ten mock high-resolution quasar spectra. Applying this method to real quasar spectra and comparing to a more realistic IGM model from hydrodynamical simulations would enable precise measurements of the thermal and cosmological parameters governing the IGM, albeit with somewhat larger uncertainties given the increased flexibility of the model.Comment: 21 pages (+ Appendix), accepted at Ap

Protocluster Discovery in Tomographic Lyα\alpha Forest Flux Maps

We present a new method of finding protoclusters using tomographic maps of Ly$\alpha$ Forest flux. We review our method of creating tomographic flux maps and discuss our new high performance implementation, which makes large reconstructions computationally feasible. Using a large N-body simulation, we illustrate how protoclusters create large-scale flux decrements, roughly 10 $h^{-1}$Mpc across, and how we can use this signal to find them in flux maps. We test the performance of our protocluster finding method by running it on the ideal, noiseless map and tomographic reconstructions from mock surveys, and comparing to the halo catalog. Using the noiseless map, we find protocluster candidates with about 90% purity, and recover about 75% of the protoclusters that form massive clusters ($> 3 \times 10^{14} \, h^{-1} M_{\odot}$). We construct mock surveys similar to the ongoing COSMOS Lyman-Alpha Mapping And Tomography Observations (CLAMATO) survey. While the existing data has an average sightline separation of 2.3 $h^{-1}$Mpc, we test separations of 2 - 6 $h^{-1}$Mpc to see what can be tolerated for our application. Using reconstructed maps from small separation mock surveys, the protocluster candidate purity and completeness are very close what was found in the noiseless case. As the sightline separation increases, the purity and completeness decrease, although they remain much higher than we initially expected. We extended our test cases to mock surveys with an average separation of 15 $h^{-1}$Mpc, meant to reproduce high source density areas of the BOSS survey. We find that even with such a large sightline separation, the method can still be used to find some of the largest protoclusters.Comment: 18 pages, 12 figure

Observational Requirements for Lyman-alpha Forest Tomographic Mapping of Large-Scale Structure at z ~ 2

The z > 2 Lyman-alpha (Lya) forest traces the underlying dark-matter distribution on large scales and, given sufficient sightlines, can be used to create 3D maps of large-scale structure. We examine the observational requirements to construct such maps and estimate the signal-to-noise as a function of exposure time and sightline density. Sightline densities at z = 2.25 are n_los = [360, 1200,3300] deg^{-2} at limiting magnitudes of g =[24.0, 24.5,25.0], resulting in transverse sightline separations of d_perp = [3.6, 1.9, 1.2] h^{-1} Mpc, which roughly sets the reconstruction scale. We simulate these reconstructions using mock spectra with realistic noise properties, and find that spectra with S/N = 4 per angstrom can be used to generate maps that clearly trace the underlying dark-matter at overdensities of rho/ ~ 1. For the VLT/VIMOS spectrograph, exposure times t_exp = [4, 6, 10] hrs are sufficient for maps with spatial resolution epsilon_3d = [5.0, 3.2, 2.3] h^{-1} Mpc. Assuming ~ 250 h^{-1} Mpc is probed along the line-of-sight, 1 deg^2 of survey area would cover a comoving volume of ~ 10^6 h^{-3} Mpc^3 at =2.3, enabling efficient mapping of large volumes with 8-10m telescopes. These maps could be used to study galaxy environments, detect proto-clusters, and study the topology of large-scale structure at high-z.Comment: 18 pages, 10 figures. Accepted by Ap

TARDIS Paper II: Synergistic Density Reconstruction from Lyman-alpha Forest and Spectroscopic Galaxy Surveys with Applications to Protoclusters and the Cosmic Web

In this work we expand upon the Tomographic Absorption Reconstruction and Density Inference Scheme (TARDIS) in order to include multiple tracers while reconstructing matter density fields at Cosmic Noon (z ~ 2-3). In particular, we jointly reconstruct the underlying density field from simulated Lyman-Alpha forest observations at $z\sim 2.5$ and an overlapping galaxy survey. We find that these data are synergistic, with the Lyman Alpha forest providing reconstruction of low density regions and galaxy surveys tracing the density peaks. We find a more accurate power spectra reconstruction going to higher scales when fitting these two data-sets simultaneously than if using either one individually. When applied to cosmic web analysis, we find performing the joint analysis is equivalent to a Lyman Alpha survey with significantly increased sight-line spacing. Since we reconstruct the velocity field and matter field jointly, we demonstrate the ability to evolve the mock observed volume further to z=0, allowing us to create a rigorous definition of "proto-cluster" as regions which will evolve into clusters. We apply our reconstructions to study protocluster structure and evolution, finding for realistic survey parameters we can provide accurate mass estimates of the z \approx 2 structures and their z = 0 fate.Comment: 16 pages, 12 figure

A New Precision Measurement of the Small-Scale Line-of-Sight Power Spectrum of the Lyα{\alpha} Forest

We present a new measurement of the Ly${\alpha}$ forest power spectrum at $1.8 < z < 3.4$ using 74 Keck/HIRES and VLT/UVES high-resolution, high-S/N quasar spectra. We developed a custom pipeline to measure the power spectrum and its uncertainty, which fully accounts for finite resolution and noise, and corrects for the bias induced by masking missing data, DLAs, and metal absorption lines. Our measurement results in unprecedented precision on the small-scale modes $k > 0.02\,\mathrm{s\,km^{-1}}$, unaccessible to previous SDSS/BOSS analyses. It is well known that these high-$k$ modes are highly sensitive to the thermal state of the intergalactic medium, however contamination by narrow metal lines is a significant concern. We quantify the effect of metals on the small-scale power, and find a modest effect on modes with $k < 0.1\,\mathrm{s\,km^{-1}}$. As a result, by masking metals and restricting to $k < 0.1\,\mathrm{s\,km^{-1}}$ their impact is completely mitigated. We present an end-to-end Bayesian forward modeling framework whereby mock spectra with the same noise, resolution, and masking as our data are generated from Ly${\alpha}$ forest simulations. These mocks are used to build a custom emulator, enabling us to interpolate between a sparse grid of models and perform MCMC fits. Our results agree well with BOSS on scales $k < 0.02\,\mathrm{s\,km^{-1}}$ where the measurements overlap. The combination of BOSS' percent level low-$k$ precision with our $5-15\%$ high-$k$ measurements, results in a powerful new dataset for precisely constraining the thermal history of the intergalactic medium, cosmological parameters, and the nature of dark matter. The power spectra and their covariance matrices are provided as electronic tables.Comment: 24 pages, 12 figures, accepted for publication in ApJ, machine readable tables will be made available after publication in the journa