MF-Box: Multi-fidelity and multi-scale emulation for the matter power spectrum
Authors:
Ming-Feng Ho,
Simeon Bird,
Martin A. Fernandez,
Christian R. Shelton
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 suites (L1 and L2) at $256 \,\mathrm{Mpc}/h$ and $100 \,\mathrm{Mpc}/h$, each with $128^3$ particles, and a high-fidelity…
▽ More
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 suites (L1 and L2) at $256 \,\mathrm{Mpc}/h$ and $100 \,\mathrm{Mpc}/h$, each with $128^3$ particles, and a high-fidelity suite (HF) with $512^3$ particles at $256 \,\mathrm{Mpc}/h$, representing a higher particle load compared to the low-fidelity 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 3 HF simulations using a Gaussian process. MF-Box successfully emulates power spectra from our HF testing set with a relative error of $< 3\%$ up to $k \simeq 7 \,h/\mathrm{Mpc}$ at $z \in [0, 3]$, while maintaining a cost similar to our previous multi-fidelity approach, which was accurate only up to $z = 1$. The addition of an extra low-fidelity node in a smaller box significantly improves emulation accuracy for MF-Box at $k > 2 \,h/\mathrm{Mpc}$, 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.
△ Less
Submitted 4 October, 2023; v1 submitted 5 June, 2023;
originally announced June 2023.
Multi-Fidelity Emulation for the Matter Power Spectrum using Gaussian Processes
Authors:
Ming-Feng Ho,
Simeon Bird,
Christian R. Shelton
Abstract:
We present methods for emulating the matter power spectrum by combining information from cosmological $N$-body simulations at different resolutions. An emulator allows estimation of simulation output by interpolating across the parameter space of a limited number of simulations. We present the first implementation in cosmology of multi-fidelity emulation, where many low-resolution simulations are…
▽ More
We present methods for emulating the matter power spectrum by combining information from cosmological $N$-body simulations at different resolutions. An emulator allows estimation of simulation output by interpolating across the parameter space of a limited number of simulations. We present the first implementation in cosmology of multi-fidelity emulation, where many low-resolution simulations are combined with a few high-resolution simulations to achieve an increased emulation accuracy. The power spectrum's dependence on cosmology is learned from the low-resolution simulations, which are in turn calibrated using high-resolution simulations. We show that our multi-fidelity emulator predicts high-fidelity counterparts to percent-level relative accuracy when using only $3$ high-fidelity simulations and outperforms a single-fidelity emulator that uses $11$ simulations, although we do not attempt to produce a converged emulator with high absolute accuracy. With a fixed number of high-fidelity training simulations, we show that our multi-fidelity emulator is $\simeq 100$ times better than a single-fidelity emulator at $k \leq 2 h \mathrm{Mpc}^{-1}$, and $\simeq 20$ times better at $3 \leq k < 6.4 h \mathrm{Mpc}^{-1}$. Multi-fidelity emulation is fast to train, using only a simple modification to standard Gaussian processes. Our proposed emulator shows a new way to predict non-linear scales by fusing simulations from different fidelities.
△ Less
Submitted 27 October, 2021; v1 submitted 3 May, 2021;
originally announced May 2021.
Automated Measurement of Quasar Redshift with a Gaussian Process
Authors:
Leah Fauber,
Ming-Feng Ho,
Simeon Bird,
Christian R. Shelton,
Roman Garnett,
Ishita Korde
Abstract:
We develop an automated technique to measure quasar redshifts in the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey (SDSS). Our technique is an extension of an earlier Gaussian process method for detecting damped Lyman-alpha absorbers (DLAs) in quasar spectra with known redshifts. We apply this technique to a subsample of SDSS DR12 with BAL quasars removed and redsh…
▽ More
We develop an automated technique to measure quasar redshifts in the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey (SDSS). Our technique is an extension of an earlier Gaussian process method for detecting damped Lyman-alpha absorbers (DLAs) in quasar spectra with known redshifts. We apply this technique to a subsample of SDSS DR12 with BAL quasars removed and redshift larger than 2.15. We show that we are broadly competitive to existing quasar redshift estimators, disagreeing with the PCA redshift by more than 0.5 in only 0.38% of spectra. Our method produces a probabilistic density function for the quasar redshift, allowing quasar redshift uncertainty to be propagated to downstream users. We apply this method to detecting DLAs, accounting in a Bayesian fashion for redshift uncertainty. Compared to our earlier method with a known quasar redshift, we have a moderate decrease in our ability to detect DLAs, predominantly in the noisiest spectra. The area under curve drops from 0.96 to 0.91. Our code is publicly available.
△ Less
Submitted 15 September, 2020; v1 submitted 12 June, 2020;
originally announced June 2020.
