Reconstruction of the early Universe as a convex optimization problem

Abstract

We show that the deterministic past history of the Universe can be uniquely reconstructed from the knowledge of the present mass density field, the latter being inferred from the 3D distribution of luminous matter, assumed to be tracing the distribution of dark matter up to a known bias. Reconstruction ceases to be unique below those scales -- a few Mpc -- where multi-streaming becomes significant. Above 6 Mpc/h we propose and implement an effective Monge-Ampere-Kantorovich method of unique reconstruction. At such scales the Zel'dovich approximation is well satisfied and reconstruction becomes an instance of optimal mass transportation, a problem which goes back to Monge (1781). After discretization into N point masses one obtains an assignment problem that can be handled by effective algorithms with not more than cubic time complexity in N and reasonable CPU time requirements. Testing against N-body cosmological simulations gives over 60% of exactly reconstructed points. We apply several interrelated tools from optimization theory that were not used in cosmological reconstruction before, such as the Monge-Ampere equation, its relation to the mass transportation problem, the Kantorovich duality and the auction algorithm for optimal assignment. Self-contained discussion of relevant notions and techniques is provided.

26 pages, 14 figures; accepted to MNRAS. Version 2: numerous minour clarifications in the text, additional material on the history of the Monge-Ampere equation, improved description of the auction algorithm, updated bibliography. Version 3: several misprints corrected

References in corpus (24)

• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Determination of Cosmological Parameters
• Large-Scale Structure of the Universe and Cosmological Perturbation Theory
• Hydra: An Adaptive--Mesh Implementation of PPPM--SPH
• Precision Cosmology? Not Just Yet
• A reconstruction of the initial conditions of the Universe by optimal mass transportation
• Lagrangian Dynamics in Non-Flat Universes and Non-Linear Gravitational Evolution
• Non-Linear Approximations to Gravitational Instability: A Comparison in the Quasi-Linear Regime
• On the least action principle in cosmology
• Reconstruction of cosmological density and velocity fields in the Lagrangian Zel'dovich Approximation
• Simulating Our Cosmological Neighborhood: Mock Catalogs for Velocity Analysis
• Asymptotic behavior of a stratified perturbation in a three dimensional expanding Universe
• Kicked Burgers Turbulence
• Reconstruction of cosmological initial conditions from galaxy redshift catalogues
• Peculiar Velocity Reconstruction with Fast Action Method: Tests on Mock Redshift Surveys
• Eulerian Perturbation Theory in Non-Flat Universes: Second-Order Approximation
• Reconstruction of the primordial Universe by a Monge--Ampere--Kantorovich optimisation scheme
• Filaments and Pancakes in the IRAS 1.2 Jy Redshift Catalogue
• Recovering the Primordial Density Fluctuations: A comparison of methods
• Reconstruction Analysis of Galaxy Redshift Surveys: A Hybrid Reconstruction Method
• The Principle of Least Action and Clustering in Cosmology
• The Inverse Redshift-Space Operator: Reconstructing Cosmological Density and Velocity Fields
• Modeling Gravitational Clustering without Computing Gravitational Force
• The Least-Action Principle: Theory of Cosmological Solutions and the Radial-Velocity Action
• The Inverse Problem for the Euler-Poisson system in Cosmology

Cited by in corpus (96)

• Improving Cosmological Distance Measurements by Reconstruction of the Baryon Acoustic Peak
• Evolution of the cosmic web
• Burgers Turbulence
• Bayesian physical reconstruction of initial conditions from large scale structure surveys
• Optimal protocols and optimal transport in stochastic thermodynamics
• Refined Second Law of Thermodynamics for fast random processes
• Cosmic flow from 2MASS redshift survey: The origin of CMB dipole and implications for LCDM cosmology
• Information field theory for cosmological perturbation reconstruction and non-linear signal analysis
• ELUCID - Exploring the Local Universe with reConstructed Initial Density field I: Hamiltonian Markov Chain Monte Carlo Method with Particle Mesh Dynamics
• Precision cosmology with voids: definition, methods, dynamics
• Physical Bayesian modelling of the non-linear matter distribution: new insights into the Nearby Universe
• ELUCID - Exploring the Local Universe with reConstructed Initial Density field III: Constrained Simulation in the SDSS Volume
• Iterative initial condition reconstruction
• Sparse sampling, galaxy bias, and voids
• Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors
• Adhesive Gravitational Clustering
• Time-analyticity of Lagrangian particle trajectories in ideal fluid flow
• Reconstructing cosmological initial conditions from galaxy peculiar velocities. I. Reverse Zeldovich Approximation
• Reconstruction of primordial density fields
• Cauchy's almost forgotten Lagrangian formulation of the Euler equation for 3D incompressible flow
• Simulating Structure Formation of the Local Universe
• Rejuvenating Power Spectra II: the Gaussianized galaxy density field
• The Maslov dequantization, idempotent and tropical mathematics: A brief introduction
• Cosmic Structure and Dynamics of the Local Universe
• An optimal transport formulation of the Einstein equations of general relativity
• COSMIC BIRTH: Efficient Bayesian Inference of the Evolving Cosmic Web from Galaxy Surveys
• Matter power spectrum from a Lagrangian-space regularization of perturbation theory
• Impact of shell crossing and scope of perturbative approaches in real and redshift space
• How smooth are particle trajectories in a $\Lambda$CDM Universe?
• New method for initial density reconstruction
• Observational biases in Lagrangian reconstructions of cosmic velocity fields
• Isobaric Reconstruction of the Baryonic Acoustic Oscillation
• Higher-order initial conditions with massive neutrinos
• Precision constrained simulation of the Local Universe
• Cosmic voids detection without density measurements
• Shell-crossing in quasi-one-dimensional flow
• Causality for nonlocal phenomena
• A fast semi-discrete optimal transport algorithm for a unique reconstruction of the early Universe
• Reconstructing Cosmological Initial Conditions from Late-Time Structure with Convolutional Neural Networks
• BAO reconstruction: a swift numerical action method for massive spectroscopic surveys
• Cosmological perturbations for two cold fluids in $\Lambda$CDM
• The Monge-Ampere equation: various forms and numerical methods
• Euclid: The search for primordial features
• On nonlinear partial differential equations with an infinite-dimensional conditional symmetry
• The Cauchy-Lagrangian method for numerical analysis of Euler flow
• A review of Lorentzian synthetic theory of timelike Ricci curvature bounds
• Transport-based analysis, modeling, and learning from signal and data distributions
• Reconstructing the baryon acoustic oscillations using biased tracers
• Cosmological Vlasov-Poisson equations for dark matter: Recent developments and connections to selected plasma problems
• Bayesian cosmic density field inference from redshift space dark matter maps
• Planar distribution of the galaxies in the Local Group: a statistical and dynamical analysis
• Optimal Transport Reconstruction of Baryon Acoustic Oscillations
• The non-perturbative regime of cosmic structure formation
• Accurate Baryon Acoustic Oscillations reconstruction via semi-discrete optimal transport
• ELUCID -- Exploring the Local Universe with the reConstructed Initial Density field. II: Reconstruction diagnostics, applied to numerical halo catalogs
• On the efficiency of heat engines at the micro-scale and below
• Straightening the Density-Displacement Relation with a Logarithmic Transform
• An application of Pontryagin's principle to Brownian particle engineered equilibration
• Iterative removal of redshift space distortions from galaxy clustering
• Biased tracer reconstruction with halo mass information
• The Monge-Ampère-Kantorovich approach to reconstruction in cosmology
• Polish spaces of causal curves
• Quasi-spherical collapse of matter in $\Lambda$CDM
• Two is better than one: joint statistics of density and velocity in concentric spheres as a cosmological probe
• Neural network reconstruction of density and velocity fields from the 2MASS Redshift Survey
• Analytical growth functions for cosmic structures in a $\Lambda$CDM Universe
• Optimal Control of Underdamped Systems: An Analytic Approach
• Solution of the optimal assignment problem by diagonal scaling algorithms
• Nonlinear Reconstruction of the Velocity Field
• One-point probability distribution function from spherical collapse: Early Dark Energy (EDE) vs. $\Lambda$CDM
• BullFrog: Multi-step perturbation theory as a time integrator for cosmological simulations
• On general-relativistic Lagrangian perturbation theory and its non-perturbative generalization
• Portfolio optimization with a prescribed terminal wealth distribution
• Lagrangian reconstruction of cosmic velocity fields
• Zero tension Kardar-Parisi-Zhang equation in (d+1)- Dimensions
• Generally covariant $N$-particle dynamics
• Perturbation-theory informed integrators for cosmological simulations
• Exact Results for Propagators in the Geometrical Adhesion Model
• Non-linear density-velocity dynamics in $f(R)$ gravity from spherical collapse
• Fast and accurate collapse-time predictions for collisionless matter
• A double large deviation principle for monge-ampere gravitation
• Starting Cosmological Simulations from the Big Bang
• The Inverse Problem for the Euler-Poisson system in Cosmology
• A Transportation $L^p$ Distance for Signal Analysis
• Optimal transport with discrete long range mean field interactions
• Partial Optimal Transport for a Constant-Volume Lagrangian Mesh with Free Boundaries
• Optimal transport by omni-potential flow and cosmological reconstruction
• Evolution of Cosmic Voids in the Schrodinger-Poisson Formalism
• Using Local Volume data to constrain Dark Matter dynamics
• SHAM-OT: Rapid Subhalo Abundance Matching with Optimal Transport
• Relaxed solutions for incompressible inviscid flows: A variational and gravitational approximation to the initial value problem
• Causal evolution of probability measures and continuity equation
• Simulating fluids with a computer: Introduction and recent advances
• Large-scale semi-discrete optimal transport with distributed Voronoi diagrams
• Fast Sampling of Cosmological Initial Conditions with Gaussian Neural Posterior Estimation
• Rigidly breaking potential flows and a countable Alexandrov theorem for polytopes