Effective field theory of a single scalar pion field for large scale structure in the Universe

This paper develops a next-to-leading order effective field theory for large-scale structure in a $\Lambda$CDM universe by treating the matter fluid's velocity potential as a single scalar "pion" field, validating the approach through consistency relations and N-body simulations to extract EFT coefficients and propose new variables for analyzing cosmic structure.

The Gist

Imagine the Universe as a giant, expanding ocean. For decades, cosmologists have tried to map the waves and currents of this cosmic ocean, which is made mostly of invisible "dark matter." The problem is that as the universe ages, these waves crash into each other, swirl, and break in chaotic ways. Trying to predict exactly how they behave using standard math is like trying to forecast the weather by solving equations for every single water molecule—it's too messy and computationally impossible.

This paper proposes a clever new way to look at the problem. Instead of tracking the water (matter), the density of the water, and the pressure all separately, the authors suggest we track just one thing: a "cosmic surfboard" called the Pion Field (or simply, the "pion").

Here is a breakdown of their ideas using everyday analogies:

1. The "Surfboard" Analogy (The Pion Field)

In standard physics, describing the universe's structure is like trying to describe a storm by tracking wind speed, water pressure, and rain intensity separately. It's complicated.

The authors say: "Let's just track the surfboard."

• The Surfboard ($\pi$): This is a single mathematical field that represents the "velocity potential" of the cosmic fluid. Think of it as a map showing where the water wants to go.
• The Magic: If you know the shape of this surfboard, you can mathematically figure out the water speed, the pressure, and how dense the water is. It collapses three complex variables into one simple one.
• Why "Pion"? In particle physics, a "pion" is a particle that appears when a symmetry is broken. The authors realized that the way the universe's fluid moves is mathematically similar to how pions behave. They are the "Goldstone bosons" (a fancy term for the ripples left behind when a perfect symmetry is broken) of the expanding universe.

2. The "Traffic Jam" Problem (Non-Linearity)

In the early universe, matter was spread out smoothly, like cars driving in perfect lanes on a highway. Physics was easy to predict here.

But as the universe aged, gravity pulled matter together. Cars started swerving, merging, and crashing. This is the non-linear regime.

• The Old Way: Standard math breaks down here because the interactions become too complex. It's like trying to predict the exact path of every car in a massive traffic jam.
• The New Way (EFT): The authors use a technique called Effective Field Theory (EFT). Imagine you are a traffic engineer who doesn't care about every individual car. Instead, you treat the traffic jam as a "fluid" with an average speed and a "viscosity" (stickiness).
  • They calculate the smooth, predictable part of the traffic flow.
  • Then, they add "correction terms" to account for the chaos of the traffic jam (the crashes and swerves) without having to simulate every single car.

3. The "Sound and Viscosity" of the Universe

When matter collapses into galaxies, it doesn't just sit there; it heats up and swirls. In their "fluid" model, this creates two new properties:

• Sound Speed: How fast a "ripple" of density moves through the collapsing matter.
• Viscosity: How "sticky" the universe gets.
  • Fun Fact from the paper: When they calculated the "stickiness" of their simulated universe, they found it was roughly as thick as peanut butter or Crisco (shortening). It's not a vacuum; it has a texture!

4. The "Shockwaves"

The authors ran computer simulations (using a code they named PLASTIC) to see how this "surfboard" evolves.

• They found that eventually, the smooth waves on the surfboard get so steep that they break, forming a shockwave.
• In the real universe, this is when galaxies form and the "single-file" flow of matter breaks down into a chaotic mess.
• Their model works perfectly until that shockwave forms. Once the shock happens, the simple "surfboard" model needs to be upgraded to include more complex physics (like the "stickiness" mentioned above).

5. Why This Matters

The universe is about to get a lot more data. New telescopes will map billions of galaxies. To make sense of this data, we need better math.

• The Benefit: By using this "Pion" language, the authors found that the math becomes much cleaner. It reveals hidden symmetries (rules that the universe follows) that were previously hidden in the messy equations.
• The Consistency Check: They proved that their new math obeys the fundamental "laws of the universe" (consistency relations). If you tweak the universe slightly, their model predicts exactly how the ripples should change, and it matches what we see in supercomputer simulations.

Summary

The authors are essentially saying: "Stop trying to track every drop of water in the cosmic ocean. Instead, track the shape of the wave itself. It's simpler, it respects the universe's hidden rules, and it helps us understand how the 'traffic jams' of galaxies form without needing a supercomputer to solve every single collision."

They have built a new, simpler language to describe the messy, beautiful chaos of the universe, and they've shown that this language works just as well as the old, complicated one—but with much less headache.

Technical Summary

1. Problem Statement

Large Scale Structure (LSS) surveys are reaching unprecedented precision, offering a pathway to constrain cosmological parameters (e.g., dark energy, primordial non-Gaussianity, modifications to General Relativity). However, extracting this information is hindered by the nonlinear evolution of matter inhomogeneities.

• Limitations of Current Methods: Standard $N$-body simulations are computationally expensive on large scales. Perturbative methods (Standard Perturbation Theory) work well in the weakly nonlinear regime but fail as nonlinear interactions dominate.
• Theoretical Gap: While the Effective Field Theory of Large Scale Structure (EFTofLSS) has been successful in treating matter as an effective fluid with density, velocity, and stress tensors, the formalism involves multiple coupled variables ($\delta, \vec{v}, \Phi$). The authors seek a more unified, symmetry-based description that organizes perturbation theory systematically while keeping spacetime symmetries manifest.

2. Methodology

The authors develop an EFT based on a single scalar degree of freedom, the "pion" field ($\pi$), which corresponds to the velocity potential of the pressureless matter fluid in a $\Lambda$CDM universe ($\vec{v} = \vec{\nabla}\pi$).

A. Theoretical Framework

• Symmetry Origin: The $\pi$ field is identified as the Goldstone boson arising from the spontaneous breaking of spacetime diffeomorphism symmetries (specifically time-dependent shifts and translations) by the matter fluid.
• Derivation:
  1. Fluid Equations: Starting from the Navier-Stokes, continuity, and Poisson equations for a pressureless fluid, they derive a single equation of motion for $\pi$ (Eq. 2.11). This equation is nonlinear and non-local (due to the inverse Laplacian).
  2. Relativistic Action: They derive the same field from a relativistic action for an irrotational fluid, confirming that $\pi$ is the natural variable for organizing the EFT.
• EFT Expansion: They construct the effective action by adding terms invariant under the broken symmetries.
  • Leading Order (LO): The standard nonlinear equation of motion.
  • Next-to-Leading Order (NLO): They introduce EFT operators involving sound speed ($c_s$) and viscosity ($c_v$) parameters, represented by terms like $(\nabla^2 \pi)^2$ and $(\nabla^2 \pi')^2$. These account for the effects of small-scale physics (shell crossing, virialization) integrated out from the large-scale description.

B. Analytical Calculations

• Perturbation Theory: The authors solve the equation of motion perturbatively ($\pi = \pi_1 + \pi_2 + \pi_3 + \dots$) in the matter-dominated and $\Lambda$CDM regimes.
• Power Spectrum: They calculate the NLO corrections to the matter power spectrum ($P(k)$) using Feynman diagrams (one-loop contributions).
• Consistency Relations: They verify that their perturbative calculations satisfy the consistency relations derived from the spontaneously broken spacetime symmetries (specifically, that certain $k^2$ corrections vanish in the soft limit).

C. Numerical Simulations

• PLASTIC Code: The authors developed a custom code (Pion LAgrangian for STructure In Cosmology) to numerically evolve the $\pi$ field in 1D and 3D. They utilize Fast Fourier Transforms (FFT) to handle the non-local inverse Laplacian operator efficiently ($O(N^3 \log N)$).
• N-body Comparison: They compare their analytical and $\pi$-field simulation results against standard $N$-body simulations (using GADGET-4) to measure the EFT coefficients ($c_s, c_v$) directly from the data.

3. Key Contributions

1. Unified Variable Formalism: Demonstrated that the complex interplay of density ($\delta$), velocity ($\vec{v}$), and gravitational potential ($\Phi$) can be reduced to a single scalar field $\pi$ without loss of information (in the irrotational limit), simplifying the EFT structure.
2. Systematic NLO Derivation: Provided a systematic derivation of NLO corrections to the power spectrum within the pion field framework, explicitly showing how gravitational collapse and EFT counterterms contribute.
3. Consistency Relation Verification: Proved that the perturbative expansion of the pion field naturally satisfies the consistency relations for spontaneously broken spacetime symmetries, validating the theoretical consistency of the approach.
4. Numerical Implementation: Created and released the PLASTIC simulation code, demonstrating that the nonlinear evolution of the pion field can be simulated efficiently and reproduces the formation of shocks (shell crossing) consistent with standard fluid dynamics.
5. Parameter Extraction: Successfully measured EFT coefficients (sound speed and viscosity) from $N$-body simulations, finding values consistent with previous literature but derived within the pion field formalism.

4. Results

• Power Spectrum Corrections:
  • In the matter-dominated regime, the NLO corrections scale as $(k/k_{NL})^4$ for $k \ll k_{NL}$, consistent with mass and momentum conservation.
  • For $k \gg k_{mr}$ (matter-radiation equality scale), the corrections include logarithmic factors and scale as $(k_{mr}/k)^4$.
  • The inclusion of EFT terms (sound speed and viscosity) significantly improves the fit between perturbation theory and simulations up to $k \sim 0.15 \, h/\text{Mpc}$.
• Shock Formation: Numerical simulations show that the $\pi$ field develops shocks (where $\nabla \pi$ transitions from positive to negative) at the same scale where the single-stream approximation breaks down ($\delta \sim 1$). Beyond this point, the single-field description fails, requiring additional degrees of freedom (vorticity).
• EFT Coefficients:
  • Measured $c_s^2$ is negative and $c_v$ is positive, indicating that small-scale nonlinearities induce instability and momentum diffusion on large scales.
  • The time dependence of the combined EFT coefficient $\tilde{c}^2_{comb}$ scales approximately as the square of the linear growth factor ($D^2$), validating the ansatz used in the perturbative calculations.
• Vector Potential: The ratio of the vector potential power spectrum to the scalar potential scales as $1/k^5$ at late times, consistent with the stochastic generation of vorticity in nonlinear halos.

5. Significance

• Theoretical Clarity: The pion field formalism offers a more elegant and symmetry-motivated organization of LSS perturbation theory, making the scaling of nonlinear terms and the realization of consistency relations more transparent than in the standard multi-variable fluid approach.
• New Observables: By expressing observables in terms of $\pi$, the paper suggests new variables for analyzing LSS data that might be more robust or sensitive to specific primordial physics (e.g., non-Gaussianity).
• Bridge Between Theory and Data: The work provides a concrete link between high-energy physics-inspired EFT techniques and cosmological simulations. It demonstrates that EFT coefficients can be reliably extracted from simulations and used to correct perturbative predictions, extending the validity of analytical methods into the weakly nonlinear regime.
• Future Directions: The paper highlights that while the single-field picture breaks down at small scales (shell crossing), it serves as a robust foundation. Future work could generalize this to include vorticity and velocity dispersion, potentially leading to a more complete effective theory for the full nonlinear regime.

In summary, this paper establishes the pion field as a powerful, symmetry-based framework for modeling Large Scale Structure, successfully deriving NLO corrections, validating them against simulations, and providing a new toolset for analyzing the nonlinear universe.