Post-collapse Lagrangian perturbation theory in three dimensions

This paper introduces a novel three-dimensional post-collapse Lagrangian perturbation theory (PCPT) that overcomes shell-crossing singularities by combining high-order perturbation evolution with one-dimensional caustic analysis to accurately model nonlinear matter dynamics after collapse, as validated against high-resolution Vlasov-Poisson simulations.

The Gist

The Big Picture: When the Universe Gets "Sticky"

Imagine the early Universe as a giant, invisible ocean made of Cold Dark Matter. Unlike water, this "dark ocean" has no friction and no pressure; it's just a swarm of invisible particles drifting through space.

For a long time, these particles flow smoothly, like a calm river. Scientists have excellent maps (mathematical theories) to predict where the river will go. But eventually, the river hits a waterfall. The water crashes, swirls, and piles up on top of itself. In physics, this moment is called Shell Crossing.

Once this happens, the smooth river turns into a chaotic mess of overlapping streams. The old maps stop working because they assumed the water was a single, smooth sheet. They can't handle the chaos of particles crashing into each other and bouncing back.

The Problem: We need a new map that works after the crash, when the dark matter has collapsed into a dense, tangled knot.

The Solution: This paper introduces a new mathematical tool called PCPT (Post-Collapse Perturbation Theory). It's like a "repair kit" for our maps, allowing us to predict what happens to the dark matter after it has already collapsed and started to tangle.

---

The Analogy: The Pancake Drop

To understand how the scientists solved this, imagine you are dropping a giant, flat pancake onto a table.

1. The Fall (Pre-Collapse): As the pancake falls, it stays flat and smooth. We can easily predict its path using standard physics. This is the "Lagrangian Perturbation Theory" (LPT) part. It works perfectly until the pancake hits the table.
2. The Crash (Shell Crossing): The pancake hits the table. It doesn't just stop; it squishes, spreads out, and the edges start to fold over each other. Now, the pancake is no longer a single layer; it's a stack of layers. This is the "Multi-stream" region.
3. The Mess (The Challenge): Standard physics gets confused here. It tries to calculate the path of the pancake as if it were still a single sheet, which leads to nonsense results.

How PCPT Works: The "One-Dimensional Shortcut"

The authors realized that even though the pancake is a 3D object, the way it crashes is very specific. It crashes primarily in one direction (straight down), while spreading out sideways.

Think of it like a stack of papers being slammed onto a desk.

• The Crash: The papers compress vertically (up and down).
• The Spread: They slide sideways, but they don't really fold over each other sideways; they just get wider.

The scientists used a clever trick:

1. Flatten the Problem: They realized that because the pancake is so thin vertically, they could treat the messy "stack of papers" problem as if it were a 1D problem (just up and down).
2. Use the Old Rules: They took a known solution for 1D pancake crashes (which they had already figured out in previous papers) and applied it to this 3D situation.
3. The "Back-Reaction": When the pancake hits the table, the table pushes back. This push changes how the pancake moves. The new theory calculates this "push back" (gravitational backreaction) and adds it to the original smooth path.

In simple terms: They didn't try to solve the whole 3D mess at once. They said, "Let's pretend this is a 1D pancake crash, solve that, and then just add a little bit of 3D flavor to the sides."

The Results: A Better Map

The team tested their new "PCPT" map against super-computer simulations (which act like a virtual universe where they can watch the dark matter crash in real-time).

• Old Map (Standard LPT): After the crash, the old map went haywire. It predicted the pancake would fly off in impossible directions or disappear.
• New Map (PCPT): The new map stayed accurate. It correctly predicted how the pancake would fold, how dense the center would get, and how the particles would swirl around.

The Takeaway:
This paper gives us the first fully 3D mathematical tool that can accurately describe the chaotic moments immediately after dark matter collapses. It bridges the gap between the smooth, predictable early Universe and the messy, clumpy Universe we see today.

Why Should You Care?

Dark matter is the invisible skeleton that holds galaxies together. To understand how galaxies form, we need to understand how dark matter collapses.

• Before this paper: We had great theories for the "before" and "after," but a huge gap in the middle where the crash happens.
• Now: We have a better understanding of that crash zone. This helps cosmologists understand how the first "seeds" of galaxies (proto-halos) form, which eventually grow into the massive galaxies like our own Milky Way.

Summary Metaphor:
If the Universe is a giant dance, standard theories tell us how the dancers move in a smooth line. When they crash into a pile, the old theories say, "We don't know what happens next!" This paper says, "Actually, we do. Even in a pile, the dancers follow a specific rhythm, and here is the sheet music for that rhythm."

Technical Summary

1. Problem Statement

The gravitational collapse of collisionless Cold Dark Matter (CDM) leads to shell-crossing singularities, where matter trajectories intersect, creating multi-stream regions with non-zero velocity dispersion.

• Limitation of Standard Theories: Standard Perturbation Theory (SPT) and standard Lagrangian Perturbation Theory (LPT) rely on single-stream assumptions. Once shell-crossing occurs, their perturbative expansions diverge rapidly, losing predictive power in the highly nonlinear, multi-stream regime.
• Limitation of Existing PCPT: Previous Post-Collapse Perturbation Theory (PCPT) formalisms were restricted to one-dimensional (1D) dynamics (infinite massive sheets). Extending this to three dimensions (3D) is challenging due to the complex geometry of multi-stream regions and the fact that shell-crossing occurs at different times along different axes.
• Goal: To develop a fully perturbative, analytical 3D framework that accurately describes the dynamics of CDM immediately after the first shell-crossing, bridging the gap between the pre-collapse Zel'dovich flow and the fully nonlinear multi-stream regime.

2. Methodology

The authors extend the 1D PCPT formalism to 3D by focusing on symmetric pancake collapse (the formation of dark matter proto-halos). The methodology relies on three core ingredients:

A. System Setup and Background Flow

• Symmetry Assumption: The system is modeled as an axisymmetric proto-pancake collapsing primarily along one axis (the $x$-axis). The shell-crossing point is set at the origin ($q_{sc}=0$).
• Taylor Expansion: The background flow (pre-collapse motion) is expanded in Lagrangian coordinates ($q$) around the shell-crossing point up to cubic order.
  • The collapse direction ($x$) includes cubic terms ($q_x^3$) to capture the topology of the caustic.
  • Transverse directions ($y, z$) are treated with linear terms, justified by the "thinness" of the pancake immediately after collapse.
• Ballistic Approximation: The background flow is approximated as ballistic (linear in time) shortly after the shell-crossing time $\tau_{sc}$. This simplifies the integration of corrections.

B. Dimensional Reduction of the Poisson Equation

• Asymptotic Reduction: Because the pancake is extremely thin along the collapse axis ($x$) compared to its transverse extent ($y, z$) immediately after shell-crossing, the 3D Poisson equation is asymptotically reduced to a 1D Poisson equation along the $x$-axis.
• Green's Function: The gravitational force in the multi-stream region is calculated using the 1D Green's function method. The transverse coordinates ($q_y, q_z$) act as parameters in this 1D reduction.
• Force Decomposition: The total gravitational force is split into a background term and a perturbation term arising from the multi-stream density contrast. The authors derive explicit analytical expressions for the force field ($F_x$) inside the multi-stream region (caustic) and in the single-stream region outside it.

C. Perturbative Integration (PCPT)

• Correction Calculation: The derived gravitational backreaction force is integrated perturbatively to correct the background displacement and velocity fields.
• Domain Partitioning: The integration is split into three Lagrangian domains based on the caustic structure:
  1. Outer Multi-stream: Particles that have crossed the shell but are outside the inner caustic ($Q_c < |q_x| \leq \hat{Q}_c$).
  2. Inner Multi-stream: Particles deep within the caustic ($|q_x| \leq Q_c$).
  3. Single-stream: Particles that have not yet crossed.
• Analytical Solutions: The authors derive closed-form analytical expressions for the position ($x$) and velocity ($u$) corrections ($\delta x, \delta u$) for these regions, valid asymptotically shortly after $\tau_{sc}$.

3. Key Contributions

1. First 3D PCPT Formalism: This is the first fully perturbative approach in 3D that captures the highly nonlinear nature of matter evolution after the first shell-crossing without resorting to numerical simulations for the core dynamics.
2. Analytical Backreaction: The paper provides explicit analytical formulas for the gravitational backreaction of multi-stream flows on the background flow in 3D, generalizing previous 1D results.
3. Caustic Structure Modeling: The method accurately models the ellipsoidal caustic structure in Lagrangian space and the resulting three-stream flow topology.
4. Convergence Improvement: The framework demonstrates that including the multi-stream backreaction significantly improves the convergence of the perturbative series compared to standard high-order LPT in the multi-stream regime.

4. Results and Validation

The theoretical predictions were validated against high-resolution Vlasov-Poisson simulations using the ColDICE code (which tracks the phase-space sheet directly).

• Initial Conditions: Tests were performed on systems seeded with crossed sine waves (2-sine and 3-sine) in both Quasi-1D (Q1D) and Anisotropic (ANI) configurations.
• Phase-Space Accuracy:
  • PCPT accurately reproduces the Lagrangian phase-space slices ($x$ vs. $v_x$) and the folding structure of the caustics, even well after the first shell-crossing.
  • It qualitatively captures the second shell-crossing structure along the collapse axis.
  • In contrast, standard 15th-order LPT (15LPT) fails to reproduce the correct phase-space structure, showing significant deviations due to the lack of backreaction.
• Density Profiles: PCPT correctly predicts the density enhancement and the narrowing of the pancake along the collapse axis, matching simulation results much better than standard LPT.
• Limitations Observed:
  • Transverse Directions: While the collapse axis is well-described, the transverse ($y, z$) motions are still driven by the uncorrected LPT background. Consequently, PCPT (like LPT) predicts shell-crossing in the transverse directions slightly earlier than simulations, as the transverse acceleration is overestimated at later times.
  • Boundary Effects: In Q1D cases with large transverse amplitudes, the Taylor expansion and ballistic approximations break down near the periodic box boundaries, leading to deviations.

5. Significance and Future Directions

• Bridging Regimes: PCPT provides a crucial analytical link between the linear/pre-collapse regime and the fully nonlinear, self-similar regime of dark matter halo formation.
• Computational Efficiency: As an analytical framework, it offers a computationally cheap alternative to full N-body or Vlasov simulations for studying early post-collapse dynamics and proto-halo profiles.
• Foundation for Future Work: The authors identify that future improvements should focus on:
  • Extending the formalism to non-symmetric and random Gaussian initial conditions (more realistic for the universe).
  • Improving the treatment of the background flow (e.g., using UV completion or renormalization group techniques) to handle the breakdown of LPT convergence.
  • Iteratively solving for higher-order corrections in time and space to extend the validity of the theory to later times and larger transverse scales.

In summary, this paper successfully generalizes post-collapse perturbation theory to three dimensions, offering a robust, analytical tool to describe the complex dynamics of dark matter immediately following shell-crossing, a regime where standard perturbation theories fail.