Shocks in the Early Universe

TL;DR

A spectrum of Gaussian, linear, adiabatic, scalar, growing mode perturbations not only creates acoustic oscillations of the kind observed on very large scales today, it also leads to the production of shocks in the radiation fluid of the very early Universe.

Abstract

We point out a surprising consequence of the usually assumed initial conditions for cosmological perturbations. Namely, a spectrum of Gaussian, linear, adiabatic, scalar, growing mode perturbations not only creates acoustic oscillations of the kind observed on very large scales today, it also leads to the production of shocks in the radiation fluid of the very early universe. Shocks cause departures from local thermal equilibrium as well as creating vorticity and gravitational waves. For a scale-invariant spectrum and standard model physics, shocks form for temperatures $1$ GeV$<T<10^{7}$ GeV. For more general power spectra, such as have been invoked to form primordial black holes, shock formation and the consequent gravitational wave emission provides a signal detectable by current and planned gravitational wave experiments, allowing them to strongly constrain conditions present in the primordial universe as early as $10^{-30}$ seconds after the big bang.

Figures (3)

• Figure 1: Simulation showing cosmological initial conditions (left) evolving into shocks (right). The magnitude of the gradient of the energy density is shown in greyscale. The initial spectrum is scale-invariant and cut off at ${1\over 128}$ the box size, with rms amplitude $\epsilon=.02$. Movie available at weblink.
• Figure 2: The growing mode perturbation, in a radiation-dominated universe, in conformal Newtonian gauge. The density perturbation $\delta_k(t)$ (black), the Newtonian potential $\Phi$ (red) and the flat spacetime approximation to $\delta_k(t)$ (blue) are plotted against $kc_s t$.
• Figure 3: Entropy, in units of its equilibrium value, versus the time $t$, in units of the sound-crossing time, for $512^3$ simulations of a perfect radiation fluid with cosmological initial conditions as in (\ref{['eq1']}). The red dashed curve is a fit to the $\epsilon=.05$ curve using (\ref{['eq9']}) with $C={1\over 4}$. For $\epsilon=0.1$, $t$ has been doubled and the entropy deficit rescaled by a quarter to verify good agreement with (\ref{['eq9']}).