-
The large scale velocity field from the Cosmicflows-4 data
Authors:
Yehuda Hoffman,
Aurélien Valade,
Noam I. Libeskind,
Jenny G. Sorce,
R. Brent Tully,
Simon Pfeifer,
Stefan Gottlöber,
Daniel Pomaréde
Abstract:
The reconstruction of the large scale velocity field from the grouped Cosmicflows-4 (CF4) database is presented. The lognormal bias of the inferred distances and velocities data is corrected by the Bias Gaussianization correction (BGc) scheme, and the linear density and velocity fields are reconstructed by means of the Wiener filter (WF) and constrained realizations (CRs) algorithm. These tools ar…
▽ More
The reconstruction of the large scale velocity field from the grouped Cosmicflows-4 (CF4) database is presented. The lognormal bias of the inferred distances and velocities data is corrected by the Bias Gaussianization correction (BGc) scheme, and the linear density and velocity fields are reconstructed by means of the Wiener filter (WF) and constrained realizations (CRs) algorithm. These tools are tested against a suite of random and constrained Cosmicflows-3-like mock data. The CF4 data consists of 3 main subsamples - the 6dFGS and the SDSS data - and the `others'. The individual contributions of the subsamples have been studied. The quantitative analysis of the velocity field is done mostly by the mean overdensity ($Δ_L(R)$) and the bulk velocity ($V_{\mathrm{bulk}}(R)$) profiles of the velocity field out to $300\, h^{-1}{\rm Mpc}$.
The $V_{\mathrm{bulk}}(R)$ and $Δ_{\mathrm L}(R)$ profiles of the CF4 data without its 6dFGS component are consistent with the cosmic variance to within $1σ$. The 6dFGS sample dominates the $V_{\mathrm{bulk}}$ ($Δ_{\mathrm L}$) profile beyond $\sim120\, h^{-1}{\rm Mpc}$, and drives it to roughly a $3.4σ$ ($-1.9σ$) excess (deficiency) relative to the cosmic variance at $R\sim250\ (190)\ \, h^{-1}{\rm Mpc}$. The excess in the amplitude of $V_{\mathrm{bulk}}$ is dominated by its Supergalactic X component, roughly in the direction of the Shapley Concentration. The amplitude and alignment of the inferred velocity field from the CF4 data is at $\sim(2\,-\,3)\,σ$ discrepancy with respect to the $Λ$CDM model. Namely, it is somewhat atypical but yet there is no compelling tension with the model.
△ Less
Submitted 2 November, 2023;
originally announced November 2023.
-
A Local Universe model for constrained simulations
Authors:
Simon Pfeifer,
Aurélien Valade,
Stefan Gottlöber,
Yehuda Hoffman,
Noam I. Libeskind,
Wojciech A. Hellwing
Abstract:
The aim of cosmological simulations is to reproduce the properties of the observed Universe, serving as tools to test structure and galaxy formation models. Constrained simulations of our local cosmological region up to a few hundred Mpc/h , the local Universe, are designed to reproduce the actual cosmic web of structures as observed. A question that often arises is how to judge the quality of con…
▽ More
The aim of cosmological simulations is to reproduce the properties of the observed Universe, serving as tools to test structure and galaxy formation models. Constrained simulations of our local cosmological region up to a few hundred Mpc/h , the local Universe, are designed to reproduce the actual cosmic web of structures as observed. A question that often arises is how to judge the quality of constrained simulations against the observations of the Local Universe. Here we introduce the Local Universe model (LUM), a new methodology, whereby many constrained simulations can be judged and the ''best'' initial conditions can be identified. By characterising the Local Universe as a set of rich clusters, the model identifies haloes that serve as simulated counterparts to the observed clusters. Their merit is determined against a null hypothesis, the probability that such a counterpart could be identified in a random, unconstrained simulation. This model is applied to 100 constrained simulations using the Cosmicflows-3 data. Cluster counterparts are found for all constrained simulations, their distribution of separation from the true observed cluster position and their mass distribution are investigated. Lastly, the ''best'' constrained simulation is selected using the LUM and discussed in more detail.
△ Less
Submitted 9 May, 2023;
originally announced May 2023.
-
Testing Bayesian reconstruction methods from peculiar velocities
Authors:
Aurélien Valade,
Noam I Libeskind,
Yehuda Hoffman,
Simon Pfeifer
Abstract:
Reconstructing the large scale density and velocity fields from surveys of galaxy distances, is a major challenge for cosmography. The data is very noisy and sparse. Estimated distances, and thereby peculiar velocities, are strongly affected by the Malmquist-like lognormal bias. Two algorithms have been recently introduced to perform reconstructions from such data: the Bias Gaussian correction cou…
▽ More
Reconstructing the large scale density and velocity fields from surveys of galaxy distances, is a major challenge for cosmography. The data is very noisy and sparse. Estimated distances, and thereby peculiar velocities, are strongly affected by the Malmquist-like lognormal bias. Two algorithms have been recently introduced to perform reconstructions from such data: the Bias Gaussian correction coupled with the Wiener filter (BGc/WF) and the HAMLET implementation of the Hamiltonian Monte Carlo forward modelling. The two methods are tested here against mock catalogs that mimic the Cosmicflows-3 data. Specifically the reconstructed cosmography and moments of the velocity field (monopole, dipole) are examined. A comparison is made to the ``exact'' wiener filter as well - namely the Wiener Filter in the unrealistic case of zero observational errors. This is to understand the limits of the WF method. The following is found. In the nearby regime ($d \lesssim 40 {\rm Mpc}/h$) the two methods perform roughly equally well. HAMLET does slightly better in the intermediate regime ($ 40 \lesssim d \lesssim 120 {\rm Mpc}/h$). The main differences between the two appear in the most distant regime ($d \gtrsim 120 {\rm Mpc}/h$), close to the edge of the data. The HAMLET outperforms the BGc/WF in terms of better and tighter correlations, yet in the distant regime the HAMLET yields a somewhat biased reconstruction. Such biases are missing from the BGc/WF reconstruction. In sum, both methods perform well and create reliable reconstructions with significant differences apparent when details are examined.
△ Less
Submitted 13 September, 2022;
originally announced September 2022.
-
Hamiltonian Monte Carlo reconstruction from peculiar velocities
Authors:
Aurélien Valade,
Yehuda Hoffman,
Noam I Libeskind,
Romain Graziani
Abstract:
The problem of the reconstruction of the large scale density and velocity fields from peculiar velocities surveys is addressed here within a Bayesian framework by means of Hamiltonian Monte Carlo (HMC) sampling. The HAmiltonian Monte carlo reconstruction of the Local EnvironmenT (Hamlet) algorithm is designed to reconstruct the linear large scale density and velocity fields in conjunction with the…
▽ More
The problem of the reconstruction of the large scale density and velocity fields from peculiar velocities surveys is addressed here within a Bayesian framework by means of Hamiltonian Monte Carlo (HMC) sampling. The HAmiltonian Monte carlo reconstruction of the Local EnvironmenT (Hamlet) algorithm is designed to reconstruct the linear large scale density and velocity fields in conjunction with the undoing of lognormal bias in the derived distances and velocities of peculiar velocities surveys such as the Cosmicflows data. The Hamlet code has been tested against Cosmicflows mock catalogs consisting of up to 30 000 data points with mock errors akin to those of the Cosmicflows-3 data, within the framework of the LCDM standard model of cosmology.
The Hamlet code outperforms previous applications of Gibbs sampling MCMC reconstruction from the Cosmicflows-3 data by two to four orders of magnitude in CPU time. The gain in performance is due to the inherent higher efficiency of the HMC algorithm and due to parallel computing on GPUs rather than CPUs. This gain will enable an increase in the reconstruction of the large scale structure from the upcoming Cosmicfows-4 data and the setting of constrained initial conditions for cosmological high resolution simulations.
△ Less
Submitted 4 March, 2022;
originally announced March 2022.
-
From Cosmicflows distance moduli to unbiased distances and peculiar velocities
Authors:
Yehuda Hoffman,
Adi Nusser,
Aurelien Valade,
Noam I. Libeskind,
R. Brent Tully
Abstract:
Surveys of galaxy distances and radial peculiar velocities can be used to reconstruct the large scale structure. Other than systematic errors in the zero-point calibration of the galaxy distances the main source of uncertainties of such data are errors on the distance moduli, assumed here to be Gaussian and thus turn into lognormal errors on distances and velocities. Naively treated, it leads to s…
▽ More
Surveys of galaxy distances and radial peculiar velocities can be used to reconstruct the large scale structure. Other than systematic errors in the zero-point calibration of the galaxy distances the main source of uncertainties of such data are errors on the distance moduli, assumed here to be Gaussian and thus turn into lognormal errors on distances and velocities. Naively treated, it leads to spurious nearby outflow and strong infall at larger distances. The lognormal bias is corrected here and tested against mock data extracted from a $Λ$CDM simulation, designed to statistically follow the grouped Cosmicflows-3 (CF3) data. Considering a subsample of data points, all of which have the same true distances or same redshifts, the lognormal bias arises because the means of the distributions of observed distances and velocities are skewed off the means of the true distances and velocities. Yet, the medians are invariant under the lognormal transformation. That invariance allows the Gaussianization of the distances and velocities and the removal of the lognormal bias. This Bias Gaussianization correction (BGc) algorithm is tested against mock CF3 catalogs. The test consists of a comparison of the BGC estimated with the simulated distances and velocities and of an examination of the Wiener filter reconstruction from the BGc data. Indeed, the BGc eliminates the lognormal bias. The estimation of Hubble's ($H_{0}$) constant is also tested. The residual of the BGc estimated $H_{0}$ from the simulated values is $0.6 \pm 0.7 {\rm kms}^{-1}{\rm Mpc}^{-1}$ and is dominated by the cosmic variance. The BGc correction of the actual CF3 data yields $H_{0} = 75.8 \pm 1.1 {\rm kms}^{-1}{\rm Mpc}^{-1}$ .
△ Less
Submitted 19 May, 2021;
originally announced May 2021.
-
Waves in the radiative zones of rotating, magnetized stars
Authors:
Aurélien Valade,
Vincent Prat,
Stéphane Mathis,
Kyle Augustson
Abstract:
Asteroseismology has reached a level of accuracy that may allow us to detect the effect of a deep magnetic field on oscillation modes. We thus aim to develop an asymptotic theory for short-wavelength waves in radiative zones of rotating stars with a general magnetic field topology (toroidal and poloidal). A dispersion relation is derived for a uniformly rotating star with a low-amplitude, axisymme…
▽ More
Asteroseismology has reached a level of accuracy that may allow us to detect the effect of a deep magnetic field on oscillation modes. We thus aim to develop an asymptotic theory for short-wavelength waves in radiative zones of rotating stars with a general magnetic field topology (toroidal and poloidal). A dispersion relation is derived for a uniformly rotating star with a low-amplitude, axisymmetric magnetic field. The parameter space of this model is explored with a Hamiltonian ray-tracing method. The features of the gravito-inertial waves modified by the magnetic field are studied with different magnetic topologies.
△ Less
Submitted 12 October, 2018; v1 submitted 4 October, 2018;
originally announced October 2018.
