Showing 1–50 of 101 results for author: Lavaux, G

PineTree: A generative, fast, and differentiable halo model for wide-field galaxy surveys

Authors: Simon Ding, Guilhem Lavaux, Jens Jasche

Abstract: Mock halo catalogues are indispensable data products for developing and validating cosmological inference pipelines. A major challenge in generating mock catalogues is modelling the halo or galaxy bias, which is the mapping from matter density to dark matter halos or observable galaxies. To this end, N-body codes produce state-of-the-art catalogues. However, generating large numbers of these N-bod… ▽ More Mock halo catalogues are indispensable data products for developing and validating cosmological inference pipelines. A major challenge in generating mock catalogues is modelling the halo or galaxy bias, which is the mapping from matter density to dark matter halos or observable galaxies. To this end, N-body codes produce state-of-the-art catalogues. However, generating large numbers of these N-body simulations for big volumes, requires significant computational time. We introduce and benchmark a differentiable and physics-informed neural network that can generate mock halo catalogues of comparable quality to those obtained from full N-body codes. The model design is computationally efficient for the training procedure and the production of large mock suites. We present a neural network, relying only on 18 to 34 trainable parameters, that produces halo catalogues from dark matter overdensity fields. The reduction of network weights is realised through incorporating symmetries motivated by first principles into our model architecture. We train our model using dark matter only N-body simulations across different resolutions, redshifts, and mass bins. We validate the final mock catalogues by comparing them to N-body halo catalogues using different N-point correlation functions. Our model produces mock halo catalogues consistent with the reference simulations, showing that this novel network is a promising way to generate mock data for upcoming wide-field surveys due to its computational efficiency. Moreover, we find that the network can be trained on approximate overdensity fields to reduce the computational cost further. We also present how the trained network parameters can be interpreted to give insights into the physics of structure formation. Finally, we discuss the current limitations of our model as well as more general requirements and pitfalls for approximate halo mock generation. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 17 pages, 17 figures, 3 tables. Abstract abridged for arXiv

Constrained cosmological simulations of the Local Group using Bayesian hierarchical field-level inference

Authors: Ewoud Wempe, Guilhem Lavaux, Simon D. M. White, Amina Helmi, Jens Jasche, Stephen Stopyra

Abstract: We present a novel approach based on Bayesian field-level inference capable of resolving individual galaxies within the Local Group (LG), enabling detailed studies of its structure and formation via posterior simulations. We extend the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm with a multi-resolution approach, allowing us to reach smaller mass scales and apply observational con… ▽ More We present a novel approach based on Bayesian field-level inference capable of resolving individual galaxies within the Local Group (LG), enabling detailed studies of its structure and formation via posterior simulations. We extend the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm with a multi-resolution approach, allowing us to reach smaller mass scales and apply observational constraints based on LG galaxies. Our updated data model simultaneously accounts for observations of mass tracers within the dark haloes of the Milky Way (MW) and M31, their observed separation and relative velocity, and the quiet surrounding Hubble flow represented through the positions and velocities of galaxies at distances from one to four Mpc. Our approach delivers representative posterior samples of $Λ$CDM realisations that are statistically and simultaneously consistent with all these observations, leading to significantly tighter mass constraints than found if the individual datasets are considered separately. In particular, we estimate the virial masses of the MW and M31 to be $\log_{10}(M_{200c}/M_\odot) = 12.07\pm0.08$ and $12.33\pm0.10$, respectively, their sum to be $\log_{10}(ΣM_{200c}/M_\odot)= 12.52\pm0.07$, and the enclosed mass within spheres of radius $R$ to be $\log_{10}(M(R)/M_\odot)= 12.71\pm0.06$ and $12.96\pm0.08$ for $R=1$ Mpc and 3 Mpc, respectively. The M31-MW orbit is nearly radial for most of our $Λ$CDM LG's, and most lie in a dark matter sheet that aligns approximately with the Supergalactic Plane, even though the surrounding density field was not used explicitly as a constraint. The approximate simulations employed in our inference are accurately reproduced by high-fidelity structure formation simulations, demonstrating the potential for future high-resolution, full-physics $Λ$CDM posterior simulations of LG look-alikes. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 23 pages, 17 figures, 4 tables, submitted to A&A. Comments welcome!

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Bayesian Inference of Initial Conditions from Non-Linear Cosmic Structures using Field-Level Emulators

Authors: Ludvig Doeser, Drew Jamieson, Stephen Stopyra, Guilhem Lavaux, Florent Leclercq, Jens Jasche

Abstract: Analysing next-generation cosmological data requires balancing accurate modeling of non-linear gravitational structure formation and computational demands. We propose a solution by introducing a machine learning-based field-level emulator, within the Hamiltonian Monte Carlo-based Bayesian Origin Reconstruction from Galaxies (BORG) inference algorithm. Built on a V-net neural network architecture,… ▽ More Analysing next-generation cosmological data requires balancing accurate modeling of non-linear gravitational structure formation and computational demands. We propose a solution by introducing a machine learning-based field-level emulator, within the Hamiltonian Monte Carlo-based Bayesian Origin Reconstruction from Galaxies (BORG) inference algorithm. Built on a V-net neural network architecture, the emulator enhances the predictions by first-order Lagrangian perturbation theory to be accurately aligned with full $N$-body simulations while significantly reducing evaluation time. We test its incorporation in BORG for sampling cosmic initial conditions using mock data based on non-linear large-scale structures from $N$-body simulations and Gaussian noise. The method efficiently and accurately explores the high-dimensional parameter space of initial conditions, fully extracting the cross-correlation information of the data field binned at a resolution of $1.95h^{-1}$ Mpc. Percent-level agreement with the ground truth in the power spectrum and bispectrum is achieved up to the Nyquist frequency $k_\mathrm{N} \approx 2.79h \; \mathrm{Mpc}^{-1}$. Posterior resimulations - using the inferred initial conditions for $N$-body simulations - show that the recovery of information in the initial conditions is sufficient to accurately reproduce halo properties. In particular, we show highly accurate $M_{200\mathrm{c}}$ halo mass function and stacked density profiles of haloes in different mass bins $[0.853,16]\times 10^{14}M_{\odot}h^{-1}$. As all available cross-correlation information is extracted, we acknowledge that limitations in recovering the initial conditions stem from the noise level and data grid resolution. This is promising as it underscores the significance of accurate non-linear modeling, indicating the potential for extracting additional information at smaller scales. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: 16 pages, 11 figures

An Anti-halo Void Catalogue of the Local Super-Volume

Authors: Stephen Stopyra, Hiranya V. Peiris, Andrew Pontzen, Jens Jasche, Guilhem Lavaux

Abstract: We construct an anti-halo void catalogue of $150$ voids with radii $R > 10\,h^{-1}\mathrm{\,Mpc}$ in the Local Super-Volume ($<135\,h^{-1}\mathrm{\,Mpc}$ from the Milky Way), using posterior resimulation of initial conditions inferred by field-level inference with Bayesian Origin Reconstruction from Galaxies (\codefont{BORG}). We describe and make use of a new algorithm for creating a single, unif… ▽ More We construct an anti-halo void catalogue of $150$ voids with radii $R > 10\,h^{-1}\mathrm{\,Mpc}$ in the Local Super-Volume ($<135\,h^{-1}\mathrm{\,Mpc}$ from the Milky Way), using posterior resimulation of initial conditions inferred by field-level inference with Bayesian Origin Reconstruction from Galaxies (\codefont{BORG}). We describe and make use of a new algorithm for creating a single, unified void catalogue by combining different samples from the posterior. The catalogue is complete out to $135\,h^{-1}\mathrm{\,Mpc}$, with void abundances matching theoretical predictions. Finally, we compute stacked density profiles of those voids which are reliably identified across posterior samples, and show that these are compatible with $Λ$CDM expectations once environmental selection (e.g., the estimated $\sim 4\%$ under-density of the Local Super-Volume) is accounted for. △ Less

Submitted 16 May, 2024; v1 submitted 21 November, 2023; originally announced November 2023.

Comments: 11 pages, 4 figures. Updated to match version accepted by MNRAS. Minor clarifications and no changes to conclusion or catalogue

Distinct distributions of elliptical and disk galaxies across the Local Supercluster as a $Λ$CDM prediction

Authors: Till Sawala, Carlos Frenk, Jens Jasche, Peter H. Johansson, Guilhem Lavaux

Abstract: Galaxies of different types are not equally distributed in the Local Universe. In particular, the supergalactic plane is prominent among the brightest ellipticals, but inconspicuous among the brightest disk galaxies. This striking difference provides a unique test for our understanding of galaxy and structure formation. Here we use the SIBELIUS DARK constrained simulation to confront the predictio… ▽ More Galaxies of different types are not equally distributed in the Local Universe. In particular, the supergalactic plane is prominent among the brightest ellipticals, but inconspicuous among the brightest disk galaxies. This striking difference provides a unique test for our understanding of galaxy and structure formation. Here we use the SIBELIUS DARK constrained simulation to confront the predictions of the standard Lambda Cold Dark Matter ($Λ$CDM) model and standard galaxy formation theory with these observations. We find that SIBELIUS DARK reproduces the spatial distributions of disks and ellipticals and, in particular, the observed excess of massive ellipticals near the supergalactic equator. We show that this follows directly from the local large-scale structure and from the standard galaxy formation paradigm, wherein disk galaxies evolve mostly in isolation, while giant ellipticals congregate in the massive clusters that define the supergalactic plane. Rather than being anomalous as earlier works have suggested, the distributions of giant ellipticals and disks in the Local Universe and in relation to the supergalactic plane are key predictions of the $Λ$CDM model. △ Less

Submitted 5 December, 2023; v1 submitted 1 November, 2023; originally announced November 2023.

Comments: 21 pages, 9 figures, published in Nature Astronomy

Euclid: The search for primordial features

Authors: M. Ballardini, Y. Akrami, F. Finelli, D. Karagiannis, B. Li, Y. Li, Z. Sakr, D. Sapone, A. Achúcarro, M. Baldi, N. Bartolo, G. Cañas-Herrera, S. Casas, R. Murgia, H. A. Winther, M. Viel, A. Andrews, J. Jasche, G. Lavaux, D. K. Hazra, D. Paoletti, J. Valiviita, A. Amara, S. Andreon, N. Auricchio , et al. (104 additional authors not shown)

Abstract: Primordial features, in particular oscillatory signals, imprinted in the primordial power spectrum of density perturbations represent a clear window of opportunity for detecting new physics at high-energy scales. Future spectroscopic and photometric measurements from the $Euclid$ space mission will provide unique constraints on the primordial power spectrum, thanks to the redshift coverage and hig… ▽ More Primordial features, in particular oscillatory signals, imprinted in the primordial power spectrum of density perturbations represent a clear window of opportunity for detecting new physics at high-energy scales. Future spectroscopic and photometric measurements from the $Euclid$ space mission will provide unique constraints on the primordial power spectrum, thanks to the redshift coverage and high-accuracy measurement of nonlinear scales, thus allowing us to investigate deviations from the standard power-law primordial power spectrum. We consider two models with primordial undamped oscillations superimposed on the matter power spectrum, one linearly spaced in $k$-space the other logarithmically spaced in $k$-space. We forecast uncertainties applying a Fisher matrix method to spectroscopic galaxy clustering, weak lensing, photometric galaxy clustering, cross correlation between photometric probes, spectroscopic galaxy clustering bispectrum, CMB temperature and $E$-mode polarization, temperature-polarization cross correlation, and CMB weak lensing. We also study a nonlinear density reconstruction method to retrieve the oscillatory signals in the primordial power spectrum. We find the following percentage relative errors in the feature amplitude with $Euclid$ primary probes for the linear (logarithmic) feature model: 21% (22%) in the pessimistic settings and 18% (18%) in the optimistic settings at 68.3% confidence level (CL) using GC$_{\rm sp}$+WL+GC$_{\rm ph}$+XC. Combining all the sources of information explored expected from $Euclid$ in combination with future SO-like CMB experiment, we forecast ${\cal A}_{\rm lin} \simeq 0.010 \pm 0.001$ at 68.3% CL and ${\cal A}_{\rm log} \simeq 0.010 \pm 0.001$ for GC$_{\rm sp}$(PS rec + BS)+WL+GC$_{\rm ph}$+XC+SO-like both for the optimistic and pessimistic settings over the frequency range $(1,\,10^{2.1})$. △ Less

Submitted 29 March, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 23 pages, 9 figures, 4 tables

Journal ref: A&A 683, A220 (2024)

Where shadows lie: reconstruction of anisotropies in the neutrino sky

Authors: Willem Elbers, Carlos S. Frenk, Adrian Jenkins, Baojiu Li, Silvia Pascoli, Jens Jasche, Guilhem Lavaux, Volker Springel

Abstract: The Cosmic Neutrino Background (CNB) encodes a wealth of information, but has not yet been observed directly. To determine the prospects of detection and to study its information content, we reconstruct the phase-space distribution of local relic neutrinos from the three-dimensional distribution of matter within 200 Mpc/h of the Milky Way. Our analysis relies on constrained realization simulations… ▽ More The Cosmic Neutrino Background (CNB) encodes a wealth of information, but has not yet been observed directly. To determine the prospects of detection and to study its information content, we reconstruct the phase-space distribution of local relic neutrinos from the three-dimensional distribution of matter within 200 Mpc/h of the Milky Way. Our analysis relies on constrained realization simulations and forward modelling of the 2M++ galaxy catalogue. We find that the angular distribution of neutrinos is anti-correlated with the projected matter density, due to the capture and deflection of neutrinos by massive structures along the line of sight. Of relevance to tritium capture experiments, we find that the gravitational clustering effect of the large-scale structure on the local number density of neutrinos is more important than that of the Milky Way for neutrino masses less than 0.1 eV. Nevertheless, we predict that the density of relic neutrinos is close to the cosmic average, with a suppression or enhancement over the mean of (-0.3%, +7%, +27%) for masses of (0.01, 0.05, 0.1) eV. This implies no more than a marginal increase in the event rate for tritium capture experiments like PTOLEMY. We also predict that the CNB and CMB rest frames coincide for 0.01 eV neutrinos, but that neutrino velocities are significantly perturbed for masses larger than 0.05 eV. Regardless of mass, we find that the angle between the neutrino dipole and the ecliptic plane is small, implying a near-maximal annual modulation in the bulk velocity. Along with this paper, we publicly release our simulation data, comprising more than 100 simulations for six different neutrino masses. △ Less

Submitted 28 December, 2023; v1 submitted 6 July, 2023; originally announced July 2023.

Comments: 27 pages, 8 figures, JCAP, accepted with minor revisions

Towards Accurate Field-Level Inference of Massive Cosmic Structures

Authors: Stephen Stopyra, Hiranya V. Peiris, Andrew Pontzen, Jens Jasche, Guilhem Lavaux

Abstract: We investigate the accuracy requirements for field-level inference of cluster and void masses using data from galaxy surveys. We introduce a two-step framework that takes advantage of the fact that cluster masses are determined by flows on larger scales than the clusters themselves. First, we determine the integration accuracy required to perform field-level inference of cosmic initial conditions… ▽ More We investigate the accuracy requirements for field-level inference of cluster and void masses using data from galaxy surveys. We introduce a two-step framework that takes advantage of the fact that cluster masses are determined by flows on larger scales than the clusters themselves. First, we determine the integration accuracy required to perform field-level inference of cosmic initial conditions on these large scales, by fitting to late-time galaxy counts using the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm. A 20-step COLA integrator is able to accurately describe the density field surrounding the most massive clusters in the Local Super-Volume ($<135\,h^{-1}\,\mathrm{Mpc}$), but does not by itself lead to converged virial mass estimates. Therefore we carry out `posterior resimulations', using full $N$-body dynamics while sampling from the inferred initial conditions, and thereby obtain estimates of masses for nearby massive clusters. We show that these are in broad agreement with existing estimates, and find that mass functions in the Local Super-Volume are compatible with $Λ$CDM. △ Less

Submitted 8 January, 2024; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 14 pages, 7 figures. Updated to match version accepted by MNRAS (published 2023/10/18). Minor changes and 1 new figure

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 1, January 2024, Pages 1244-1256

Field-level inference of cosmic shear with intrinsic alignments and baryons

Authors: Natalia Porqueres, Alan Heavens, Daniel Mortlock, Guilhem Lavaux, T. Lucas Makinen

Abstract: We construct a field-based Bayesian Hierarchical Model for cosmic shear that includes, for the first time, the important astrophysical systematics of intrinsic alignments and baryon feedback, in addition to a gravity model. We add to the BORG-WL framework the tidal alignment and tidal torquing model (TATT) for intrinsic alignments and compare them with the non-linear alignment (NLA) model. With sy… ▽ More We construct a field-based Bayesian Hierarchical Model for cosmic shear that includes, for the first time, the important astrophysical systematics of intrinsic alignments and baryon feedback, in addition to a gravity model. We add to the BORG-WL framework the tidal alignment and tidal torquing model (TATT) for intrinsic alignments and compare them with the non-linear alignment (NLA) model. With synthetic data, we have shown that adding intrinsic alignments and sampling the TATT parameters does not reduce the constraining power of the method and the field-based approach lifts the weak lensing degeneracy. We add baryon effects at the field level using the enthalpy gradient descent (EGD) model. This model displaces the dark matter particles without knowing whether they belong to a halo and allows for self-calibration of the model parameters, which are inferred from the data. We have also illustrated the effects of model misspecification for the baryons. The resulting model now contains the most important physical effects and is suitable for application to data. △ Less

Submitted 10 April, 2023; originally announced April 2023.

LyAl-Net: A high-efficiency Lyman-$α$ forest simulation with a neural network

Authors: Chotipan Boonkongkird, Guilhem Lavaux, Sebastien Peirani, Yohan Dubois, Natalia Porqueres, Eleni Tsaprazi

Abstract: The inference of cosmological quantities requires accurate and large hydrodynamical cosmological simulations. Unfortunately, their computational time can take millions of CPU hours for a modest coverage in cosmological scales ($\approx (100 {h^{-1}}\,\text{Mpc})^3)$). The possibility to generate large quantities of mock Lyman-$α$ observations opens up the possibility of much better control on cova… ▽ More The inference of cosmological quantities requires accurate and large hydrodynamical cosmological simulations. Unfortunately, their computational time can take millions of CPU hours for a modest coverage in cosmological scales ($\approx (100 {h^{-1}}\,\text{Mpc})^3)$). The possibility to generate large quantities of mock Lyman-$α$ observations opens up the possibility of much better control on covariance matrices estimate for cosmological parameters inference, and on the impact of systematics due to baryonic effects. We present a machine learning approach to emulate the hydrodynamical simulation of intergalactic medium physics for the Lyman-$α$ forest called LyAl-Net. The main goal of this work is to provide highly efficient and cheap simulations retaining interpretation abilities about the gas field level, and as a tool for other cosmological exploration. We use a neural network based on the U-net architecture, a variant of convolutional neural networks, to predict the neutral hydrogen physical properties, density, and temperature. We train the LyAl-Net model with the Horizon-noAGN simulation, though using only 9% of the volume. We also explore the resilience of the model through tests of a transfer learning framework using cosmological simulations containing different baryonic feedback. We test our results by analysing one and two-point statistics of emulated fields in different scenarios, as well as their stochastic properties. The ensemble average of the emulated Lyman-$α$ forest absorption as a function of redshift lies within 2.5% of one derived from the full hydrodynamical simulation. The computation of individual fields from the dark matter density agrees well with regular physical regimes of cosmological fields. The results tested on IllustrisTNG100 showed a drastic improvement in the Lyman-$α$ forest flux without arbitrary rescaling. △ Less

Submitted 31 March, 2023; originally announced March 2023.

Fixed-kinetic Neural Hamiltonian Flows for enhanced interpretability and reduced complexity

Authors: Vincent Souveton, Arnaud Guillin, Jens Jasche, Guilhem Lavaux, Manon Michel

Abstract: Normalizing Flows (NF) are Generative models which are particularly robust and allow for exact sampling of the learned distribution. They however require the design of an invertible mapping, whose Jacobian determinant has to be computable. Recently introduced, Neural Hamiltonian Flows (NHF) are based on Hamiltonian dynamics-based Flows, which are continuous, volume-preserving and invertible and th… ▽ More Normalizing Flows (NF) are Generative models which are particularly robust and allow for exact sampling of the learned distribution. They however require the design of an invertible mapping, whose Jacobian determinant has to be computable. Recently introduced, Neural Hamiltonian Flows (NHF) are based on Hamiltonian dynamics-based Flows, which are continuous, volume-preserving and invertible and thus make for natural candidates for robust NF architectures. In particular, their similarity to classical Mechanics could lead to easier interpretability of the learned mapping. However, despite being Physics-inspired architectures, the originally introduced NHF architecture still poses a challenge to interpretability. For this reason, in this work, we introduce a fixed kinetic energy version of the NHF model. Inspired by physics, our approach improves interpretability and requires less parameters than previously proposed architectures. We then study the robustness of the NHF architectures to the choice of hyperparameters. We analyze the impact of the number of leapfrog steps, the integration time and the number of neurons per hidden layer, as well as the choice of prior distribution, on sampling a multimodal 2D mixture. The NHF architecture is robust to these choices, especially the fixed-kinetic energy model. Finally, we adapt NHF to the context of Bayesian inference and illustrate our method on sampling the posterior distribution of two cosmological parameters knowing type Ia supernovae observations. △ Less

Submitted 3 February, 2023; originally announced February 2023.

Higher-order statistics of the large-scale structure from photometric redshifts

Authors: Eleni Tsaprazi, Jens Jasche, Guilhem Lavaux, Florent Leclercq

Abstract: The large-scale structure is a major source of cosmological information. However, next-generation photometric galaxy surveys will only provide a distorted view of cosmic structures due to large redshift uncertainties. To address the need for accurate reconstructions of the large-scale structure in presence of photometric uncertainties, we present a framework that constrains the three-dimensional d… ▽ More The large-scale structure is a major source of cosmological information. However, next-generation photometric galaxy surveys will only provide a distorted view of cosmic structures due to large redshift uncertainties. To address the need for accurate reconstructions of the large-scale structure in presence of photometric uncertainties, we present a framework that constrains the three-dimensional dark matter density jointly with galaxy photometric redshift probability density functions (PDFs), exploiting information from galaxy clustering. Our forward model provides Markov Chain Monte Carlo realizations of the primordial and present-day dark matter density, inferred jointly from data. Our method goes beyond 2-point statistics via field-level inference. It accounts for all observational uncertainties and the survey geometry. We showcase our method using mock catalogs that emulate next-generation surveys with a worst-case redshift uncertainty, equivalent to ${\sim}300$ Mpc. On scales $150$ Mpc, we improve the cross-correlation of the photometric galaxy positions with the ground truth from $28\%$ to $86\%$. The improvement is significant down to $13$ Mpc. On scales $150$ Mpc, we achieve a cross-correlation of $80-90\%$ with the ground truth for the dark matter density, radial peculiar velocities, tidal shear and gravitational potential. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: 16 pages, 12 figures, submitted to A&A

Constraints on dark matter annihilation and decay from the large-scale structure of the nearby universe

Authors: Deaglan J. Bartlett, Andrija Kostić, Harry Desmond, Jens Jasche, Guilhem Lavaux

Abstract: Decaying or annihilating dark matter particles could be detected through gamma-ray emission from the species they decay or annihilate into. This is usually done by modelling the flux from specific dark matter-rich objects such as the Milky Way halo, Local Group dwarfs, and nearby groups. However, these objects are expected to have significant emission from baryonic processes as well, and the analy… ▽ More Decaying or annihilating dark matter particles could be detected through gamma-ray emission from the species they decay or annihilate into. This is usually done by modelling the flux from specific dark matter-rich objects such as the Milky Way halo, Local Group dwarfs, and nearby groups. However, these objects are expected to have significant emission from baryonic processes as well, and the analyses discard gamma-ray data over most of the sky. Here we construct full-sky templates for gamma-ray flux from the large-scale structure within $\sim$200 Mpc by means of a suite of constrained $N$-body simulations (CSiBORG) produced using the Bayesian Origin Reconstruction from Galaxies algorithm. Marginalising over uncertainties in this reconstruction, small-scale structure, and parameters describing astrophysical contributions to the observed gamma-ray sky, we compare to observations from the Fermi Large Area Telescope to constrain dark matter annihilation cross sections and decay rates through a Markov Chain Monte Carlo analysis. We rule out the thermal relic cross section for $s$-wave annihilation for all $m_χ\lesssim 7 {\rm \, GeV}/c^2$ at 95\% confidence if the annihilation produces gluons or quarks less massive than the bottom quark. We infer a contribution to the gamma-ray sky with the same spatial distribution as dark matter decay at $3.3σ$. Although this could be due to dark matter decay via these channels with a decay rate $Γ\approx 6 \times 10^{-28} {\rm \, s^{-1}}$, we find that a power-law spectrum of index $p=-2.75^{+0.71}_{-0.46}$, likely of baryonic origin, is preferred by the data. △ Less

Submitted 30 November, 2022; v1 submitted 25 May, 2022; originally announced May 2022.

Comments: 24 pages, 9 figures, 1 table. Accepted for publication in Physical Review D

Journal ref: Phys. Rev. D 106, 103526 (2022)

Euclid: Cosmological forecasts from the void size function

Authors: S. Contarini, G. Verza, A. Pisani, N. Hamaus, M. Sahlén, C. Carbone, S. Dusini, F. Marulli, L. Moscardini, A. Renzi, C. Sirignano, L. Stanco, M. Aubert, M. Bonici, G. Castignani, H. M. Courtois, S. Escoffier, D. Guinet, A. Kovacs, G. Lavaux, E. Massara, S. Nadathur, G. Pollina, T. Ronconi, F. Ruppin , et al. (101 additional authors not shown)

Abstract: The Euclid mission $-$ with its spectroscopic galaxy survey covering a sky area over $15\,000 \ \mathrm{deg}^2$ in the redshift range $0.9<z<1.8\ -$ will provide a sample of tens of thousands of cosmic voids. This paper explores for the first time the constraining power of the void size function on the properties of dark energy (DE) from a survey mock catalogue, the official Euclid Flagship simula… ▽ More The Euclid mission $-$ with its spectroscopic galaxy survey covering a sky area over $15\,000 \ \mathrm{deg}^2$ in the redshift range $0.9<z<1.8\ -$ will provide a sample of tens of thousands of cosmic voids. This paper explores for the first time the constraining power of the void size function on the properties of dark energy (DE) from a survey mock catalogue, the official Euclid Flagship simulation. We identify voids in the Flagship light-cone, which closely matches the features of the upcoming Euclid spectroscopic data set. We model the void size function considering a state-of-the art methodology: we rely on the volume conserving (Vdn) model, a modification of the popular Sheth & van de Weygaert model for void number counts, extended by means of a linear function of the large-scale galaxy bias. We find an excellent agreement between model predictions and measured mock void number counts. We compute updated forecasts for the Euclid mission on DE from the void size function and provide reliable void number estimates to serve as a basis for further forecasts of cosmological applications using voids. We analyse two different cosmological models for DE: the first described by a constant DE equation of state parameter, $w$, and the second by a dynamic equation of state with coefficients $w_0$ and $w_a$. We forecast $1σ$ errors on $w$ lower than $10\%$, and we estimate an expected figure of merit (FoM) for the dynamical DE scenario $\mathrm{FoM}_{w_0,w_a} = 17$ when considering only the neutrino mass as additional free parameter of the model. The analysis is based on conservative assumptions to ensure full robustness, and is a pathfinder for future enhancements of the technique. Our results showcase the impressive constraining power of the void size function from the Euclid spectroscopic sample, both as a stand-alone probe, and to be combined with other Euclid cosmological probes. △ Less

Submitted 25 November, 2022; v1 submitted 23 May, 2022; originally announced May 2022.

Comments: 19 pages, 7 figures, 4 tables - published in A&A

Journal ref: A&A 667, A162 (2022)

The Milky Way's plane of satellites: consistent with $Λ$CDM

Authors: Till Sawala, Marius Cautun, Carlos S. Frenk, John Helly, Jens Jasche, Adrian Jenkins, Peter H. Johansson, Guilhem Lavaux, Stuart McAlpine, Matthieu Schaller

Abstract: The "plane of satellites problem" describes the arrangement of the Milky Way's 11 brightest satellite galaxies in a remarkably thin plane, possibly supported by rotation. This is in apparent contradiction to the standard cosmological model, wherein the Galaxy is surrounded by a dispersion-supported dark matter halo. Here, we show that the reported exceptional anisotropy of the satellite system is… ▽ More The "plane of satellites problem" describes the arrangement of the Milky Way's 11 brightest satellite galaxies in a remarkably thin plane, possibly supported by rotation. This is in apparent contradiction to the standard cosmological model, wherein the Galaxy is surrounded by a dispersion-supported dark matter halo. Here, we show that the reported exceptional anisotropy of the satellite system is strongly contingent on a lopsided radial distribution, which earlier simulations have failed to reproduce, combined with the close but fleeting conjunction of the two most distant satellites, Leo I and Leo II. Using Gaia proper motions, we show that the orbital pole alignment is much more common than previously reported, and reveal the plane of satellites to be transient rather than rotationally supported. Comparing to new simulations, where such short-lived planes are common, we find the Milky Way satellites to be compatible with standard model expectations. △ Less

Submitted 5 May, 2022; originally announced May 2022.

LyMAS reloaded: improving the predictions of the large-scale Lyman-α forest statistics from dark matter density and velocity fields

Authors: S. Peirani, S. Prunet, S. Colombi, C. Pichon, D. H. Weinberg, C. Laigle, G. Lavaux, Y. Dubois, J. Devriendt

Abstract: We present LyMAS2, an improved version of the "Lyman-α Mass Association Scheme" aiming at predicting the large-scale 3d clustering statistics of the Lyman-α forest (Ly-α) from moderate resolution simulations of the dark matter (DM) distribution, with prior calibrations from high resolution hydrodynamical simulations of smaller volumes. In this study, calibrations are derived from the Horizon-AGN s… ▽ More We present LyMAS2, an improved version of the "Lyman-α Mass Association Scheme" aiming at predicting the large-scale 3d clustering statistics of the Lyman-α forest (Ly-α) from moderate resolution simulations of the dark matter (DM) distribution, with prior calibrations from high resolution hydrodynamical simulations of smaller volumes. In this study, calibrations are derived from the Horizon-AGN suite simulations, (100 Mpc/h)^3 comoving volume, using Wiener filtering, combining information from dark matter density and velocity fields (i.e. velocity dispersion, vorticity, line of sight 1d-divergence and 3d-divergence). All new predictions have been done at z=2.5 in redshift-space, while considering the spectral resolution of the SDSS-III BOSS Survey and different dark matter smoothing (0.3, 0.5 and 1.0 Mpc/h comoving). We have tried different combinations of dark matter fields and found that LyMAS2, applied to the Horizon-noAGN dark matter fields, significantly improves the predictions of the Ly-α 3d clustering statistics, especially when the DM overdensity is associated with the velocity dispersion or the vorticity fields. Compared to the hydrodynamical simulation trends, the 2-point correlation functions of pseudo-spectra generated with LyMAS2 can be recovered with relative differences of ~5% even for high angles, the flux 1d power spectrum (along the light of sight) with ~2% and the flux 1d probability distribution function exactly. Finally, we have produced several large mock BOSS spectra (1.0 and 1.5 Gpc/h) expected to lead to much more reliable and accurate theoretical predictions. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: Submitted to MNRAS (one iteration with the referee) - 29 pages

Field-Based Physical Inference From Peculiar Velocity Tracers

Authors: James Prideaux-Ghee, Florent Leclercq, Guilhem Lavaux, Alan Heavens, Jens Jasche

Abstract: We present a Bayesian hierarchical modelling approach to reconstruct the initial cosmic matter density field constrained by peculiar velocity observations. As our approach features a model for the gravitational evolution of dark matter to connect the initial conditions to late-time observations, it reconstructs the final density and velocity fields as natural byproducts. We implement this field-ba… ▽ More We present a Bayesian hierarchical modelling approach to reconstruct the initial cosmic matter density field constrained by peculiar velocity observations. As our approach features a model for the gravitational evolution of dark matter to connect the initial conditions to late-time observations, it reconstructs the final density and velocity fields as natural byproducts. We implement this field-based physical inference approach by adapting the Bayesian Origin Reconstruction from Galaxies (BORG) algorithm, which explores the high-dimensional posterior through the use of Hamiltonian Monte Carlo sampling. We test the self-consistency of the method using random sets of mock tracers, and assess its accuracy in a more complex scenario where peculiar velocity tracers are non-linearly evolved mock haloes. We find that our framework self-consistently infers the initial conditions, density and velocity fields, and shows some robustness to model mis-specification. As compared to the state-of-the-art approach of constrained Gaussian random fields/Wiener filtering, our method produces more accurate final density and velocity field reconstructions. It also allows us to constrain the initial conditions by peculiar velocity observations, complementing in this aspect previous field-based approaches based on other cosmological observables. △ Less

Submitted 7 December, 2022; v1 submitted 31 March, 2022; originally announced April 2022.

Comments: 23 pages, 15 figures. Accepted for publication in MNRAS

Bayesian field-level inference of primordial non-Gaussianity using next-generation galaxy surveys

Authors: Adam Andrews, Jens Jasche, Guilhem Lavaux, Fabian Schmidt

Abstract: Detecting and measuring a non-Gaussian signature of primordial origin in the density field is a major science goal of next-generation galaxy surveys. The signal will permit us to determine primordial physics processes and constrain models of cosmic inflation. While traditional approaches utilise a limited set of statistical summaries of the galaxy distribution to constrain primordial non-Gaussiani… ▽ More Detecting and measuring a non-Gaussian signature of primordial origin in the density field is a major science goal of next-generation galaxy surveys. The signal will permit us to determine primordial physics processes and constrain models of cosmic inflation. While traditional approaches utilise a limited set of statistical summaries of the galaxy distribution to constrain primordial non-Gaussianity, we present a field-level approach by Bayesian forward-modelling the entire three-dimensional galaxy survey. Our method naturally and fully self-consistently exploits the entirety of the large-scale structure, e.g., higher-order statistics, peculiar velocity fields, and scale-dependent galaxy bias, to extract information on the local non-Gaussianity parameter, $\fnl$. We demonstrate the performance of our approach through various tests with mock galaxy data emulating relevant features of the \sdssiii{}-like survey, and additional tests with a \textit{Stage IV} mock data set. These tests reveal that the method infers unbiased values of $\fnl$ by accurately handling survey geometries, noise, and unknown galaxy biases. We demonstrate that our method can achieve constraints of $σ_{\fnl} \approx 8.78$ for \sdssiii{}-like data, an improvement of a factor $\sim 2.5$ over currently published constraints. Tests with next-generation mock data show that significant further improvements are feasible with sufficiently high resolution. Furthermore, the results demonstrate that our method can consistently marginalise all nuisance parameters of the data model. The method further provides an inference of the three-dimensional primordial density field, providing opportunities to explore additional signatures of primordial physics. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: 17 pages, 16 figures. Submitted to MNRAS. Comments are welcome

SIBELIUS-DARK: a galaxy catalogue of the Local Volume from a constrained realisation simulation

Authors: Stuart McAlpine, John C. Helly, Matthieu Schaller, Till Sawala, Guilhem Lavaux, Jens Jasche, Carlos S. Frenk, Adrian Jenkins, John R. Lucey, Peter H. Johansson

Abstract: We present SIBELIUS-DARK, a constrained realisation simulation of the local volume to a distance of 200~Mpc from the Milky Way. SIBELIUS-DARK is the first study of the \textit{Simulations Beyond The Local Universe} (SIBELIUS) project, which has the goal of embedding a model Local Group-like system within the correct cosmic environment. The simulation is dark-matter-only, with the galaxy population… ▽ More We present SIBELIUS-DARK, a constrained realisation simulation of the local volume to a distance of 200~Mpc from the Milky Way. SIBELIUS-DARK is the first study of the \textit{Simulations Beyond The Local Universe} (SIBELIUS) project, which has the goal of embedding a model Local Group-like system within the correct cosmic environment. The simulation is dark-matter-only, with the galaxy population calculated using the semi-analytic model of galaxy formation, GALFORM. We demonstrate that the large-scale structure that emerges from the SIBELIUS constrained initial conditions matches well the observational data. The inferred galaxy population of SIBELIUS-DARK also match well the observational data, both statistically for the whole volume and on an object-by-object basis for the most massive clusters. For example, the $K$-band number counts across the whole sky, and when divided between the northern and southern Galactic hemispheres, are well reproduced by SIBELIUS-DARK. We find that the local volume is somewhat unusual in the wider context of $Λ$CDM: it contains an abnormally high number of supermassive clusters, as well as an overall large-scale underdensity at the level of $\approx 5$\% relative to the cosmic mean. However, whilst rare, the extent of these peculiarities does not significantly challenge the $Λ$CDM model. SIBELIUS-DARK is the most comprehensive constrained realisation simulation of the local volume to date, and with this paper we publicly release the halo and galaxy catalogues at $z=0$, which we hope will be useful to the wider astronomy community. △ Less

Submitted 8 February, 2022; originally announced February 2022.

Comments: 26 pages, 19 figures. Paper as been accepted for publication in MNRAS. Comments welcome

Rubin-Euclid Derived Data Products: Initial Recommendations

Authors: Leanne P. Guy, Jean-Charles Cuillandre, Etienne Bachelet, Manda Banerji, Franz E. Bauer, Thomas Collett, Christopher J. Conselice, Siegfried Eggl, Annette Ferguson, Adriano Fontana, Catherine Heymans, Isobel M. Hook, Éric Aubourg, Hervé Aussel, James Bosch, Benoit Carry, Henk Hoekstra, Konrad Kuijken, Francois Lanusse, Peter Melchior, Joseph Mohr, Michele Moresco, Reiko Nakajima, Stéphane Paltani, Michael Troxel , et al. (95 additional authors not shown)

Abstract: This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum… ▽ More This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum and a series of virtual meetings. Strong interest in enhancing science with joint DDPs emerged from across a wide range of astrophysical domains: Solar System, the Galaxy, the Local Volume, from the nearby to the primaeval Universe, and cosmology. △ Less

Submitted 13 October, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: Report of the Rubin-Euclid Derived Data Products Working Group, 78 pages, 11 figures

Field-level inference of galaxy intrinsic alignment from the SDSS-III BOSS survey

Authors: Eleni Tsaprazi, Nhat-Minh Nguyen, Jens Jasche, Fabian Schmidt, Guilhem Lavaux

Abstract: As a large-scale overdensity collapses, it affects the orientation and shape of galaxies that form, by exerting tidal shear along their axes. Therefore, the shapes of elliptical galaxies align with the tidal field of cosmic structures. This intrinsic alignment provides insights into galaxy formation and the primordial universe, complements late-time cosmological probes and constitutes a significan… ▽ More As a large-scale overdensity collapses, it affects the orientation and shape of galaxies that form, by exerting tidal shear along their axes. Therefore, the shapes of elliptical galaxies align with the tidal field of cosmic structures. This intrinsic alignment provides insights into galaxy formation and the primordial universe, complements late-time cosmological probes and constitutes a significant systematic effect for weak gravitational lensing observations. In the present study, we provide constraints on the linear alignment model using a fully Bayesian field-level approach, using galaxy shape measurements from the SDSS-III BOSS LOWZ sample and three-dimensional tidal fields constrained with the LOWZ and CMASS galaxy samples of the SDSS-III BOSS survey. We find 4$σ$ evidence of intrinsic alignment, with an amplitude of $A_I=2.9 \pm 0.7$ at 20$h^{-1}\;\mathrm{Mpc}$. △ Less

Submitted 29 June, 2022; v1 submitted 8 December, 2021; originally announced December 2021.

Comments: 26 pages, 8 figures, accepted by JCAP

Reconstructing dark matter distribution with peculiar velocities: Bayesian forward modelling with corrections for inhomogeneous Malmquist bias

Authors: Supranta S. Boruah, Guilhem Lavaux, Michael J. Hudson

Abstract: We present a forward-modelled velocity field reconstruction algorithm that performs the reconstruction of the mass density field using only peculiar velocity data. Our method consistently accounts for the inhomogeneous Malmquist bias using analytic integration along the line-of-sight. By testing our method on a simulation, we show that our method gives an unbiased reconstruction of the velocity fi… ▽ More We present a forward-modelled velocity field reconstruction algorithm that performs the reconstruction of the mass density field using only peculiar velocity data. Our method consistently accounts for the inhomogeneous Malmquist bias using analytic integration along the line-of-sight. By testing our method on a simulation, we show that our method gives an unbiased reconstruction of the velocity field. We show that not accounting for the inhomogeneous Malmquist bias can lead to significant biases in the forward-modelled reconstructions. We applied our method to a peculiar velocity data set consisting of the SFI++ and 2MTF Tully-Fisher catalogues and the A2 supernovae compilation, thus obtaining a novel velocity reconstruction in the local Universe. Our velocity reconstructions have a cosmological power spectrum consistent with the theoretical expectation. Furthermore, we obtain a full description of the uncertainties on reconstruction through samples of the posterior distribution. We validate our velocity reconstruction of the local Universe by comparing it to an independent reconstruction using the 2M++ galaxy catalogue, obtaining good agreement between the two reconstructions. Using Bayesian model comparison, we find that our velocity model performs better than the adaptive kernel smoothed velocity with the same peculiar velocity data. However, our velocity model does not perform as well as the velocity reconstruction from the 2M++ galaxy catalogue, due to the sparse and noisy nature of the peculiar velocity tracer samples. The method presented here provides a way to include peculiar velocity data in initial condition reconstruction frameworks. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Comments: 15 pages, 12 figures, To be submitted to MNRAS

Detectability of large correlation length inflationary magnetic field with Cherenkov telescopes

Authors: Alexander Korochkin, Andrii Neronov, Guilhem Lavaux, Marius Ramsoy, Dmitri Semikoz

Abstract: Magnetic fields occupying the voids of the large scale structure may be a relic from the Early Universe originating from either Inflation or from cosmological phase transitions. We explore the possibility of identifying the inflationary origin of the void magnetic fields and measuring its parameters with gamma-ray astronomy methods. The large correlation length inflationary field is expected to im… ▽ More Magnetic fields occupying the voids of the large scale structure may be a relic from the Early Universe originating from either Inflation or from cosmological phase transitions. We explore the possibility of identifying the inflationary origin of the void magnetic fields and measuring its parameters with gamma-ray astronomy methods. The large correlation length inflationary field is expected to impose a characteristic asymmetry of extended gamma-ray emission that is correlated between different sources on the sky. We show that a set of nearby blazars for which the extended emission is observable in the 0.1-1 TeV band with CTA can be used for the test of inflationary origin of the void magnetic fields. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: 8 pages, 6 figures, accepted to JETP

GLADE+: An Extended Galaxy Catalogue for Multimessenger Searches with Advanced Gravitational-wave Detectors

Authors: G. Dálya, R. Díaz, F. R. Bouchet, Z. Frei, J. Jasche, G. Lavaux, R. Macas, S. Mukherjee, M. Pálfi, R. S. de Souza, B. D. Wandelt, M. Bilicki, P. Raffai

Abstract: We present GLADE+, an extended version of the GLADE galaxy catalogue introduced in our previous paper for multimessenger searches with advanced gravitational-wave detectors. GLADE+ combines data from six separate but not independent astronomical catalogues: the GWGC, 2MPZ, 2MASS XSC, HyperLEDA, and WISExSCOSPZ galaxy catalogues, and the SDSS-DR16Q quasar catalogue. To allow corrections of CMB-fram… ▽ More We present GLADE+, an extended version of the GLADE galaxy catalogue introduced in our previous paper for multimessenger searches with advanced gravitational-wave detectors. GLADE+ combines data from six separate but not independent astronomical catalogues: the GWGC, 2MPZ, 2MASS XSC, HyperLEDA, and WISExSCOSPZ galaxy catalogues, and the SDSS-DR16Q quasar catalogue. To allow corrections of CMB-frame redshifts for peculiar motions, we calculated peculiar velocities along with their standard deviations of all galaxies having $B$-band magnitude data within redshift $z=0.05$ using the "Bayesian Origin Reconstruction from Galaxies" formalism. GLADE+ is complete up to luminosity distance $d_L=47^{+4}_{-2}$ Mpc in terms of the total expected $B$-band luminosity of galaxies, and contains all of the brightest galaxies giving 90\% of the total $B$-band and $K$-band luminosity up to $d_L\simeq 130$ Mpc. We include estimations of stellar masses and individual binary neutron star merger rates for galaxies with $W1$ magnitudes. These parameters can help in ranking galaxies in a given gravitational wave localization volume in terms of their likelihood of being hosts, thereby possibly reducing the number of pointings and total integration time needed to find the electromagnetic counterpart. △ Less

Submitted 2 June, 2022; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: 9 pages, 4 figures, accepted for publication in MNRAS

Euclid: Forecasts from redshift-space distortions and the Alcock-Paczynski test with cosmic voids

Authors: N. Hamaus, M. Aubert, A. Pisani, S. Contarini, G. Verza, M. -C. Cousinou, S. Escoffier, A. Hawken, G. Lavaux, G. Pollina, B. D. Wandelt, J. Weller, M. Bonici, C. Carbone, L. Guzzo, A. Kovacs, F. Marulli, E. Massara, L. Moscardini, P. Ntelis, W. J. Percival, S. Radinović, M. Sahlén, Z. Sakr, A. G. Sánchez , et al. (105 additional authors not shown)

Abstract: Euclid is poised to survey galaxies across a cosmological volume of unprecedented size, providing observations of more than a billion objects distributed over a third of the full sky. Approximately 20 million of these galaxies will have their spectroscopy available, allowing us to map the 3D large-scale structure of the Universe in great detail. This paper investigates prospects for the detection… ▽ More Euclid is poised to survey galaxies across a cosmological volume of unprecedented size, providing observations of more than a billion objects distributed over a third of the full sky. Approximately 20 million of these galaxies will have their spectroscopy available, allowing us to map the 3D large-scale structure of the Universe in great detail. This paper investigates prospects for the detection of cosmic voids therein and the unique benefit they provide for cosmology. In particular, we study the imprints of dynamic and geometric distortions of average void shapes and their constraining power on the growth of structure and cosmological distance ratios. To this end, we made use of the Flagship mock catalog, a state-of-the-art simulation of the data expected to be observed with Euclid. We arranged the data into four adjacent redshift bins, each of which contains about 11000 voids and estimated the stacked void-galaxy cross-correlation function in every bin. Fitting a linear-theory model to the data, we obtained constraints on $f/b$ and $D_M H$, where $f$ is the linear growth rate of density fluctuations, $b$ the galaxy bias, $D_M$ the comoving angular diameter distance, and $H$ the Hubble rate. In addition, we marginalized over two nuisance parameters included in our model to account for unknown systematic effects. With this approach, Euclid will be able to reach a relative precision of about 4% on measurements of $f/b$ and 0.5% on $D_M H$ in each redshift bin. Better modeling or calibration of the nuisance parameters may further increase this precision to 1% and 0.4%, respectively. Our results show that the exploitation of cosmic voids in Euclid will provide competitive constraints on cosmology even as a stand-alone probe. For example, the equation-of-state parameter $w$ for dark energy will be measured with a precision of about 10%, consistent with previous more approximate forecasts. △ Less

Submitted 2 December, 2021; v1 submitted 23 August, 2021; originally announced August 2021.

Comments: 15 pages, 7 figures. Accepted by A&A (Oct. 31)

Journal ref: A&A 658, A20 (2022)

Lifting weak lensing degeneracies with a field-based likelihood

Authors: Natalia Porqueres, Alan Heavens, Daniel Mortlock, Guilhem Lavaux

Abstract: We present a field-based approach to the analysis of cosmic shear data to infer jointly cosmological parameters and the dark matter distribution. This forward modelling approach samples the cosmological parameters and the initial matter fluctuations, using a physical gravity model to link the primordial fluctuations to the non-linear matter distribution. Cosmological parameters are sampled and upd… ▽ More We present a field-based approach to the analysis of cosmic shear data to infer jointly cosmological parameters and the dark matter distribution. This forward modelling approach samples the cosmological parameters and the initial matter fluctuations, using a physical gravity model to link the primordial fluctuations to the non-linear matter distribution. Cosmological parameters are sampled and updated consistently through the forward model, varying (1) the initial matter power spectrum, (2) the geometry through the distance-redshift relationship, and (3) the growth of structure and light-cone effects. Our approach extracts more information from the data than methods based on two-point statistics. We find that this field-based approach lifts the strong degeneracy between the cosmological matter density, $Ω_\mathrm{m}$, and the fluctuation amplitude, $σ_8$, providing tight constraints on these parameters from weak lensing data alone. In the simulated four-bin tomographic experiment we consider, the field-based likelihood yields marginal uncertainties on $σ_8$ and $Ω_\mathrm{m}$ that are, respectively, a factor of 3 and 5 smaller than those from a two-point power spectrum analysis applied to the same underlying data. △ Less

Submitted 3 November, 2021; v1 submitted 10 August, 2021; originally announced August 2021.

Comments: Accepted in MNRAS

Optimal machine-driven acquisition of future cosmological data

Authors: Andrija Kostić, Jens Jasche, Doogesh Kodi Ramanah, Guilhem Lavaux

Abstract: We present maps classifying regions of the sky according to their information gain potential as quantified by the Fisher information. These maps can guide the optimal retrieval of relevant physical information with targeted cosmological searches. Specifically, we calculate the response of observed cosmic structures to perturbative changes in the cosmological model and chart their respective contri… ▽ More We present maps classifying regions of the sky according to their information gain potential as quantified by the Fisher information. These maps can guide the optimal retrieval of relevant physical information with targeted cosmological searches. Specifically, we calculate the response of observed cosmic structures to perturbative changes in the cosmological model and chart their respective contributions to the Fisher information. Our physical forward modeling machinery transcends the limitations of contemporary analyses based on statistical summaries to yield detailed characterizations of individual 3D structures. We demonstrate this using galaxy counts data and showcase the potential of our approach by studying the information gain of the Coma cluster. We find that regions in the vicinity of the filaments and cluster core, where mass accretion ensues from gravitational infall, are the most informative about our physical model of structure formation in the Universe. Hence, collecting data in those regions would be most optimal for testing our model predictions. The results presented in this work are the first of their kind and elucidate the inhomogeneous distribution of cosmological information in the Universe. This study paves a new way forward to perform efficient targeted searches for the fundamental physics of the Universe, where search strategies are progressively refined with new cosmological data sets within an active learning framework. △ Less

Submitted 5 December, 2021; v1 submitted 1 July, 2021; originally announced July 2021.

Comments: 6 + 3 pages, 5 figures. Accepted to A&A Letters

Journal ref: A&A 657, L17 (2022)

Bayesian cosmological inference through implicit cross-correlation statistics

Authors: Guilhem Lavaux, Jens Jasche

Abstract: Analyzes of next-generation galaxy data require accurate treatment of systematic effects such as the bias between observed galaxies and the underlying matter density field. However, proposed models of the phenomenon are either numerically expensive or too inaccurate to achieve unbiased inferences of cosmological parameters even at mildly-nonlinear scales of the data. As an alternative to construct… ▽ More Analyzes of next-generation galaxy data require accurate treatment of systematic effects such as the bias between observed galaxies and the underlying matter density field. However, proposed models of the phenomenon are either numerically expensive or too inaccurate to achieve unbiased inferences of cosmological parameters even at mildly-nonlinear scales of the data. As an alternative to constructing accurate galaxy bias models, requiring understanding galaxy formation, we propose to construct likelihood distributions for Bayesian forward modeling approaches that are insensitive to linear, scale-dependent bias and provide robustness against model misspecification. We use maximum entropy arguments to construct likelihood distributions designed to account only for correlations between data and inferred quantities. By design these correlations are insensitive to linear galaxy biasing relations, providing the desired robustness. The method is implemented and tested within a Markov Chain Monte Carlo approach. The method is assessed using a halo mock catalog based on standard full, cosmological, N-body simulations. We obtain unbiased and tight constraints on cosmological parameters exploiting only linear cross-correlation rates for $k\le 0.10$ Mpc/h. Tests for halos of masses ~10$^{12}$ M$_\odot$ to ~10$^{13}$ M$_\odot$ indicate that it is possible to ignore all details of the linear, scale dependent, bias function while obtaining robust constraints on cosmology. Our results provide a promising path forward to analyzes of galaxy surveys without the requirement of having to accurately model the details of galaxy biasing but by designing robust likelihoods for the inference. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: 11 pages, 9 figures, comments welcome

The SIBELIUS Project: E Pluribus Unum

Authors: Till Sawala, Stuart McAlpine, Jens Jasche, Guilhem Lavaux, Adrian Jenkins, Peter H. Johansson, Carlos S. Frenk

Abstract: We introduce "Simulations Beyond The Local Universe" (SIBELIUS) that connect the Local Group to its cosmic environment. We show that introducing hierarchical small-scale perturbations to a density field constrained on large scales by observations provides an efficient way to explore the sample space of Local Group analogues. From more than 60 000 simulations, we identify a hierarchy of Local Group… ▽ More We introduce "Simulations Beyond The Local Universe" (SIBELIUS) that connect the Local Group to its cosmic environment. We show that introducing hierarchical small-scale perturbations to a density field constrained on large scales by observations provides an efficient way to explore the sample space of Local Group analogues. From more than 60 000 simulations, we identify a hierarchy of Local Group characteristics emanating from different scales: the total mass, orientation, orbital energy and the angular momentum are largely determined by modes above $λ$ = 1.6 comoving Mpc (cMpc) in the primordial density field. Smaller scale variations are mostly manifest as perturbations to the MW-M31 orbit, and we find that the observables commonly used to describe the Local Group -- the MW-M31 separation and radial velocity -- are transient and depend on specifying scales down to 0.2 cMpc in the primordial density field. We further find that the presence of M33/LMC analogues significantly affects the MW-M31 orbit and its sensitivity to small-scale perturbations. We construct initial conditions that lead to the formation of a Local Group whose primary observables precisely match the current observations. △ Less

Submitted 22 March, 2021; originally announced March 2021.

Comments: 14 pages, 15 figures. Submitted to MNRAS, comments welcome!

Gravitation And the Universe from large Scale-Structures: The GAUSS mission concept

Authors: Alain Blanchard, Aubourg Éric, Philippe Brax, Francisco J. Castender, Sandrine Codis, Stéphanie Escoffier, Fabien Dournac, Agnès Ferté, Fabio Finelli, Pablo Fosalba, Emmanuel Gangler, A Gontcho Satya Gontcho, Adam Hawken, Stéphane Ilić, Jean-Paul Kneib, Martin Kunz, Guilhem Lavaux, Olivier Le Fèvre, Julien Lesgourgues, Yannick Mellier, Jérémy Neveu, Yann Rasera, Cécile Renault, Marina Ricci, Ziad Sakr , et al. (3 additional authors not shown)

Abstract: Today, thanks in particular to the results of the ESA Planck mission, the concordance cosmological model appears to be the most robust to describe the evolution and content of the Universe from its early to late times. It summarizes the evolution of matter, made mainly of dark matter, from the primordial fluctuations generated by inflation around $10^{-30}$ second after the Big-Bang to galaxies an… ▽ More Today, thanks in particular to the results of the ESA Planck mission, the concordance cosmological model appears to be the most robust to describe the evolution and content of the Universe from its early to late times. It summarizes the evolution of matter, made mainly of dark matter, from the primordial fluctuations generated by inflation around $10^{-30}$ second after the Big-Bang to galaxies and clusters of galaxies, 13.8 billion years later, and the evolution of the expansion of space, with a relative slowdown in the matter-dominated era and, since a few billion years, an acceleration powered by dark energy. But we are far from knowing the pillars of this model which are inflation, dark matter and dark energy. Comprehending these fundamental questions requires a detailed mapping of our observable Universe over the whole of cosmic time. The relic radiation provides the starting point and galaxies draw the cosmic web. JAXA's LiteBIRD mission will map the beginning of our Universe with a crucial test for inflation (its primordial gravity waves), and the ESA Euclid mission will map the most recent half part, crucial for dark energy. The mission concept, described in this White Paper, GAUSS, aims at being a mission to fully map the cosmic web up to the reionization era, linking early and late evolution, to tackle and disentangle the crucial degeneracies persisting after the Euclid era between dark matter and inflation properties, dark energy, structure growth and gravitation at large scale. △ Less

Submitted 7 February, 2021; originally announced February 2021.

Comments: 24 pages, 3 figures, ESA Voyage2050 White Paper

Journal ref: Experimental Astronomy, Volume 51, Issue 3, p.1623-1640 (2021)

Bayesian forward modelling of cosmic shear data

Authors: Natalia Porqueres, Alan Heavens, Daniel Mortlock, Guilhem Lavaux

Abstract: We present a Bayesian hierarchical modelling approach to infer the cosmic matter density field, and the lensing and the matter power spectra, from cosmic shear data. This method uses a physical model of cosmic structure formation to infer physically plausible cosmic structures, which accounts for the non-Gaussian features of the gravitationally evolved matter distribution and light-cone effects. W… ▽ More We present a Bayesian hierarchical modelling approach to infer the cosmic matter density field, and the lensing and the matter power spectra, from cosmic shear data. This method uses a physical model of cosmic structure formation to infer physically plausible cosmic structures, which accounts for the non-Gaussian features of the gravitationally evolved matter distribution and light-cone effects. We test and validate our framework with realistic simulated shear data, demonstrating that the method recovers the unbiased matter distribution and the correct lensing and matter power spectrum. While the cosmology is fixed in this test, and the method employs a prior power spectrum, we demonstrate that the lensing results are sensitive to the true power spectrum when this differs from the prior. In this case, the density field samples are generated with a power spectrum that deviates from the prior, and the method recovers the true lensing power spectrum. The method also recovers the matter power spectrum across the sky, but as currently implemented, it cannot determine the radial power since isotropy is not imposed. In summary, our method provides physically plausible inference of the dark matter distribution from cosmic shear data, allowing us to extract information beyond the two-point statistics and exploiting the full information content of the cosmological fields. △ Less

Submitted 21 January, 2021; v1 submitted 13 November, 2020; originally announced November 2020.

Comments: 10 pages, 13 figures. Accepted for publication in MNRAS

Impacts of the physical data model on the forward inference of initial conditions from biased tracers

Authors: Nhat-Minh Nguyen, Fabian Schmidt, Guilhem Lavaux, Jens Jasche

Abstract: We investigate the impact of each ingredient in the employed physical data model on the Bayesian forward inference of initial conditions from biased tracers at the field level. Specifically, we use dark matter halos in a given cosmological simulation volume as tracers of the underlying matter density field. We study the effect of tracer density, grid resolution, gravity model, bias model and likel… ▽ More We investigate the impact of each ingredient in the employed physical data model on the Bayesian forward inference of initial conditions from biased tracers at the field level. Specifically, we use dark matter halos in a given cosmological simulation volume as tracers of the underlying matter density field. We study the effect of tracer density, grid resolution, gravity model, bias model and likelihood on the inferred initial conditions. We find that the cross-correlation coefficient between true and inferred phases reacts weakly to all ingredients above, and is well predicted by the theoretical expectation derived from a Gaussian model on a broad range of scales. The bias in the amplitude of the inferred initial conditions, on the other hand, depends strongly on the bias model and the likelihood. We conclude that the bias model and likelihood hold the key to an unbiased cosmological inference. Together they must keep the systematics -- which arise from the sub-grid physics that are marginalized over -- under control in order to obtain an unbiased inference. △ Less

Submitted 16 March, 2021; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 29 pages, 11 figures, 1 table; main results in Sec 4.2; comments welcome. v2: Match JCAP accepted version (no major change), update references. v3: Fix links, typos and further update reference info. v4: Match JCAP published version, fix an incorrect link and line spacing

Journal ref: JCAP03(2021)058

Peculiar velocities in the local Universe: comparison of different models and the implications for $H_0$ and dark matter

Authors: Supranta S. Boruah, Michael J. Hudson, Guilhem Lavaux

Abstract: When measuring the value of the Hubble parameter, $H_0$, it is necessary to know the recession velocity free of the effects of peculiar velocities. In this work, we study different models of peculiar velocity in the local Universe. In particular, we compare models based on density reconstruction from galaxy redshift surveys and kernel smoothing of peculiar velocity data. The velocity field from th… ▽ More When measuring the value of the Hubble parameter, $H_0$, it is necessary to know the recession velocity free of the effects of peculiar velocities. In this work, we study different models of peculiar velocity in the local Universe. In particular, we compare models based on density reconstruction from galaxy redshift surveys and kernel smoothing of peculiar velocity data. The velocity field from the density reconstruction is obtained using the 2M++ galaxy redshift compilation, which is compared to two adaptive kernel-smoothed velocity fields: the first obtained from the 6dF Fundamental Plane sample and the other using a Tully-Fisher catalogue obtained by combining SFI++ and 2MTF. We highlight that smoothed velocity fields should be rescaled to obtain unbiased velocity estimates. Comparing the predictions of these models to the observations from a few test sets of peculiar velocity data, obtained from the Second Amendment Supernovae catalogue and the Tully-Fisher catalogues, we find that 2M++ reconstruction provides a better model of the peculiar velocity in the local Universe than the kernel-smoothed peculiar velocity models. We study the impact of peculiar velocities on the measurement of $H_0$ from gravitational waves and megamasers. In doing so, we introduce a probabilistic framework to marginalize over the peculiar velocity corrections along the line-of-sight. For the megamasers, we find $H_0 = 69^{+2.9}_{-2.8}$ km s^{-1} Mpc^{-1} using the 2M++ velocity field. We also study the peculiar velocity of the the galaxy NGC1052-DF2, concluding that a short $\sim$ 13 Mpc distance is not a likely explanation of the anomalously low dark matter fraction of that galaxy. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: Submitted to MNRAS

Taking measurements of the kinematic Sunyaev-Zel'dovich effect forward: including uncertainties from velocity reconstruction with forward modeling

Authors: Nhat-Minh Nguyen, Jens Jasche, Guilhem Lavaux, Fabian Schmidt

Abstract: We measure the kinematic Sunyaev-Zel'dovich (kSZ) effect, imprinted by maxBCG clusters, on the Planck SMICA map of the Cosmic Microwave Background (CMB). Our measurement, for the first time, directly accounts for uncertainties in the velocity reconstruction step through the process of Bayesian forward modeling. We show that this often neglected uncertainty budget typically increases the final unce… ▽ More We measure the kinematic Sunyaev-Zel'dovich (kSZ) effect, imprinted by maxBCG clusters, on the Planck SMICA map of the Cosmic Microwave Background (CMB). Our measurement, for the first time, directly accounts for uncertainties in the velocity reconstruction step through the process of Bayesian forward modeling. We show that this often neglected uncertainty budget typically increases the final uncertainty on the measured kSZ signal amplitude by $\simeq15\%$ at cluster scale. We observe evidence for the kSZ effect, at a significance of $\simeq2σ$. Our analysis, when applied to future higher-resolution CMB data, together with minor improvements in map-filtering and signal-modeling methods, should yield both significant and unbiased measurements of the kSZ signal, which can then be used to probe and constrain baryonic content of galaxy clusters and galaxy groups. △ Less

Submitted 6 October, 2020; v1 submitted 27 July, 2020; originally announced July 2020.

Comments: 29 pages, 15 figures, 2 tables; main results in Figs. 4, 8, 10; comments welcome. v2: Add arguments and a simplified treatment for photo-z clusters (sec. 3.3), update numerical results accordingly, main results and conclusions remain unchanged. v3: Add statements regarding data and code availability, add citations, update acknowledgement, match JCAP version

Journal ref: JCAP12(2020)011

Precision cosmology with voids in the final BOSS data

Authors: Nico Hamaus, Alice Pisani, Jin-Ah Choi, Guilhem Lavaux, Benjamin D. Wandelt, Jochen Weller

Abstract: We report novel cosmological constraints obtained from cosmic voids in the final BOSS DR12 dataset. They arise from the joint analysis of geometric and dynamic distortions of average void shapes (i.e., the stacked void-galaxy cross-correlation function) in redshift space. Our model uses tomographic deprojection to infer real-space void profiles and self-consistently accounts for the Alcock-Paczyns… ▽ More We report novel cosmological constraints obtained from cosmic voids in the final BOSS DR12 dataset. They arise from the joint analysis of geometric and dynamic distortions of average void shapes (i.e., the stacked void-galaxy cross-correlation function) in redshift space. Our model uses tomographic deprojection to infer real-space void profiles and self-consistently accounts for the Alcock-Paczynski (AP) effect and redshift-space distortions (RSD) without any prior assumptions on cosmology or structure formation. It is derived from first physical principles and provides an extremely good description of the data at linear perturbation order. We validate this model with the help of mock catalogs and apply it to the final BOSS data to constrain the RSD and AP parameters $f/b$ and $D_AH/c$, where $f$ is the linear growth rate, $b$ the linear galaxy bias, $D_A$ the comoving angular diameter distance, $H$ the Hubble rate, and $c$ the speed of light. In addition, we include two nuisance parameters in our analysis to marginalize over potential systematics. We obtain $f/b=0.540\pm0.091$ and $D_AH/c=0.588\pm0.004$ from the full void sample at a mean redshift of $z=0.51$. In a flat $Λ$CDM cosmology, this implies $Ω_\mathrm{m}=0.312\pm0.020$ for the present-day matter density parameter. When we use additional information from the survey mocks to calibrate our model, these constraints improve to $f/b=0.347\pm0.023$, $D_AH/c=0.588\pm0.003$, and $Ω_\mathrm{m}=0.310\pm0.017$. However, we emphasize that the calibration depends on the specific model of cosmology and structure formation assumed in the mocks, so the calibrated results should be considered less robust. Nevertheless, our calibration-independent constraints are among the tightest of their kind to date, demonstrating the immense potential of using cosmic voids for cosmology in current and future data. △ Less

Submitted 19 November, 2020; v1 submitted 15 July, 2020; originally announced July 2020.

Comments: 46 pages, 18 figures, 1 table. JCAP accepted version after peer review, typos corrected

Journal ref: Journal of Cosmology and Astroparticle Physics, Issue 12, article id. 023 (2020)

A hierarchical field-level inference approach to reconstruction from sparse Lyman-$α$ forest data

Authors: Natalia Porqueres, Oliver Hahn, Jens Jasche, Guilhem Lavaux

Abstract: We address the problem of inferring the three-dimensional matter distribution from a sparse set of one-dimensional quasar absorption spectra of the Lyman-$α$ forest. Using a Bayesian forward modelling approach, we focus on extending the dynamical model to a fully self-consistent hierarchical field-level prediction of redshift-space quasar absorption sightlines. Our field-level approach rests on a… ▽ More We address the problem of inferring the three-dimensional matter distribution from a sparse set of one-dimensional quasar absorption spectra of the Lyman-$α$ forest. Using a Bayesian forward modelling approach, we focus on extending the dynamical model to a fully self-consistent hierarchical field-level prediction of redshift-space quasar absorption sightlines. Our field-level approach rests on a recently developed semiclassical analogue to Lagrangian perturbation theory (LPT), which improves over noise problems and interpolation requirements of LPT. It furthermore allows for a manifestly conservative mapping of the optical depth to redshift space. In addition, this new dynamical model naturally introduces a coarse-graining scale, which we exploited to accelerate the Markov chain Monte-Carlo (MCMC) sampler using simulated annealing. By gradually reducing the effective temperature of the forward model, we were able to allow it to first converge on large spatial scales before the sampler became sensitive to the increasingly larger space of smaller scales. We demonstrate the advantages, in terms of speed and noise properties, of this field-level approach over using LPT as a forward model, and, using mock data, we validated its performance to reconstruct three-dimensional primordial perturbations and matter distribution from sparse quasar sightlines. △ Less

Submitted 18 August, 2020; v1 submitted 26 May, 2020; originally announced May 2020.

Journal ref: A&A 642, A139 (2020)

Unbiased Cosmology Inference from Biased Tracers using the EFT Likelihood

Authors: Fabian Schmidt, Giovanni Cabass, Jens Jasche, Guilhem Lavaux

Abstract: We present updates on the cosmology inference using the effective field theory (EFT) likelihood presented previously in Schmidt et al., 2018, Elsner et al., 2019 [1,2]. Specifically, we add a cutoff to the initial conditions that serve as starting point for the matter forward model. We show that this cutoff, which was not employed in any previous related work, is important to regularize loop integ… ▽ More We present updates on the cosmology inference using the effective field theory (EFT) likelihood presented previously in Schmidt et al., 2018, Elsner et al., 2019 [1,2]. Specifically, we add a cutoff to the initial conditions that serve as starting point for the matter forward model. We show that this cutoff, which was not employed in any previous related work, is important to regularize loop integrals that otherwise involve small-scale, non-perturbative modes. We then present results on the inferred value of the linear power spectrum normalization $σ_{8}$ from rest-frame halo catalogs using both second- and third-order bias expansions, imposing uniform priors on all bias parameters. Due to the perfect bias-$σ_{8}$ degeneracy at linear order, constraints on $σ_{8}$ rely entirely on nonlinear information. The results show the expected convergence behavior when lowering the cutoff in wavenumber, $Λ$. When including modes up to $k \leq Λ= 0.1\,h\,{\rm Mpc}^{-1}$ in the second-order case, $σ_{8}$ is recovered to within $\lesssim 6\,\%$ for a range of halo masses and redshifts. The systematic bias shrinks to $4\,\%$ or less for the third-order bias expansion on the same range of scales. Together with additional evidence we provide, this shows that the residual mismatch in $σ_{8}$ can be attributed to higher-order bias contributions. We conclude that the EFT likelihood is able to infer unbiased cosmological constraints, within expected theoretical systematic errors, from physical biased tracers on quasilinear scales △ Less

Submitted 5 November, 2020; v1 submitted 14 April, 2020; originally announced April 2020.

Comments: 36 pages, 10 figures, 2 tables; main results in Figs. 2, 4, 6; v2: Sec. 3 clarified, Sec. 5.4 expanded; matches published JCAP version

Journal ref: JCAP 11 (2020), 008

Perfectly parallel cosmological simulations using spatial comoving Lagrangian acceleration

Authors: Florent Leclercq, Baptiste Faure, Guilhem Lavaux, Benjamin D. Wandelt, Andrew H. Jaffe, Alan F. Heavens, Will J. Percival, Camille Noûs

Abstract: Existing cosmological simulation methods lack a high degree of parallelism due to the long-range nature of the gravitational force, which limits the size of simulations that can be run at high resolution. To solve this problem, we propose a new, perfectly parallel approach to simulate cosmic structure formation, which is based on the spatial COmoving Lagrangian Acceleration (sCOLA) framework. Buil… ▽ More Existing cosmological simulation methods lack a high degree of parallelism due to the long-range nature of the gravitational force, which limits the size of simulations that can be run at high resolution. To solve this problem, we propose a new, perfectly parallel approach to simulate cosmic structure formation, which is based on the spatial COmoving Lagrangian Acceleration (sCOLA) framework. Building upon a hybrid analytical and numerical description of particles' trajectories, our algorithm allows for an efficient tiling of a cosmological volume, where the dynamics within each tile is computed independently. As a consequence, the degree of parallelism is equal to the number of tiles. We optimised the accuracy of sCOLA through the use of a buffer region around tiles and of appropriate Dirichlet boundary conditions around sCOLA boxes. As a result, we show that cosmological simulations at the degree of accuracy required for the analysis of the next generation of surveys can be run in drastically reduced wall-clock times and with very low memory requirements. The perfect scalability of our algorithm unlocks profoundly new possibilities for computing larger cosmological simulations at high resolution, taking advantage of a variety of hardware architectures. △ Less

Submitted 16 September, 2022; v1 submitted 10 March, 2020; originally announced March 2020.

Comments: 24 pages, 8 figures, 2 tables. Typos corrected with respect to A&A published version. The code is publicly available at https://bitbucket.org/florent-leclercq/simbelmyne/

Journal ref: A&A 639, A91 (2020)

Setting the Stage: Structures from Gaussian Random Fields

Authors: Till Sawala, Adrian Jenkins, Stuart McAlpine, Jens Jasche, Guilhem Lavaux, Peter H. Johansson, Carlos S. Frenk

Abstract: We study structure formation in a set of cosmological simulations to uncover the scales in the initial density field that gave rise to the formation of present-day structures. Our simulations share a common primordial power spectrum (here Lambda-CDM), but the introduction of hierarchical variations of the phase information allows us to systematically study the scales that determine the formation o… ▽ More We study structure formation in a set of cosmological simulations to uncover the scales in the initial density field that gave rise to the formation of present-day structures. Our simulations share a common primordial power spectrum (here Lambda-CDM), but the introduction of hierarchical variations of the phase information allows us to systematically study the scales that determine the formation of structure at later times. We consider the variance in z=0 statistics such as the matter power spectrum and halo mass function. We also define a criterion for the existence of individual haloes across simulations, and determine what scales in the initial density field contain sufficient information for the non-linear formation of unique haloes. We study how the characteristics of individual haloes such as the mass and concentration, as well as the position and velocity, are affected by variations on different scales, and give scaling relations for haloes of different mass. Finally, we use the example of a cluster-mass halo to show how our hierarchical parametrisation of the initial density field can be used to create variants of particular objects. With properties such as mass, concentration, kinematics and substructure of haloes set on distinct and well-determined scales, and its unique ability to introduce variations localised in real space, our method is a powerful tool to study structure formation in cosmological simulations. △ Less

Submitted 9 March, 2020; originally announced March 2020.

Comments: 18 pages, 22 figures, submitted to MNRAS

Cosmic flows in the nearby Universe: new peculiar velocities from SNe and cosmological constraints

Authors: Supranta S. Boruah, Michael J. Hudson, Guilhem Lavaux

Abstract: The peculiar velocity field offers a unique way to probe dark matter density field on large scales at low redshifts. In this work, we have compiled a new sample of 465 peculiar velocities from low redshift $(z < 0.067)$ Type Ia supernovae. We compare the reconstructed velocity field derived from the 2M++ galaxy redshift compilation to the supernovae, the SFI++ and the 2MTF Tully-Fisher distance ca… ▽ More The peculiar velocity field offers a unique way to probe dark matter density field on large scales at low redshifts. In this work, we have compiled a new sample of 465 peculiar velocities from low redshift $(z < 0.067)$ Type Ia supernovae. We compare the reconstructed velocity field derived from the 2M++ galaxy redshift compilation to the supernovae, the SFI++ and the 2MTF Tully-Fisher distance catalogues. We used a forward method to jointly infer the distances and the velocities of distance indicators by comparing the observations to the reconstruction. Comparison of the reconstructed peculiar velocity fields to observations allows us to infer the cosmological parameter combination $fσ_8$, and the bulk flow velocity arising from outside the survey volume. The residual bulk flow arising from outside the 2M++ volume is inferred to be $171^{+11}_{-11}$ km s$^{-1}$ in the direction $l=301^{\circ} \pm 4^{\circ}$ and $b=0^{\circ} \pm 3^{\circ}$. We obtain $fσ_{8} = 0.400 \pm 0.017$, equivalent to $S_8 \approx σ_8(Ω_m/0.3)^{0.55}=0.776\pm0.033$, which corresponds to an approximately $ 4\%\,$ statistical uncertainty on the value of $fσ_8$. Our inferred value is consistent with other low redshift results in the literature. △ Less

Submitted 2 October, 2020; v1 submitted 19 December, 2019; originally announced December 2019.

Comments: MNRAS accepted version

Velocity correction for Hubble constant measurements from standard sirens

Authors: Suvodip Mukherjee, Guilhem Lavaux, François R. Bouchet, Jens Jasche, Benjamin D. Wandelt, Samaya M. Nissanke, Florent Leclercq, Kenta Hotokezaka

Abstract: Gravitational wave (GW) sources are an excellent probe of the luminosity distance and offer a novel measure of the Hubble constant, $H_0$. This estimation of $H_0$ from standard sirens requires an accurate estimation of the cosmological redshift of the host galaxy of the GW source, after correcting for its peculiar velocity. Absence of an accurate peculiar velocity correction affects both the prec… ▽ More Gravitational wave (GW) sources are an excellent probe of the luminosity distance and offer a novel measure of the Hubble constant, $H_0$. This estimation of $H_0$ from standard sirens requires an accurate estimation of the cosmological redshift of the host galaxy of the GW source, after correcting for its peculiar velocity. Absence of an accurate peculiar velocity correction affects both the precision and accuracy of the measurement of $H_0$, particularly for nearby sources. We propose a framework to incorporate such a peculiar velocity correction for GW sources. A first implementation of our method to the event GW170817 combined with the Very Large Baseline Interferometry (VLBI) observation leads to a revised value of $H_0= 68.3^{+ 4.6}_{-4.5}$ km/s/Mpc. While this revision is minor, it demonstrates that our method makes it possible for obtaining an unbiased and accurate measurements of $H_0$ at the precision required for the standard siren cosmology. △ Less

Submitted 24 November, 2020; v1 submitted 18 September, 2019; originally announced September 2019.

Comments: 12 pages, 10 figures. Matches the version accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 646, A65 (2021)

Systematic-free inference of the cosmic matter density field from SDSS3-BOSS data

Authors: Guilhem Lavaux, Jens Jasche, Florent Leclercq

Abstract: We perform an analysis of the three-dimensional cosmic matter density field traced by galaxies of the SDSS-III/BOSS galaxy sample. The systematic-free nature of this analysis is confirmed by two elements: the successful cross-correlation with the gravitational lensing observations derived from Planck 2018 data and the absence of bias at scales $k \simeq 10^{-3}-10^{-2}h$ Mpc$^{-1}$ in the a poster… ▽ More We perform an analysis of the three-dimensional cosmic matter density field traced by galaxies of the SDSS-III/BOSS galaxy sample. The systematic-free nature of this analysis is confirmed by two elements: the successful cross-correlation with the gravitational lensing observations derived from Planck 2018 data and the absence of bias at scales $k \simeq 10^{-3}-10^{-2}h$ Mpc$^{-1}$ in the a posteriori power spectrum of recovered initial conditions. Our analysis builds upon our algorithm for Bayesian Origin Reconstruction from Galaxies (BORG) and uses a physical model of cosmic structure formation to infer physically meaningful cosmic structures and their corresponding dynamics from deep galaxy observations. Our approach accounts for redshift-space distortions and light-cone effects inherent to deep observations. We also apply detailed corrections to account for known and unknown foreground contaminations, selection effects and galaxy biases. We obtain maps of residual, so far unexplained, systematic effects in the spectroscopic data of SDSS-III/BOSS. Our results show that unbiased and physically plausible models of the cosmic large scale structure can be obtained from present and next-generation galaxy surveys. △ Less

Submitted 13 September, 2019; originally announced September 2019.

Comments: 19 pages, 15 figures, submitted to MNRAS

Neural physical engines for inferring the halo mass distribution function

Authors: Tom Charnock, Guilhem Lavaux, Benjamin D. Wandelt, Supranta Sarma Boruah, Jens Jasche, Michael J. Hudson

Abstract: An ambitious goal in cosmology is to forward-model the observed distribution of galaxies in the nearby Universe today from the initial conditions of large-scale structures. For practical reasons, the spatial resolution at which this can be done is necessarily limited. Consequently, one needs a mapping between the density of dark matter averaged over ~Mpc scales, and the distribution of dark matter… ▽ More An ambitious goal in cosmology is to forward-model the observed distribution of galaxies in the nearby Universe today from the initial conditions of large-scale structures. For practical reasons, the spatial resolution at which this can be done is necessarily limited. Consequently, one needs a mapping between the density of dark matter averaged over ~Mpc scales, and the distribution of dark matter halos (used as a proxy for galaxies) in the same region. Here we demonstrate a method for determining the halo mass distribution function by learning the tracer bias between density fields and halo catalogues using a neural bias model. The method is based on the Bayesian analysis of simple, physically motivated, neural network-like architectures, which we denote as neural physical engines, and neural density estimation. As a result, we are able to sample the initial phases of the dark matter density field whilst inferring the parameters describing the halo mass distribution function, providing a fully Bayesian interpretation of both the initial dark matter density distribution and the neural bias model. We successfully run an upgraded BORG inference using our new likelihood and neural bias model with halo catalogues derived from full N-body simulations. We notice orders of magnitude improvement in modelling compared to classical biasing techniques. △ Less

Submitted 13 September, 2019; originally announced September 2019.

Comments: 12 pages, 5 figures

Inferring high redshift large-scale structure dynamics from the Lyman-alpha forest

Authors: Natalia Porqueres, Jens Jasche, Guilhem Lavaux, Torsten Enßlin

Abstract: One of the major science goals over the coming decade is to test fundamental physics with probes of the cosmic large-scale structure out to high redshift. Here we present a fully Bayesian approach to infer the three-dimensional cosmic matter distribution and its dynamics at $z>2$ from observations of the Lyman-$α$ forest. We demonstrate that the method recovers the unbiased mass distribution and t… ▽ More One of the major science goals over the coming decade is to test fundamental physics with probes of the cosmic large-scale structure out to high redshift. Here we present a fully Bayesian approach to infer the three-dimensional cosmic matter distribution and its dynamics at $z>2$ from observations of the Lyman-$α$ forest. We demonstrate that the method recovers the unbiased mass distribution and the correct matter power spectrum at all scales. Our method infers the three-dimensional density field from a set of one-dimensional spectra, interpolating the information between the lines of sight. We show that our algorithm provides unbiased mass profiles of clusters, becoming an alternative for estimating cluster masses complementary to weak lensing or X-ray observations. The algorithm employs a Hamiltonian Monte Carlo method to generate realizations of initial and evolved density fields and the three-dimensional large-scale flow, revealing the cosmic dynamics at high redshift. The method correctly handles multi-modal parameter distributions, which allow constraining the physics of the intergalactic medium (IGM) with high accuracy. We performed several tests using realistic simulated quasar spectra to test and validate our method. Our results show that detailed and physically plausible inference of three-dimensional large-scale structures at high redshift has become feasible. △ Less

Submitted 17 September, 2019; v1 submitted 5 July, 2019; originally announced July 2019.

Journal ref: A&A 630, A151 (2019)

Wiener filtering and pure E/B decomposition of CMB maps with anisotropic correlated noise

Authors: Doogesh Kodi Ramanah, Guilhem Lavaux, Benjamin D. Wandelt

Abstract: We present an augmented version of our dual messenger algorithm for spin field reconstruction on the sphere, while accounting for highly non-trivial and realistic noise models such as modulated correlated noise. We also describe an optimization method for the estimation of noise covariance from Monte Carlo simulations. Using simulated Planck polarized cosmic microwave background (CMB) maps as a sh… ▽ More We present an augmented version of our dual messenger algorithm for spin field reconstruction on the sphere, while accounting for highly non-trivial and realistic noise models such as modulated correlated noise. We also describe an optimization method for the estimation of noise covariance from Monte Carlo simulations. Using simulated Planck polarized cosmic microwave background (CMB) maps as a showcase, we demonstrate the capabilities of the algorithm in reconstructing pure E and B maps, guaranteed to be free from ambiguous modes resulting from the leakage or coupling issue that plagues conventional methods of E/B separation. Due to its high speed execution, coupled with lenient memory requirements, the algorithm can be optimized in exact global Bayesian analyses of state-of-the-art CMB data for a statistically optimal separation of pure E and B modes. Our algorithm, therefore, has a potentially key role in the data analysis of high-resolution and high-sensitivity CMB data, especially with the range of upcoming CMB experiments tailored for the detection of the elusive primordial B-mode signal. △ Less

Submitted 30 September, 2019; v1 submitted 25 June, 2019; originally announced June 2019.

Comments: 16 pages, 7 figures. Accepted for publication in MNRAS

Cosmology Inference from Biased Tracers using the EFT-based Likelihood

Authors: Franz Elsner, Fabian Schmidt, Jens Jasche, Guilhem Lavaux, Nhat-Minh Nguyen

Abstract: The effective-field-theory (EFT) approach to the clustering of galaxies and other biased tracers allows for an isolation of the cosmological information that is protected by symmetries, in particular the equivalence principle, and thus is robust to the complicated dynamics of dark matter, gas, and stars on small scales. All existing implementations proceed by making predictions for the lowest-orde… ▽ More The effective-field-theory (EFT) approach to the clustering of galaxies and other biased tracers allows for an isolation of the cosmological information that is protected by symmetries, in particular the equivalence principle, and thus is robust to the complicated dynamics of dark matter, gas, and stars on small scales. All existing implementations proceed by making predictions for the lowest-order $n$-point functions of biased tracers, as well as their covariance, and comparing with measurements. Recently, we presented an EFT-based expression for the conditional probability of the density field of a biased tracer given the matter density field, which in principle combines information from arbitrarily high order $n$-point functions. Here, we report results based on this likelihood by applying it to halo catalogs in real space, specifically on the inference of the power spectrum normalization $σ_8$. We include bias terms up to second order as well as the leading higher-derivative term. For a cutoff value of $Λ= 0.1 h\,{\rm Mpc}^{-1}$, we recover the ground-truth value of $σ_8$ to within 95% CL for different halo samples and redshifts. We discuss possible sources for the remaining systematic bias in $σ_8$ as well as future developments. △ Less

Submitted 31 January, 2020; v1 submitted 17 June, 2019; originally announced June 2019.

Comments: 23 pages, 5 figures; v2: some numerical results updated and clarifications added, matches published version

Journal ref: JCAP 01,029 (2020)

Painting halos from cosmic density fields of dark matter with physically motivated neural networks

Authors: Doogesh Kodi Ramanah, Tom Charnock, Guilhem Lavaux

Abstract: We present a novel halo painting network that learns to map approximate 3D dark matter fields to realistic halo distributions. This map is provided via a physically motivated network with which we can learn the non-trivial local relation between dark matter density field and halo distributions without relying on a physical model. Unlike other generative or regressive models, a well motivated prior… ▽ More We present a novel halo painting network that learns to map approximate 3D dark matter fields to realistic halo distributions. This map is provided via a physically motivated network with which we can learn the non-trivial local relation between dark matter density field and halo distributions without relying on a physical model. Unlike other generative or regressive models, a well motivated prior and simple physical principles allow us to train the mapping network quickly and with relatively little data. In learning to paint halo distributions from computationally cheap, analytical and non-linear density fields, we bypass the need for full particle mesh simulations and halo finding algorithms. Furthermore, by design, our halo painting network needs only local patches of dark matter density to predict the halos, and as such, it can predict the 3D halo distribution for any arbitrary simulation box size. Our neural network can be trained using small simulations and used to predict large halo distributions, as long as the resolutions are equivalent. We evaluate our model's ability to generate 3D halo count distributions which reproduce, to a high degree, summary statistics such as the power spectrum and bispectrum, of the input or reference realizations. △ Less

Submitted 25 July, 2019; v1 submitted 25 March, 2019; originally announced March 2019.

Comments: 18 pages, 14 figures. Accepted for publication in PRD

Journal ref: Phys. Rev. D 100, 043515 (2019)

The peculiar velocity field up to $z \sim 0.05$ by forward-modeling Cosmicflows-3 data

Authors: R. Graziani, H. M. Courtois, G. Lavaux, Y. Hoffman, R. B. Tully, Y. Copin, D. Pomarède

Abstract: A hierarchical Bayesian model is applied to the Cosmicflows-3 catalog of galaxy distances in order to derive the peculiar velocity field and distribution of matter within $z \sim 0.054$. The model assumes the $Λ$CDM model within the linear regime and includes the fit of the galaxy distances together with the underlying density field. By forward modeling the data, the method is able to mitigate bia… ▽ More A hierarchical Bayesian model is applied to the Cosmicflows-3 catalog of galaxy distances in order to derive the peculiar velocity field and distribution of matter within $z \sim 0.054$. The model assumes the $Λ$CDM model within the linear regime and includes the fit of the galaxy distances together with the underlying density field. By forward modeling the data, the method is able to mitigate biases inherent to peculiar velocity analyses, such as the Homogeneous Malmquist bias or the log-normal distribution of peculiar velocities. The statistical uncertainty on the recovered velocity field is about 150 km/s depending on the location, and we study systematics coming from the selection function and calibration of distance indicators. The resulting velocity field and related density fields recover the cosmography of the Local Universe which is presented in an unprecedented volume of universe 10 times larger than previously reached. This methodology open the doors to reconstruction of initial conditions for larger and more accurate constrained cosmological simulations. This work is also preparatory to larger peculiar velocity datasets coming from Wallaby, TAIPAN or LSST. △ Less

Submitted 7 January, 2019; originally announced January 2019.

Comments: 15 pages, 14 figures, 1 table. Accepted by MNRAS

Explicit Bayesian treatment of unknown foreground contaminations in galaxy surveys

Authors: Natalia Porqueres, Doogesh Kodi Ramanah, Jens Jasche, Guilhem Lavaux

Abstract: The treatment of unknown foreground contaminations will be one of the major challenges for galaxy clustering analyses of coming decadal surveys. These data contaminations introduce erroneous large-scale effects in recovered power spectra and inferred dark matter density fields. In this work, we present an effective solution to this problem in the form of a robust likelihood designed to account for… ▽ More The treatment of unknown foreground contaminations will be one of the major challenges for galaxy clustering analyses of coming decadal surveys. These data contaminations introduce erroneous large-scale effects in recovered power spectra and inferred dark matter density fields. In this work, we present an effective solution to this problem in the form of a robust likelihood designed to account for effects due to unknown foreground and target contaminations. Conceptually, this robust likelihood marginalizes over the unknown large-scale contamination amplitudes. We showcase the effectiveness of this novel likelihood via an application to a mock SDSS-III data set subject to dust extinction contamination. In order to illustrate the performance of our proposed likelihood, we infer the underlying dark-matter density field and reconstruct the matter power spectrum, being maximally agnostic about the foregrounds. The results are compared to those of an analysis with a standard Poissonian likelihood, as typically used in modern large-scale structure analyses. While the standard Poissonian analysis yields excessive power for large-scale modes and introduces an overall bias in the power spectrum, our likelihood provides unbiased estimates of the matter power spectrum over the entire range of Fourier modes considered in this work. Further, we demonstrate that our approach accurately accounts for and corrects the effects of unknown foreground contaminations when inferring three-dimensional density fields. Robust likelihood approaches, as presented in this work, will be crucial to control unknown systematic error and maximize the outcome of the decadal surveys. △ Less

Submitted 20 March, 2019; v1 submitted 12 December, 2018; originally announced December 2018.

Journal ref: A&A 624, A115 (2019)