Search NASASearch

SEARCH · Search NASA

Results for “INTERGALACTIC MEDIUM”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Measurements of the thermal and ionization state of the intergalactic medium during the cosmic afternoon

We perform the first measurement of the thermal and ionization state of the intergalactic medium (IGM) across 0.9 < z < 1.5 using 301 Ly α absorption lines fitted from 12 archival Hubble Space Telescope Space Telescope Imaging Spectrograph quasar spectra. We employ the machine-learning-based inference method that uses joint Doppler parameter–column density (⁠b-N HI ⁠) distributions obtained from Ly α forest decomposition. Our results show that the Γ HI photoionization rates, ⁠, agree with recent ultraviolet background synthesis models, with log(Γ HI /s -1 ) = $-11.79^{+0.18}_{-0.15}$, $-11.98^{+0.09}_{-0.09}$⁠, and $-12.32^{+0.10}_{-0.12}$⁠, at z = 1.4, 1.2, and 1, respectively. We obtain the IGM temperature at the mean density, T 0 ⁠, and the adiabatic index, γ⁠, as [log(T 0 /K), γ] = $[4.13^{+0.12}_{-0.10}, 1.34^{+0.10}_{-0.15}]$, $[3.79^{+0.11}_{-0.11}, 1.70^{+0.09}_{-0.09}]$, and $[4.12^{+0.15}_{-0.25}, 1.34^{+0.21}_{-0.26}]$ at z = 1.4⁠, 1.2, and 1. Our measurements of T 0 at z = 1.4 and 1.2 are consistent with the trend predicted from previous z < 3 temperature measurements and theoretical expectations, where the IGM cools down after $He\tiny{II}$ reionization in the absence of any non-standard heating. However, our T 0 measurement at z = 1 unexpectedly high IGM temperature. Given the relatively large uncertainty in these measurements, where σ T$_0$ ~ 5000 K, mostly emanating from the limited size of our data set, we cannot conclude whether the IGM cools down as expected. Lastly, we generate mock data sets to test the constraining power of future measurement with larger data sets. The results demonstrate that, with redshift path-length Δz ~ 2 for each redshift bin, three times the current data set, we can constrain the T 0 of IGM within 1500 K, which would be sufficient to constrain the IGM thermal history at z < 1.5 conclusively.

79 ASTRONOMY AND ASTROPHYSICS

Cosmic Reionization on Computers: Statistical Properties of the Distributions of Mean Opacities

Quasar absorption lines provide a unique window to the relationship between galaxies and the intergalactic medium during the Epoch of Reionization. In particular, high redshift quasars enable measurements of the neutral hydrogen content of the universe. However, the limited sample size of observed quasar spectra, particularly at the highest redshifts, hampers our ability to fully characterize the intergalactic medium during this epoch from observations alone. In this work, we characterize the distributions of mean opacities of the intergalactic medium in simulations from the Cosmic Reionization on Computers (CROC) project. We find that the distribution of mean opacities along sightlines follows a non-trivial distribution that cannot be easily approximated by a known distribution. When comparing the cumulative distribution function of mean opacities measurements in subsamples of sample sizes similar to observational measurements from the literature, we find consistency between CROC and observations at redshifts $z\lesssim 5.7$. However, at higher redshifts ($z\gtrsim5.7$), the cumulative distribution function of mean opacities from CROC is notably narrower than those from observed quasar sightlines implying that observations probe a systematically more opaque intergalactic medium at higher redshifts than the intergalactic medium in CROC boxes at these same redshifts. This is consistent with previous analyses that indicate that the universe is reionized too early in CROC simulations.

79 ASTRONOMY AND ASTROPHYSICS

Tracing the Evolution of the Cool Gas in CGM and IGM Environments through Mg II Absorption from Redshift z = 0.75 to z = 1.65 Using DESI-Y1 Data

We present a measurement of the mean absorption of cool gas traced by Mg II (λλ2796, 2803) around emission line galaxies (ELGs), spanning spatial scales from 20 kpc to 10 Mpc. The measurement is based on crossmatching the positions ELGs at z = 0.75–1.65 and the metal absorption in the spectra of background quasars with data provided by the Year 1 sample of the Dark Energy Spectroscopic Instrument. The ELGs are divided into two redshift intervals: 0.75 < z < 1.0 and 1.0 < z < 1.65. We find that the composite quasar spectra constructed by stacking the ELG-QSO pairs evolve with redshift, with z > 1 having a systematically higher signal of Mg II absorption. Within 1 Mpc, the covering fraction of the cool gas at z > 1 is higher than that of z < 1. The enhancement becomes less apparent especially if the projected distance r p > 1 Mpc. ELGs with higher stellar mass and star formation rate (SFR) yield higher clustering of Mg ii absorbers at z < 1. For z > 1, the covering fractions with different SFRs show little difference. The higher Mg II absorption at higher redshift supports the observations of higher star formation at cosmic noon. Converting the Mg II absorbers to unsaturated Si II , our estimate indicates that the metal abundance of Si II ranges from 0.7 to 1.2 × 10 −6 from z = 0.9 to 1.3. The growth of low-ionization metal abundance strongly suggests a metal-enriched circumgalactic medium and an increased presence of cool gas in the intergalactic medium toward higher redshifts.

circumgalactic medium

The impact of varying inhomogeneous reionization histories on metrics of Ly α opacity

The epoch of hydrogen reionization is complete by z = 5⁠, but its progression at higher redshifts is uncertain. Measurements of Ly α forest opacity show large scatter at z < 6⁠, suggestive of spatial fluctuations in neutral fraction, temperature, or ionizing background, either individually or in combination. However, there are degeneracies in the impact of such fluctuations, necessitating careful modelling. We develop a framework for modelling the reionization history and associated temperature fluctuations, with the intention of incorporating ionizing background fluctuations at a later time. We generate several reionization histories using seminumerical code AMBER, and implement them in the Nyx cosmological hydrodynamics code to examine the impact on the evolution of gas within the simulation and the associated metrics of the Ly α forest opacity. We find that the pressure smoothing scale within the intergalactic medium is strongly correlated with the adiabatic index of the temperature–density relation. We find that while models with 20 000 K photoheating at reionization are better able to reproduce the shape of the observed z = 5 1D flux power spectrum than colder ones, they fail to match the highest wavenumbers. The simulated autocorrelation function and optical depth distributions are systematically low and narrow, respectively, compared to the observed values, but are in better agreement when the reionization history is longer in duration, more symmetric in its distribution of reionization redshifts, or if there are remaining neutral regions at z < 6.

79 ASTRONOMY AND ASTROPHYSICS

Constraints on dark photon dark matter from Lyman- α forest simulations and an ultrahigh signal-to-noise quasar spectrum

The ultralight dark photon is a well-motivated, hypothetical dark matter candidate. In a dilute plasma, they can resonantly convert into photons, and heat up the intergalactic medium between galaxies. In this work, we explore the dark photon dark matter parameter space by comparing synthetic Lyman- α forest data from cosmological hydrodynamical simulations to observational data from VLT/UVES of the quasar HE0940-1050 ( z em = 3.09 ). We use a novel flux normalization technique that targets underdense gas, reshaping the flux probability distribution. Not only do we place robust constraints on the kinetic mixing parameter of dark photon dark matter, but notably our findings suggest that this model can still reconcile simulated and observed Doppler parameter distributions of z ∼ 0 Lyman- α lines, as seen by HST/COS. This work opens new pathways for the use of the Lyman- α forest to explore new physics, and can be extended to other scenarios such as primordial black hole evaporation, dark matter decay, and annihilation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

The Cosmic Evolution of C IV Absorbers at 1.4 < z < 4.5: Insights from 100,000 Systems in DESI Quasars

We present the largest catalog to date of triply ionized carbon (C IV ) absorbers detected in quasar spectra from the Dark Energy Spectroscopic Instrument. Using an automated matched-kernel convolution method with adaptive signal-to-noise thresholds, we identify 101,487 C IV systems in the redshift range 1.4 < z < 4.5 from 300,637 quasar spectra. Completeness is estimated via Monte Carlo simulations, and the catalog is 50% complete at EW C IV ≥ 0.4 Å. The differential equivalent width frequency distribution declines exponentially and shows weak redshift evolution. The absorber incidence per unit comoving path increases by a factor of 2–5 from z ≈ 4.5 to z ≈ 1.4, with stronger redshift evolution for strong systems. Using column densities derived from the apparent optical depth method, we constrain the cosmic mass density of C IV , Ω C IV , which increases by a factor of ∼3.8 from (0.82 ± 0.05) × 10 −8 at z ≈ 4.5 to (3.16 ± 0.2) × 10 −8 at z ≈ 1.4. From Ω C IV , we estimate a lower limit on intergalactic medium metallicity ${\mathrm{log}}({Z}_{{\rm{IGM}}}/{Z}_{\odot })\gtrsim -3.25$ at z ∼ 2.3, with a smooth decline at higher redshifts. These trends trace the cosmic star formation history and He II photoheating rate, suggesting a link between C IV enrichment, star formation, and UV background over ∼3 Gyr. The catalog also provides a critical resource for future studies connecting circumgalactic metals to galaxy evolution, especially near cosmic noon.

79 ASTRONOMY AND ASTROPHYSICS

Suppression of pair beam instabilities in a laboratory analogue of blazar pair cascades

The generation of dense electron–positron pair beams in the laboratory can enable direct tests of theoretical models of γ-ray bursts and active galactic nuclei. We have successfully achieved this using ultrarelativistic protons accelerated by the Super Proton Synchrotron at (CERN). In the first application of this experimental platform, the stability of the pair beam is studied as it propagates through a meter-length plasma, analogous to TeV γ-ray-induced pair cascades in the intergalactic medium. It has been argued that pair beam instabilities disrupt the cascade, thus accounting for the observed lack of reprocessed GeV emission from TeV blazars. If true, this would remove the need for a moderate strength intergalactic magnetic field to explain the observations. We find that the pair beam instability is suppressed if the beam is not perfectly collimated or monochromatic, hence the lower limit to the intergalactic magnetic field inferred from γ-ray observations of blazars is robust.

blazar jets

Modeling of the high column density systems in the Lyman-Alpha forest

The Lyman-α forests observed in the spectra of high-redshift quasars can be used as a tracer of the cosmological matter density to study baryon acoustic oscillations (BAO) and the Alcock-Paczynski effect. Extraction of cosmological information from these studies requires modeling of the forest correlations. While the models depend most importantly on the bias parameters of the intergalactic medium (IGM), they also depend on the numbers and characteristics of high-column-density systems (HCDs) ranging from Lyman-limit systems with column densities log N HI /1cm -2 > 17 to damped Lyman-α systems (DLAs) with log N HI /1cm -2 > 20.2. These HCDs introduce broad damped absorption characteristic of a Voigt profile. Consequently they imprint a component on the power spectrum whose modes in the radial direction are suppressed, leading to a scale-dependent bias. Using mock data sets of known HCD content, we test a model that describes this effect in terms of the distribution of column densities of HCDs, the Fourier transforms of their Voigt profiles and the bias of the halos containing the HCDs. Our results show that this physically well-motivated model describes the effects of HCDs with an accuracy comparable to that of the ad-hoc models used in published forest analyses. We also discuss the problems of applying the model to real data, where the HCD content and their bias is uncertain.

Lyman alpha forest

ForestFlow: predicting the Lyman-α forest clustering from linear to nonlinear scales

On large scales, the Lyman-α forest provides insights into the expansion history of the Universe, while on small scales, it imposes strict constraints on the growth history, the nature of dark matter, and the sum of neutrino masses. This work introduces ForestFlow, a novel framework that bridges the gap between large- and small-scale analyses, which have traditionally relied on distinct modeling approaches. Using conditional normalizing flows, ForestFlow predicts the two Lyman-α linear biases (b δ and b η ) and six parameters describing small-scale deviations of the three-dimensional flux power spectrum (P 3D ) from linear theory as a function of cosmology and intergalactic medium physics. These are then combined with a Boltzmann solver to make consistent predictions, from arbitrarily large scales down to the nonlinear regime, for P 3D and any other statistics derived from it. Trained on a suite of 30 fixed-and-paired cosmological hydrodynamical simulations spanning redshifts from z = 2 to 4.5, ForestFlow achieves 3 and 1.5% precision in describing P 3D and the one-dimensional flux power spectrum (P 1D ) from linear scales to k = 5 Mpc −1 and k ∥ = 4 Mpc −1 , respectively. Thanks to its conditional parameterization, ForestFlow shows similar performance for ionization histories and two ΛCDM model extensions – massive neutrinos and curvature – even though none of these are included in the training set. This framework will enable full-scale cosmological analyses of Lyman-α forest measurements from the DESI survey.

79 ASTRONOMY AND ASTROPHYSICS

Characterization of contaminants in the Lyman-alpha forest auto-correlation with DESI

Baryon Acoustic Oscillations can be measured with sub-percent precision above redshift two with the Lyman-α (Lyα) forest auto-correlation and its cross-correlation with quasar positions. This is one of the key goals of the Dark Energy Spectroscopic Instrument (DESI) which started its main survey in May 2021. We present in this paper a study of the contaminants to the Lyα forest which are mainly caused by correlated signals introduced by the spectroscopic data processing pipeline as well as astrophysical contaminants due to foreground absorption in the intergalactic medium. Notably, an excess signal caused by the sky background subtraction noise is present in the Lyα auto-correlation in the first line-of-sight separation bin. We use synthetic data to isolate this contribution, we also characterize the effect of spectro-photometric calibration noise, and propose a simple model to account for both effects in the analysis of the Lyα forest. We then measure the auto-correlation of the quasar flux transmission fraction of low redshift quasars, where there is no Lyα forest absorption but only its contaminants. We demonstrate that we can interpret the data with a two-component model: data processing noise and triply ionized Silicon and Carbon auto-correlations. This result can be used to improve the modeling of the Lyα auto-correlation function measured with DESI.

79 ASTRONOMY AND ASTROPHYSICS

DESI DR1 Lyα 1D power spectrum: the optimal estimator measurement

The one-dimensional power spectrum P 1D of Lyα forest offers rich insights into cosmological and astrophysical parameters, including constraints on the sum of neutrino masses, warm dark matter models, and the thermal state of the intergalactic medium. We present the measurement of P 1D using the optimal quadratic maximum likelihood estimator applied to over 300,000 Lyα quasars from Data Release 1 (DR1) of the Dark Energy Spectroscopic Instrument (DESI) survey. This sample represents the largest to date for P 1D measurements and is larger than the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) by a factor of 1.7. We conduct a meticulous investigation of instrumental and analysis systematics and quantify their impact on P 1D . This includes the development of a cross-exposure estimator that eliminates the need to model the pipeline noise and has strong potential for future P 1D measurements. We also present new insights into metal contamination through the 1D correlation function. Using a fitting function we measure the evolution of the Lyα forest bias with high precision: b F (z) = (-0.218 ± 0.002) × ((1 + z)/4) 2.96±0.06 . In a companion validation paper, we substantially extend our previous suite of CCD image simulations to quantify the pipeline's exquisite performance accurately. In another companion paper, we present DR1 P 1D measurements using the Fast Fourier Transform (FFT) approach to power spectrum estimation. These two measurements produce a forest bias parameter that differs by 2.2 sigma. However, our model is simplistic, so this disagreement will be investigated in future work.

Lyman alpha forest

Generation of lognormal synthetic Lyman-α forest spectra for P 1 D analysis

The one-dimensional flux power spectrum (P 1 D ) of the Lyman-α forest probes small-scale structure in the intergalactic medium (IGM) and is therefore sensitive to a variety of cosmological and astrophysical parameters. These include the amplitude and shape of the matter power spectrum, the thermal history of the IGM, the sum of neutrino masses, and potential small-scale fluctuations due to the nature of dark matter. However, P 1 D is also highly sensitive to observational and instrumental systematics, making accurate synthetic spectra essential for validating analyses and quantifying these effects, especially in high-volume surveys like the Dark Energy Spectroscopic Instrument (DESI). We present an efficient lognormal mock framework for generating one-dimensional Lyman-α forest spectra tailored for P 1 D analysis. Our method captures the redshift evolution of the mean transmitted flux and the scale-dependent shape and amplitude of the one-dimensional flux power spectrum by tuning Gaussian field correlations and transformation parameters. Across the DESI Early Data Release (EDR) redshift range (2.0 ≤ z ≤ 3.8), and a wide range of scales (10 -4 s km -1 ≤ k ≤ 1.0 s km -1 ), our mocks recover the mean flux evolution with redshift to sub-percent accuracy, and the P 1 D at the percent level. Additionally, we discuss potential extensions of this framework, such as the incorporation of astrophysical contaminants, continuum uncertainties, and instrumental effects. Such improvements would expand its utility in ongoing and upcoming surveys and enable a broader range of validation efforts and systematics studies for P 1 D inference and precision cosmology.

Lyman alpha forest

Effects of primordial black holes on IGM history

Currently the asteroid mass window (mass ∼ 10$^{17}$- 10$^{21}$ grams) remains unconstrained for Primordial Black Holes (PBHs) to make up all of the dark matter content of the universe. Given these PBHs have very small masses, their Hawking temperature can be up to hundreds of keV. This study investigates the potential impacts of PBH Hawking radiation on the intergalactic medium from z ∼ 800-25, namely studying the ionization history, kinetic gas temperature, and ultimately the 21 cm signature. We find that for masses on the low edge of the asteroid mass window, there are up two orders of magnitude increases in the ionization fraction and kinetic gas temperature by redshift 25, and the 21 cm spin temperature can differ from non-PBH cosmology by factors of a few. This analysis results in maximum differential brightness temperatures of +17 mK for our lightest PBH masses of 2.12 × 10$^{16}$g. We also show maximal 53 mK discrepancies in differential brightness temperatures between our PBH and non-PBH cosmologies for our lightest PBH mass, while our heaviest PBH mass of 1.65 × 10$^{17}$g shows only 0.5 mK variations. We find the Hawking-radiated electrons and positrons are instrumental in driving these IGM modifications. This study shows the necessity for a rigorous treatment of Hawking radiation in PBH cosmological observables from the dark ages through cosmic dawn.

dark matter theory

Cosmological analysis of the DESI DR1 Lyα 1D power spectrum

We present the cosmological analysis of the one-dimensional Lyman-α flux power spectrum from the first data release of the Dark Energy Spectroscopic Instrument (DESI). We capture the dependence of the signal on cosmology and intergalactic medium physics using an emulator trained on a cosmological suite of hydrodynamical simulations, and we correct its predictions for the impact of astrophysical contaminants and systematics, many of these not considered in previous analyses. We employ this framework to constrain the amplitude and logarithmic slope of the linear matter power spectrum at k ★ = 0.009 km -1 s and redshift z = 3, obtaining Δ 2 ★ = 0.379 ± 0.032 and n ★ = -2.309 ± 0.019 https://github.com/igmhub/cobaya_lya_p1d. The robustness of these constraints is validated through the analysis of mocks and a large number of alternative data analysis variations, with cosmological parameters kept blinded throughout the validation process. We then combine our results with constraints from DESI BAO and temperature, polarization, and lensing measurements from Planck, ACT, and SPT-3G to set constraints on ΛCDM extensions. While our measurements do not significantly tighten the limits on the sum of neutrino masses from the combination of these probes, they sharpen the constraints on the effective number of relativistic species, N eff = 3.02 ± 0.10, the running of the spectral index, α s = 0.0014 ± 0.0041, and the running of the running, β s = -0.0006 ± 0.0048, by a factor of 1.18, 1.27, and 1.90, respectively. We conclude by outlining the improvements needed to fully reach the level of confidence implied by these uncertainties.

79 ASTRONOMY AND ASTROPHYSICS

Reconstructing Quasar Spectra and Measuring the Lyα Forest with SpenderQ

Quasar spectra carry the imprint of foreground intergalactic medium (IGM) through absorption features. In particular, absorption caused by neutral hydrogen gas, the "Lyα forest," is a key spectroscopic tracer for cosmological analyses used to measure cosmic expansion and test physics beyond the standard model. Despite their importance, current methods for measuring Lyα absorption cannot directly derive the intrinsic quasar continuum and make strong assumptions on its shape, thus distorting the measured Lyα clustering. We present SpenderQ , a ML-based approach for directly reconstructing the intrinsic quasar spectra and measuring the Lyα forest from observations. SpenderQ uses the Spender spectrum autoencoder to learn a compact and redshift-invariant latent encoding of quasar spectra, combined with an iterative procedure to identify and mask absorption regions. To demonstrate its performance, we apply SpenderQ to 400,000 synthetic quasar spectra created to validate the Dark Energy Spectroscopic Instrument Year 1 Lyα cosmological analyses. SpenderQ accurately reconstructs the true intrinsic quasar spectra, including the broad Lyβ, Lyα, SiIV, CIV, and CIII emission lines. Redward of Lyα, SpenderQ provides percent-level reconstructions of the true quasar spectra. Blueward of Lyα, SpenderQ reconstructs the true spectra to < 5%. SpenderQ reproduces the shapes of individual quasar spectra more robustly than the current state-of-the-art. We, thus, expect it will significantly reduce biases in Lyα clustering measurements and enable studies of quasars and their physical properties. SpenderQ also provides informative latent variable encodings that can be used to, e.g., classify quasars with Broad Absorption Lines. Overall, SpenderQ provides a new data-driven approach for unbiased Lyα forest measurements in cosmological, quasar, and IGM studies.

Hahn, ChangHoon [Arizona U., Astron. Dept. - Stewa

Multiwavelength investigations of PKS 2300–18: S-shaped radio quasar with precessing jets and double-peaked broad emission-line spectrum

ABSTRACT S-shaped radio galaxy jets are prime sources for investigating the dynamic interplay between the central active galactic nucleus (AGN), the jets, and the ambient intergalactic medium. These sources are excellent candidates for studying jet precession, as their S-shaped inversion symmetry strongly indicates underlying precession. We present a multiwavelength analysis of the giant inversion-symmetric S-shaped radio galaxy PKS 2300$-$18, which spans 0.76 Mpc. The host is a quasar at a redshift of 0.128, displaying disturbed optical morphology due to an ongoing merger with a companion galaxy. We conducted a broad-band radio spectral study using multifrequency data ranging from 183 MHz to 6 GHz, incorporating dedicated observations with the upgraded Giant Metrewave Radio Telescope (uGMRT) and Karl G. Jansky Very Large Array (JVLA) alongside archival radio data. A particle injection model was fitted to the spectra of different regions of the source to perform ageing analysis, which was supplemented with a kinematic jet precession model. The ageing analysis revealed a maximum plasma age of $\sim$ 40 Myr, while the jet precession model indicated a precession period of $\sim$ 12 Myr. ROentgen SATellite (ROSAT) data revealed an X-ray halo of Mpc size, and from Chandra the AGN X-ray spectrum was modelled using thermal and power-law components. The optical spectrum displaying double-peaked broad emission lines was modelled, indicating complex broad-line region kinematics at the centre with the possibility of a binary SMBH. We present the results of our multiwavelength analysis of the source, spanning scales from a few light-days to a few Mpc, and discuss its potential evolutionary path.

Misra, Arpita

Forecasting constraints on the high-z IGM thermal state from the Lyman-α forest flux autocorrelation function

ABSTRACT The autocorrelation function of the Lyman-$\alpha$ (Ly $\alpha$) forest flux from high-z quasars probes the small-scale structure of the intergalactic medium (IGM). The thermal state of the IGM, determined by the physics of reionization, sets the small-scale power observed in the Ly $\alpha$ forest. To explore the sensitivity of the autocorrelation function to the IGM’s thermal state, we compute the autocorrelation function from a cosmological hydrodynamical simulation with an instantaneous reionization model and 135 post-processed thermal states. Using mock data sets of 20 quasars, we forecast constraints on $T_0$ and $\gamma$, which characterize the post-processed IGM thermal state, at $5.4 \le z \le 6$. While this model simplifies the IGM’s thermal state, it serves as a key first step in assessing future observational prospects. We also perform an inference test on mocks and re-weight out posterior distributions to guarantee that they exhibit statistically correct behaviour. At $z = 5.4$, we find that an idealized data set constrains $T_0$ to 59 per cent and $\gamma$ to 16 per cent at the 1$\sigma$ equivalent confidence level. To explore more realistic, non-instantaneous reionization scenarios, we analyse four models combining temperature and ultraviolet background (UVB) fluctuations at $z = 5.8$. We find that mock data generated from a model with both temperature and UVB fluctuations can rule out a model with only temperature fluctuations at the $> 1\sigma$ level 73.9 per cent of the time.

Wolfson, Molly

Modeling the Cosmological Lyman-𝛼 Forest at the Field Level

The distribution of absorption lines in the spectra of distant quasars, called the Lyman-𝛼 (Ly-𝛼) forest, is a unique probe of cosmology and the intergalactic medium at high redshifts and small scales. The statistical power of ongoing redshift surveys demands precise theoretical tools to model the Ly-𝛼 forest. We address this challenge by developing an analytic, perturbative forward model to predict the Ly-𝛼 forest at the field level for a given set of cosmological initial conditions. Our model shows a remarkable performance when compared with the Sherwood hydrodynamic simulations: it reproduces the Ly-𝛼 forest flux power spectrum, its cross-correlation with dark matter halos, and the one-point probability distribution function of both fields at the percent level down to scales of a few Mpc. Our work provides crucial tools that bridge analytic modeling on large scales with simulations on small scales, enabling field-level inference from Ly-𝛼 forest data and simulation-based priors for cosmological analyses. Furthermore, this is especially timely for realizing the full scientific potential of the Ly-𝛼 forest measurements by the dark energy spectroscopic instrument.

Cosmological parameters