-
Cosmological Constraints from the eBOSS Lyman-$α$ Forest using the PRIYA Simulations
Authors:
M. A. Fernandez,
Simeon Bird,
Ming-Feng Ho
Abstract:
We present new cosmological parameter constraints from the eBOSS Lyman-$α$ forest survey. We use a new theoretical model and likelihood based on the PRIYA simulation suite. PRIYA is the first suite to resolve the Lyman-$α$ forest in a ($120$ Mpc/h)$^3$ volume, using a multi-fidelity emulation technique. We use PRIYA to predict Lyman-$α$ forest observables with $\lesssim 1\%$ interpolation error ov…
▽ More
We present new cosmological parameter constraints from the eBOSS Lyman-$α$ forest survey. We use a new theoretical model and likelihood based on the PRIYA simulation suite. PRIYA is the first suite to resolve the Lyman-$α$ forest in a ($120$ Mpc/h)$^3$ volume, using a multi-fidelity emulation technique. We use PRIYA to predict Lyman-$α$ forest observables with $\lesssim 1\%$ interpolation error over an $11$ dimensional ($9$ simulated, $2$ in post-processing) parameter space. We identify an internal tension within the flux power spectrum data. Once the discrepant data is removed, we find the scalar spectral index at $k = 0.78$ h/Mpc to be $n_P = 0.97 - 0.995$ at $68\%$ confidence from the Lyman-$α$ forest flux power spectrum alone, in good agreement with Planck. The amplitude of matter fluctuations is $σ_8 = 0.733 \pm 0.026$ at $68\%$ confidence, in agreement with Dark Energy Survey weak lensing measurements and other small-scale structure probes and in tension with CMB measurements from Planck and ACT. The effective optical depth to Lyman-$α$ photons from our pipeline is in good agreement with earlier measurements. We add measurements of the mean temperature of the intergalactic gas from $z=3.8 - 2.2$ and use them to constrain the duration and heating rate of helium reionization, finding a preference for an early, hot, helium reionization event, as suggested by measurements from the helium Lyman-$α$ forest. Adding the mean IGM temperature data also increases the significance of the $σ_8$ tension. In the near future we will use our pipeline to infer cosmological parameters from the DESI Lyman-$α$ data.
△ Less
Submitted 7 September, 2023;
originally announced September 2023.
-
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.
-
A Multi-Fidelity Emulator for the Lyman-$α$ Forest Flux Power Spectrum
Authors:
M. A. Fernandez,
Ming-Feng Ho,
Simeon Bird
Abstract:
In this work we extend our recently developed multi-fidelity emulation technique to the simulated Lyman-$α$ forest flux power spectrum. Multi-fidelity 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) simulations, we show that multi-…
▽ More
In this work we extend our recently developed multi-fidelity emulation technique to the simulated Lyman-$α$ forest flux power spectrum. Multi-fidelity 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) simulations, we show that multi-fidelity 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\%$. 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.
△ Less
Submitted 24 October, 2022; v1 submitted 13 July, 2022;
originally announced July 2022.
-
Effect of Separate Initial Conditions on the Lyman-$α$ Forest in Simulations
Authors:
M. A. Fernandez,
Simeon Bird,
Phoebe Upton Sanderbeck
Abstract:
Using a set of high resolution simulations, we quantify the effect of species specific initial transfer functions on probes of the IGM via the Lyman-$α$ forest. We focus on redshifts $2-6$, after H~{\sc i} reionization. We explore the effect of these initial conditions on measures of the thermal state of the low density IGM: the curvature, Doppler width cutoff, and Doppler width distribution. We a…
▽ More
Using a set of high resolution simulations, we quantify the effect of species specific initial transfer functions on probes of the IGM via the Lyman-$α$ forest. We focus on redshifts $2-6$, after H~{\sc i} reionization. We explore the effect of these initial conditions on measures of the thermal state of the low density IGM: the curvature, Doppler width cutoff, and Doppler width distribution. We also examine the matter and flux power spectrum, and potential consequences for constraints on warm dark matter models. We find that the curvature statistic is at most affected at the $\approx2\%$ level at $z=6$. The Doppler width cutoff parameters are affected by $\approx5\%$ for the intercept, and $\approx8\%$ for the fit slope, though this is subdominant to sample variation. The Doppler width distribution shows a $\approx30\%$ effect at $z=3$, however the distribution is not fully converged with simulation box size and resolution. The flux power spectrum is at most affected by $\approx5\%$ at high redshift and small scales. We discuss numerical convergence with simulation parameters.
△ Less
Submitted 12 March, 2021; v1 submitted 18 September, 2020;
originally announced September 2020.
-
Cosmic Filaments from Cosmic Strings
Authors:
M. A. Fernandez,
Simeon Bird,
Yanou Cui
Abstract:
Cosmic strings are generically predicted in many extensions of the Standard Model of particle physics. We propose a new avenue for detecting cosmic strings through their effect on the filamentary structure in the cosmic web. Using cosmological simulations of the density wake from a cosmic string, we examine a variety of filament structure probes. We show that the largest effect of the cosmic strin…
▽ More
Cosmic strings are generically predicted in many extensions of the Standard Model of particle physics. We propose a new avenue for detecting cosmic strings through their effect on the filamentary structure in the cosmic web. Using cosmological simulations of the density wake from a cosmic string, we examine a variety of filament structure probes. We show that the largest effect of the cosmic string is an overdensity in the filament distribution around the string wake. The signal from the overdensity is stronger at higher redshift, and more robust with a wider field. We analyze the spatial distribution of filaments from a publicly available catalog of filaments built from SDSS galaxies. With existing data, we find no evidence for the presence of a cosmic string wake with string tension parameter $Gμ$ above $5\times 10^{-6}$. However, we project WFIRST will be able to detect a signal from such a wake at the $99\%$ confidence level at redshift $z=2$, with significantly higher confidence and the possibility of probing lower tensions ($Gμ\sim 10^{-6}$), at $z=10$. The sensitivity of this method is not competitive with constraints derived from the CMB. However, it provides an independent discovery channel at low redshift, which could be a smoking-gun in scenarios where the CMB bound can be weakened.
△ Less
Submitted 7 August, 2020; v1 submitted 28 April, 2020;
originally announced April 2020.
-
IN-SYNC VI. Identification and Radial Velocity Extraction for 100+ Double-Lined Spectroscopic Binaries in the APOGEE/IN-SYNC Fields
Authors:
M. A. Fernandez,
Kevin R. Covey,
Nathan De Lee,
S. Drew Chojnowski,
David Nidever,
Richard Ballantyne,
Michiel Cottaar,
Nicola Da Rio,
Jonathan B. Foster,
Steven R. Majewski,
Michael R. Meyer,
A. M. Reyna,
G. W. Roberts,
Jacob Skinner,
Keivan Stassun,
Jonathan C. Tan,
Nicholas Troup,
Gail Zasowski
Abstract:
We present radial velocity measurements for 70 high confidence, and 34 potential binary systems in fields containing the Perseus Molecular Cloud, Pleiades, NGC 2264, and the Orion A star forming region. 18 of these systems have been previously identified as binaries in the literature. Candidate double-lined spectroscopic binaries (SB2s) are identified by analyzing the cross-correlation functions (…
▽ More
We present radial velocity measurements for 70 high confidence, and 34 potential binary systems in fields containing the Perseus Molecular Cloud, Pleiades, NGC 2264, and the Orion A star forming region. 18 of these systems have been previously identified as binaries in the literature. Candidate double-lined spectroscopic binaries (SB2s) are identified by analyzing the cross-correlation functions (CCFs) computed during the reduction of each APOGEE spectrum. We identify sources whose CCFs are well fit as the sum of two Lorentzians as likely binaries, and provide an initial characterization of the system based on the radial velocities indicated by that dual fit. For systems observed over several epochs, we present mass ratios and systemic velocities; for two systems with observations on eight or more epochs, and which meet our criteria for robust orbital coverage, we derive initial orbital parameters. The distribution of mass ratios for multi-epoch sources in our sample peaks at q=1, but with a significant tail toward lower q values. Tables reporting radial velocities, systemic velocities, and mass ratios are provided online. We discuss future improvements to the radial velocity extraction method we employ, as well as limitations imposed by the number of epochs currently available in the APOGEE database. The Appendix contains brief notes from the literature on each system in the sample, and more extensive notes for select sources of interest.
△ Less
Submitted 12 June, 2017; v1 submitted 4 June, 2017;
originally announced June 2017.
-
Young Stellar Clusters Containing Massive Young Stellar Objects in the VVV Survey
Authors:
J. Borissova,
S. Ramírez Alegría,
J. Alonso,
P. W. Lucas,
R. Kurtev,
N. Medina,
C. Navarro,
M. Kuhn,
M. Gromadzki,
G. Retamales,
M. A. Fernandez,
C. Agurto-Gangas,
A. -N. Chené,
D. Minniti,
C. Contreras Peña,
M. Catelan,
I. Decany,
M. A. Thompson,
E. F. E. Morales,
P. Amigo
Abstract:
The purpose of this research is to study the connection of global properties of eight young stellar clusters projected in the Vista Variables in the Via Lactea (VVV) ESO Large Public Survey disk area and their young stellar object population. The analysis in based on the combination of spectroscopic parallax-based reddening and distance determinations with main sequence and pre-main sequence ishoc…
▽ More
The purpose of this research is to study the connection of global properties of eight young stellar clusters projected in the Vista Variables in the Via Lactea (VVV) ESO Large Public Survey disk area and their young stellar object population. The analysis in based on the combination of spectroscopic parallax-based reddening and distance determinations with main sequence and pre-main sequence ishochrone fitting to determine the basic parameters (reddening, age, distance) of the sample clusters. The lower mass limit estimations show that all clusters are low or intermediate mass (between 110 and 1800 Mo), the slope Gamma of the obtained present-day mass functions of the clusters is close to the Kroupa initial mass function. On the other hand, the young stellar objects in the surrounding cluster's fields are classified by low resolution spectra, spectral energy distribution fit with theoretical predictions, and variability, taking advantage of multi-epoch VVV observations. All spectroscopically confirmed young stellar objects (except one) are found to be massive (more than 8 Mo). Using VVV and GLIMPSE color-color cuts we have selected a large number of new young stellar object candidates, which are checked for variability and 57% are found to show at least low-amplitude variations. In few cases it was possible to distinguish between YSO and AGB classification on the basis of the light curves.
△ Less
Submitted 26 October, 2016;
originally announced October 2016.