Magrathea-Pathfinder: A 3D adaptive-mesh code for geodesic ray tracing in $N$-body simulations

Abstract

We introduce Magrathea-Pathfinder, a relativistic ray-tracing framework that can reconstruct the past light cone of observers in cosmological simulations. The code directly computes the 3D trajectory of light rays through the null geodesic equations, with the weak-field limit as its only approximation. This approach offers high levels of versatility while removing the need for many of the standard ray-tracing approximations such as plane-parallel, Born, or multiple-lens. Moreover, the use of adaptive integration steps and interpolation strategies based on adaptive-mesh refinement (AMR) grids allows Magrathea-Pathfinder to accurately account for the non-linear regime of structure formation and fully take advantage of the small-scale gravitational clustering. To handle very large N-body simulations, the framework has been designed as a high-performance computing post-processing tool relying on a hybrid parallelization that combines MPI tasks with C++11 std::threads. In this paper, we describe how realistic cosmological observables can be computed from numerical simulation using ray-tracing techniques. We discuss in particular the production of simulated catalogues and sky maps that account for all the observational effects considering first-order metric perturbations (such as peculiar velocities, gravitational potential, integrated Sachs-Wolfe, time-delay, and gravitational lensing). We perform convergence tests of our gravitational lensing algorithms and conduct performance benchmarks of the null geodesic integration procedures. Magrathea-Pathfinder introduces sophisticated ray-tracing tools to make the link between the space of N-body simulations and light-cone observables. This should provide new ways of exploring existing cosmological probes and building new ones beyond standard assumptions in order to prepare for the next generation of large-scale structure surveys.

Fix typos, code available at https://github.com/vreverdy/magrathea-pathfinder

References in corpus (65)

• Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation
• The cosmological simulation code GADGET-2
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• The Complete Light-curve Sample of Spectroscopically Confirmed Type Ia Supernovae from Pan-STARRS1 and Cosmological Constraints from The Combined Pantheon Sample
• Weak Gravitational Lensing
• Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES
• Stable clustering, the halo model and nonlinear cosmological power spectra
• Large Magellanic Cloud Cepheid Standards Provide a 1% Foundation for the Determination of the Hubble Constant and Stronger Evidence for Physics Beyond LambdaCDM
• The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Cosmological Implications from two Decades of Spectroscopic Surveys at the Apache Point observatory
• The DESI Experiment Part I: Science,Targeting, and Survey Design
• H0LiCOW XIII. A 2.4% measurement of $H_0$ from lensed quasars: $5.3\sigma$ tension between early and late-Universe probes
• Euclid Definition Study Report
• Enzo: An Adaptive Mesh Refinement Code for Astrophysics
• The Carnegie-Chicago Hubble Program. VIII. An Independent Determination of the Hubble Constant Based on the Tip of the Red Giant Branch
• Using gravitational-wave standard sirens
• The linear power spectrum of observed source number counts
• The Cosmic Linear Anisotropy Solving System (CLASS) I: Overview
• What galaxy surveys really measure
• New Perspective on Galaxy Clustering as a Cosmological Probe: General Relativistic Effects
• Ray-tracing through the Millennium Simulation: Born corrections and lens-lens coupling in cosmic shear and galaxy-galaxy lensing
• First Cosmology Results using Type Ia Supernovae from the Dark Energy Survey: Constraints on Cosmological Parameters
• The Uchuu Simulations: Data Release 1 and Dark Matter Halo Concentrations
• Fast estimation of polarization power spectra using correlation functions
• The onion universe: all sky lightcone simulations in spherical shells
• MLAPM - a C code for cosmological simulations
• Precision calculations of the cosmic shear power spectrum projection
• The $\text{Abacus}$ Cosmological $N$-body Code
• Asymmetric galaxy correlation functions
• Probability Distribution Functions of Cosmological Lensing: Convergence, Shear, and Magnification
• The Outer Rim Simulation: A Path to Many-Core Supercomputers
• Constraining early and interacting dark energy with gravitational wave standard sirens: the potential of the eLISA mission
• Full-sky map of the ISW and Rees-Sciama effect from Gpc simulations
• CMB weak-lensing beyond the Born approximation: a numerical approach
• GLAMER Part I: A Code for Gravitational Lensing Simulations with Adaptive Mesh Refinement
• Testing the equivalence principle on cosmological scales
• Imprints of relativistic effects on the asymmetry of the halo cross-correlation function: from linear to non-linear scales
• Observable Deviations from Homogeneity in an Inhomogeneous Universe
• Bias and scatter in the Hubble diagram from cosmological large-scale structure
• Modeling relativistic contributions to the halo power spectrum dipole
• On the Bias of the Distance-Redshift Relation from Gravitational Lensing
• EMMA: an AMR cosmological simulation code with radiative transfer
• Average magnification effect of clumping of matter
• LIGER: mock relativistic light-cones from Newtonian simulations
• Cosmic variance limited Baryon Acoustic Oscillations from the DEUS-FUR $\Lambda$CDM simulation
• The Accuracy of Weak Lensing Simulations
• Weak-lensing observables in relativistic N-body simulations
• The relativistic dipole and gravitational redshift on LSS
• Cosmological ensemble and directional averages of observables
• Weak gravitational lensing of finite beams
• Disentangling dark sector models using weak lensing statistics
• Weak lensing in the Horizon-AGN simulation lightcone. Small scale baryonic effects
• DEUS Full Observable ΛCDM Universe Simulation: the numerical challenge
• Probing dark energy models with extreme pairwise velocities of galaxy clusters from the DEUS-FUR simulations
• Methods for studying the accuracy of light propagation in N-body simulations
• Wide-angle redshift-space distortions at quasi-linear scales: cross-correlation functions from Zel'dovich approximation
• Cosmic convergence and shear with extended sources
• Gravitational Lensing with Three-Dimensional Ray Tracing
• Detectability of the gravitational redshift effect from the asymmetric galaxy clustering
• Modelling the asymmetry of the halo cross-correlation function with relativistic effects at quasi-linear scales
• Weak lensing distortions beyond shear
• Theoretical and numerical perspectives on cosmic distance averages
• Shape, shear and flexion II - Quantifying the flexion formalism for extended sources with the ray-bundle method
• The RayGalGroupSims cosmological simulation suite for the study of relativistic effects: an application to lensing-matter clustering statistics
• Ray tracing the integrated Sachs-Wolfe effect through the light cones of the Dark Energy Universe Simulation -- Full Universe Runs
• Impact of lensing magnification on the analysis of galaxy clustering in redshift space

Cited by in corpus (6)

• The RayGalGroupSims cosmological simulation suite for the study of relativistic effects: an application to lensing-matter clustering statistics
• Cosmological distances with general-relativistic ray tracing: framework and comparison to cosmographic predictions
• Impact of lensing magnification on the analysis of galaxy clustering in redshift space
• Relativistic matter bispectrum of cosmic structures on the light cone
• Linear Relativistic Corrections in the Spherical Fourier-Bessel Power Spectrum
• The spherical Fast Multipole Method (sFMM) for Gravitational Lensing Simulation