An official website of the United States government
ABSTRACT In this work, we extend our recently developed multifidelity emulation technique to the simulated Lyman-α forest flux power spectrum. Multifidelity emulation allows interpolation of simulation outputs between cosmological parameters using many cheap low-fidelity simulations and a few expensive high-fidelity simulations. Using a test suite of small-box (30 Mpc h−1) simulations, we show that multifidelity emulation is able to reproduce the Lyman-α forest flux power spectrum well, achieving an average accuracy when compared to a test suite of $$0.8\, {\rm {per\ cent}}$$. We further show that it has a substantially increased accuracy over single-fidelity emulators, constructed using either the high- or low-fidelity simulations only. In particular, it allows the extension of an existing simulation suite to smaller scales and higher redshifts.
more »
« less
MF-Box: multifidelity and multiscale emulation for the matter power spectrum
ABSTRACT We introduce MF-Box, an extended version of MFEmulator, designed as a fast surrogate for power spectra, trained using N-body simulation suites from various box sizes and particle loads. To demonstrate MF-Box’s effectiveness, we design simulation suites that include low-fidelity (LF) suites (L1 and L2) at 256 and $$100 \, \rm {Mpc\, ~}h^{-1}$$, each with 1283 particles, and a high-fidelity (HF) suite with 5123 particles at $$256 \, \rm {Mpc\, ~}h^{-1}$$, representing a higher particle load compared to the LF suites. MF-Box acts as a probabilistic resolution correction function, learning most of the cosmological dependencies from L1 and L2 simulations and rectifying resolution differences with just three HF simulations using a Gaussian process. MF-Box successfully emulates power spectra from our HF testing set with a relative error of $$\lt 3~{{\ \rm per\ cent}}$$ up to $$k \simeq 7 \, h\rm {Mpc}{^{-1}}$$ at z ∈ [0, 3], while maintaining a cost similar to our previous multifidelity approach, which was accurate only up to z = 1. The addition of an extra LF node in a smaller box significantly improves emulation accuracy for MF-Box at $$k \gt 2 \, h\rm {Mpc}{^{-1}}$$, increasing it by a factor of 10. We conduct an error analysis of MF-Box based on computational budget, providing guidance for optimizing budget allocation per fidelity node. Our proposed MF-Box enables future surveys to efficiently combine simulation suites of varying quality, effectively expanding the range of emulation capabilities while ensuring cost efficiency.
more »
« less
- Award ID(s):
- 2215705
- PAR ID:
- 10468250
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 526
- Issue:
- 2
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 2903-2919
- Size(s):
- p. 2903-2919
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Bayesian optimization (BO) is a sequential optimization strategy that is increasingly employed in a wide range of areas including materials design. In real world applications, acquiring high-fidelity (HF) data through physical experiments or HF simulations is the major cost component of BO. To alleviate this bottleneck, multi-fidelity (MF) methods are increasingly used to forgo the sole reliance on the expensive HF data and reduce the sampling costs by querying inexpensive low-fidelity (LF) sources whose data are correlated with HF samples. Existing multi-fidelity BO (MFBO) methods operate under the following two assumptions: (1) Leveraging global (rather than local) correlation between HF and LF sources, and (2) Associating all the data sources with the same noise process. These assumptions dramatically reduce the performance of MFBO when LF sources are only locally correlated with the HF source or when the noise variance varies across the data sources. To dispense with these incorrect assumptions, we propose an MF emulation method that (1) learns a noise model for each data source, and (2) enables BO to leverage highly biased LF sources which are only locally correlated with the HF source. We illustrate the performance of our method through analytical examples and engineering problems on materials design.more » « less
-
ABSTRACT We implement a model for the two-point statistics of biased tracers that combines dark matter dynamics from N-body simulations with an analytic Lagrangian bias expansion. Using Aemulus, a suite of N-body simulations built for emulation of cosmological observables, we emulate the cosmology dependence of these non-linear spectra from redshifts z = 0 to z = 2. We quantify the accuracy of our emulation procedure, which is sub-per cent at $$k=1\, h \,{\rm Mpc}^{-1}$$ for the redshifts probed by upcoming surveys and improves at higher redshifts. We demonstrate its ability to describe the statistics of complex tracer samples, including those with assembly bias and baryonic effects, reliably fitting the clustering and lensing statistics of such samples at redshift z ≃ 0.4 to scales of $$k_{\rm max} \approx 0.6\, h\,\mathrm{Mpc}^{-1}$$. We show that the emulator can be used for unbiased cosmological parameter inference in simulated joint clustering and galaxy–galaxy lensing analyses with data drawn from an independent N-body simulation. These results indicate that our emulator is a promising tool that can be readily applied to the analysis of current and upcoming data sets from galaxy surveys.more » « less
-
ABSTRACT In this work, we establish and test methods for implementing dynamical friction (DF) for massive black hole pairs that form in large volume cosmological hydrodynamical simulations that include galaxy formation and black hole growth. We verify our models and parameters both for individual black hole dynamics and for the black hole population in cosmological volumes. Using our model of DF from collisionless particles, black holes can effectively sink close to the galaxy centre, provided that the black hole’s dynamical mass is at least twice that of the lowest mass resolution particles in the simulation. Gas drag also plays a role in assisting the black holes’ orbital decay, but it is typically less effective than that from collisionless particles, especially after the first billion years of the black hole’s evolution. DF from gas becomes less than $$1{{\ \rm per\ cent}}$$ of DF from collisionless particles for BH masses >107 M⊙. Using our best DF model, we calculate the merger rate down to z = 1.1 using an Lbox = 35 Mpc h−1 simulation box. We predict ∼2 mergers per year for z > 1.1 peaking at z ∼ 2. These merger rates are within the range obtained in previous work using similar resolution hydrodynamical simulations. We show that the rate is enhanced by factor of ∼2 when DF is taken into account in the simulations compared to the no-DF run. This is due to $${\gt}40{{\ \rm per\ cent}}$$ more black holes reaching the centre of their host halo when DF is added.more » « less
-
null (Ed.)ABSTRACT Galaxy intrinsic alignments (IAs) have long been recognized as a significant contaminant to weak lensing-based cosmological inference. In this paper we seek to quantify the impact of a common modelling assumption in analytic descriptions of IAs: that of spherically symmetric dark matter haloes. Understanding such effects is important as the current generation of IA models are known to be limited, particularly on small scales, and building an accurate theoretical description will be essential for fully exploiting the information in future lensing data. Our analysis is based on a catalogue of 113 560 galaxies between z = 0.06 and 1.00 from massiveblack-ii, a hydrodynamical simulation of box length $$100 \, h^{-1}$$ Mpc. We find satellite anisotropy contributes at the level of $$\ge 30\!-\!40{{\ \rm per\ cent}}$$ to the small-scale alignment correlation functions. At separations larger than $$1 \, h^{-1}$$ Mpc the impact is roughly scale independent, inducing a shift in the amplitude of the IA power spectra of $$\sim 20{{\ \rm per\ cent}}$$. These conclusions are consistent across the redshift range and between the massiveblack-ii and the illustris simulations. The cosmological implications of these results are tested using a simulated likelihood analysis. Synthetic cosmic shear data are constructed with the expected characteristics (depth, area, and number density) of a future LSST-like survey. Our results suggest that modelling alignments using a halo model based upon spherical symmetry could potentially induce cosmological parameter biases at the ∼1.5σ level for S8 and w.more » « less
