Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cosmological hydrodynamical simulation”

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 91 records · Page 5

Demonstrating the viability of Lagrangian in situ reduction on supercomputers

Performing exploratory analysis and visualization of large-scale time-varying computational science applications is challenging due to inaccuracies that arise from under-resolved data. In recent years, Lagrangian representations of the vector field computed using in situ processing are being increasingly researched and have emerged as a potential solution to enable exploration. However, prior works have offered limited estimates of the encumbrance on the simulation code as they consider “theoretical” in situ environments. Further, the effectiveness of this approach varies based on the nature of the vector field, benefitting from an in-depth investigation for each application area. With this study, an extended version of Sane et al. (2021), we contribute an evaluation of Lagrangian analysis viability and efficacy for simulation codes executing at scale on a supercomputer. We investigated previously unexplored cosmology and seismology applications as well as conducted a performance benchmarking study by using a hydrodynamics mini-application targeting exascale computing. Here, to inform encumbrance, we integrated in situ infrastructure with simulation codes, and evaluated Lagrangian in situ reduction in representative homogeneous and heterogeneous HPC environments. To inform post hoc accuracy, we conducted a statistical analysis across a range of spatiotemporal configurations as well as a qualitative evaluation. Additionally, our study contributes cost estimates for distributed-memory post hoc reconstruction. In all, we demonstrate viability for each application — data reduction to less than 1% of the total data via Lagrangian representations, while maintaining accurate reconstruction and requiring under 10% of total execution time in over 90% of our experiments.

97 MATHEMATICS AND COMPUTING↗

The aluminium-26 distribution in a cosmological simulation of a Milky Way-type Galaxy

Context . The 1.8 MeV γ -rays corresponding to the decay of the radioactive isotope 26 Al (with a half-life of 0.72 Myr ) have been observed by the SPI detector on the INTEGRAL spacecraft and extensively used as a tracer of star formation and current nucleosynthetic activity in the Milky Way Galaxy. Further information is encoded in the observation related to the higher 26 Al content found in regions of the Galaxy with the highest line-of-sight (LoS) velocity relative to an observer located in the Solar System. However, this feature remains unexplained. Aims . We ran a cosmological “zoom-in” chemodynamical simulation of a Milky Way-type galaxy, including the production and decays of radioactive nuclei in a fully self-consistent way. We then analyzed the results to follow the evolution of 26 Al throughout the lifetime of the simulated galaxy to provide a new method for interpreting the 26 Al observations. Methods . We included the massive star sources of 26 Al in the Galaxy and its radioactive decay into a state-of-the-art galactic chemical evolution model, coupled with cosmological growth and hydrodynamics. This approach allowed us to follow the spatial and temporal evolution of the 26 Al content in the simulated galaxy. Results . Our results are in agreement with the observations with respect to the fact that gas particles in the simulation with relatively higher 26 Al content also have the highest LoS velocities. On the other hand, gas particles with relatively lower 26 Al content (i.e., not bright enough to be observed) generally display the lowest LoS velocities. However, this result is not conclusive because the overall rotational velocity of our simulated galaxy is higher than that observed for cold CO gas in the Milky Way Galaxy. Furthermore, we found no significant correlation between gas temperature, rotational velocity, and 26 Al content at any given radius. We also found the presence of transient 26 Al-rich spots at low LoS velocities and we show that one such spot had been captured by the INTEGRAL/SPI data. Based on our model, we present a prediction for the detection of 1.8 MeV γ -rays by the future COSI mission. We find that according to our model, the new instrument will be able to observe similar 26 Al-emission patterns to those seen by INTEGRAL/SPI.

Wehmeyer, B. (ORCID:0009000254149292)↗

The AGORA High-resolution Galaxy Simulations Comparison Project. IV. Halo and Galaxy Mass Assembly in a Cosmological Zoom-in Simulation at z ≤ 2

In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift z = 2 and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at z = 4 and 3, and before the last major merger, focusing on the formation of well-defined rotationally supported disks, the mass–metallicity relation, the specific star formation rate, the gas metallicity gradients, and the nonaxisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semianalytic models at z ~ 2. We see that almost all the hydro codes develop rotationally supported structures at low redshifts. Most agree within 0.5 dex with the observed mass–metallicity relation at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the intercode differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to z = 2 and the auxiliary data will be made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Combining Planck and SPT Cluster Catalogs: Cosmological Analysis and Impact on the Planck Scaling Relation Calibration

We provide the first combined cosmological analysis of the South Pole Telescope (SPT) and Planck cluster catalogs. The aim is to provide an independent calibration for Planck scaling relations, exploiting the cosmological constraining power of the SPT-SZ cluster catalog and its dedicated weak lensing (WL) and X-ray follow-up observations. We build a new version of the Planck cluster likelihood. In the νΛ CDM scenario, focusing on the mass slope and mass bias of Planck scaling relations, we find ${\alpha }_{\mathrm{SZ}}={1.49}_{-0.10}^{+0.07}$ and ${\left(1-b\right)}_{\mathrm{SZ}}={0.69}_{-0.14}^{+0.07}$, respectively. The results for the mass slope show a ~4 σ departure from the self-similar evolution, α SZ ~ 1.8. This shift is mainly driven by the matter density value preferred by SPT data, Ωm = 0.30 ± 0.03, lower than the one obtained by Planck data alone, ${{\rm{\Omega }}}_{m}={0.37}_{-0.06}^{+0.02}$. The mass bias constraints are consistent both with outcomes of hydrodynamical simulations and external WL calibrations, (1 – b) ~ 0.8, and with results required by the Planck cosmic microwave background cosmology, (1 – b) ~ 0.6. From this analysis, we obtain a new catalog of Planck cluster masses M 500 . We estimate the ratio between the published Planck M SZ masses and our derived masses M 500 , as a "measured mass bias," ${\left(1-b\right)}_{M}$. We analyze the mass, redshift, and detection noise dependence of ${\left(1-b\right)}_{M}$, finding an increasing trend toward high redshift and low mass. These results mimic the effect of departure from self-similarity in cluster evolution, showing different dependencies for the low-mass, high-mass, low-z, and high-z regimes.

79 ASTRONOMY AND ASTROPHYSICS↗

Dwarf Galaxies at Cosmic Noon: New JWST Constraints on Satellite Models and Subhalo Tidal Evolution

The advent of JWST has revolutionized the study of faint satellite galaxies at z ≳ 1, enabling statistical constraints on galaxy evolution and the galaxy–halo connection in a previously unexplored mass and redshift regime. We compare satellite abundances at 1 < z < 3.5 from recent JWST observations with predictions from cosmological dark-matter-only zoom-in simulations. We identify and quantify several sources of biases that can impact theoretical satellite counts, finding that assumptions about subhalo tidal evolution introduce the largest uncertainty in predictions for the satellite mass function. Using a flexible galaxy disruption model, we explore a range of disruption scenarios, spanning hydrodynamically motivated and idealized prescriptions, to bracket plausible physical outcomes. We show that varying galaxy durability can change the predicted satellite mass functions by a factor of ∼3.5. The JWST data and our fiducial model are consistent within 1σ–2σ across the full redshift (1 < z < 3.5) and stellar mass (M⋆ > 10 7 M ⊙ ) range probed. We find evidence that subhalos are at least as long-lived as predicted by hydrodynamic simulations. Our framework will enable robust constraints on the tidal evolution of subhalos with future observations. This work presents the first direct comparison between cosmological models and observations of the high-redshift satellite population in this low-mass regime. These results showcase JWST’s emerging power to test structure formation in the first half of the Universe in a new domain and to constrain the physical processes driving the evolution of low-mass galaxies across cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗

The R -process Alliance: Enrichment of r -process Elements in a Simulated Milky Way–like Galaxy

We study the formation of stars with varying amounts of heavy elements synthesized by the rapid neutron-capture process (r-process) based on our detailed cosmological zoom-in simulation of a Milky Way–like galaxy with an N-body/smoothed particle hydrodynamics code, asura. Most stars with no overabundance in r-process elements, as well as the strongly r-process-enhanced (RPE) r-II stars ([Eu/Fe] > +0.7), are formed in dwarf galaxies accreted by the Milky Way within the 6 Gyr after the Big Bang. In contrast, over half of the moderately enhanced r-I stars (+0.3 < [Eu/Fe] ≤ +0.7) are formed in the main in situ disk after 6 Gyr. Our results suggest that the fraction of r-I and r-II stars formed in disrupted dwarf galaxies is larger the higher their [Eu/Fe] is. Accordingly, the most strongly enhanced r-III stars ([Eu/Fe] > +2.0) are formed in accreted components. These results suggest that non-r-process-enhanced stars and r-II stars are mainly formed in low-mass dwarf galaxies that hosted either none or a single neutron star merger, while the r-I stars tend to form in the well-mixed in situ disk. We compare our findings with high-resolution spectroscopic observations of RPE metal-poor stars in the halo and dwarf galaxies, including those collected by the R-Process Alliance. We conclude that observed [Eu/Fe] and [Eu/Mg] ratios can be employed in chemical tagging of the Milky Way’s accretion history.

Chemical enrichment↗

Evidence for large baryonic feedback at low and intermediate redshifts from kinematic Sunyaev-Zel’dovich observations with ACT and DESI photometric galaxies

Recent advances in cosmological observations have provided an unprecedented opportunity to investigate the distribution of baryons relative to the underlying matter. In this work, we show that the gas is more extended than the dark matter, and the amount of baryonic feedback at $z \lesssim 1$ disfavors low-feedback models such as that of state-of-the-art hydrodynamical simulation IllustrisTNG compared with high-feedback models such as that of the original Illustris simulation. This has important implications for bridging the gap between theory and observations and understanding galaxy formation and evolution. Furthermore, a better grasp of the baryon-dark matter link is critical to future cosmological analyses, which are currently impeded by our limited knowledge of baryonic feedback. Here, we measure the kinematic Sunyaev-Zel'dovich (kSZ) effect from the Atacama Cosmology Telescope (ACT), stacked on the luminous red galaxy (LRG) sample of the Dark Energy Spectroscopic Instrument (DESI) imaging survey. This is the first analysis to use photometric redshifts for reconstructing galaxy velocities. Due to the large number of galaxies comprising the DESI imaging survey, this is the highest signal-to-noise stacked kSZ measurement to date: we detect the signal at 13$σ$, finding strong evidence that the gas is more spread out than the dark matter, as well as a preference for larger feedback compared to some commonly used state-of-the-art hydrodynamical simulations. Here, our work opens up the possibility of recalibrating large hydrodynamical simulations using the kSZ effect. In addition, our findings highlight the importance of properly accounting for baryonic feedback with future surveys such as LSST through direct probes such as the kSZ, and shed light on long-standing enigmas in astrophysics, such as the “missing baryon” problem.

79 ASTRONOMY AND ASTROPHYSICS↗

Pericentric passage-driven star formation in satellite galaxies and their hosts: CLUES from local group simulations

Local Group satellite galaxies show a wide diversity of star formation histories (SFHs) whose origin is yet to be fully understood. Using hydrodynamical simulations from the Constrained Local UniversE project, we study the SFHs of satellites of Milky Way-like galaxies in a cosmological context: while in the majority of the cases the accretion on to their host galaxy causes the satellites to lose their gas, with a subsequent suppression in SF, in about 25 per cent of our sample we observe a clear enhancement of SF after infall. Peaks in SF clearly correlate with the satellite pericentric passage around its host and, in one case, with a satellite–satellite interaction. We identify two key ingredients that result in enhanced SF after infall: galaxies must enter the host’s virial radius with a reservoir of cold gas M gas , inf/M vir , inf ≳ 10 -2 and with a minimum pericentric distance ≳10 kpc (mean distance ~50 kpc for the full sample), in order to form new stars due to compression of cold gas at pericentric passage. On the other hand, satellites that infall with little gas or whose pericentric distance is too small, have their gas ram-pressure stripped and subsequent SF quenched. The pericentric passage of satellites likewise correlates with SF peaks in their hosts, suggesting that this mechanism induces bursts of SF in satellites and central galaxies alike, in agreement with recent studies of our Galaxy’s SFH. Our research findings can explain the recently reported multiple stellar populations observed in dwarf galaxies such as Carina and Fornax, and should be taken into account in semi-analytic models of galaxy formation and satellite quenching.

79 ASTRONOMY AND ASTROPHYSICS↗

Mitigating baryonic effects with a theoretical error covariance

ABSTRACT One of the primary sources of uncertainties in modelling the cosmic-shear power spectrum on small scales is the effect of baryonic physics. Accurate cosmology for stage-IV surveys requires knowledge of the matter power spectrum deep in the non-linear regime at the per cent level. Therefore, it is important to develop reliable mitigation techniques to take into account baryonic uncertainties if information from small scales is to be considered in the cosmological analysis. In this work, we develop a new mitigation method for dealing with baryonic physics for the case of the shear angular power spectrum. The method is based on an augmented covariance matrix that incorporates baryonic uncertainties informed by hydrodynamical simulations. We use the results from 13 hydrodynamical simulations and the residual errors arising from a fit to a ΛCDM model using the extended halo model code HMCode to account for baryonic physics. These residual errors are used to model a so-called theoretical error covariance matrix that is added to the original covariance matrix. In order to assess the performance of the method, we use the 2D tomographic shear from four hydrodynamical simulations that have different extremes of baryonic parameters as mock data and run a likelihood analysis comparing the residual bias on Ωm and σ8 of our method and the HMCode for an LSST-like survey. We use different modelling of the theoretical error covariance matrix to test the robustness of the method. We show that it is possible to reduce the bias in the determination of the tested cosmological parameters at the price of a modest decrease in the precision.

79 ASTRONOMY AND ASTROPHYSICS↗

Fast, High-fidelity Lyα Forests with Convolutional Neural Networks

Full-physics cosmological simulations are powerful tools for studying the formation and evolution of structure in the universe but require extreme computational resources. Here, we train a convolutional neural network to use a cheaper N-body-only simulation to reconstruct the baryon hydrodynamic variables (density, temperature, and velocity) on scales relevant to the Lyα forest, using data from Nyx simulations. We show that our method enables rapid estimation of these fields at a resolution of ~20 kpc, and captures the statistics of the Lyα forest with much greater accuracy than existing approximations. Because our model is fully convolutional, we can train on smaller simulation boxes and deploy on much larger ones, enabling substantial computational savings. Furthermore, as our method produces an approximation for the hydrodynamic fields instead of Lyα flux directly, it is not limited to a particular choice of ionizing background or mean transmitted flux.

79 ASTRONOMY AND ASTROPHYSICS↗

Effective bias expansion for 21-cm cosmology in redshift space

A near-future detection of the 21-cm signal from the epoch of reionization will provide unique opportunities to probe the underlying cosmology, provided that such cosmological information can be extracted with precision. To this end, we further develop effective field theory (EFT) inspired techniques for the 21-cm brightness temperature field during the epoch of reionization, incorporating renormalized bias and a treatment of redshift space distortions. Notably, we confirm that in redshift space, measures of the 21-cm brightness, e.g., the power spectrum, should have irreducible contributions that lack a bias coefficient and therefore contain direct, astrophysics-free information about the cosmological density field; in this work, we study this effect beyond linear order. To validate our theoretical treatment, we fit the predicted EFT Fourier-space shapes to the thesan suite of hydrodynamical simulations of reionization at the field level, where the considerable number of modes prevents overfitting. We find agreement at the level of a few percent between the 21-cm power spectrum from the EFT fits and simulations over the wave number range k ≲ 0.8 h/Mpc and neutral fraction x HI ≳ 0.4, which is imminently measurable by the Hydrogen Epoch of Reionization Array and future experiments. Furthermore, the ability of the EFT to describe the 21-cm signal extends to simulations that have different astrophysical prescriptions for reionization as well as simulations with interacting dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗

Galaxy cluster profiles: a Gaussian mixture model approach to halo miscentering

Measurements of the galaxy density and weak-lensing profiles of galaxy clusters typically rely on an assumed cluster center, which is taken to be the brightest cluster galaxy or other proxies for the true halo center defined as the minimum in the potential well. Departure of the assumed cluster center from the true halo center bias the resultant profile measurements, an effect known as miscentering bias. Currently, miscentering is typically modeled in stacked profiles of clusters with a two parameter model. We use an alternate approach in which the profiles of individual clusters are used with the corresponding likelihood computed using a Gaussian mixture model. We test the approach using halos and the corresponding subhalo profiles from the IllustrisTNG hydrodynamic simulations. We obtain significantly improved estimates of the miscentering parameters for both 3D and projected 2D profiles relevant for imaging surveys. We discuss applications to upcoming cosmological surveys. Our Python package for the Gaussian mixture model is publicly available at https://github.com/KyleMiller1/Halo-Miscentering-Mixture-Model.

Bayesian reasoning↗

Probing cosmic velocities with the pairwise kinematic Sunyaev-Zel’dovich signal in DESI Bright Galaxy Sample DR1 and ACT DR6

We present a measurement of the pairwise kinematic Sunyaev-Zel’dovich (kSZ) signal using the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Sample (BGS) Data Release 1 (DR1) galaxy sample overlapping with the Atacama Cosmology Telescope (ACT) CMB temperature map. Our analysis makes use of 1.6 million galaxies with stellar masses log⁡ 𝑀 ⋆ /𝑀 ⊙ >10, and we explore measurements across a range of aperture sizes (2.1′ <𝜃 ap <3.5′) and stellar mass selections. This statistic directly probes the velocity field of the large-scale structure, a unique observable of cosmic dynamics and modified gravity. In particular, at low redshifts, this quantity is especially interesting, as deviations from General Relativity are expected to be largest. Notably, our result represents the highest-significance low-redshift (𝑧 ∼ 0.3) detection of the kSZ pairwise effect yet. In our most optimal configuration (𝜃 ap =3.3′, log⁡ 𝑀 ⋆ >11), we achieve a 5⁢𝜎 detection. Assuming that an estimate of the optical depth and galaxy bias of the sample exists via e.g., external observables, this measurement constrains the fundamental cosmological combination 𝐻 0 ⁡𝑓⁡𝜎$^2_8$. A key challenge is the degeneracy with the galaxy optical depth. We address this by combining CMB lensing, which allows us to infer the halo mass and galaxy population properties, with hydrodynamical simulation estimates of the mean optical depth, $\bar{𝜏}$ . We stress that this is a proof-of-concept analysis; with BGS DR2 data we expect to improve the statistical precision by roughly a factor of two, paving the way toward robust tests of modified gravity with kSZ-informed velocity-field measurements at low redshift.

Hadzhiyska, Boryana [Institute of Astronomy; Kavli↗

Cosmology with second- and third-order shear statistics for the Dark Energy Survey: Methods and simulated analysis

We present a new pipeline designed for the robust inference of cosmological parameters using both second- and third-order shear statistics. We build a theoretical model for rapid evaluation of three-point correlations using our fastnc code and integrate it into the cosmosis framework. We measure the two-point functions 𝜉 ± and the full configuration-dependent three-point shear correlation functions across all auto- and cross-redshift bins. We compress the three-point functions into the mass aperture statistic ⟨ℳ$^{3}_{ap}$⟩ for a set of 796 simulated shear maps designed to model the Dark Energy Survey Year 3 data. We estimate from it the full covariance matrix and model the effects of intrinsic alignments, shear calibration biases and photometric redshift uncertainties. We apply scale cuts to minimize the contamination from the baryonic signal as modeled through hydrodynamical simulations. We find a significant improvement of 83% on the figure of merit in the Ω m − 𝑆 8 plane when we add the ⟨ℳ$^{3}_{ap}$⟩ data to 𝜉 ± . Here, we present our findings for all relevant cosmological and systematic uncertainty parameters and discuss the complementarity of third-order and second-order statistics.

79 ASTRONOMY AND ASTROPHYSICS↗

A non-linear solution to the S 8 tension – II. Analysis of DES Year 3 cosmic shear

ABSTRACT Weak galaxy lensing surveys have consistently reported low values of the S8 parameter compared to the Planck lambda cold dark matter (ΛCDM) cosmology. Amon & Efstathiou used KiDS-1000 cosmic shear measurements to propose that this tension can be reconciled if the matter fluctuation spectrum is suppressed more strongly on non-linear scales than assumed in state-of-the-art hydrodynamical simulations. In this paper, we investigate cosmic shear data from the Dark Energy Survey (DES) Year 3. The non-linear suppression of the matter power spectrum required to resolve the S8 tension between DES and the Planck ΛCDM model is not as strong as inferred using KiDS data, but is still more extreme than predictions from recent numerical simulations. An alternative possibility is that non-standard dark matter contributes to the required suppression. We investigate the redshift and scale dependence of the suppression of the matter power spectrum. If our proposed explanation of the S8 tension is correct, the required suppression must extend into the mildly non-linear regime to wavenumbers $k\sim 0.2 \, h\, {\rm Mpc}^{-1}$. In addition, all measures of S8 using linear scales should agree with the Planck ΛCDM cosmology, an expectation that will be testable to high precision in the near future.

Preston, Calvin (ORCID:0000000281136916)↗

On the possibility of baryon acoustic oscillation measurements at redshift z > 7.6 with the Roman space telescope

ABSTRACT The Nancy Grace Roman Space Telescope (RST), with its field of view and high sensitivity will make surveys of cosmological large-scale structure possible at high redshifts. We investigate the possibility of detecting baryon acoustic oscillations (BAO) at redshifts z > 7.6 for use as a standard ruler. We use data from the hydrodynamic simulation bluetides in conjunction with the gigaparsec-scale Outer Rim simulation and a model for patchy reionization to create mock RST High Latitude Survey grism data for Lyman α emission line selected galaxies at redshifts z = 7.4 to z = 10, covering 2280 deg2. We measure the monopoles of galaxies in the mock catalogues and fit the BAO features. We find that for a line flux of $L = 7\times 10^{-17} \ {\rm erg\, s^{-1}\, cm}^{-2}$, the 5σ detection limit for the current design, the BAO feature is partially detectable (measured in three out of four survey quadrants analysed independently). The resulting root mean square error on the angular diameter distance to z = 7.7 is 7.9 ${{\ \rm per\ cent}}$. If we improve the detection sensitivity by a factor of two (i.e. $L = 3.5\times 10^{-17} \ {\rm erg\, s^{-1}\, cm}^{-2}$), the distance error reduces to $1.4{{\ \rm per\ cent}}$. We caution that many more factors are yet to be modelled, including dust obscuration, the damping wing due to the intergalactic medium, and low redshift interlopers. If these issues do not strongly affect the results, or different observational techniques (such as use of multiple lines) can mitigate them, RST, or similar instruments may be able to constrain the angular diameter distance to the high redshift universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Accurate predictions from small boxes: variance suppression via the Zel'dovich approximation

Simulations have become an indispensable tool for accurate modelling of observables measured in galaxy surveys, but can be expensive if very large dynamic range in scale is required. We describe how to combine Lagrangian perturbation theory models with N-body simulations to reduce the effects of finite computational volume in the prediction of ensemble average properties in the simulations within the context of control variates. In particular we use the fact that Zel'dovich displacements, computed during initial condition generation for any simulation, correlate strongly with the final density field. Since all the correlators of biased tracers can be computed with arbitrary precision for these displacements, pairing the Zel'dovich `simulation' with the N-body realization allows hundredfold reductions in sample variance for power spectrum or correlation function estimation. Additionally, Zel'dovich control variates can accurately extend matter or tracer field emulators to larger scales than previously possible, as well as improving measurements of statistics in simulations which are inherently limited to small volumes, such as hydrodynamical simulations of galaxy formation and reionization.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-linearities in the Lyman-α forest and in its cross-correlation with dark matter halos

Here three-dimensional correlations of the Lyman-α (Lyα) forest and cross correlations between the Lyα forest and quasars have been measured on large scales, allowing a precise measurement of the baryon acoustic oscillation (BAO) feature at redshifts z > 2. These 3D correlations are often modelled using linear perturbation theory, but full-shape analyses to extract cosmological information beyond BAO will require more realistic models capable of describing non-linearities present at smaller scales. We present a measurement of the Lyα forest flux power spectrum from large hydrodynamic simulations — the Sherwood simulations — and compare it to different models describing the small-scale deviations from linear theory. We confirm that the model presented in Arinyo-i-Prats et al. (2015) fits the measured 3D power up to k = 10 h Mpc -1 with an accuracy better than 5%, and show that the same model can also describe the 1D correlations with similar precision. We also present, for the first time, an equivalent study for the cross-power spectrum of halos with the Lyα forest, and we discuss different challenges we face when modelling the cross-power spectrum beyond linear scales. We make all our measured power spectra public in https://github.com/andreufont/sherwoo_p3d. This study is a step towards joint analyses of 1D and 3D flux correlations, and towards using the quasar-Lyα cross-correlation beyond BAO analyses.

79 ASTRONOMY AND ASTROPHYSICS↗