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:
- 10534751
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 526
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 2903 to 2919
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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 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
-
Accurate assessment of driver visibility is crucial in automotive design and safety enhancement, particularly in situations where A-pillars obstruct the driver’s field of view. To address this challenge, this research develops a multi-fidelity Gaussian Process (MF-GP) modeling framework to enhance visibility prediction by integrating low-fidelity (LF) image segmentation data with high-fidelity digital human modeling (DHM) simulations. By leveraging a limited set of high-fidelity samples, the proposed MF-GP framework systematically calibrates low-fidelity data to improve predictive accuracy while reducing computational costs. Two A-pillar cutout designs (3.75 cm and 5 cm) were analyzed under varying HF sampling densities of 3%, 7%, and 10%. Results indicate that the 3.75 cm cutout is more sensitive to sparse HF sampling, requiring a denser HF dataset to achieve stable calibration. In contrast, the 5 cm cutout, benefiting from improved LF-HF alignment, achieves comparable accuracy with fewer HF samples. Model validation using root mean square error (RMSE) and coefficient of determination (R2) confirms that increasing HF sampling enhances surrogate model accuracy, with the effect being more pronounced in cases where model performance is susceptible to high-fidelity data. The proposed framework provides a computationally efficient methodology for driver visibility prediction and human-in-the-loop design applications. Future research could explore adaptive HF sampling strategies and ensemble surrogate modeling techniques to further enhance multi-fidelity learning efficiency.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
An official website of the United States government

