SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2022

An intergalactic medium temperature from a giant radio galaxy

ABSTRACT The warm–hot intergalactic medium (warm–hot IGM, or WHIM) pervades the filaments of the Cosmic Web and harbours half of the Universe’s baryons. The WHIM’s thermodynamic properties are notoriously hard to measure. Here we estimate a galaxy group–WHIM boundary temperature using a new method. In particular, we use a radio image of the giant radio galaxy (giant RG, or GRG) created by NGC 6185, a massive nearby spiral. We analyse this extraordinary object with a Bayesian 3D lobe model and deduce an equipartition pressure Peq = 6 × 10−16 Pa – among the lowest found in RGs yet. Using an X-ray-based statistical conversion for Fanaroff–Riley II RGs, we find a true lobe pressure $P = 1.5\substack{+1.7\\ -0.4}\times 10^{-15}\ \mathrm{Pa}$. Cosmic Web reconstructions, group catalogues, and MHD simulations furthermore imply an Mpc–scale IGM density $1 + \delta _\mathrm{IGM} = 40\substack{+30\\ -10}$. The buoyantly rising lobes are crushed by the IGM at their inner side, where an approximate balance between IGM and lobe pressure occurs: PIGM ≈ P. The ideal gas law then suggests an IGM temperature $T_\mathrm{IGM} = 11\substack{+12\\-5} \times 10^6\ \mathrm{K}$, or $k_\mathrm{B}T_\mathrm{IGM} = 0.9\substack{+1.0\\-0.4}\ \mathrm{keV}$, at the virial radius – consistent with X-ray-derived temperatures of similarly massive groups. Interestingly, the method is not performing at its limit: in principle, estimates TIGM ∼ 4 × 106 K are already possible – rivalling the lowest X-ray measurements available. The technique’s future scope extends from galaxy group outskirts to the WHIM. In conclusion, we demonstrate that observations of GRGs in Cosmic Web filaments are finally sensitive enough to probe the thermodynamics of galaxy groups and beyond.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Emergence of the temperature–density relation in the low-density intergalactic medium

We examine the evolution of the phase diagram of the low-density intergalactic medium during the Epoch of Reionization in simulation boxes with varying reionization histories from the Cosmic Reionization on Computers project. The probability density function (PDF) of gas temperature at fixed density exhibits two clear modes: a warm and a cold temperature mode, corresponding to the gas inside and outside of ionized bubbles. We find that the transition between the two modes is ‘universal’ in the sense that its timing is accurately parametrized by the value of the volume-weighted neutral fraction for any reionization history. This ‘universality’ is more complex than just a reflection of the fact that ionized gas is warm and neutral gas is cold: it holds for the transition at a fixed value of gas density, and gas at different densities transitions from the cold to the warm mode at different values of the neutral fraction, reflecting a non-trivial relationship between the ionization history and the evolving gas density PDF. Furthermore, the ‘emergence’ of the tight temperature–density relation in the warm mode is also approximately ‘universally’ controlled by the volume-weighted neutral fraction for any reionization history. In particular, the ‘emergence’ of the temperature–density relation (as quantified by the rapid decrease in its width) occurs when the neutral fraction is 10 –4 ≲ X HI ≲ 10 –3 for any reionization history. Our results indicate that the neutral fraction is a primary quantity controlling the various properties of the temperature–density relation, regardless of reionization history.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

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↗

DOE OSTI2020

Effective Opacity of the Intergalactic Medium from Galaxy Spectra Analysis

We measure the effective opacity (τ{sub eff}) of the intergalactic medium from the composite spectra of 281 Lyman-break galaxies in the redshift range 2 ≲ z ≲ 3. Our spectra are taken from the COSMOS Lyα Mapping And Tomographic Observations survey derived from the Low Resolution Imaging Spectrometer on the W.M. Keck I telescope. We generate composite spectra in two redshift intervals and fit them with spectral energy distribution (SED) models composed of simple stellar populations. Extrapolating these SED models into the Lyα forest, we measure the effective Lyα opacity (τ{sub eff}) in the 2.02 ≤ z ≤ 2.44 range. At z = 2.22, we estimate τ{sub eff} =0.159±0.001 from a power-law fit to the data. These measurements are consistent with estimates from quasar analyses at z < 2.5 indicating that the systematic errors associated with normalizing quasar continua are not substantial. We provide a Gaussian processes model of our results and previous τ{sub eff} measurements that describes the steep redshift evolution in τ{sub eff} from z = 1.5–4.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Inferring the Thermal History of the Intergalactic Medium from the Properties of the Hydrogen and Helium Lyα Forest

The filamentary network of intergalactic medium (IGM) gas that gives origin to the Lyα forest in the spectra of distant quasars encodes information on the physics of structure formation and the early thermodynamics of diffuse baryonic material. Here we use a massive suite of more than 400 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit–accelerated code Cholla to study the IGM at high spatial resolution maintained over the entire computational volume. The simulations capture a wide range of possible IGM thermal histories by varying the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei. A statistical comparison of synthetic spectra with the observed 1D flux power spectra of hydrogen at redshifts 2.2 ≤ z ≤ 5.0 and with the helium Lyα opacity at redshifts 2.4 < z < 2.9 tightly constrains the photoionization and photoheating history of the IGM. By leveraging the constraining power of the available Lyα forest data to break model degeneracies, we find that the IGM experienced two main reheating events over 1.2 Gyr of cosmic time. For our best-fit model, hydrogen reionization completes by z R ≈ 6.0 with a first IGM temperature peak of T 0 ≃ 1.3 × 10 4 K and is followed by the reionization of He ii that completes by z R ≈ 3.0 and yields a second temperature peak of T 0 ≃ 1.4 × 10 4 K. We discuss how our results can be used to obtain information on the timing and the sources of hydrogen and helium reionization.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Lyman-α transmission properties of the intergalactic medium in the CoDaII simulation

ABSTRACT The decline in abundance of Lyman-α (Lyα) emitting galaxies at z ≳ 6 is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of ∼94 comoving Mpc side length, to compute the Lyα transmission properties of the intergalactic medium (IGM) at z ∼ 5.8 to 7. Our results mainly confirm previous studies, i.e. we find a declining Lyα transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Lyα peaks at high redshift, we also analyse the IGM transmission on the blue side, which shows a rapid decline at z ≳ 6 of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as $n_{\rm HI}\sim 10^{-9}\, \mathrm{cm}^{-3}$ (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that $\sim 1{{\ \rm per\ cent}}$ of sightlines towards M1600AB ∼ −21 galaxies at z ∼ 7 are transparent enough to allow a transmission of a blue Lyα peak. We discuss our results in the context of the interpretation of observations.

Gronke, Max (ORCID:000000032491060X)↗

DOE OSTI2025

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↗

DOE OSTI2022

Revision of Faraday rotation measure constraints on the primordial magnetic field using the IllustrisTNG simulation

ABSTRACT Previously derived Faraday rotation constraints on the volume-filling intergalactic magnetic field (IGMF) have used analytical models that made a range of simplifying assumptions about magnetic field evolution in the intergalactic medium and did not consider the effect of baryonic feedback on large-scale structures. In this work, we revise existing Faraday rotation constraints on the IGMF using a numerical model of the intergalactic medium from the IllustrisTNG cosmological simulation that includes a sophisticated model of the baryonic feedback. We use the IllustrisTNG model to calculate the rotation measure and compare the resulting mean and median of the absolute value of the rotation measure with data from the NRAO VLA Sky Survey (NVSS). The numerical model of the intergalactic medium includes a full magnetohydrodynamic model of the compressed primordial magnetic field as well as a model of the regions where the magnetic field is not primordial, but is rather produced by the process of baryonic feedback. Separating these two types of regions, we are able to assess the influence of the primordial magnetic field on the Faraday rotation signal. We find that by correcting for regions of compressed primordial field and accounting for the fact that part of the intergalactic medium is occupied by magnetic fields spread by baryonic feedback processes rather than by the primordial field relaxes the Faraday rotation bound by a factor of ≃3. This results in B0 < 1.8 × 10−9 G for large correlation length IGMFs.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Probing large-scale UV background inhomogeneity associated with quasars using metal absorption

ABSTRACT We study large-scale UV background inhomogeneity in three-dimensions associated with the observed quasar populations at high redshift. We do this by measuring metal absorption through quasar absorption spectrum stacking as a function distance to closest quasar in SDSS-IV/eBOSS on 10 s of comoving megaparsec scales. We study both intergalactic medium absorbers and mixed circumgalactic medium absorbers and probe absorption in O vi, C iv, and Si iv, and Si iii. Overall, stronger high ionization species absorption is seen closer to quasars at 2.4 < z < 3.1. O vi absorption shows a particularly strong change, with effects in C iv evident in some cases, and more marginal effects in Si iii and Si iv. We further study 2.05 < z < 2.4 (with weak signs of increasing homogeneity with time) and explore the study of metal absorption as a function of integrated SDSS-r band flux quasar flux (yielding consistent but less significant results). While the metal absorption does show sensitivity to large-scale 3D quasar proximity, the current incomplete quasar samples limit detailed interpretation. This work does, however, demonstrate that UV background inhomogeneities exist on scales of several 10 s of comoving megaparsecs associated with quasars and that they can be measured with precision by examining metal absorption in the intergalactic medium.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

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↗

DOE OSTI2022

LAGER Ly α Luminosity Function at z ~ 7: Implications for Reionization

We present a new measurement of the Lyα luminosity function (LF) at redshift z = 6.9, finding moderate evolution from z = 5.7 that is consistent with a fully or largely ionized z ~ 7 intergalactic medium. Our result is based on four fields of the LAGER (Lyman Alpha Galaxies in the Epoch of Reionization) project. Our survey volume of 6.1 × 10 6 Mpc 3 is double that of the next largest z ~ 7 survey. We combine two new LAGER fields (WIDE12 and GAMA15A) with two previously reported LAGER fields (COSMOS and CDFS). In the new fields, we identify N = 95 new z = 6.9 Lyα emitter (LAEs) candidates, characterize our survey's completeness and reliability, and compute Lyα LFs. The best-fit Schechter LF parameters for all four LAGER fields are in good general agreement. Two fields (COSMOS and WIDE12) show evidence for a bright-end excess above the Schechter function fit. We find that the Lyα luminosity density declines at the same rate as the UV continuum LF from z = 5.7 to 6.9. This is consistent with an intergalactic medium that was fully ionized as early as redshift z ~ 7 or with a volume-averaged neutral hydrogen fraction of x H I < 0.33 at 1σ.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Outside the Lyman-break box: detecting Lyman continuum emitters at 3.5 < z < 5.1 with CLAUDS

ABSTRACT Identifying non-contaminated sample of high-redshift galaxies with escaping Lyman continuum (LyC) flux is important for understanding the sources and evolution of cosmic reionization. We present CLAUDS (CFHT Large Area u-band deep survey) u-band photometry of the COSMOS field to probe LyC radiation from spectroscopically confirmed galaxies at $z$ ≥ 3.5 and outside the standard Lyman-break galaxy colour-selection expectations. Complementary to the CLAUDS data, we use Subaru multifilter photometry, Hubble Space Telescope (HST) multifilter imaging, and the spectroscopic surveys D10K, VUDS, and 3D-HST. We present a sample of Lyman continuum galaxy (LCG) candidates in the redshift range 3.5 ≲ $z$ ≲ 5.1. Here, we introduce 5 LCG candidates, where two are flagged quality 1 and three quality 2. The estimated $f_{\rm esc}^{\rm abs}$ for quality 1 candidates are in the range $\sim 5 - 73{{\ \rm per\ cent}}$ and $\sim 30 - 93{{\ \rm per\ cent}}$. These estimates are based on our derived parameters from individual galaxies as inputs to a range of BPASS models as well as mean intergalactic medium (IGM) and maximal intergalactic and circumgalactic media (IGM+CGM) transmission. We conclude that our search for LCGs is most likely biased to lines of sight with low H i densities or free from Lyman limit systems. Our two best LCG candidates have EW (Lyα) ≤ 50 Å and we find no correlation or anticorrelation between EW (Lyα), $f_{\rm esc}^{\rm abs}$, and Robs, the ratio of ionizing to non-ionizing observed flux in the measured passbands. Stacking candidates without solid LyC detections (S/N < 3) results in an estimated $f_{\rm esc}^{\rm abs}$ from galaxies not greater than $1{{\ \rm per\ cent}}$.

Meštrić, U.↗

DOE OSTI2023

The Bimodal Absorption System Imaging Campaign (BASIC). I. A Dual Population of Low-metallicity Absorbers at z < 1

The bimodal absorption system imaging campaign (BASIC) aims to characterize the galaxy environments of a sample of 36 H I-selected partial Lyman limit systems (pLLSs) and Lyman limit systems (LLSs) in 23 QSO fields at z ≲ 1. These pLLSs/LLSs provide a unique sample of absorbers with unbiased and well-constrained metallicities, allowing us to explore the origins of metal-rich and low-metallicity circumgalactic medium (CGM) at z < 1. Here we present Keck/KCWI and Very Large Telescope/MUSE observations of 11 of these QSO fields (19 pLLSs) that we combine with Hubble Space Telescope/Advanced Camera for Surveys imaging to identify and characterize the absorber-associated galaxies at 0.16 ≲ z ≲ 0.84. We find 23 unique absorber-associated galaxies, with an average of one associated galaxy per absorber. For seven absorbers, all with <10% solar metallicities, we find no associated galaxies with logM * ≳ 9.0 within ρ/R vir and |Δv|/v esc ≤ 1.5 with respect to the absorber. We do not find any strong correlations between the metallicities or H I column densities of the gas and most of the galaxy properties, except for the stellar mass of the galaxies: the low-metallicity ([X/H] ≤ –1.4) systems have a probability of $0.39_{–0.15}^{+0.16}$ for having a host galaxy with logM * ≥ 9.0 within ρ/R vir ≤ 1.5, while the higher metallicity absorbers have a probability of $0.78_{–0.13}^{+0.10}$. This implies metal-enriched pLLSs/LLSs at z < 1 are typically associated with the CGM of galaxies with logM * > 9.0, whereas low-metallicity pLLSs/LLSs are found in more diverse locations, with one population arising in the CGM of galaxies and another more broadly distributed in overdense regions of the universe. Using absorbers not associated with galaxies, we estimate the unweighted geometric mean metallicity of the intergalactic medium to be [X/H] ≲ –2.1 at z < 1, which is lower than previously estimated.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

A Hydro-particle-mesh Code for Efficient and Rapid Simulations of the Intracluster Medium

We introduce the cosmological HYPER code based on an innovative hydro-particle-mesh (HPM) algorithm for efficient and rapid simulations of gas and dark matter. For the HPM algorithm, we update the approach of Gnedin & Hui to expand the scope of its application from the lower-density intergalactic medium (IGM) to the higher-density intracluster medium (ICM). While the original algorithm tracks only one effective particle species, the updated version separately tracks the gas and dark matter particles, as they do not exactly trace each other on small scales. For the approximate hydrodynamics solver, the pressure term in the gas equations of motion is calculated using robust physical models. In particular, we use a dark matter halo model, ICM pressure profile, and IGM temperature–density relation, all of which can be systematically varied for parameter-space studies. We show that the HYPER simulation results are in good agreement with the halo model expectations for the density, temperature, and pressure radial profiles. Simulated galaxy cluster scaling relations for Sunyaev–Zel’dovich (SZ) and X-ray observables are also in good agreement with mean predictions, with scatter comparable to that found in hydrodynamic simulations. HYPER also produces lightcone catalogs of dark matter halos and full-sky tomographic maps of the lensing convergence, SZ effect, and X-ray emission. These simulation products are useful for testing data analysis pipelines, generating training data for machine learning, understanding selection and systematic effects, and for interpreting astrophysical and cosmological constraints.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Intergalactic magnetic field studies by means of γ -ray emission from GRB 190114C

The presence of delayed GeV emission after a strong transient, such as a gamma-ray burst (GRB), in the very-high energy (VHE) E > 100 GeV band can be the signature of a nonzero magnetic field in the intergalactic medium. Here we used a synchrotron self-Compton multiwavelength model to infer an analytical description of the intrinsic VHE spectrum (corrected for absorption by the extragalactic background light) of GRB 190114C to predict the light curves and spectral energy distributions of the delayed emission with Monte Carlo simulations for different intergalactic magnetic field (IGMF) configurations (strengths B = 8 × 10 -21 G, 10 -20 G, 3 × 10 -20 G and correlation length $\lambda$ > 1 Mpc), and compared them with the Fermi Large Area Telescope (Fermi LAT) limits computed for several exposure times. We found that the Fermi LAT is not sensitive enough to constrain any IGMF strengths using GRB 190114C.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

MUSE Analysis of Gas around Galaxies (MAGG) – IV. The gaseous environment of z ∼ 3–4 Ly α emitting galaxies

ABSTRACT We study the link between galaxies and H i-selected absorption systems at z ∼ 3–4 in the MUSE Analysis of Gas around Galaxies (MAGG) survey, an ESO large programme consisting of integral field spectroscopic observations of 28 quasar fields hosting 61 strong absorbers with $\rm N_{\rm H\,{\small I}}\gtrsim 10^{16.5}~\rm cm^{-2}$. We identify 127 Ly α emitting galaxies (LAEs) around the absorbers, corresponding to a detection rate of 82 ± 16 per cent. The luminosity function of these LAEs is ≈5 times higher in normalization than the field population and we detect a significant clustering of galaxies with respect to the gas, confirming that high column density absorbers and LAEs trace each other. Between 30 and 40 per cent of the absorbers are associated with multiple LAEs, which lie preferentially along filaments. Galaxies in groups also exhibit a three times higher covering factor of optically thick gas compared to isolated systems. No significant correlations are identified between the emission properties of LAEs and the absorption properties of optically thick gas clouds, except for a weak preference of brighter and multiple galaxies to reside near broad absorbers. Based on the measured impact parameters and the covering factor, we conclude that the near totality of optically thick gas in the Universe can be found in the outer circumgalactic medium (CGM) of LAEs or in the intergalactic medium (IGM) in proximity to these galaxies. Thus, LAEs act as tracers of larger scale structures within which both galaxies and optically thick clouds are embedded. The patchy and inhomogeneous nature of the CGM and IGM explains the lack of correlations between absorption and emission properties. This implies that very large samples are needed to unveil the trends that encode the properties of the baryon cycle.

Astronomy & Astrophysics↗

DOE OSTI2024

Cosmic Reionization on Computers: Statistics, Physical Properties, and Environments of Lyman Limit Systems at z ~ 6

Lyman limit systems (LLSs) are dense hydrogen clouds with high enough H i column densities to absorb Lyman continuum photons emitted from distant quasars. Their high column densities imply an origin in dense environments; however, the statistics and distribution of LLSs at high redshifts still remain uncertain. In this paper, we use self-consistent radiative transfer cosmological simulations from the Cosmic Reionization on Computers (CROC) project to study the physical properties of LLSs at the tail end of cosmic reionization at z ~ 6. We generate 3000 synthetic quasar sight lines to obtain a large number of LLS samples in the simulations. In addition, with the high physical fidelity and resolution of CROC, we are able to quantify the association between these LLS samples and nearby galaxies. Our results show that the fraction of LLSs spatially associated with nearby galaxies increases with H i column density. Moreover, we find that LLSs that are not near any galaxy typically reside in filamentary structures connecting neighboring galaxies in the intergalactic medium (IGM). This quantification of the distribution and association of LLSs to large-scale structure informs our understanding of the IGM–galaxy connection during the "Epoch of Reionization," and provides a theoretical basis for interpreting future observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Modelling the Lyman-α forest with Eulerian and SPH hydrodynamical methods

ABSTRACT We compare two state-of-the-art numerical codes to study the overall accuracy in modelling the intergalactic medium and reproducing Lyman-α forest observables for DESI and high-resolution data sets. The codes employ different approaches to solving both gravity and modelling the gas hydrodynamics. The first code, Nyx, solves the Poisson equation using the Particle-Mesh (PM) method and the Euler equations using a finite-volume method. The second code, CRK-HACC , uses a Tree-PM method to solve for gravity, and an improved Lagrangian smoothed particle hydrodynamics (SPH) technique, where fluid elements are modelled with particles, to treat the intergalactic gas. We compare the convergence behaviour of the codes in flux statistics as well as the degree to which the codes agree in the converged limit. We find good agreement overall with differences being less than observational uncertainties, and a particularly notable ≲1 per cent agreement in the 1D flux power spectrum. This agreement was achieved by applying a tessellation methodology for reconstructing the density in CRK-HACC instead of using an SPH kernel as is standard practice. We show that use of the SPH kernel can lead to significant and unnecessary biases in flux statistics; this is especially prominent at high redshifts, z ∼ 5, as the Lyman-α forest mostly comes from lower-density regions that are intrinsically poorly sampled by SPH particles.

79 ASTRONOMY AND ASTROPHYSICS↗