Showing 1–50 of 84 results for author: Jasche, J

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

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)

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.

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

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

Constraints on quantum gravity and the photon mass from gamma ray bursts

Authors: Deaglan J. Bartlett, Harry Desmond, Pedro G. Ferreira, Jens Jasche

Abstract: Lorentz invariance violation in quantum gravity (QG) models or a nonzero photon mass, $m_γ$, would lead to an energy-dependent propagation speed for photons, such that photons of different energies from a distant source would arrive at different times, even if they were emitted simultaneously. By developing source-by-source, Monte Carlo-based forward models for such time delays from gamma ray burs… ▽ More Lorentz invariance violation in quantum gravity (QG) models or a nonzero photon mass, $m_γ$, would lead to an energy-dependent propagation speed for photons, such that photons of different energies from a distant source would arrive at different times, even if they were emitted simultaneously. By developing source-by-source, Monte Carlo-based forward models for such time delays from gamma ray bursts, and marginalising over empirical noise models describing other contributions to the time delay, we derive constraints on $m_γ$ and the QG length scale, $\ell_{\rm QG}$, using spectral lag data from the BATSE satellite. We find $m_γ< 4.0 \times 10^{-5} \, h \, {\rm eV}/c^2$ and $\ell_{\rm QG} < 5.3 \times 10^{-18} \, h \, {\rm \, GeV^{-1}}$ at 95% confidence, and demonstrate that these constraints are robust to the choice of noise model. The QG constraint is among the tightest from studies which consider multiple gamma ray bursts and the constraint on $m_γ$, although weaker than from using radio data, provides an independent constraint which is less sensitive to the effects of dispersion by electrons. △ Less

Submitted 18 November, 2021; v1 submitted 16 September, 2021; originally announced September 2021.

Comments: 6 pages, 1 figure, 1 table, accepted for publication in Physical Review D

Journal ref: Phys. Rev. D 104, 103516 (2021)

The large-scale environment of thermonuclear and core-collapse supernovae

Authors: Eleni Tsaprazi, Jens Jasche, Ariel Goobar, Hiranya V. Peiris, Igor Andreoni, Michael W. Coughlin, Christoffer U. Fremling, Matthew J. Graham, Mansi Kasliwal, Shri R. Kulkarni, Ashish A. Mahabal, Reed Riddle, Jesper Sollerman, Anastasios Tzanidakis

Abstract: The new generation of wide-field time-domain surveys has made it feasible to study the clustering of supernova (SN) host galaxies in the large-scale structure (LSS) for the first time. We investigate the LSS environment of SN populations, using 106 dark matter density realisations with a resolution of $\sim$ 3.8 Mpc, constrained by the 2M++ galaxy survey. We limit our analysis to redshift… ▽ More The new generation of wide-field time-domain surveys has made it feasible to study the clustering of supernova (SN) host galaxies in the large-scale structure (LSS) for the first time. We investigate the LSS environment of SN populations, using 106 dark matter density realisations with a resolution of $\sim$ 3.8 Mpc, constrained by the 2M++ galaxy survey. We limit our analysis to redshift $z<0.036$, using samples of 498 thermonuclear and 782 core-collapse SNe from the Zwicky Transient Facility's Bright Transient Survey and Census of the Local Universe catalogues. We detect clustering of SNe with high significance; the observed clustering of the two SNe populations is consistent with each other. Further, the clustering of SN hosts is consistent with that of the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey (BOSS) DR12 spectroscopic galaxy sample in the same redshift range. Using a tidal shear classifier, we classify the LSS into voids, sheets, filaments and knots. We find that both SNe and SDSS galaxies are predominantly found in sheets and filaments. SNe are significantly under-represented in voids and over-represented in knots compared to the volume fraction in these structures. This work opens the potential for using forthcoming wide-field deep SN surveys as a complementary LSS probe. △ Less

Submitted 1 December, 2021; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: 8 pages, 4 figures, accepted by MNRAS

Quantifying the Rarity of the Local Super-Volume

Authors: Stephen Stopyra, Hiranya V. Peiris, Andrew Pontzen, Jens Jasche, Priyamvada Natarajan

Abstract: We investigate the extent to which the number of clusters of mass exceeding $10^{15}\,M_{\odot}\,h^{-1}$ within the local super-volume ($<135\mathrm{\,Mpc}h^{-1}$) is compatible with the standard $Λ$CDM cosmological model. Depending on the mass estimator used, we find that the observed number $N$ of such massive structures can vary between $0$ and $5$. Adopting $N=5$ yields $Λ$CDM likelihoods as l… ▽ More We investigate the extent to which the number of clusters of mass exceeding $10^{15}\,M_{\odot}\,h^{-1}$ within the local super-volume ($<135\mathrm{\,Mpc}h^{-1}$) is compatible with the standard $Λ$CDM cosmological model. Depending on the mass estimator used, we find that the observed number $N$ of such massive structures can vary between $0$ and $5$. Adopting $N=5$ yields $Λ$CDM likelihoods as low as $2.4\times 10^{-3}$ (with $σ_8=0.81$) or $3.8\times 10^{-5}$ (with $σ_8=0.74$). However, at the other extreme ($N=0$), the likelihood is of order unity. Thus, while potentially very powerful, this method is currently limited by systematic uncertainties in cluster mass estimates. This motivates efforts to reduce these systematics with additional observations and improved modelling. △ Less

Submitted 14 July, 2021; originally announced July 2021.

Comments: 7 pages, 3 figures, 1 table

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)

Constraints on Equivalence Principle Violation from Gamma Ray Bursts

Authors: Deaglan J. Bartlett, Dexter Bergsdal, Harry Desmond, Pedro G. Ferreira, Jens Jasche

Abstract: Theories of gravity that obey the Weak Equivalence Principle have the same Parametrised Post-Newtonian parameter $γ$ for all particles at all energies. The large Shapiro time delays of extragalactic sources allow us to put tight constraints on differences in $γ$ between photons of different frequencies from spectral lag data, since a non-zero $Δγ$ would result in a frequency-dependent arrival time… ▽ More Theories of gravity that obey the Weak Equivalence Principle have the same Parametrised Post-Newtonian parameter $γ$ for all particles at all energies. The large Shapiro time delays of extragalactic sources allow us to put tight constraints on differences in $γ$ between photons of different frequencies from spectral lag data, since a non-zero $Δγ$ would result in a frequency-dependent arrival time. The majority of previous constraints have assumed that the Shapiro time delay is dominated by a few local massive objects, although this is a poor approximation for distant sources. In this work we consider the cosmological context of these sources by developing a source-by-source, Monte Carlo-based forward model for the Shapiro time delays by combining constrained realisations of the local density field using the Bayesian origin reconstruction from galaxies algorithm with unconstrained large-scale modes. Propagating uncertainties in the density field reconstruction and marginalising over an empirical model describing other contributions to the time delay, we use spectral lag data of Gamma Ray Bursts from the BATSE satellite to constrain $Δγ< 2.1 \times 10^{-15}$ at $1 σ$ confidence between photon energies of $25 {\rm \, keV}$ and $325 {\rm \, keV}$. △ Less

Submitted 5 October, 2021; v1 submitted 29 June, 2021; originally announced June 2021.

Comments: 14 pages, 7 figures, 3 tables, accepted for publication in Physical Review D

Journal ref: Phys. Rev. D 104, 084025 (2021)

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!

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

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

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

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

Targeting Earth: CRPropa learns to aim

Authors: Jens Jasche, Arjen van Vliet, Jörg P. Rachen

Abstract: Realistic predictions for the arrival directions of ultra-high-energy cosmic rays require extensive simulations of UHECR propagation through 3D space, potentially even including cosmological evolution and timing effects. Such 3D or 4D simulations of cosmic-ray propagation suffer from the fact that a relatively small target - the observer sphere - needs to be hit. If particles are ejected in any di… ▽ More Realistic predictions for the arrival directions of ultra-high-energy cosmic rays require extensive simulations of UHECR propagation through 3D space, potentially even including cosmological evolution and timing effects. Such 3D or 4D simulations of cosmic-ray propagation suffer from the fact that a relatively small target - the observer sphere - needs to be hit. If particles are ejected in any direction from the source according to the source emission geometry, such simulations are tremendously inefficient. We present here a targeting mechanism which finds an optimal emission geometry to maximize the number of hits while remaining unbiased in the arrival-direction distribution. This can lead to speedups by many of orders of magnitude, depending on the simulation setup. We present the basic mathematics to produce unbiased results from targeted simulations, demonstrate its effectiveness with the simulation package CRPropa 3 for various propagation scenarios, and discuss prospects to include this mechanism as a standard part of CRPropa in the future. △ Less

Submitted 12 November, 2019; originally announced November 2019.

Comments: Presented at the 36th International Cosmic Ray Conference (ICRC 2019)

Journal ref: PoS(ICRC2019)447

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)

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)

IMAGINE: modeling the Galactic magnetic field

Authors: Marijke Haverkorn, François Boulanger, Torsten Enßlin, Jörg R. Hörandel, Tess Jaffe, Jens Jasche, Jörg P. Rachen, Anvar Shukurov

Abstract: The IMAGINE Consortium aims to bring modeling of the magnetic field of the Milky Way to a next level, by using Bayesian inference. IMAGINE includes an open-source modular software pipeline that optimizes parameters in a user-defined Galactic magnetic field model against various selected observational datasets. Bayesian priors can be added as external probabilistic constraints of the model paramete… ▽ More The IMAGINE Consortium aims to bring modeling of the magnetic field of the Milky Way to a next level, by using Bayesian inference. IMAGINE includes an open-source modular software pipeline that optimizes parameters in a user-defined Galactic magnetic field model against various selected observational datasets. Bayesian priors can be added as external probabilistic constraints of the model parameters. These conference proceedings describe the science goals of the IMAGINE Consortium, the software pipeline and its inputs, viz observational data sets, Galactic magnetic field models, and Bayesian priors. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Accepted to Galaxies Special Issue "The Power of Faraday Tomography", proceedings of workshop of the same name in Miyazaki, Japan, May 2018

Journal ref: Galaxies 2019, 7(1), 17

Primordial power spectrum and cosmology from black-box galaxy surveys

Authors: Florent Leclercq, Wolfgang Enzi, Jens Jasche, Alan Heavens

Abstract: We propose a new, likelihood-free approach to inferring the primordial matter power spectrum and cosmological parameters from arbitrarily complex forward models of galaxy surveys where all relevant statistics can be determined from numerical simulations, i.e. black-boxes. Our approach, which we call simulator expansion for likelihood-free inference (SELFI), builds upon approximate Bayesian computa… ▽ More We propose a new, likelihood-free approach to inferring the primordial matter power spectrum and cosmological parameters from arbitrarily complex forward models of galaxy surveys where all relevant statistics can be determined from numerical simulations, i.e. black-boxes. Our approach, which we call simulator expansion for likelihood-free inference (SELFI), builds upon approximate Bayesian computation using a novel effective likelihood, and upon the linearisation of black-box models around an expansion point. Consequently, we obtain simple "filter equations" for an effective posterior of the primordial power spectrum, and a straightforward scheme for cosmological parameter inference. We demonstrate that the workload is computationally tractable, fixed a priori, and perfectly parallel. As a proof of concept, we apply our framework to a realistic synthetic galaxy survey, with a data model accounting for physical structure formation and incomplete and noisy galaxy observations. In doing so, we show that the use of non-linear numerical models allows the galaxy power spectrum to be safely fitted up to at least $k_\mathrm{max} = 0.5$ $h$/Mpc, outperforming state-of-the-art backward-modelling techniques by a factor of $\sim 5$ in the number of modes used. The result is an unbiased inference of the primordial matter power spectrum across the entire range of scales considered, including a high-fidelity reconstruction of baryon acoustic oscillations. It translates into an unbiased and robust inference of cosmological parameters. Our results pave the path towards easy applications of likelihood-free simulation-based inference in cosmology. We have made our code pySELFI and our data products publicly available at http://pyselfi.florent-leclercq.eu. △ Less

Submitted 8 October, 2019; v1 submitted 26 February, 2019; originally announced February 2019.

Comments: 22 pages, 12 figures. Matches MNRAS published version. The main results are equations (25), (26), and figure 10. The pySELFI code is publicly available on the first author's website (currently http://pyselfi.florent-leclercq.eu) and on GitHub (https://github.com/florent-leclercq/pyselfi/)

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)

Cosmological inference from Bayesian forward modelling of deep galaxy redshift surveys

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

Abstract: We present a large-scale Bayesian inference framework to constrain cosmological parameters using galaxy redshift surveys, via an application of the Alcock-Paczyński (AP) test. Our physical model of the non-linearly evolved density field, as probed by galaxy surveys, employs Lagrangian perturbation theory (LPT) to connect Gaussian initial conditions to the final density field, followed by a coordin… ▽ More We present a large-scale Bayesian inference framework to constrain cosmological parameters using galaxy redshift surveys, via an application of the Alcock-Paczyński (AP) test. Our physical model of the non-linearly evolved density field, as probed by galaxy surveys, employs Lagrangian perturbation theory (LPT) to connect Gaussian initial conditions to the final density field, followed by a coordinate transformation to obtain the redshift space representation for comparison with data. We generate realizations of primordial and present-day matter fluctuations given a set of observations. This hierarchical approach encodes a novel AP test, extracting several orders of magnitude more information from the cosmological expansion compared to classical approaches, to infer cosmological parameters and jointly reconstruct the underlying 3D dark matter density field. The novelty of this AP test lies in constraining the comoving-redshift transformation to infer the appropriate cosmology which yields isotropic correlations of the galaxy density field, with the underlying assumption relying purely on the cosmological principle. Such an AP test does not rely explicitly on modelling the full statistics of the field. We verify in depth via simulations that this renders our test robust to model misspecification. This leads to another crucial advantage, namely that the cosmological parameters exhibit extremely weak dependence on the currently unresolved phenomenon of galaxy bias, thereby circumventing a potentially key limitation. This is consequently among the first methods to extract a large fraction of information from statistics other than that of direct density contrast correlations, without being sensitive to the amplitude of density fluctuations. We perform several statistical efficiency and consistency tests on a mock galaxy catalogue, using the SDSS-III survey as template. △ Less

Submitted 7 November, 2018; v1 submitted 22 August, 2018; originally announced August 2018.

Comments: 11 + 7 pages, 13 figures. Accepted for publication in A&A

Journal ref: A&A 621, A69 (2019)

A rigorous EFT-based forward model for large-scale structure

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

Abstract: Conventional approaches to cosmology inference from galaxy redshift surveys are based on n-point functions, which are under rigorous perturbative control on sufficiently large scales. Here, we present an alternative approach, which employs a likelihood at the level of the galaxy density field. By integrating out small-scale modes based on effective-field theory arguments, we prove that this likeli… ▽ More Conventional approaches to cosmology inference from galaxy redshift surveys are based on n-point functions, which are under rigorous perturbative control on sufficiently large scales. Here, we present an alternative approach, which employs a likelihood at the level of the galaxy density field. By integrating out small-scale modes based on effective-field theory arguments, we prove that this likelihood is under perturbative control if certain specific conditions are met. We further show that the information captured by this likelihood is equivalent to the combination of the next-to-leading order galaxy power spectrum, leading-order bispectrum, and BAO reconstruction. Combined with MCMC sampling and MAP optimization techniques, our results allow for fully Bayesian cosmology inference from large-scale structure that is under perturbative control. We illustrate this via a first demonstration of unbiased cosmology inference from nonlinear large-scale structure using this likelihood. In particular, we show unbiased estimates of the power spectrum normalization $σ_{8}$ from a catalog of simulated dark matter halos, where nonlinear information is crucial in breaking the $b_{1} - σ_{8}$ degeneracy. △ Less

Submitted 10 January, 2019; v1 submitted 6 August, 2018; originally announced August 2018.

Comments: 38 pages, 1 figure; v2: numerous clarifications, matches JCAP accepted version

Journal ref: JCAP 01,042 (2019)

Fifth force constraints from galaxy warps

Authors: Harry Desmond, Pedro G Ferreira, Guilhem Lavaux, Jens Jasche

Abstract: Intra-galaxy signals contain a wealth of information on fundamental physics, both the dark sector and the nature of gravity. While so far largely unexplored, such probes are set to rise dramatically in importance as upcoming surveys provide data of unprecedented quantity and quality on galaxy structure and dynamics. In this paper, we use warping of stellar disks to test the chameleon- or symmetron… ▽ More Intra-galaxy signals contain a wealth of information on fundamental physics, both the dark sector and the nature of gravity. While so far largely unexplored, such probes are set to rise dramatically in importance as upcoming surveys provide data of unprecedented quantity and quality on galaxy structure and dynamics. In this paper, we use warping of stellar disks to test the chameleon- or symmetron-screened fifth forces which generically arise when new fields couple to matter. We take r-band images of mostly late-type galaxies from the Nasa Sloan Atlas and develop an automated algorithm to quantify the degree of U-shaped warping they exhibit. We then forward-model the warp signal as a function of fifth-force strength $ΔG/G_N$ and range $λ_C$, and the gravitational environments and internal properties of the galaxies, including full propagation of the non-Gaussian uncertainties. Convolving this fifth-force likelihood function with a Gaussian describing astrophysical and observational noise and then constraining $ΔG/G_N$ and $λ_C$ by Markov Chain Monte Carlo, we find the overall likelihood to be significant increased ($Δ\log(\mathcal{L}) \simeq 20$) by adding a screened fifth force with $λ_C \simeq 2$ Mpc, $ΔG/G_N \simeq 0.01$. The variation of $Δ\log(\mathcal{L})$ with $λ_C$ is quantitatively as expected from the correlation of the magnitude of the fifth-force field with the force's range, and a similar model without screening achieves no increase in likelihood over the General Relativistic case $ΔG=0$. Although these results are in good agreement with a previous analysis of the same model using offsets between galaxies' stellar and gas mass centroids (Desmond et al. 2018), we caution that the effects of confounding baryonic and dark matter physics must be thoroughly investigated for the results of the inference to be unambiguous. △ Less

Submitted 29 September, 2018; v1 submitted 31 July, 2018; originally announced July 2018.

Comments: 19 pages, 15 figures; minor revision to match PRD published version

Journal ref: Phys. Rev. D 98, 083010 (2018)

Fifth force constraints from the separation of galaxy mass components

Authors: Harry Desmond, Pedro G Ferreira, Guilhem Lavaux, Jens Jasche

Abstract: One of the most common consequences of extensions to the standard models of particle physics or cosmology is the emergence of a fifth force. While generic fifth forces are tightly constrained at Solar System scales and below, they may escape detection by means of a screening mechanism which effectively removes them in dense environments. We constrain the strength $ΔG/G_N$ and range $λ_C$ of a cham… ▽ More One of the most common consequences of extensions to the standard models of particle physics or cosmology is the emergence of a fifth force. While generic fifth forces are tightly constrained at Solar System scales and below, they may escape detection by means of a screening mechanism which effectively removes them in dense environments. We constrain the strength $ΔG/G_N$ and range $λ_C$ of a chameleon- or symmetron-screened fifth force with Yukawa coupling -- as well as an unscreened fifth force with differential coupling to galactic mass components -- by searching for the displacement it predicts between galaxies' stellar and gas mass centroids. Taking data from the Alfalfa HI survey, identifying galaxies' gravitational environments with the maps of Desmond et al. (2018a) and forward-modelling with a Bayesian likelihood framework, we find $6.6σ$ evidence for $ΔG>0$ at $λ_C \simeq 2$ Mpc, with $ΔG/G_N = 0.025$ at maximum-likelihood. A similar fifth-force model without screening gives no increase in likelihood over the case $ΔG = 0$ for any $λ_C$. Although we validate these results by several methods, we do not claim screened modified gravity to provide the only possible explanation for the signal: this would require knowing that "galaxy formation" physics could not be responsible. We show also the results of a more conservative -- though less well motivated -- noise model which yields only upper limits on $ΔG/G_N$, ranging from $\sim10^{-1}$ for $λ_C \simeq 0.5$ Mpc to $\sim \: \text{few} \times 10^{-4}$ at $λ_C \simeq 50$ Mpc. We show how these constraints may be improved by future galaxy surveys and identify the key features of an observational programme for directly constraining fifth forces on galactic scales. This paper provides a complete description of the analysis summarised in Desmond et al. (2018b). △ Less

Submitted 8 September, 2018; v1 submitted 4 July, 2018; originally announced July 2018.

Comments: 29 pages, 16 figures; minor revision to match PRD published version

Journal ref: Phys. Rev. D 98, 064015 (2018)

Physical Bayesian modelling of the non-linear matter distribution: new insights into the Nearby Universe

Authors: Jens Jasche, Guilhem Lavaux

Abstract: Accurate analyses of present and next-generation galaxy surveys require new ways to handle effects of non-linear gravitational structure formation in data. To address these needs we present an extension of our previously developed algorithm for Bayesian Origin Reconstruction from Galaxies to analyse matter clustering at non-linear scales in observations. This is achieved by incorporating a numeric… ▽ More Accurate analyses of present and next-generation galaxy surveys require new ways to handle effects of non-linear gravitational structure formation in data. To address these needs we present an extension of our previously developed algorithm for Bayesian Origin Reconstruction from Galaxies to analyse matter clustering at non-linear scales in observations. This is achieved by incorporating a numerical particle mesh model of structure formation into our Bayesian inference framework. The algorithm simultaneously infers the 3D primordial matter fluctuations from which present non-linear observations formed and provides reconstructions of velocity fields and structure formation histories. The physical forward modelling approach automatically accounts for non-Gaussian features in evolved matter density fields and addresses the redshift space distortion problem associated with peculiar motions of galaxies. Our algorithm employs a hierarchical Bayes approach to jointly account for observational effects, such as galaxy biases, selection effects, and noise. Corresponding parameters are marginalized out via a sophisticated Markov Chain Monte Carlo approach relying on a combination of a multiple block sampling framework and a Hamiltonian Monte Carlo sampler. We demonstrate the performance of the method by applying it to the 2M++ galaxy compilation, tracing the matter distribution of the Nearby Universe. We show accurate and detailed inferences of the 3D non-linear dark matter distribution of the Nearby Universe. As exemplified in the case of the Coma cluster, we provide mass estimates that are compatible with those obtained from weak lensing and X-ray observations. For the first time, we reconstruct the vorticity of the non-linear velocity field from observations. In summary, our method provides plausible and detailed inferences of dark matter and velocity fields of our cosmic neighbourhood. △ Less

Submitted 28 June, 2018; originally announced June 2018.

Comments: 33 pages, 17 figures

Journal ref: A&A 625, A64 (2019)

IMAGINE: A comprehensive view of the interstellar medium, Galactic magnetic fields and cosmic rays

Authors: Francois Boulanger, Torsten Ensslin, Andrew Fletcher, Philipp Girichides, Stefan Hackstein, Marijke Haverkorn, Joerg R. Hoerandel, Tess R. Jaffe, Jens Jasche, Michael Kachelriess, Kumiko Kotera, Christoph Pfrommer, Jorg P. Rachen, Luiz F. S. Rodrigues, Beatriz Ruiz-Granados, Amit Seta, Anvar Shukurov, Gunter Sigl, Theo Steininger, Valentina Vacca, Ellert van der Velden, Arjen van Vliet, Jiaxin Wang

Abstract: In this white paper we introduce the IMAGINE Consortium and its scientific background, goals and structure. Our purpose is to coordinate and facilitate the efforts of a diverse group of researchers in the broad areas of the interstellar medium, Galactic magnetic fields and cosmic rays, and our goal is to develop more comprehensive insights into the structures and roles of interstellar magnetic fie… ▽ More In this white paper we introduce the IMAGINE Consortium and its scientific background, goals and structure. Our purpose is to coordinate and facilitate the efforts of a diverse group of researchers in the broad areas of the interstellar medium, Galactic magnetic fields and cosmic rays, and our goal is to develop more comprehensive insights into the structures and roles of interstellar magnetic fields and their interactions with cosmic rays. To achieve a higher level of self-consistency, depth and rigour can only be achieved by the coordinated efforts of experts in diverse areas of astrophysics involved in observational, theoretical and numerical work. We present our view of the present status of this topic, identify its key unsolved problems and suggest a strategy that will underpin our work. The backbone of the consortium is the Interstellar MAGnetic field INference Engine, a publicly available Bayesian platform that employs robust statistical methods to explore the multi-dimensional likelihood space using any number of modular inputs. It provides an interpretation and modelling framework that has the power and flexibility to include a variety of observational, theoretical and numerical lines of evidence into a self-consistent and comprehensive picture of the thermal and non-thermal interstellar media. An important innovation is that a consistent understanding of the phenomena that are directly or indirectly influenced by the Galactic magnetic field, such as the deflection of ultra-high energy cosmic rays or extragalactic backgrounds, is made an integral part of the modelling. The IMAGINE Consortium, which is informal by nature and open to new participants, hereby presents a methodological framework for the modelling and understanding of Galactic magnetic fields that is available to all communities whose research relies on a state-of-the-art solution to this problem. (Abridged.) △ Less

Submitted 6 September, 2018; v1 submitted 7 May, 2018; originally announced May 2018.

Comments: 56 pages, 15 figures, minor updates to match published JCAP version

Journal ref: JCAP08(2018)049

The primordial magnetic field in our cosmic backyard

Authors: Sebastian Hutschenreuter, Sebastian Dorn, Jens Jasche, Franco Vazza, Daniela Paoletti, Guilhem Lavaux, Torsten A. Enßlin

Abstract: We reconstruct the 3D structure of magnetic fields, which were seeded by density perturbations during the radiation dominated epoch of the Universe and later on were evolved by structure formation. To achieve this goal, we rely on three dimensional initial density fields inferred from the 2M++ galaxy compilation via the Bayesian $\texttt{BORG}$ algorithm. Using those, we estimate the magnetogenesi… ▽ More We reconstruct the 3D structure of magnetic fields, which were seeded by density perturbations during the radiation dominated epoch of the Universe and later on were evolved by structure formation. To achieve this goal, we rely on three dimensional initial density fields inferred from the 2M++ galaxy compilation via the Bayesian $\texttt{BORG}$ algorithm. Using those, we estimate the magnetogenesis by the so called Harrison mechanism. This effect produced magnetic fields exploiting the different photon drag on electrons and ions in vortical motions, which are exited due to second order perturbation effects in the Early Universe. Subsequently we study the evolution of these seed fields through the non-linear cosmic structure formation by virtue of a MHD simulation to obtain a 3D estimate for the structure of this primordial magnetic field component today. At recombination we obtain a reliable lower limit on the large scale magnetic field strength around $10^{-23} \mathrm{G}$, with a power spectrum peaking at about $ 2\, \mathrm{Mpc}^{-1}h$ in comoving scales. At present we expect this evolved primordial field to have strengthts above $\approx 10^{-27}\, \mathrm{G}$ and $\approx 10^{-29}\, \mathrm{G}$ in clusters of galaxies and voids, respectively. We also calculate the corresponding Faraday rotation measure map and show the magnetic field morphology and strength for specific objects of the Local Universe. △ Less

Submitted 18 October, 2019; v1 submitted 7 March, 2018; originally announced March 2018.

Comments: 14 pages, 19 figures. Revised version after peer review. Scientific and computational results are unchanged. This is an author-created, un-copyedited version of an article published in Classical and Quantum Gravity. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. This article is published under a CC BY licence

Journal ref: Class. Quantum Grav. 35 154001 (2018)

The Fifth Force in the Local Cosmic Web

Authors: Harry Desmond, Pedro G Ferreira, Guilhem Lavaux, Jens Jasche

Abstract: Extensions of the standard models of particle physics and cosmology often lead to long-range fifth forces with properties dependent on gravitational environment. Fifth forces on astrophysical scales are best studied in the cosmic web where perturbation theory breaks down. We present constraints on chameleon- and symmetron-screened fifth forces with Yukawa coupling and megaparsec range -- as well a… ▽ More Extensions of the standard models of particle physics and cosmology often lead to long-range fifth forces with properties dependent on gravitational environment. Fifth forces on astrophysical scales are best studied in the cosmic web where perturbation theory breaks down. We present constraints on chameleon- and symmetron-screened fifth forces with Yukawa coupling and megaparsec range -- as well as unscreened fifth forces with differential coupling to galactic mass components -- by searching for the displacements they predict between galaxies' stars and gas. Taking data from the Alfalfa HI survey, identifying galaxies' gravitational environments with the maps of Desmond et al. (2018a) and forward-modelling with a Bayesian likelihood framework, we set upper bounds on fifth-force strength relative to Newtonian gravity from $ΔG/G_N < \text{few} \: \times 10^{-4}$ for range $λ_C = 50$ Mpc, to $ΔG/G_N \lesssim 0.1$ for $λ_C = 500$ kpc. In $f(R)$ gravity this requires $f_{R0} < \text{few} \: \times \: 10^{-8}$. The analogous bounds without screening are $ΔG/G_N < \text{few} \: \times 10^{-4}$ and $ΔG/G_N < \text{few} \times 10^{-3}$. These are the tightest and among the only fifth-force constraints on galaxy scales. We show how our results may be strengthened with future survey data and identify the key features of an observational programme for furthering fifth-force tests beyond the Solar System. △ Less

Submitted 23 November, 2018; v1 submitted 20 February, 2018; originally announced February 2018.

Comments: 6 pages, 2 figures; revised to match MNRASL published version

Inferring Galactic magnetic field model parameters using IMAGINE - An Interstellar MAGnetic field INference Engine

Authors: Theo Steininger, Torsten A. Enßlin, Maksim Greiner, Tess Jaffe, Ellert van der Velden, Jiaxin Wang, Marijke Haverkorn, Jörg R. Hörandel, Jens Jasche, Jörg P. Rachen

Abstract: Context. The Galactic magnetic field (GMF) has a huge impact on the evolution of the Milky Way. Yet currently there exists no standard model for it, as its structure is not fully understood. In the past many parametric GMF models of varying complexity have been developed that all have been fitted to an individual set of observational data complicating comparability. Aims. Our goal is to systematiz… ▽ More Context. The Galactic magnetic field (GMF) has a huge impact on the evolution of the Milky Way. Yet currently there exists no standard model for it, as its structure is not fully understood. In the past many parametric GMF models of varying complexity have been developed that all have been fitted to an individual set of observational data complicating comparability. Aims. Our goal is to systematize parameter inference of GMF models. We want to enable a statistical comparison of different models in the future, allow for simple refitting with respect to newly available data sets and thereby increase the research area's transparency. We aim to make state-of-the-art Bayesian methods easily available and in particular to treat the statistics related to the random components of the GMF correctly. Methods. To achieve our goals, we built IMAGINE, the Interstellar Magnetic Field Inference Engine. It is a modular open source framework for doing inference on generic parametric models of the Galaxy. We combine highly optimized tools and technology such as the MultiNest sampler and the information field theory framework NIFTy in order to leverage existing expertise. Results. We demonstrate the steps needed for robust parameter inference and model comparison. Our results show how important the combination of complementary observables like synchrotron emission and Faraday depth is while building a model and fitting its parameters to data. IMAGINE is open-source software available under the GNU General Public License v3 (GPL-3) at: https://gitlab.mpcdf.mpg.de/ift/IMAGINE △ Less

Submitted 12 January, 2018; originally announced January 2018.

Imprints of the large-scale structure on AGN formation and evolution

Authors: Natàlia Porqueres, Jens Jasche, Torsten A. Enßlin, Guilhem Lavaux

Abstract: Black hole masses are found to correlate with several global properties of their host galaxies, suggesting that black holes and galaxies have an intertwined evolution and that active galactic nuclei (AGN) have a significant impact on galaxy evolution. Since the large-scale environment can also affect AGN, this work studies how their formation and properties depend on the environment. We have used… ▽ More Black hole masses are found to correlate with several global properties of their host galaxies, suggesting that black holes and galaxies have an intertwined evolution and that active galactic nuclei (AGN) have a significant impact on galaxy evolution. Since the large-scale environment can also affect AGN, this work studies how their formation and properties depend on the environment. We have used a reconstructed three-dimensional high-resolution density field obtained from a Bayesian large-scale structure reconstruction method applied to the 2M++ galaxy sample. A web-type classification relying on the shear tensor is used to identify different structures on the cosmic web, defining voids, sheets, filaments, and clusters. We confirm that the environmental density affects the AGN formation and their properties. We found that the AGN abundance is equivalent to the galaxy abundance, indicating that active and inactive galaxies reside in similar dark matter halos. However, occurrence rates are different for each spectral type and accretion rate. These differences are consistent with the AGN evolutionary sequence suggested by previous authors, Seyferts and Transition objects transforming into LINERs (Low-Ionization Nuclear Emission Line Regions), the weaker counterpart of Seyferts. We conclud that AGN properties depend on the environmental density more than on the web-type. More powerful starbursts and younger stellar populations are found in high densities, where interactions and mergers are more likely. AGN hosts show smaller masses in clusters for Seyferts and Transition objects, which might be due to gas stripping. In voids, the AGN population is dominated by the most massive galaxy hosts. △ Less

Submitted 15 January, 2018; v1 submitted 20 October, 2017; originally announced October 2017.

Comments: 10 pages, 8 figures

Journal ref: A&A 612, A31 (2018)