Cosmological long-wavelength solutions in non-adiabatic… — Plain-Language Explanation

The Big Picture: The Universe as a Busy Kitchen

Imagine the early universe not as a single, smooth soup, but as a busy kitchen with several different chefs working at once. Some chefs are making soup (radiation), others are baking bread (matter), and maybe some are frying eggs (dark matter or other fluids).

Usually, when scientists study the early universe, they pretend all these ingredients are mixed into one perfect, uniform batter. They assume that if you stir the batter, everything moves together perfectly. This is called an adiabatic system (like a smoothie where everything is blended).

However, this paper argues that in the real early universe, the "chefs" didn't always blend perfectly. Sometimes the soup was hot while the bread was cold, or the eggs were overcooked while the soup was underdone. This mismatch is called non-adiabaticity. The paper asks: What happens to the universe's shape and density when these different fluids don't move in perfect sync?

The Problem: The Universe is Too Big to Measure Directly

The scientists are looking at the universe on scales so huge that they are larger than the distance light could have traveled since the Big Bang (called "superhorizon" scales). It's like trying to understand the shape of the entire Earth while standing on a tiny island; you can't see the curvature directly.

To solve this, they use a mathematical trick called a gradient expansion. Imagine looking at a bumpy road. If you stand very close, the bumps look huge. But if you zoom out far enough, the road looks almost flat. The scientists zoom out so far that the "bumps" (density fluctuations) look very gentle. They treat these gentle slopes as a small parameter (a tiny number, $\epsilon$ ) and solve the equations step-by-step, starting with the flattest, simplest version and adding the bumps back in.

The Main Discovery: The "Separate Universes"

The paper uses a framework called ADM formalism (a way of slicing spacetime like a loaf of bread to study it layer by layer). They found that on these giant scales, the universe behaves like a collection of "separate universes."

Imagine a giant field of independent gardens. In each garden, the sun rises and sets, and the plants grow, but they don't talk to each other.

In a single-fluid universe (one type of plant), if you know how one garden is growing, you know how all of them are growing. They are all in sync.
In this multi-fluid universe (different plants), each garden can grow at its own pace. One garden might be full of fast-growing vines (radiation), while another has slow-growing trees (matter). Because they grow at different rates, the "shape" of the garden (the curvature) changes over time in a way that depends on the specific mix of plants in that spot.

The Two Key Ingredients: Adiabatic vs. Entropy

The authors break down the chaos in the kitchen into two types of "noise":

Adiabatic Perturbations (The "Volume" Knob): This is when the whole kitchen gets louder or quieter at the same time. If you turn up the volume, the soup gets louder, the bread gets louder, and the eggs get louder. The ratio between them stays the same. This is the "standard" way the universe expands.
Entropy Perturbations (The "Recipe" Knob): This is when the recipe changes from one spot to another. In one garden, you have too much soup and not enough bread. In another, it's the opposite. The total volume might be the same, but the mix is different. This is called entropy (or isocurvature) perturbation.

The Big Twist: In a universe with only one fluid, the "Recipe Knob" doesn't exist. But in a multi-fluid universe, the "Recipe Knob" is real and powerful. The paper shows that this "Recipe Knob" can actually change the shape of the universe (the curvature) over time, even on the largest scales. This is a surprise because, in simpler models, the shape of the universe was thought to be frozen once it formed.

The "Geodesic Slice": The Observer's View

To make sense of this, the authors had to choose a specific way to watch the universe evolve, which they call the geodesic slice.

Imagine you are a tiny ant walking on a rubber sheet (spacetime). If the sheet stretches, you move with it. This is the "geodesic" view.
The paper shows that if you watch the universe from this specific "ant's-eye view," the "Recipe Knob" (entropy) causes the curvature of the sheet to wiggle and change as the different fluids (radiation vs. matter) take turns dominating the kitchen.

The Demonstration: Matter vs. Radiation

The authors tested their theory with a specific scenario: a universe filled with Radiation (hot, fast-moving particles) and Matter (slower, clumpy stuff).

Early Times: Radiation dominates. The universe acts like it has one fluid. The "Recipe Knob" is barely noticeable.
The Transition: As the universe expands, Radiation cools down faster than Matter. Eventually, Matter takes over.
The Result: During this transition, the "Recipe Knob" spins wildly. The curvature of space (how much the universe bends) changes significantly. It's not constant. The density of matter and radiation fluctuates in complex, non-linear ways that simple math couldn't predict before.

Why This Matters (According to the Paper)

The authors built this mathematical "engine" to create initial conditions for computer simulations.

If you want to simulate how Primordial Black Holes (tiny black holes formed right after the Big Bang) are born, you need to start the simulation with the right "bumps" in the universe.
Previous models assumed the universe was a smooth, single fluid. This paper says, "No, it's a mix of fluids, and the mix matters."
By using their new formulas, scientists can now feed more realistic starting data into their supercomputers to see if these "Recipe Knob" fluctuations are strong enough to crush matter into black holes.

Summary in One Sentence

This paper provides a new mathematical toolkit to describe how a universe made of different "fluids" (like radiation and matter) evolves when they don't move in perfect sync, revealing that the "mix" of these fluids can actively change the shape of space over time, which is crucial for understanding how the first black holes might have formed.

Technical Summary: Cosmological long-wavelength solutions in non-adiabatic multi-fluid systems

Problem Statement
The formation of primordial black holes (PBHs) from large-amplitude density fluctuations in the early universe requires a rigorous treatment of nonlinear gravitational collapse. While linear perturbation theory suffices for small fluctuations, the collapse process is inherently nonlinear, necessitating the solution of the full Einstein equations. Existing long-wavelength approximation methods, particularly the gradient expansion scheme, have been extensively developed for single-fluid or effectively single-fluid systems where the equation of state is barotropic and perturbations are adiabatic. However, the early universe likely underwent transitional epochs (e.g., reheating, QCD phase transitions) involving multiple fluid components with distinct thermodynamic properties. In such non-adiabatic multi-fluid systems, entropy (isocurvature) perturbations play a crucial role. Previous formulations often assumed one component dominated the energy budget or treated the system as effectively adiabatic. There is a lack of explicit, fully nonlinear long-wavelength solutions for general multi-fluid systems where multiple components contribute comparably to the energy density, and where intrinsic non-adiabaticity affects the leading-order behavior of curvature and density perturbations. This gap hinders the construction of consistent physical initial conditions for numerical simulations of PBH formation in multi-component scenarios.

Methodology
The authors develop a formulation of nonlinear cosmological perturbations on superhorizon scales using the Arnowitt–Deser–Misner (ADM) formalism combined with a spatial gradient expansion. The expansion is characterized by a small parameter $\epsilon \equiv k/(aH)$ , where $k$ is the comoving wavenumber, $a(t)$ is the scale factor of a fiducial flat Friedmann-Lemaître-Robertson-Walker (FLRW) spacetime, and $H$ is the Hubble parameter. The analysis assumes $\epsilon \ll 1$ , focusing on scales larger than the Hubble radius.

Key methodological steps include:

ADM Decomposition: The spacetime metric is decomposed into a lapse function, shift vector, and induced metric. The energy-momentum tensor for $N$ collisionless perfect fluids is decomposed into energy density, momentum density, and stress tensor components relative to the normal vector of the spatial hypersurfaces.
Conformal Decomposition: The induced metric is decomposed into a scale factor, a conformal factor $\psi$ , and a traceless part. The extrinsic curvature is similarly decomposed into its trace (mean curvature $K$ ) and a traceless part.
Gradient Expansion: The field equations (Hamiltonian constraint, momentum constraint, and evolution equations) are expanded in powers of $\epsilon$ . The authors derive the leading-order ( $O(\epsilon^0)$ ) equations, assuming the spacetime is locally homogeneous and isotropic in the long-wavelength limit.
Definition of Perturbations: The authors define adiabatic and entropy perturbations in a fully nonlinear manner for multi-fluid systems. An adiabatic perturbation is defined by the existence of a common spatially dependent time shift (or e-folding number shift $\delta N$ ) that maps the perturbed energy densities to the background densities for all fluid components. Entropy perturbations are defined as the deviations from this condition.
Slicing Conditions: The study investigates the behavior of solutions under various slicing conditions, including the uniform e-folding ( $N$ ) slice, constant mean curvature (CMC) slice, and the geodesic slice (where the lapse function $\alpha=1$ ).

Key Contributions and Results

Nonlinear Definition of Adiabatic and Entropy Modes: The paper provides a fully nonlinear definition of adiabatic and entropy perturbations applicable to general multi-fluid cosmologies. It demonstrates that in multi-fluid systems, even at the leading order of the gradient expansion, the system admits both adiabatic and entropy modes due to intrinsic non-adiabaticity. The authors show that the spatial integration functions arising from the conservation equations of individual fluids represent genuine physical degrees of freedom that cannot be removed by gauge transformations.
Explicit Long-Wavelength Solutions: The authors explicitly construct nonlinear long-wavelength solutions for general multi-fluid systems. They express the leading-order energy density of each fluid component in terms of integral functions $C^{(\alpha)}(x^k)$ and the conformal factor $\Psi$ .
Time Dependence of Curvature Perturbations: A critical finding is that in non-adiabatic multi-fluid systems, the curvature perturbation (represented by the conformal factor $\Psi$ or the local e-folding number) is not necessarily conserved on superhorizon scales, even at the leading order. Unlike single-fluid adiabatic systems where $\zeta$ is constant, the curvature perturbation in multi-fluid systems evolves in time, driven by the entropy perturbations.
Non-Uniqueness of Pure Entropy Perturbations: The paper highlights that defining a "pure" entropy perturbation (where the adiabatic mode is set to zero) is non-unique. The authors propose three distinct methods to isolate the adiabatic mode and define pure entropy initial conditions:
1. Primeval Isocurvature Condition: Setting the entropy perturbation to zero at an early time when one component dominates.
2. Zero-Average Condition: Defining the adiabatic mode as the average of the individual fluid e-folding shifts.
3. Initially Isocurvature Condition: Fixing the conformal factor $\Psi=1$ at an initial time $t_0$ .
  The authors demonstrate that these different choices lead to different time evolutions for the curvature perturbation, although the physical density perturbations of individual components remain invariant under these definitions.
Two-Fluid System Analysis: The framework is applied to two-fluid systems. The authors derive explicit solutions showing that while the density perturbations of individual components depend only on the entropy mode (specifically the ratio of the integral functions), the curvature perturbation depends on the specific definition of the adiabatic mode chosen.
Matter-Radiation Demonstration: As a concrete example, the authors analyze a matter-radiation system. They show that for a pure entropy initial condition, the curvature perturbation evolves significantly around the equality epoch (when matter and radiation densities are comparable) but asymptotically approaches a constant value in both the early (radiation-dominated) and late (matter-dominated) limits, where the system behaves effectively as a single adiabatic fluid. The density perturbations are shown to grow nonlinearly and eventually settle into time-independent profiles.

Significance and Claims
The paper claims to establish a consistent theoretical foundation for studying nonlinear perturbations in non-adiabatic multi-component universes. By constructing explicit long-wavelength solutions that account for non-negligible contributions from multiple fluids, the authors provide the necessary physical initial conditions for numerical simulations of PBH formation in such systems. The work extends the conventional $\delta N$ formalism (which relies on the separate universe picture) to scenarios where the equation of state is not barotropic and entropy perturbations are significant at leading order. The authors emphasize that the time evolution of the curvature perturbation in these systems is a direct consequence of intrinsic non-adiabaticity, a feature that distinguishes multi-fluid dynamics from single-fluid adiabatic scenarios. The solutions derived are intended to be used as initial data for numerical relativity simulations to investigate the threshold and dynamics of PBH formation in the early universe.

Cosmological long-wavelength solutions in non-adiabatic multi-fluid systems