Gravitational collapse of cold dark matter leads to infinite-density caustics that seed the primordial dark-matter halos in the large-scale structure. The development of these caustics begins, generically, as an almost one-dimensional phenomenon with the formation of pancakes. Focusing on the one-dimensional case, we identify a landscape of so far unknown singularities in the particle acceleration that emerge after the first crossing of particle trajectories. We complement our fully analytical studies by high-resolution N-body simulations and find outstanding agreement, particularly shortly after the first crossing. We develop the methods in 1D but outline briefly the necessary steps for the 3D case.

Blow-up of solutions for the cosmological fluid equations, often dubbed shell-crossing or orbit crossing, denotes the breakdown of the single-stream regime of the cold-dark-matter fluid. At this instant, the velocity becomes multi-valued and the density singular. Shell-crossing is well understood in one dimension (1D), but not in higher dimensions. This paper is about quasi-one-dimensional (Q1D) flow that depends on all three coordinates but differs only slightly from a strictly 1D flow, thereby allowing a perturbative treatment of shell-crossing using the Euler--Poisson equations written in Lagrangian coordinates. The signature of shell-crossing is then just the vanishing of the Jacobian of the Lagrangian map, a regular perturbation problem. In essence the problem of the first shell-crossing, which is highly singular in Eulerian coordinates, has been desingularized by switching to Lagrangian coordinates, and can then be handled by perturbation theory. Here, all-order recursion relations are obtained for the time-Taylor coefficients of the displacement field, and it is shown that the Taylor series has an infinite radius of convergence. This allows the determination of the time and location of the first shell-crossing, which is generically shown to be taking place earlier than for the unperturbed 1D flow. The time variable used for these statements is not the cosmic time $t$ but the linear growth time $\tau \sim t^{2/3}$. For simplicity, calculations are restricted to an Einstein--de Sitter universe in the Newtonian approximation, and tailored initial data are used. However it is straightforward to relax these limitations, if needed.

It is shown here that in a flat, cold dark matter (CDM) dominated Universe with positive cosmological constant ($\Lambda$), modelled in terms of a Newtonian and collisionless fluid, particle trajectories are analytical in time (representable by a convergent Taylor series) until at least a finite time after decoupling. The time variable used for this statement is the cosmic scale factor, i.e., the "$a$-time", and not the cosmic time. For this, a Lagrangian-coordinates formulation of the Euler-Poisson equations is employed, originally used by Cauchy for 3-D incompressible flow. Temporal analyticity for $\Lambda$CDM is found to be a consequence of novel explicit all-order recursion relations for the $a$-time Taylor coefficients of the Lagrangian displacement field, from which we derive the convergence of the $a$-time Taylor series. A lower bound for the $a$-time where analyticity is guaranteed and shell-crossing is ruled out is obtained, whose value depends only on $\Lambda$ and on the initial spatial smoothness of the density field. The largest time interval is achieved when $\Lambda$ vanishes, i.e., for an Einstein-de Sitter universe. Analyticity holds also if, instead of the $a$-time, one uses the linear structure growth $D$-time, but no simple recursion relations are then obtained. The analyticity result also holds when a curvature term is included in the Friedmann equation for the background, but inclusion of a radiation term arising from the primordial era spoils analyticity.

It is known that the Eulerian and Lagrangian structures of fluid flow can be drastically different; for example, ideal fluid flow can have a trivial (static) Eulerian structure, while displaying chaotic streamlines. Here we show that ideal flow with limited spatial smoothness (an initial vorticity that is just a little better than continuous), nevertheless has time-analytic Lagrangian trajectories before the initial limited smoothness is lost. For proving such results we use a little-known Lagrangian formulation of ideal fluid flow derived by Cauchy in 1815 in a manuscript submitted for a prize of the French Academy. This formulation leads to simple recurrence relations among the time-Taylor coefficients of the Lagrangian map from initial to current fluid particle positions; the coefficients can then be bounded using elementary methods. We first consider various classes of incompressible fluid flow, governed by the Euler equations, and then turn to a case of compressible flow of cosmological relevance, governed by the Euler-Poisson equations.

One of the simplest models used in studying the dynamics of large-scale structure in cosmology, known as the Zeldovich approximation, is equivalent to the three-dimensional inviscid Burgers equation for potential flow. For smooth initial data and sufficiently short times it has the property that the mapping of the positions of fluid particles at any time $t_1$ to their positions at any time $t_2\ge t_1$ is the gradient of a convex potential, a property we call omni-potentiality. Are there other flows with this property, that are not straightforward generalizations of Zeldovich flows? This is answered in the affirmative in both two and three dimensions. How general are such flows? Using a WKB technique we show that in two dimensions, for sufficiently short times, there are omni-potential flows with arbitrary smooth initial velocity. Mappings with a convex potential are known to be associated with the quadratic-cost optimal transport problem. This has important implications for the problem of reconstructing the dynamical history of the Universe from the knowledge of the present mass distribution.

A reconstruction method for recovering the initial conditions of the Universe starting from the present galaxy distribution is presented which guarantees uniqueness of solutions. We show how our method can be used to obtain the peculiar velocities of a large number of galaxies, hence trace galaxies orbits back in time and obtain the entire past dynamical history of the Universe above scales where multi-streaming has not occurred. When tested against a 128^3 LambdaCDM simulation in a box of 200 Mpc/h length, we obtain 60% exact reconstruction on scales above 6 Mpc/h. We apply our method to a real galaxy redshift catalogue, the updated NBG (Nearby Galaxies), containing 1483 galaxies, groups and clusters in a radius of 30 Mpc/h, and reconstruct the peculiar velocity fields in the local neighbourhood. Our reconstructed distances are well-matched to the observed values outside the collapsed regions if Omega_m(t)= 0.20 exp(-0.26 (t-13)) where t is the age of the Universe in Gyrs.

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.

A method for the reconstruction of the primordial density fluctuation field is presented. Various previous approaches to this problem rendered {\it non-unique} solutions. Here, it is demonstrated that the initial positions of dark matter fluid elements, under the hypothesis that their displacement is the gradient of a convex potential, can be reconstructed uniquely. In our approach, the cosmological reconstruction problem is reformulated as an assignment problem in optimisation theory. When tested against numerical simulations, our scheme yields excellent reconstruction on scales larger than a few megaparsecs.

Reconstructing the density fluctuations in the early Universe that evolved into the distribution of galaxies we see today is a challenge of modern cosmology [ref.]. An accurate reconstruction would allow us to test cosmological models by simulating the evolution starting from the reconstructed state and comparing it to the observations. Several reconstruction techniques have been proposed [8 refs.], but they all suffer from lack of uniqueness because the velocities of galaxies are usually not known. Here we show that reconstruction can be reduced to a well-determined problem of optimisation, and present a specific algorithm that provides excellent agreement when tested against data from N-body simulations. By applying our algorithm to the new redshift surveys now under way [ref.], we will be able to recover reliably the properties of the primeval fluctuation field of the local Universe and to determine accurately the peculiar velocities (deviations from the Hubble expansion) and the true positions of many more galaxies than is feasible by any other method. A version of the paper with higher-quality figures is available at http://www.obs-nice.fr/etc7/nature.pdf

This is a review of selected work on the one- and multi-dimensional random Burgers equation (burgulence) with emphasis on questions generally asked for incompressible Navier--Stokes turbulence, such as the law of decay of the energy and the pdf of velocity gradients. Most of the material is devoted to decaying (unforced) burgulence. For more details see the Table of Contents.

Extending work of E, Khanin, Mazel and Sinai (1997 PRL 78:1904-1907) on the one-dimensional Burgers equation, we show that density pdf's have universal power-law tails with exponent -7/2. This behavior stems from singularities, other than shocks, whose nature is quite different in one and several dimensions. We briefly discuss the possibility of detecting singularities of Navier-Stokes turbulence using pdf tails.

We are interested in the tail behavior of the pdf of mass density within the one and $d$-dimensional Burgers/adhesion model used, e.g., to model the formation of large-scale structures in the Universe after baryon-photon decoupling. We show that large densities are localized near ``kurtoparabolic'' singularities residing on space-time manifolds of codimension two ($d \le 2$) or higher ($d \ge 3$). For smooth initial conditions, such singularities are obtained from the convex hull of the Lagrangian potential (the initial velocity potential minus a parabolic term). The singularities contribute {\em \hbox{universal} power-law tails} to the density pdf when the initial conditions are random. In one dimension the singularities are preshocks (nascent shocks), whereas in two and three dimensions they persist in time and correspond to boundaries of shocks; in all cases the corresponding density pdf has the exponent -7/2, originally proposed by E, Khanin, Mazel and Sinai (1997 Phys. Rev. Lett. 78, 1904) for the pdf of velocity gradients in one-dimensional forced Burgers turbulence. We also briefly consider models permitting particle crossings and thus multi-stream solutions, such as the Zel'dovich approximation and the (Jeans)--Vlasov--Poisson equation with single-stream initial data: they have singularities of codimension one, yielding power-law tails with exponent -3.