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.

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.

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.