Liu, Jia | Engineering Papers

DOE OSTI2025

The HalfDome multi-survey cosmological simulations: N-body simulations

Upcoming cosmological surveys have the potential to reach groundbreaking discoveries on multiple fronts, including the neutrino mass, dark energy, and inflation. Most of the key science goals require the joint analysis of datasets from multiple surveys to break parameter degeneracies and calibrate systematics. To realize such analyses, a large set of mock simulations that realistically model correlated observables is required. In this paper we present the N-body component of the HalfDome cosmological simulations, designed for the joint analysis of Stage-IV cosmological surveys, such as Rubin LSST, Euclid, SPHEREx, Roman, DESI, PFS, Simons Observatory, CMB-S4, and LiteBIRD. Our 300TB initial data release includes full-sky lightcones and halo catalogs between z = 0–4 for 11 fixed cosmology realizations, as well as an additional run with local primordial non-Gaussianity (f NL = 20). The simulations evolve 61443 particles in a 3.75 h -1 Gpc box, reaching a minimum halo mass of ∼6 × 1012 h -1 M ⊙ and maximum scale of k ∼ h Mpc-1. Our data is publicly available: instructions to access the data and plans for future data releases can be found at https://halfdomesims.github.io.

Bayer, Adrian E↗

DOE OSTI2024

Impact of baryonic feedback on HSC-Y1 weak lensing non-Gaussian statistics

Baryonic feedback is a major systematic in weak lensing cosmology. Its most studied effect is the suppression of the lensing power spectrum, a second-order statistic, on small scales. Motivated by the growing interest in statistics beyond the second order, we investigate the effect of baryons on lensing non-Gaussian statistics and the resulting biases in the matter clustering amplitude S 8 = σ 8 Ω m / 0.3 . We focus on the Subaru Hyper Suprime-Cam Year 1 (HSC-Y1) data which, with its high source number density, closely resembles those expected from the upcoming Euclid and Rubin Legacy Survey of Space and Time. We study four non-Gaussian statistics of convergence maps—peak counts, minimum counts, the probability distribution function, and the scattering transform coefficients—in addition to the usual power spectrum. We first estimate the biases in S 8 using mock observations built from the IllustrisTNG and BAHAMAS hydrodynamical simulations and theoretical models built from dark-matter-only simulations. We find up to 1 σ bias in S 8 when the smallest scales (2 arcmin) and the highest feedback level are considered. We then analyze the HSC-Y1 data and compare the S 8 obtained for each statistic with different smoothing scales or scale cuts. As we expect that baryons mostly affect the small scales, comparing the results obtained from including and excluding small scales can indicate the level of impact from baryons. With HSC data, we find only minor ( ≤ 0.5 σ ) differences in S 8 for all statistics, even when considering very small scales (2 arcmin). Our results suggest that the effect of baryons is insignificant at the level of HSC-Y1 down to 2 arcmin for all statistics examined here, or it is canceled by other scale-dependent systematics.

79 ASTRONOMY AND ASTROPHYSICS↗