🔭 astrophysics

Primordial Asymmetries, Primordial Equation of State & Primordial Black Holes

This paper investigates how primordial lepton and baryon asymmetries influence the early Universe's equation of state and thermal history, ultimately determining their effects on the primordial black hole mass spectrum and the resulting gravitational wave signals detectable by current ground-based interferometers.

Original authors: Maël Gonin, Julien Froustey, Albert Escrivà, Alberto Magaraggia, Günther Hasinger, Florian Kühnel

Published 2026-08-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Maël Gonin, Julien Froustey, Albert Escrivà, Alberto Magaraggia, Günther Hasinger, Florian Kühnel

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The early universe was a place of extreme heat and density, a time when the fundamental forces of nature were still settling into the forms we know today. In this primordial soup, particles were not just floating randomly; they carried subtle imbalances, tiny differences between the amount of matter and antimatter that would eventually shape the entire cosmos. Scientists have long known that the universe contains a slight excess of matter over antimatter, a fact that allowed stars, planets, and life to exist. However, there is another, less visible imbalance involving neutrinos, the ghostly particles that rarely interact with anything. These "lepton asymmetries" are difficult to measure directly, yet they leave a profound mark on how the universe expanded and cooled. By studying the thermal history of this ancient era, researchers can trace how these hidden imbalances influenced the formation of the first structures, including the mysterious primordial black holes that may still lurk in the dark corners of the cosmos.

A team of researchers has now mapped out how these early imbalances changed the behavior of the universe's expansion, specifically focusing on a period when the temperature dropped from ten billion degrees down to a few thousand degrees. This era covers critical moments when the universe underwent phase transitions, similar to water freezing into ice, but involving the fundamental building blocks of matter. As the universe cooled, quarks and gluons, which were previously free-moving, locked together to form protons and neutrons. Shortly after, heavy particles like muons and taus vanished, and finally, electrons and positrons annihilated each other. The researchers calculated how the presence of lepton asymmetries altered the pressure and energy density of the universe during these transitions. They found that these asymmetries acted like a hidden hand, smoothing out the sharp changes in the universe's behavior that would have occurred in a perfectly balanced scenario.

To understand the impact of these changes, the team used a combination of advanced computer codes to track the chemical potentials of different particles as the universe cooled. They simulated four different scenarios for how these asymmetries might have been distributed among the three types of neutrinos. In the standard view, these imbalances are tiny and evenly spread. However, the researchers explored more extreme possibilities where one type of neutrino carried a much larger imbalance than the others. Their simulations revealed that when a large imbalance existed, particularly in the tau neutrino sector, it significantly changed the thermodynamics of the early universe. Instead of a sharp drop in the universe's ability to resist gravity during the quark-to-hadron transition, the presence of these asymmetries created a gentler slope. This subtle shift meant that the universe remained slightly more "stiff" against gravitational collapse than previously thought.

This change in the universe's stiffness has a direct consequence for the formation of primordial black holes. These objects are thought to have formed when dense regions of the early universe collapsed under their own gravity. The likelihood of such a collapse depends heavily on the pressure of the surrounding material; if the pressure drops too low, gravity wins, and a black hole is born. The researchers found that the smoothing effect of lepton asymmetries reduced the depth of the pressure drop during the quark transition. Consequently, the peak in the number of black holes forming at a specific mass—roughly the mass of our sun—became less pronounced. In scenarios with large asymmetries, the distribution of black hole masses became flatter and broader, spreading the population across a wider range of sizes rather than concentrating them in a single, sharp peak.

The study also looked at how these different mass distributions would appear to modern detectors. By simulating the gravitational waves produced when pairs of these primordial black holes merge, the researchers showed that the shape of the signal depends on the underlying mass spectrum. If the universe had significant lepton asymmetries, the resulting gravitational wave signals would be distributed differently across mass and mass ratios compared to the standard model. This offers a potential way for current and future gravitational wave observatories, like the LIGO and Virgo detectors, to indirectly detect these ancient asymmetries. By analyzing the masses of the black holes they observe, scientists could potentially infer whether the early universe was perfectly balanced or carried these hidden, uneven charges.

The researchers were careful to note that their findings rely on complex simulations of the early universe's thermal history and the statistical theory of how black holes form. They did not claim to have proven that these asymmetries exist, but rather demonstrated how they would alter the cosmic landscape if they did. Their work highlights a deep connection between the microscopic world of particle physics and the macroscopic structure of the cosmos. The subtle differences in the number of neutrinos versus antineutrinos in the first moments of the universe could have dictated the abundance and size of the black holes that populate the universe today. As gravitational wave astronomy matures, the data from these cosmic collisions may soon provide the first direct evidence of these primordial imbalances, turning the quiet whispers of the early universe into a clear voice for the laws of physics.

Technical Summary: Primordial Asymmetries, Primordial Equation of State & Primordial Black Holes

Problem Statement
The formation of Primordial Black Holes (PBHs) is driven by the gravitational collapse of large-amplitude primordial curvature perturbations during the radiation era. The abundance and mass spectrum of these PBHs are exponentially sensitive to the equation of state (EoS) of the early Universe, specifically the sound speed and pressure support during phase transitions. While the standard cosmological model assumes vanishing lepton asymmetries, recent constraints from Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB) allow for significant primordial lepton asymmetries (LAU). These asymmetries, along with the baryon asymmetry (BAU), introduce non-zero chemical potentials that alter the thermodynamics of the primordial plasma. The paper addresses the lack of a comprehensive understanding of how these non-zero asymmetries modify the cosmic EoS across key epochs (QCD transition, muon/tau annihilation, neutrino decoupling) and how these modifications subsequently reshape the PBH mass spectrum and associated Gravitational Wave (GW) signals.

Methodology
The authors employ a multi-code approach to track the evolution of the Universe's EoS from $T = 10$ GeV down to $T = 1$ keV, spanning the QCD transition, hadron gas phase, and neutrino decoupling.

Thermodynamic Evolution:
- $T > 37$ MeV: The code CosmicEoS is used, relying on lattice QCD simulations and a generalized Beth-Uhlenbeck (GBU) approach. It solves conservation equations for baryon, electric charge, and lepton flavor numbers to determine chemical potentials ( $\mu_B, \mu_Q, \mu_{\nu_\alpha}$ ) in the presence of asymmetries.
- $10 \le T \le 37$ MeV: The Thermal-FIST hadron resonance gas model is utilized. The authors ensure a smooth connection between CosmicEoS and Thermal-FIST at $T=37$ MeV to avoid discontinuities.
- $T < 10$ MeV: The codes nudec_BSM (for standard symmetric cases) and NEVO (for asymmetric cases) are used. NEVO is critical for modeling neutrino decoupling and collective neutrino oscillations, which redistribute flavor asymmetries and affect the EoS.
- Models: Four specific LAU scenarios are tested: a "Standard" sphaleron-induced model, two "Type I" models with large muon/tau asymmetries, and a "BKOS-like" model based on previous literature.
PBH Mass Spectrum Calculation:
- Peak Theory: The PBH abundance is calculated using peak theory, accounting for the Laplacian of the curvature fluctuation ( $\Delta\zeta$ ).
- Collapse Thresholds: Critical thresholds for gravitational collapse ( $\mu_{2,c}$ and $\mathcal{C}_c$ ) are determined using fully relativistic numerical simulations (SPriBHoS codes) adapted for the specific time-dependent EoS derived above.
- Critical Collapse: The mapping from horizon mass ( $M_H$ ) to PBH mass ( $M_{PBH}$ ) incorporates critical collapse scaling relations.
- GW Signal: The paper calculates the merger rate of PBH binaries formed in the late Universe (gravitational capture in dark matter halos) and estimates the resulting GW signal detectable by ground-based interferometers (e.g., Advanced LIGO).

Key Contributions

Unified EoS with Asymmetries: The study provides the first continuous calculation of the cosmic EoS from the QCD transition down to neutrino decoupling that consistently incorporates non-zero lepton and baryon asymmetries, including the effects of neutrino oscillations on flavor redistribution.
Impact on Phase Transitions: The authors demonstrate that large lepton asymmetries (specifically $\mathcal{L}_\tau$ ) significantly alter the thermodynamic behavior during the QCD transition and particle annihilations ( $\tau^\pm, \mu^\pm, e^\pm$ ).
PBH Spectrum Modification: By combining the modified EoS with updated collapse thresholds, the paper quantifies how LAU flattens the PBH mass spectrum and shifts the characteristic "QCD peak" to lower masses.
GW Constraints: The work connects these modified mass spectra to observable GW signals, showing how current interferometer data could constrain or hint at specific LAU configurations and spectral indices of curvature perturbations.

Results

Equation of State: In the presence of large LAU, the characteristic "dip" in the EoS parameter $w$ at the QCD transition is mitigated. This is because the chemical potentials increase the thermodynamic weight of neutrinos (which behave as pure radiation, $w=1/3$ ), pulling the total EoS back toward the radiation value.
Oscillation Effects: Neutrino oscillations effectively "wash out" large primordial flavor asymmetries, but the process is irreversible and creates entropy. The final state depends on the neutrino mass ordering (Normal vs. Inverted), with Normal Ordering (NO) retaining more asymmetry and thus a stronger modification to the EoS compared to Inverted Ordering (IO).
PBH Mass Spectra:
- The introduction of LAU generally flattens the PBH mass distribution.
- For models with large $\mathcal{L}_\tau$ (Type I), the primary peak associated with the QCD transition shifts to smaller masses ( $M_{PBH} \sim 10^{-1} - 1 M_\odot$ ).
- A secondary "shoulder" in the mass spectrum around $10 - 100 M_\odot$ (associated with $\mu^\pm$ annihilation) gains significant weight, particularly in the BKOS-like model.
- The double-peak structure seen in the EoS for the Type I-0.08 model does not translate into a distinct double peak in the PBH spectrum due to the smoothing effects of the nonlinear collapse dynamics and critical scaling.
Gravitational Waves: The modified mass spectra alter the predicted distribution of PBH mergers. The authors show that the relative detection weight in the $(\log_{10} M_B, q)$ plane changes, potentially allowing current GW observations to distinguish between standard and asymmetric scenarios, provided the spectral index of curvature perturbations is favorable.

Significance and Claims
The paper claims that the thermal history of the early Universe, specifically the presence of primordial asymmetries, leaves a "remarkable feature" in the form of the PBH mass distribution. By linking the microphysics of the QCD transition and neutrino decoupling to the macroscopic formation of black holes, the study suggests that:

PBHs as Probes: The mass distribution of PBHs can serve as a probe for fundamental physics in the pre-recombination era, specifically the magnitude of lepton asymmetries and the nature of the QCD phase transition.
GW Constraints: Current ground-based GW observations (LVK) have the potential to constrain or provide positive evidence for specific LAU scenarios and curvature perturbation spectral indices, provided the PBHs constitute a significant fraction of dark matter.
Refinement of Models: The work highlights that assuming a standard, symmetric early Universe may lead to incorrect predictions for PBH abundances, and that future constraints on PBHs must account for the variety of EoS realizations allowed by BBN and CMB data.

The authors remain modest regarding definitive conclusions, noting that current constraints on PBHs often assume monochromatic mass distributions, whereas their results produce broad, flattened spectra. They emphasize that while their models are illustrative benchmarks, the interplay between thermal history and PBH formation is a robust feature of nature that requires further exploration to fully resolve observational tensions.

More like this