-
Sesame: A power spectrum emulator pipeline for beyond-$Λ$CDM models
Authors:
Renate Mauland,
Hans A. Winther,
Cheng-Zong Ruan
Abstract:
The mysterious nature of the dark sector of the $Λ$CDM model is one of the main motivators behind the study of alternative cosmological models. A central quantity of interest for these models is the matter power spectrum, which quantifies structure formation on various scales and can be cross-validated through theory, simulations, and observations. We present a tool that can be used to create emul…
▽ More
The mysterious nature of the dark sector of the $Λ$CDM model is one of the main motivators behind the study of alternative cosmological models. A central quantity of interest for these models is the matter power spectrum, which quantifies structure formation on various scales and can be cross-validated through theory, simulations, and observations. We present a tool that can be used to create emulators for the non-linear matter power spectrum, and similar global clustering statistics, for models beyond $Λ$CDM with very little computation effort and without the need for supercomputers. We use fast approximate N-body simulations to emulate the boost, $B(k,z)=P_{\text{beyond}-Λ\rm CDM}(k,z) / P_{Λ\rm CDM}(k, z)$, and then rely on existing high-quality emulators made for $Λ$CDM to isolate $P_{\text{beyond}-Λ\rm CDM}(k,z)$. Since both the $Λ$CDM and beyond-$Λ$CDM models are simulated in our approach, some of the lack of power on small scales is factored out, allowing us to extend the emulator to $k \sim 3-5\,h\,\mathrm{Mpc}^{-1}$ and still maintain good accuracy. In addition, errors from the simulation and emulation process can be estimated and factored into the covariance when using the emulator on data. As an example of using the pipeline, we create an emulator for the well-studied $f(R)$ model with massive neutrinos, using approximately 3000 CPU hours of computation time. Provided with the paper is a fully functioning pipeline that generates parameter samples, runs a Boltzmann solver to produce initial conditions, runs the simulations, and then gathers all the data and runs it through a machine learning module to develop the emulator. This tool, named Sesame, can be used by anyone to generate a power spectrum emulator for the cosmological model of their choice.
△ Less
Submitted 10 June, 2024; v1 submitted 23 September, 2023;
originally announced September 2023.
-
The void-galaxy cross-correlation function with massive neutrinos and modified gravity
Authors:
Renate Mauland,
Øystein Elgarøy,
David Fonseca Mota,
Hans Arnold Winther
Abstract:
Massive neutrinos and $f(R)$ modified gravity have degenerate observational signatures that can impact the interpretation of results in galaxy survey experiments, such as cosmological parameter estimations and gravity model tests. Because of this, it is important to investigate astrophysical observables that can break these degeneracies. Cosmic voids are sensitive to both massive neutrinos and mod…
▽ More
Massive neutrinos and $f(R)$ modified gravity have degenerate observational signatures that can impact the interpretation of results in galaxy survey experiments, such as cosmological parameter estimations and gravity model tests. Because of this, it is important to investigate astrophysical observables that can break these degeneracies. Cosmic voids are sensitive to both massive neutrinos and modifications of gravity and provide a promising ground for disentangling the above mentioned degeneracies. In order to analyse cosmic voids in the context of non-$Λ$CDM cosmologies, we must first understand how well the current theoretical framework operates in these settings. We performed a suite of simulations with the RAMSES-based N-body code ANUBISIS, including massive neutrinos and $f(R)$ modified gravity both individually and simultaneously. The data from the simulations were compared to models of the void velocity profile and the void-halo cross-correlation function (CCF). This was done both with the real space simulation data as model input and by applying a reconstruction method to the redshift space data. In addition, we ran Markov chain Monte Carlo (MCMC) fits on the data sets to assess the capability of the models to reproduce the fiducial simulation values of $fσ_8(z)$ and the Alcock-Paczynski parameter, $ε$. The void modelling applied performs similarly for all simulated cosmologies, indicating that more accurate models and higher resolution simulations are needed in order to directly observe the effects of massive neutrinos and $f(R)$ modified gravity through studies of the void-galaxy CCF. The MCMC fits show that the choice of void definition plays an important role in the recovery of the correct cosmological parameters, but otherwise, there is no clear distinction between the ability to reproduce $fσ_8$ and $ε$ for the various simulations.
△ Less
Submitted 21 June, 2023; v1 submitted 10 March, 2023;
originally announced March 2023.
-
Euclid: Modelling massive neutrinos in cosmology -- a code comparison
Authors:
J. Adamek,
R. E. Angulo,
C. Arnold,
M. Baldi,
M. Biagetti,
B. Bose,
C. Carbone,
T. Castro,
J. Dakin,
K. Dolag,
W. Elbers,
C. Fidler,
C. Giocoli,
S. Hannestad,
F. Hassani,
C. Hernández-Aguayo,
K. Koyama,
B. Li,
R. Mauland,
P. Monaco,
C. Moretti,
D. F. Mota,
C. Partmann,
G. Parimbelli,
D. Potter
, et al. (111 additional authors not shown)
Abstract:
The measurement of the absolute neutrino mass scale from cosmological large-scale clustering data is one of the key science goals of the Euclid mission. Such a measurement relies on precise modelling of the impact of neutrinos on structure formation, which can be studied with $N$-body simulations. Here we present the results from a major code comparison effort to establish the maturity and reliabi…
▽ More
The measurement of the absolute neutrino mass scale from cosmological large-scale clustering data is one of the key science goals of the Euclid mission. Such a measurement relies on precise modelling of the impact of neutrinos on structure formation, which can be studied with $N$-body simulations. Here we present the results from a major code comparison effort to establish the maturity and reliability of numerical methods for treating massive neutrinos. The comparison includes eleven full $N$-body implementations (not all of them independent), two $N$-body schemes with approximate time integration, and four additional codes that directly predict or emulate the matter power spectrum. Using a common set of initial data we quantify the relative agreement on the nonlinear power spectrum of cold dark matter and baryons and, for the $N$-body codes, also the relative agreement on the bispectrum, halo mass function, and halo bias. We find that the different numerical implementations produce fully consistent results. We can therefore be confident that we can model the impact of massive neutrinos at the sub-percent level in the most common summary statistics. We also provide a code validation pipeline for future reference.
△ Less
Submitted 8 August, 2023; v1 submitted 22 November, 2022;
originally announced November 2022.
-
Quantized Vortices in Superfluid Dark Matter
Authors:
Renate Mauland,
Øystein Elgarøy
Abstract:
In 2015 Berezhiani & Khoury proposed a Superfluid Dark Matter (SFDM) model where dark matter condenses and forms a superfluid on galactic scales. In the superfluid state phonons interact with baryons, resulting in a behavior similar to that of Modified Newtonian Dynamics (MOND). If one assumes that the DM condensate rotates along with the galaxy, a grid of vortices should form throughout the super…
▽ More
In 2015 Berezhiani & Khoury proposed a Superfluid Dark Matter (SFDM) model where dark matter condenses and forms a superfluid on galactic scales. In the superfluid state phonons interact with baryons, resulting in a behavior similar to that of Modified Newtonian Dynamics (MOND). If one assumes that the DM condensate rotates along with the galaxy, a grid of vortices should form throughout the superfluid component if the rotation is fast enough. We aim to investigate the size and impact of the vortices on surrounding baryons, and to further investigate the parameter space of the model. We also look for a possible vortex solution of the Lagrangian presented for the SFDM theory. We first take a simple approach and investigate vortex properties in a constant density DM halo, applying knowledge from condensed matter physics. We then use the zero-temperature condensate density profile as a template to vary the DM particle mass and the energy scale, $Λ$, of the SFDM model. Further, we attempt to find a vortex solution of the theory by extracting the Euler-Lagrange equation with respect to the modulus of the condensate wavefunction from the full relativistic SFDM Lagrangian. For the constant density approach we find that the vortices are on millimeter scale, and separated by distances $\sim0.002\,\rm{AU}$. The parameter space of the model is found to be substantial and a reduction in the DM particle mass leads to larger vortices with a higher energy. However, none of the parameter combinations explored here give both realistic values of $Λ$ and vortices energetic enough to have an observational impact on the galaxy as a whole. The vortex equation extracted from the Lagrangian of the model is unstable, and no solution exhibiting the standard properties of a vortex solution is found.
△ Less
Submitted 31 January, 2022; v1 submitted 2 September, 2021;
originally announced September 2021.
-
The Ubiquitous Imprint of Radiative Acceleration in the Mean Absorption Spectrum of Quasar Outflows
Authors:
Lluís Mas-Ribas,
Renate Mauland
Abstract:
Observational evidence revealing the main mechanisms that accelerate quasar outflows has proven difficult to obtain, due to the complexity of the absorption features that this gas produces in the spectra of the emission sources. We build 36 composite outflow spectra, covering a large range of outflow and quasar parameters, by stacking broad ($>450\,{\rm km\,s^{-1}}$) absorption line systems in the…
▽ More
Observational evidence revealing the main mechanisms that accelerate quasar outflows has proven difficult to obtain, due to the complexity of the absorption features that this gas produces in the spectra of the emission sources. We build 36 composite outflow spectra, covering a large range of outflow and quasar parameters, by stacking broad ($>450\,{\rm km\,s^{-1}}$) absorption line systems in the spectra of SDSS-III/BOSS DR12 quasars. The two lines of the atomic doublet of CIV, with a separation of $\approx 497\,{\rm km\,s^{-1}}$, as well as those of other species appear well resolved in most of our composites. This agrees with broad outflow troughs consisting in the superposition of narrow absorbers. We also report on the ubiquitous detection of the radiative-acceleration signature known as line locking in all our composite outflow spectra, including one spectrum strictly built from broad absorption line (BAL) systems. This is the first line-locking detection in BAL composite spectra. Line locking is driven by the CIV atomic doublet, and is visible on the blue side of most strong absorption transitions. Similar effects from the doublets of OVI, SiIV, or NV, however, seem to not be present. Our results confirm that radiation pressure is a prevalent mechanism for accelerating outflows in quasars.
△ Less
Submitted 17 October, 2019; v1 submitted 11 February, 2019;
originally announced February 2019.