Reconstruction of the early Universe as a convex optimization problem

Y. Brenier, U. Frisch, M. Hénon, G. Loeper, S. Matarrese, R. Mohayaee, A. Sobolevskiǐ

Research output: Contribution to journalArticleResearchpeer-review

71 Citations (Scopus)

Abstract

We show that the deterministic past history of the Universe can be uniquely reconstructed from knowledge of the present mass density field, the latter being inferred from the three-dimensional 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 multistreaming becomes significant. Above 6 h-1 Mpc we propose and implement an effective Monge-Ampère-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. After discretization into N point masses one obtains an assignment problem that can be handled by effective algorithms with not more than O(N3) time complexity and reasonable CPU time requirements. Testing against N-body cosmological simulations gives over 60 per cent of exactly reconstructed points. We apply several interrelated tools from optimization theory that were not used in cosmological reconstruction before, such as the Monge-Ampère equation, its relation to the mass transportation problem, the Kantorovich duality and the auction algorithm for optimal assignment. A self-contained discussion of relevant notions and techniques is provided.

Original languageEnglish
Pages (from-to)501-524
Number of pages24
JournalMonthly Notices of the Royal Astronomical Society
Volume346
Issue number2
DOIs
Publication statusPublished - 1 Dec 2003
Externally publishedYes

Keywords

  • Cosmology: theory
  • Early Universe
  • Hydrodynamics
  • Large-scale structure of Universe

Cite this

Brenier, Y. ; Frisch, U. ; Hénon, M. ; Loeper, G. ; Matarrese, S. ; Mohayaee, R. ; Sobolevskiǐ, A. / Reconstruction of the early Universe as a convex optimization problem. In: Monthly Notices of the Royal Astronomical Society. 2003 ; Vol. 346, No. 2. pp. 501-524.
@article{b8b9927d1d13400589785ef08118ddd4,
title = "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 knowledge of the present mass density field, the latter being inferred from the three-dimensional 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 multistreaming becomes significant. Above 6 h-1 Mpc we propose and implement an effective Monge-Amp{\`e}re-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. After discretization into N point masses one obtains an assignment problem that can be handled by effective algorithms with not more than O(N3) time complexity and reasonable CPU time requirements. Testing against N-body cosmological simulations gives over 60 per cent of exactly reconstructed points. We apply several interrelated tools from optimization theory that were not used in cosmological reconstruction before, such as the Monge-Amp{\`e}re equation, its relation to the mass transportation problem, the Kantorovich duality and the auction algorithm for optimal assignment. A self-contained discussion of relevant notions and techniques is provided.",
keywords = "Cosmology: theory, Early Universe, Hydrodynamics, Large-scale structure of Universe",
author = "Y. Brenier and U. Frisch and M. H{\'e}non and G. Loeper and S. Matarrese and R. Mohayaee and A. Sobolevskiǐ",
year = "2003",
month = "12",
day = "1",
doi = "10.1046/j.1365-2966.2003.07106.x",
language = "English",
volume = "346",
pages = "501--524",
journal = "Monthly Notices of the Royal Astronomical Society",
issn = "0035-8711",
publisher = "Oxford University Press, USA",
number = "2",

}

Brenier, Y, Frisch, U, Hénon, M, Loeper, G, Matarrese, S, Mohayaee, R & Sobolevskiǐ, A 2003, 'Reconstruction of the early Universe as a convex optimization problem', Monthly Notices of the Royal Astronomical Society, vol. 346, no. 2, pp. 501-524. https://doi.org/10.1046/j.1365-2966.2003.07106.x

Reconstruction of the early Universe as a convex optimization problem. / Brenier, Y.; Frisch, U.; Hénon, M.; Loeper, G.; Matarrese, S.; Mohayaee, R.; Sobolevskiǐ, A.

In: Monthly Notices of the Royal Astronomical Society, Vol. 346, No. 2, 01.12.2003, p. 501-524.

Research output: Contribution to journalArticleResearchpeer-review

TY - JOUR

T1 - Reconstruction of the early Universe as a convex optimization problem

AU - Brenier, Y.

AU - Frisch, U.

AU - Hénon, M.

AU - Loeper, G.

AU - Matarrese, S.

AU - Mohayaee, R.

AU - Sobolevskiǐ, A.

PY - 2003/12/1

Y1 - 2003/12/1

N2 - We show that the deterministic past history of the Universe can be uniquely reconstructed from knowledge of the present mass density field, the latter being inferred from the three-dimensional 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 multistreaming becomes significant. Above 6 h-1 Mpc we propose and implement an effective Monge-Ampère-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. After discretization into N point masses one obtains an assignment problem that can be handled by effective algorithms with not more than O(N3) time complexity and reasonable CPU time requirements. Testing against N-body cosmological simulations gives over 60 per cent of exactly reconstructed points. We apply several interrelated tools from optimization theory that were not used in cosmological reconstruction before, such as the Monge-Ampère equation, its relation to the mass transportation problem, the Kantorovich duality and the auction algorithm for optimal assignment. A self-contained discussion of relevant notions and techniques is provided.

AB - We show that the deterministic past history of the Universe can be uniquely reconstructed from knowledge of the present mass density field, the latter being inferred from the three-dimensional 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 multistreaming becomes significant. Above 6 h-1 Mpc we propose and implement an effective Monge-Ampère-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. After discretization into N point masses one obtains an assignment problem that can be handled by effective algorithms with not more than O(N3) time complexity and reasonable CPU time requirements. Testing against N-body cosmological simulations gives over 60 per cent of exactly reconstructed points. We apply several interrelated tools from optimization theory that were not used in cosmological reconstruction before, such as the Monge-Ampère equation, its relation to the mass transportation problem, the Kantorovich duality and the auction algorithm for optimal assignment. A self-contained discussion of relevant notions and techniques is provided.

KW - Cosmology: theory

KW - Early Universe

KW - Hydrodynamics

KW - Large-scale structure of Universe

UR - http://www.scopus.com/inward/record.url?scp=0344585375&partnerID=8YFLogxK

U2 - 10.1046/j.1365-2966.2003.07106.x

DO - 10.1046/j.1365-2966.2003.07106.x

M3 - Article

VL - 346

SP - 501

EP - 524

JO - Monthly Notices of the Royal Astronomical Society

JF - Monthly Notices of the Royal Astronomical Society

SN - 0035-8711

IS - 2

ER -